title: Observation of Wakefields and Resonances in Coherent Synchrotron Radiation

abstract: We report on high resolution measurements of resonances in the spectrum of coherent synchrotron radiation (CSR) at the Canadian Light Source (CLS). The resonances permeate the spectrum at wavenumber intervals of $0.074 ~\textrm{cm}^{-1}$, and are highly stable under changes in the machine setup (energy, bucket filling pattern, CSR in bursting or continuous mode). Analogous resonances were predicted long ago in an idealized theory as eigenmodes of a smooth toroidal vacuum chamber driven by a bunched beam moving on a circular orbit. A corollary of peaks in the spectrum is the presence of pulses in the wakefield of the bunch at well defined spatial intervals. Through experiments and further calculations we elucidate the resonance and wakefield mechanisms in the CLS vacuum chamber, which has a fluted form much different from a smooth torus. The wakefield is observed directly in the 30-110 GHz range by RF diodes, and indirectly by an interferometer in the THz range. The wake pulse sequence found by diodes is less regular than in the toroidal model, and depends on the point of observation, but is accounted for in a simulation of fields in the fluted chamber. Attention is paid to polarization of the observed fields, and possible coherence of fields produced in adjacent bending magnets. Low frequency wakefield production appears to be mainly local in a single bend, but multi-bend effects cannot be excluded entirely, and could play a role in high frequency resonances. New simulation techniques have been developed, which should be invaluable in further work.


title: Electromagnetically induced transparency for guided light in an atomic array outside an optical nanofiber

abstract: We study the propagation of guided light along an array of three-level atoms in the vicinity of an optical nanofiber under the condition of electromagnetically induced transparency. We examine two schemes of atomic levels and field polarizations where the guided probe field is quasilinearly polarized along the major or minor principal axis, which is parallel or perpendicular, respectively, to the radial direction of the atomic position. Our numerical calculations indicate that 200 cesium atoms in a linear array with a length of 100 $\mu$m at a distance of 200 nm from the surface of a nanofiber with a radius of 250 nm can slow down the speed of guided probe light by a factor of about $3.5\times 10^6$ (the corresponding group delay is about 1.17 $\mu$s). In the neighborhood of the Bragg resonance, a significant fraction of the guided probe light can be reflected back with a negative group delay. The reflectivity and the group delay of the reflected field do not depend on the propagation direction of the probe field. However, when the input guided light is quasilinearly polarized along the major principal axis, the transmittivity and the group delay of the transmitted field substantially depend on the propagation direction of the probe field. Under the Bragg resonance condition, an array of atoms prepared in an appropriate internal state can transmit guided light polarized along the major principal in one specific direction even in the limit of infinitely large atom numbers. The directionality of transmission of guided light through the array of atoms is a consequence of the existence of a longitudinal component of the guided light field as well as the ellipticity of both the field polarization and the atomic dipole vector.


title: Polynomial behavior of the Honda formal group law

abstract: This note provides the calculation of the formal group law $F(x,y)$ in modulo $p$ Morava $K$-theory at prime $p$ and $s>1$ as an element in $K(s)^*[x][[y]]$ and one application to relevant examples.


title: Resolving unoccupied electronic states with laser ARPES in bismuth-based cuprate superconductors

abstract: Angle-resolved photoemission spectroscopy (ARPES) is typically used to study only the occupied electronic band structure of a material. Here we use laser-based ARPES to observe a feature in bismuth-based superconductors that, in contrast, is related to the unoccupied states. Specifically, we observe a dispersive suppression of intensity cutting across the valence band, which, when compared with relativistic one-step calculations, can be traced to two final-state gaps in the bands 6 eV above the Fermi level. This finding opens up possibilities to bring the ultra-high momentum resolution of existing laser-ARPES instruments to the unoccupied electron states. For cases where the final-state gap is not the object of study, we find that its effects can be made to vanish under certain experimental conditions.


title: Extracting Hidden Groups and their Structure from Streaming Interaction Data

abstract: When actors in a social network interact, it usually means they have some general goal towards which they are collaborating. This could be a research collaboration in a company or a foursome planning a golf game. We call such groups \emph{planning groups}. In many social contexts, it might be possible to observe the \emph{dyadic interactions} between actors, even if the actors do not explicitly declare what groups they belong too. When groups are not explicitly declared, we call them \emph{hidden groups}. Our particular focus is hidden planning groups. By virtue of their need to further their goal, the actors within such groups must interact in a manner which differentiates their communications from random background communications. In such a case, one can infer (from these interactions) the composition and structure of the hidden planning groups. We formulate the problem of hidden group discovery from streaming interaction data, and we propose efficient algorithms for identifying the hidden group structures by isolating the hidden group's non-random, planning-related, communications from the random background communications. We validate our algorithms on real data (the Enron email corpus and Blog communication data). Analysis of the results reveals that our algorithms extract meaningful hidden group structures.


title: Measurement of J/psi and psi(2S) prompt double-differential cross sections in pp collisions at sqrt(s) = 7 TeV

abstract: The double-differential cross sections of promptly produced J/psi and psi(2S) mesons are measured in pp collisions at sqrt(s) = 7 TeV, as a function of transverse momentum pt and absolute rapidity abs(y). The analysis uses J/psi and psi(2S) dimuon samples collected by CMS, corresponding to integrated luminosities of 4.55 and 4.90 inverse femtobarns, respectively. The results are based on a two-dimensional analysis of the dimuon invariant mass and decay length, and extend to pt = 120 and 100 GeV for the J/psi and psi(2S), respectively, when integrated over the interval abs(y) < 1.2. The ratio of the psi(2S) to J/psi cross sections is also reported for abs(y) < 1.2, over the range 10 < pt < 100 GeV. These are the highest pt values for which the cross sections and ratio have been measured.


title: Towards Biologically Plausible Deep Learning

abstract: Neuroscientists have long criticised deep learning algorithms as incompatible with current knowledge of neurobiology. We explore more biologically plausible versions of deep representation learning, focusing here mostly on unsupervised learning but developing a learning mechanism that could account for supervised, unsupervised and reinforcement learning. The starting point is that the basic learning rule believed to govern synaptic weight updates (Spike-Timing-Dependent Plasticity) arises out of a simple update rule that makes a lot of sense from a machine learning point of view and can be interpreted as gradient descent on some objective function so long as the neuronal dynamics push firing rates towards better values of the objective function (be it supervised, unsupervised, or reward-driven). The second main idea is that this corresponds to a form of the variational EM algorithm, i.e., with approximate rather than exact posteriors, implemented by neural dynamics. Another contribution of this paper is that the gradients required for updating the hidden states in the above variational interpretation can be estimated using an approximation that only requires propagating activations forward and backward, with pairs of layers learning to form a denoising auto-encoder. Finally, we extend the theory about the probabilistic interpretation of auto-encoders to justify improved sampling schemes based on the generative interpretation of denoising auto-encoders, and we validate all these ideas on generative learning tasks.


title: Liouville type theorems for stable solutions of the weighted elliptic system

abstract: We examine the weighted elliptic system \begin{equation*} \begin{cases} -\Delta u=(1+|x|^2)^{\frac{\alpha}{2}} v,\\ -\Delta v=(1+|x|^2)^{\frac{\alpha}{2}} u^p, \end{cases} \quad \mbox{in}\;\ \mathbb{R}^N, \end{equation*}where $N \ge 5$, $p>1$ and $\alpha >0$. We prove Liouville type results for the classical positive (nonnegative) stable solutions in dimension $N<\ell+\dfrac{\alpha (\ell-2)}{2}$ ($N <\ell+\dfrac{\alpha (\ell-2)(p+3)}{4(p+1)}$) and $\ell \ge 5$, $p \in (1,p_*(\ell))$. In particular, for any $p>1$ and $\alpha > 0$, we obtain the nonexistence of classical positive (nonnegative) stable solutions for any $N \le 12+5 \alpha$ ($N\le 12+\dfrac{5\alpha (p+3)}{2(p+1)}$).


title: New calculation of antiproton production by cosmic ray protons and nuclei

abstract: A dramatic increase in the accuracy and statistics of space-borne cosmic ray (CR) measurements has yielded several breakthroughs over the last several years. The most puzzling is the rise in the positron fraction above ~10 GeV over the predictions of the propagation models assuming pure secondary production. The accuracy of the antiproton production cross section is critical for astrophysical applications and searches for new physics since antiprotons in CRs seem to hold the keys to many puzzles including the origin of those excess positrons. However, model calculations of antiproton production in CR interactions with interstellar gas are often employing parameterizations that are out of date or are using outdated physical concepts. That may lead to an incorrect interpretation of antiproton data which could have broad consequences for other areas of astrophysics. In this work, we calculate antiproton production in pp-, pA-, and AA-interactions using EPOS-LHC and QGSJET-II-04, two of the most advanced Monte Carlo (MC) generators tuned to numerous accelerator data including those from the Large Hadron Collider (LHC). We show that the antiproton yields obtained with these MC generators differ by up to an order of magnitude from yields of parameterizations commonly used in astrophysics.


title: The Photomultiplier Tube Calibration System of the MicroBooNE Experiment

abstract: We report on the design and construction of a LED-based fiber calibration system for large liquid argon time projection detectors. This system was developed to calibrate the optical systems of the MicroBooNE experiment. As well as detailing the materials and installation procedure, we provide technical drawings and specifications so that the system may be easily replicated in future LArTPC detectors.


title: An Exercise (?) in Fourier Analysis on the Heisenberg Group

abstract: Let H(n) be the group of 3x3 uni-uppertriangular matrices with entries in Z/nZ, the integers mod n. We show that the simple random walk converges to the uniform distribution in order n^2 steps. The argument uses Fourier analysis and is surprisingly challenging. It introduces novel techniques for bounding the spectrum which are useful for a variety of walks on a variety of groups.


title: Energy shift of magnons in a ferromagnetic spinor-dipolar Bose-Einstein condensate

abstract: Motivated by the recent experiment performed by the Berkeley group [G. E. Marti {\it et al.}, Phys. Rev. Lett. {\bf 113}, 155302 (2014)], we consider the dynamics of magnons in a spin-1 spinor-dipolar Bose-Einstein condensate, using mean-field theory. We show that the effective mass of a magnon is increased by the magnetic dipole-dipole interaction, as observed in the experiment. The magnon mass is also decreased by changing the direction of the magnetic field. The increase and decrease in the magnon mass manifest themselves in the acceleration of the magnons.


title: Transverse-momentum-dependent wave functions with Glauber gluons in $B \to \pi \pi$, $\rho \rho$ decays

abstract: We investigate the Glauber-gluon effect on the $B \to \pi \pi$ and $\rho \rho$ decays, which is introduced via a convolution of a universal Glauber phase factor with transverse-momentum-dependent(TMD) meson wave functions in the $k_T$ factorization theorem. For an appropriate parametrization of the Glauber phase, it is observed that a TMD wave function for the pion ($\rho$ meson) with a weak (strong) falloff in parton transverse momentum $k_T$ leads to significant (moderate) modification of the $B^0 \to \pi^0 \pi^0$ ($B^0 \to \rho^0 \rho^0$) branching ratio: the former (latter) is enhanced (reduced) by about a factor of 2 (15\%). This observation is consistent with the dual role of the pion as a massless Nambu-Goldstone boson and as a $q \bar q$ bound state, which requires a tighter spatial distribution of its leading Fock state relative to higher Fock states. The agreement between the theoretical predictions and the data for all the $B \to \pi \pi$ and $\rho^0 \rho^0$ branching ratios is then improved simultaneously, and it is possible to resolve the $B \to \pi \pi$ puzzle.


title: A Distributional Representation Model For Collaborative Filtering

abstract: In this paper, we propose a very concise deep learning approach for collaborative filtering that jointly models distributional representation for users and items. The proposed framework obtains better performance when compared against current state-of-art algorithms and that made the distributional representation model a promising direction for further research in the collaborative filtering.


title: Think Different: Applying the Old Macintosh Mantra to the Computability of the SUSY Auxiliary Field Problem

abstract: Starting with valise supermultiplets obtained from 0-branes plus field redefinitions, valise adinkra networks, and the "Garden Algebra," we discuss an architecture for algorithms that (starting from on-shell theories and, through a well-defined computation procedure), search for off-shell completions. We show in one dimension how to directly attack the notorious "off-shell auxiliary field" problem of supersymmetry with algorithms in the adinkra network-world formulation.


title: Coexistence of Wi-Fi and Heterogeneous Small Cell Networks Sharing Unlicensed Spectrum

abstract: As two major players in terrestrial wireless communications, Wi-Fi systems and cellular networks have different origins and have largely evolved separately. Motivated by the exponentially increasing wireless data demand, cellular networks are evolving towards a heterogeneous and small cell network architecture, wherein small cells are expected to provide very high capacity. However, due to the limited licensed spectrum for cellular networks, any effort to achieve capacity growth through network densification will face the challenge of severe inter-cell interference. In view of this, recent standardization developments have started to consider the opportunities for cellular networks to use the unlicensed spectrum bands, including the 2.4 GHz and 5 GHz bands that are currently used by Wi-Fi, Zigbee and some other communication systems. In this article, we look into the coexistence of Wi-Fi and 4G cellular networks sharing the unlicensed spectrum. We introduce a network architecture where small cells use the same unlicensed spectrum that Wi-Fi systems operate in without affecting the performance of Wi-Fi systems. We present an almost blank subframe (ABS) scheme without priority to mitigate the co-channel interference from small cells to Wi-Fi systems, and propose an interference avoidance scheme based on small cells estimating the density of nearby Wi-Fi access points to facilitate their coexistence while sharing the same unlicensed spectrum. Simulation results show that the proposed network architecture and interference avoidance schemes can significantly increase the capacity of 4G heterogeneous cellular networks while maintaining the service quality of Wi-Fi systems.


title: The extragalactic background light, the Hubble constant, and anomalies: conclusions from 20 years of TeV gamma-ray observations

abstract: Ground-based observatories have been collecting 0.2 - 20 TeV gamma rays from blazars for about twenty years. These gamma rays can experience absorption along the line of sight due to interactions with the extragalactic background light (EBL). In this paper, we show that the gamma-ray optical depth can be reduced to the convolution product of an EBL kernel with the EBL intensity, assuming a particular form for the EBL evolution. We extract the absorption signal from the most extensive set of TeV spectra from blazars collected so far and unveil a broad-band EBL spectrum from mid-ultraviolet to far infrared. This spectrum is in good agreement with the accumulated emission of galaxies, constraining unresolved populations of sources. We propose a data-driven estimate of the Hubble constant based on the comparison of local and gamma-ray measurements of the EBL. After setting stringent upper-limits on the redshift of four TeV blazars, we investigate the 106 gamma-ray spectra in our sample and find no significant evidence for anomalies. The intrinsic TeV spectra are not harder than their GeV counterpart, and no spectral upturn is visible at the highest optical depths. Finally, we investigate a modification of the pair-creation threshold due to Lorentz invariance violation. A mild excess prevents us from ruling out an effect at the Planck energy, and we constrain for the first time the energy scale of the modification to values larger than sixty percent of the Planck energy.


title: A new parameter in attractor single-field inflation

abstract: We revisit the notion of slow-roll in the context of general single-field inflation. As a generalization of slow-roll dynamics, we consider an inflaton $\phi$ in an attractor phase where the time derivative of $\phi$ is determined by a function of $\phi$, $\dot\phi=\dot\phi(\phi)$. In other words, we consider the case when the number of $e$-folds $N$ counted backward in time from the end of inflation is solely a function of $\phi$, $N=N(\phi)$. In this case, it is found that we need a new independent parameter to properly describe the dynamics of the inflaton field in general, in addition to the standard parameters conventionally denoted by $\epsilon$, $\eta$, $c_s^2$ and $s$. Two illustrative examples are presented to discuss the non-slow-roll dynamics of the inflaton field consistent with observations.


title: Nonparametric regression using needlet kernels for spherical data

abstract: Needlets have been recognized as state-of-the-art tools to tackle spherical data, due to their excellent localization properties in both spacial and frequency domains. This paper considers developing kernel methods associated with the needlet kernel for nonparametric regression problems whose predictor variables are defined on a sphere. Due to the localization property in the frequency domain, we prove that the regularization parameter of the kernel ridge regression associated with the needlet kernel can decrease arbitrarily fast. A natural consequence is that the regularization term for the kernel ridge regression is not necessary in the sense of rate optimality. Based on the excellent localization property in the spacial domain further, we also prove that all the $l^{q}$ $(01\leq q < \infty)$ kernel regularization estimates associated with the needlet kernel, including the kernel lasso estimate and the kernel bridge estimate, possess almost the same generalization capability for a large range of regularization parameters in the sense of rate optimality. This finding tentatively reveals that, if the needlet kernel is utilized, then the choice of $q$ might not have a strong impact in terms of the generalization capability in some modeling contexts. From this perspective, $q$ can be arbitrarily specified, or specified merely by other no generalization criteria like smoothness, computational complexity, sparsity, etc..


title: Computationally Tractable Algorithms for Finding a Subset of Non-defective Items from a Large Population

abstract: In the classical non-adaptive group testing setup, pools of items are tested together, and the main goal of a recovery algorithm is to identify the "complete defective set" given the outcomes of different group tests. In contrast, the main goal of a "non-defective subset recovery" algorithm is to identify a "subset" of non-defective items given the test outcomes. In this paper, we present a suite of computationally efficient and analytically tractable non-defective subset recovery algorithms. By analyzing the probability of error of the algorithms, we obtain bounds on the number of tests required for non-defective subset recovery with arbitrarily small probability of error. Our analysis accounts for the impact of both the additive noise (false positives) and dilution noise (false negatives). By comparing with the information theoretic lower bounds, we show that the upper bounds on the number of tests are order-wise tight up to a $\log^2K$ factor, where $K$ is the number of defective items. We also provide simulation results that compare the relative performance of the different algorithms and provide further insights into their practical utility. The proposed algorithms significantly outperform the straightforward approaches of testing items one-by-one, and of first identifying the defective set and then choosing the non-defective items from the complement set, in terms of the number of measurements required to ensure a given success rate.


title: Human Factors in Agile Software Development

abstract: Through our four years experiments on students' Scrum based agile software development (ASD) process, we have gained deep understanding into the human factors of agile methodology. We designed an agile project management tool - the HASE collaboration development platform to support more than 400 students self-organized into 80 teams to practice ASD. In this thesis, Based on our experiments, simulations and analysis, we contributed a series of solutions and insights in this researches, including 1) a Goal Net based method to enhance goal and requirement management for ASD process, 2) a novel Simple Multi-Agent Real-Time (SMART) approach to enhance intelligent task allocation for ASD process, 3) a Fuzzy Cognitive Maps (FCMs) based method to enhance emotion and morale management for ASD process, 4) the first large scale in-depth empirical insights on human factors in ASD process which have not yet been well studied by existing research, and 5) the first to identify ASD process as a human-computation system that exploit human efforts to perform tasks that computers are not good at solving. On the other hand, computers can assist human decision making in the ASD process.


title: A combined experimental and theoretical study of the electronic and vibrational properties of bulk and few-layer Td-WTe2

abstract: The recent discovery of non-saturating giant positive magnetoresistance in Td-WTe2 has aroused great interest in this material. We have studied the structural, electronic and vibrational properties of bulk and few-layer Td-WTe2 experimentally and theoretically. Spin-orbit coupling is found to govern the semi-metallic character of Td-WTe2. Its structural link with the metallic 1T form provides an understanding of its structural stability. We observe a metal to insulator transition and a change in the sign of the Seebeck coefficient around 373 K. Lattice vibrations in Td-WTe2 have been analyzed by first principle calculations. Out of the 33 possible zone-center Raman active modes, five distinct Raman bands are observed around 112, 118, 134, 165 and 212 cm-1 in bulk Td-WTe2. Based on symmetry analysis and the calculated Raman tensors, we assign the intense bands at 165 cm-1 and 212 cm-1 to the A_1^' and A_1^" modes respectively. We have examined the effect of temperature and the number of layers on the Raman spectrum. Most of the bands of Td-WTe2 stiffen, and the ratio of the integrated intensities of the A_1^" to A_1^' bands decreases in the few-layer sample, while all the bands soften in both bulk and few-layer samples with increasing temperature.


title: Digital confocal microscopy through a multimode fiber

abstract: Acquiring high-contrast optical images deep inside biological tissues is still a challenging problem. Confocal microscopy is an important tool for biomedical imaging since it improves image quality by rejecting background signals. However, it suffers from low sensitivity in deep tissues due to light scattering. Recently, multimode fibers have provided a new paradigm for minimally invasive endoscopic imaging by controlling light propagation through them. Here we introduce a combined imaging technique where confocal images are acquired through a multimode fiber. We achieve this by digitally engineering the excitation wavefront and then applying a virtual digital pinhole on the collected signal. In this way, we are able to acquire images through the fiber with significantly increased contrast. With a fiber of numerical aperture 0.22, we achieve a lateral resolution of 1.5um, and an axial resolution of 12.7um. The point-scanning rate is currently limited by our spatial light modulator (20Hz).


title: Elastic Coulomb-levitation: why is ice so slippery?

abstract: The elastic, less dense, polarized, and thermally stable supersolid skin lubricates ice. Molecular undercoordination shortens the H-O bond and lengthens the O:H nonbond through O-O repulsion, which is associated with low-frequency and high-magnitude of O:H vibration and a dual O-O polarization. The softer O:H springs attached with stronger molecular dipoles provide forces levitating objects sliding on ice, like Maglev or Hovercraft.


title: Probabilistic Models for High-Order Projective Dependency Parsing

abstract: This paper presents generalized probabilistic models for high-order projective dependency parsing and an algorithmic framework for learning these statistical models involving dependency trees. Partition functions and marginals for high-order dependency trees can be computed efficiently, by adapting our algorithms which extend the inside-outside algorithm to higher-order cases. To show the effectiveness of our algorithms, we perform experiments on three languages---English, Chinese and Czech, using maximum conditional likelihood estimation for model training and L-BFGS for parameter estimation. Our methods achieve competitive performance for English, and outperform all previously reported dependency parsers for Chinese and Czech.


title: A comparison of automorphic and Artin L-series of GL(2)-type agreeing at degree one primes

abstract: Let $F/k$ be a cyclic extension of number fields of prime degree. Let $\rho$ be an irreducible $2$-dimensional representation of Artin type of the absolute Galois group of $F$, and $\pi$ a cuspidal automorphic representation of GL$_2(\mathbb A_F)$, such that the $L$-functions $L(s,\rho_v)$ and $L(s,\pi_v)$ agree at all (but finitely many of) the places $v$ of degree one over $k$. We prove in this case that we have the global identity $L(s,\rho)=L(s,\pi)$, with $\rho_v \leftrightarrow \pi_v$ being given by the local Langlands correspondence at all $v$. In particular, $\pi$ is tempered and $L(s,\rho)$ is entire.


title: Negative Group Velocity Pulse Propagation Through a Left-Handed Transmission Line

abstract: In this paper, the microwave pulse propagation transferred through a left-handed transmission line using Complementary Omega-Like Structures (COLS) loaded was studied. There was a stop band in transmission from 5.6GHz to 6.1GHz, and the anomalous dispersion was causes in this band. Negative group velocity corresponds to the case in which the peak of the pulse exited before the peak of the incident pulse had entered the sample. The negative group velocity reached (-0.27c~-1.85c).


title: Dark Energy: Observational Evidence and Theoretical Models

abstract: The book elucidates the current state of the dark energy problem and presents the results of the authors, who work in this area. It describes the observational evidence for the existence of dark energy, the methods and results of constraining of its parameters, modeling of dark energy by scalar fields, the space-times with extra spatial dimensions, especially Kaluza---Klein models, the braneworld models with a single extra dimension as well as the problems of positive definition of gravitational energy in General Relativity, energy conditions and consequences of their violation in the presence of dark energy. This monograph is intended for science professionals, educators and graduate students, specializing in general relativity, cosmology, field theory and particle physics.


title: Second order concentration on the sphere

abstract: Sharpened forms of the concentration of measure phenomenon for classes of functions on the sphere are developed in terms of Hessians of these functions.


title: Microwave probes Dipole Blockade and van der Waals Forces in a Cold Rydberg Gas

abstract: We show that microwave spectroscopy of a dense Rydberg gas trapped on a superconducting atom chip in the dipole blockade regime reveals directly the dipole-dipole many-body interaction energy spectrum. We use this method to investigate the expansion of the Rydberg cloud under the effect of repulsive van der Waals forces and the breakdown of the frozen gas approximation. This study opens a promising route for quantum simulation of many-body systems and quantum information transport in chains of strongly interacting Rydberg atoms.


title: New mechanism of kinetic exchange interaction induced by strong magnetic anisotropy

abstract: It is well known that the kinetic exchange interaction between single-occupied magnetic orbitals (s-s) is always antiferromagnetic, while between single- and double-occupied orbitals (s-d) is always ferromagnetic and much weaker. Here we show that the exchange interaction between strongly anisotropic doublets of lanthanides, actinides and transition metal ions with unquenched orbital momentum contains a new s-d kinetic contribution equal in strength with the s-s one. In non-collinear magnetic systems, this s-d kinetic mechanism can cause an overall ferromagnetic exchange interaction which can become very strong for transition metal ions. These findings are fully confirmed by DFT based analysis of exchange interaction in several Ln$^{3+}$ complexes.


title: Novel Electrostatically Doped Planar Field-Effect Transistor for High Temperature Applications

abstract: In this paper, we present experimental results and simulation data of an electrostatically doped and therefore voltage-programmable, planar, CMOS-compatible field-effect transistor (FET) structure. This planar device is based on our previously published Si-nanowire (SiNW) technology. Schottky barrier source/drain (S/D) contacts and a silicon-on-insulator (SOI) technology platform are the key features of this dual-gated but single channel universal FET. The combination of two electrically independent gates, one back-gate for S/D Schottky barrier modulation as well as channel formation to establish Schottky barrier FET (SBFET) operation and one front-gate forming a junctionless FET (JLFET) for actual current control, significantly increases the temperature robustness of the device.


title: Study of quark mass dependence of binding energy for light nuclei in 2+1 flavor lattice QCD

abstract: We investigate the formation of light nuclei with the nuclear mass number less than or equal to four in 2+1 flavor QCD using a non-perturbative improved Wilson quark and Iwasaki gauge actions. The quark mass is decreased from our previous work to the one corresponding to the pion mass of 0.30 GeV. In each multi-nucleon channel, the energy shift of the ground state relative to the assembly of free nucleons is calculated on two volumes, whose spatial extents are 4.3 fm and 5.8 fm. From the volume dependence of the energy shift, we distinguish a bound state of multi nucleons from an attractive scattering state. We find that all the ground states measured in this calculation are bound states. As in the previous studies at larger $m_\pi$, our result indicates that at $m_\pi = 0.30$ GeV the effective interaction between nucleons in the light nuclei is relatively stronger than the one in nature, since the results for the binding energies are larger than the experimental values and a bound state appears in the dineutron channel, which is not observed in experiment. Possible sources of systematic error in this calculation are discussed.


title: Weitzenb\"ock's Torsion, Fermi Coordinates and Adapted Frames

abstract: We study Weitzenb\"ock's torsion and discuss its properties. Specifically, we calculate the measured components of Weitzenb\"ock's torsion tensor for a frame field adapted to static observers in a Fermi normal coordinate system that we establish along the world line of an arbitrary accelerated observer in general relativity. A similar calculation is carried out in the standard Schwarzschild-like coordinates for static observers in the exterior Kerr spacetime; we then compare our results with the corresponding curvature components. Our work supports the contention that in the extended general relativistic framework involving both the Levi-Civita and Weitzenb\"ock connections, curvature and torsion provide complementary representations of the gravitational field.


title: Promotion and resignation in employee networks

abstract: Enterprises have put more and more emphasis on data analysis so as to obtain effective management advices. Managers and researchers are trying to dig out the major factors that lead to employees' promotion and resignation. Most previous analyses were based on questionnaire survey, which usually consists of a small fraction of samples and contains biases caused by psychological defense. In this paper, we successfully collect a data set consisting of all the employees' work-related interactions (action network, AN for short) and online social connections (social network, SN for short) of a company, which inspires us to reveal the correlations between structural features and employees' career development, namely promotion and resignation. Through statistical analysis and prediction, we show that the structural features of both AN and SN are correlated and predictive to employees' promotion and resignation, and the AN has higher correlation and predictability. More specifically, the in-degree in AN is the most relevant indicator for promotion; while the k-shell index in AN and in-degree in SN are both very predictive to resignation. Our results provide a novel and actionable understanding of enterprise management and suggest that to enhance the interplays among employees, no matter work-related or social interplays, can largely improve the loyalty of employees.


title: Evidence of Collective Flow in p-p Collision at LHC

abstract: The transverse momentum distributions of hadrons produced in p-p collisions at LHC energies of Root(sNN) = 0.9 TeV, 2.76 TeV and Root(sNN) = 7.0 TeV have been studied using a unified statistical thermal freeze-out model. A good agreement is seen between the theoretical calculations and experimental data. The extracted thermal freeze-out parameters reveal the presence of significant collective flow and a well defined temperature in the system thus indicating the formation of a thermally equilibrated hydrodynamic system in p-p collisions at LHC. The transverse momentum distributions of protons and Kaons produced in p-p collisions at Root(sNN) = 200 GeV have also been reproduced successfully. The model incorporates a longitudinal as well as transverse hydrodynamic flow. The contributions from heavier decay resonances have also been taken into account.


title: Spectrum Allocation for Multi-Operator Device-to-Device Communication

abstract: In order to harvest the business potential of device-to-device (D2D) communication, direct communication between devices subscribed to different mobile operators should be supported. This would also support meeting requirements resulting from D2D relevant scenarios, like vehicle-to-vehicle communication. In this paper, we propose to allocate the multi-operator D2D communication over dedicated cellular spectral resources contributed from both operators. Ideally, the operators should negotiate about the amount of spectrum to contribute, without revealing proprietary information to each other and/or to other parties. One possible way to do that is to use the sequence of operators' best responses, i.e., the operators make offers about the amount of spectrum to contribute using a sequential updating procedure until reaching consensus. Besides spectrum allocation, we need a mode selection scheme for the multi-operator D2D users. We use a stochastic geometry framework to capture the impact of mode selection on the distribution of D2D users and assess the performance of the best response iteration algorithm. With the performance metrics considered in the paper, we show that the best response iteration has a unique Nash equilibrium that can be reached from any initial strategy. In general, asymmetric operators would contribute unequal amounts of spectrum for multi-operator D2D communication. Provided that the multi-operator D2D density is not negligible, we show that both operators may experience significant performance gains as compared to the scheme without spectrum sharing.


title: Application of Deep Neural Network in Estimation of the Weld Bead Parameters

abstract: We present a deep learning approach to estimation of the bead parameters in welding tasks. Our model is based on a four-hidden-layer neural network architecture. More specifically, the first three hidden layers of this architecture utilize Sigmoid function to produce their respective intermediate outputs. On the other hand, the last hidden layer uses a linear transformation to generate the final output of this architecture. This transforms our deep network architecture from a classifier to a non-linear regression model. We compare the performance of our deep network with a selected number of results in the literature to show a considerable improvement in reducing the errors in estimation of these values. Furthermore, we show its scalability on estimating the weld bead parameters with same level of accuracy on combination of datasets that pertain to different welding techniques. This is a nontrivial result that is counter-intuitive to the general belief in this field of research.


title: On the Jeffrey-Kirwan residue of BCD-instantons

abstract: We apply the Jeffrey-Kirwan method to compute the multiple integrals for the $BCD$ type Nekrasov partition functions of four dimensional $\mathcal{N}=2$ supersymmetric gauge theories. We construct a graphical distinction rule to determine which poles are surrounded by their integration cycles. We compute the instanton correction of the "$Sp(0)$" pure super-Yang-Mills theory and find that $Z^{Sp(0)}_{k}=(-1)^{k}(2^{k}k!\varepsilon_{1}^{k}\varepsilon_{2}^{k})^{-1}$ for $k\le 8$, which resembles the formula $Z^{U(1)}_{k}=(k!\varepsilon_{1}^{k}\varepsilon_{2}^{k})^{-1}$ for the pure super-Yang-Mills theory with gauge group $U(1)$.


title: On the probability that all eigenvalues of Gaussian, Wishart, and double Wishart random matrices lie within an interval

abstract: We derive the probability that all eigenvalues of a random matrix $\bf M$ lie within an arbitrary interval $[a,b]$, $\psi(a,b)\triangleq\Pr\{a\leq\lambda_{\min}({\bf M}), \lambda_{\max}({\bf M})\leq b\}$, when $\bf M$ is a real or complex finite dimensional Wishart, double Wishart, or Gaussian symmetric/Hermitian matrix. We give efficient recursive formulas allowing the exact evaluation of $\psi(a,b)$ for Wishart matrices, even with large number of variates and degrees of freedom. We also prove that the probability that all eigenvalues are within the limiting spectral support (given by the Mar{\v{c}}enko-Pastur or the semicircle laws) tends for large dimensions to the universal values $0.6921$ and $0.9397$ for the real and complex cases, respectively. Applications include improved bounds for the probability that a Gaussian measurement matrix has a given restricted isometry constant in compressed sensing.

categories: math.ST cs.IT math-ph math.IT math.MP stat.TH


title: An Adaptive Sampling Approach to 3D Reconstruction of Weld Joint

abstract: We present an adaptive sampling approach to 3D reconstruction of the welding joint using the point cloud that is generated by a laser sensor. We start with a randomized strategy to approximate the surface of the volume of interest through selection of a number of pivotal candidates. Furthermore, we introduce three proposal distributions over the neighborhood of each of these pivots to adaptively sample from their neighbors to refine the original randomized approximation to incrementally reconstruct this welding space. We prevent our algorithm from being trapped in a neighborhood via permanently labeling the visited samples. In addition, we accumulate the accepted candidates along with their selected neighbors in a queue structure to allow every selected sample to contribute to the evolution of the reconstructed welding space as the algorithm progresses. We analyze the performance of our adaptive sampling algorithm in contrast to the random sampling, with and without replacement, to show a significant improvement in total number of samples that are drawn to identify the region of interest, thereby expanding upon neighboring samples to extract the entire region in a fewer iterations and a shorter computation time.


title: On the $3$-representations of groups and the $2$-categorical Traces

abstract: To $2$-categorify the theory of group representations, we introduce the notions of the $3$-representation of a group in a strict $3$-category and the strict $2$-categorical action of a group on a strict $2$-category. We also $2$-categorify the concept of the trace by introducing the $2$-categorical trace of a $1$-endomorphism in a strict $3$-category. For a $3$-representation $\rho$ of a group $G$ and an element $f$ of $G$, the $2$-categorical trace $\mathbb{T}r_2 \rho_f $ is a category. Moreover, the centralizer of $f$ in $G$ acts categorically on this $2$-categorical trace. We construct the induced strict $2$-categorical action of a finite group, and show that the $2$-categorical trace $ \mathbb{T}r_2$ takes an induced strict $2$-categorical action into an induced categorical action of the initia groupoid. As a corollary, we get the $3$-character formula of the induced strict $2$-categorical action.


title: Black Hole Solutions in $R^2$ Gravity

abstract: We find static spherically symmetric solutions of scale invariant $R^2$ gravity. The latter has been shown to be equivalent to General Relativity with a positive cosmological constant and a scalar mode. Therefore, one expects that solutions of the $R^2$ theory will be identical to that of Einstein theory. Indeed, we find that the solutions of $R^2$ gravity are in one-to-one correspondence with solutions of General Relativity in the case of non-vanishing Ricci scalar. However, scalar-flat $R=0$ solutions are global minima of the $R^2$ action and they cannot in general be mapped to solutions of the Einstein theory. As we will discuss, the $R=0$ solutions arise in Einstein gravity as solutions in the tensionless, strong coupling limit $M_P\rightarrow 0$. As a further result, there is no corresponding Birkhoff theorem and the Schwarzschild black hole is by no means unique in this framework. In fact, $R^2$ gravity has a rich structure of vacuum static spherically symmetric solutions partially uncovered here. We also find charged static spherically symmetric backgrounds coupled to a $U(1)$ field. Finally, we provide the entropy and energy formulas for the $R^2$ theory and we find that entropy and energy vanish for scalar-flat backgrounds.


title: Group theory of Wannier functions providing the basis for a deeper understanding of magnetism and superconductivity

abstract: The paper presents the group theory of best localized and symmetry-adapted Wannier functions in a crystal of any given space group G or magnetic group M. Provided that the calculated band structure of the considered material is given and that the symmetry of the Bloch functions at all the points of symmetry in the Brillouin zone is known, the paper details whether or not the Bloch functions of particular energy bands can be unitarily transformed into best localized Wannier functions symmetry-adapted to the space group G, to the magnetic group M, or to a subgroup of G or M. In this context, the paper considers usual as well as spin-dependent Wannier functions, the latter representing the most general definition of Wannier functions. The presented group theory is a review of the theory published by one of the authors in several former papers and is independent of any physical model of magnetism or superconductivity. However, it is suggested to interpret the special symmetry of the best localized Wannier functions in the framework of a nonadiabatic extension of the Heisenberg model, the nonadiabatic Heisenberg model. On the basis of the symmetry of the Wannier functions, this model of strongly correlated localized electrons makes clear predictions whether or not the system can possess superconducting or magnetic eigenstates.

categories: cond-mat.str-el cond-mat.mtrl-sci cond-mat.supr-con math-ph math.MP


title: On the blow up criterion of 3D-NSE in Sobolev-Gevrey spaces

abstract: In \cite{JB1}, Benameur proved a blow-up result of the non regular solution of $(NSE)$ in the Sobolev-Gevrey spaces. In this paper we improve this result, precisely we give an exponential type explosion in Sobolev-Gevrey spaces with less regularity on the initial condition. Fourier analysis is used.


title: Wavelet-based cross-correlation analysis of structure scaling in turbulent clouds

abstract: We propose a statistical tool to compare the scaling behaviour of turbulence in pairs of molecular cloud maps. Using artificial maps with well defined spatial properties, we calibrate the method and test its limitations to ultimately apply it to a set of observed maps. We develop the wavelet-based weighted cross-correlation (WWCC) method to study the relative contribution of structures of different sizes and their degree of correlation in two maps as a function of spatial scale, and the mutual displacement of structures in the molecular cloud maps. We test the WWCC for circular structures having a single prominent scale and fractal structures showing a self-similar behavior without prominent scales. Observational noise and a finite map size limit the scales where the cross-correlation coefficients and displacement vectors can be reliably measured. For fractal maps containing many structures on all scales, the limitation from the observational noise is negligible for signal-to-noise ratios >5. (abbrev). Application of the WWCC to the observed line maps of the giant molecular cloud G333 allows to add specific scale information to the results obtained earlier using the principle component analysis. It confirms the chemical and excitation similarity of $^{13}$CO and C$^{18}$O on all scales, but shows a deviation of HCN at scales of up to 7' (~7 pc). This can be interpreted as a chemical transition scale. The largest structures also show a systematic offset along the filament, probably due to a large-scale density gradient. The WWCC can compare correlated structures in different maps of molecular clouds identifying scales that represent structural changes such as chemical and phase transitions and prominent or enhanced dimensions.


title: Long time decay for 3D-NSE in Gevrey-Sobolev spaces

abstract: In this paper we prove, if $u$ is a global solution to Navier-Stokes equations in the Sobolev-Gevrey spaces $H^1_{a,\sigma}(\mathbb R^3)$, then $\|u(t)\|_{H^1_{a,\sigma}}$ decays to zero as time goes to infinity. Fourier analysis is used.


title: Long time decay of 3D-NSE in Lei-Lin-Gevrey spaces

abstract: In this paper, we prove that there exists a unique global solution of $3D$ Navier-Stokes equation if $\exp(a|D|^{1/\sigma})u^0\in{\mathcal{X}}^{-1}(\mathbb R^3)$ and $\|u^0\|_{{\mathcal{X}}^{-1}}<\nu$. Moreover, we will show that $\|\exp(a|D|^{1/\sigma}) u(t)\|_{{\mathcal{X}}^{-1}}$ goes to zero if the time $t$ goes to infinity.


title: Effective destruction of CO by cosmic rays: implications for tracing H$_2$ gas in the Universe

abstract: We report on the effects of cosmic rays (CRs) on the abundance of CO in $\rm H_2$ clouds under conditions typical for star-forming galaxies in the Universe. We discover that this most important molecule for tracing H$_2$ gas is very effectively destroyed in ISM environments with CR energy densities $\rm U_{CR}\sim(50-10^{3})\times U_{CR,Gal}$, a range expected in numerous star-forming systems throughout the Universe. This density-dependent effect operates volumetrically rather than only on molecular cloud surfaces (i.e. unlike FUV radiation that also destroys CO), and is facilitated by: a) the direct destruction of CO by CRs, and b) a reaction channel activated by CR-produced He$^{+}$. The effect we uncover is strong enough to render Milky-Way type Giant Molecular Clouds (GMCs) very CO-poor (and thus CO-untraceable), even in ISM environments with rather modestly enhanced average CR energy densities of $\rm U_{CR}\sim(10-50)\times\rm U_{CR,Gal}$. We conclude that the CR-induced destruction of CO in molecular clouds, unhindered by dust absorption, is perhaps the single most important factor controlling the CO-visibility of molecular gas in vigorously star-forming galaxies. We anticipate that a second order effect of this CO destruction mechanism will be to make the H$_2$ distribution in the gas-rich disks of such galaxies appear much clumpier in CO $J$=1--0, 2--1 line emission than it actually is. Finally we give an analytical approximation of the CO/H$_2$ abundance ratio as a function of gas density and CR energy density for use in galaxy-size or cosmological hydrodynamical simulations, and propose some key observational tests.


title: New Barr-Zee contributions to $\mathbf{(g-2)_\mu}$ in two-Higgs-doublet models

abstract: We study the contribution of new sets of two-loop Barr-Zee type diagrams to the anomalous magnetic moment of the muon within the two-Higgs-doublet model framework. We show that some of these contributions can be quite sizeable for a large region of the parameter space and can significantly reduce, and in some cases even explain, the discrepancy between the theoretical prediction and the experimentally measured value of this observable. Analytical expressions are given for all the calculations performed in this work.


title: Gaps in the Milnor-Moore spectral sequence and the Hilali conjecture

abstract: In his study of Halperin's toral-rank conjecture, M. R. Hilali conjectured that for any simply connected rationally elliptic space $X$, one must have $dim\pi_*(X)\otimes \mathbb{Q} \leq dimH^*(X,\mathbb{Q})$. Let $(\Lambda V, d)$ denote a Sullivan minimal model of $X$ and $d_k$ the first non-zero homogeneous part of the differential $d$. In this paper, we use spectral sequence arguments to prove that if $(\Lambda V, d_k)$ is elliptic, then, there is no gaps in the $E_{\infty}$ term of the Milnor-Moore spectral sequence of $X$. Consequently, we confirm the Hilali conjecture when $V = V^{odd}$ or else when $k\geq 3$ and $(\Lambda V, d_k)$ is elliptic.


title: The function $(\cosh\sqrt{at^2+b})$ is exponentially convex

abstract: Given positive numbers $a$ and $b$, the function $\sqrt{at^2+b}$ is exponentially convex function of $t$ on the whole real axis. Three proofs of this result are presented.


title: A fast Mixed Model B-splines algorithm

abstract: A fast algorithm for B-splines in mixed models is presented. B-splines have local support and are computational attractive, because the corresponding matrices are sparse. A key element of the new algorithm is that the local character of B-splines is preserved, while in other existing methods this local character is lost. The computation time for the fast algorithm is linear in the number of B-splines, while computation time scales cubically for existing transformations.


title: Transmuted Geometric Distribution and its Prpoerties

abstract: Transmuted geometric distribution with two parameters and is proposed as a new generalization of the geometric distribution by employing the quadratic transmutation techniques of Shaw and Buckley (2007). Its important distributional and reliability properties are investigated. Parameter estimation methods is discussed.


title: Gray-Level Image Transitions Driven by Tsallis Entropic Index

abstract: The maximum entropy principle is largely used in thresholding and segmentation of images. Among the several formulations of this principle, the most effectively applied is that based on Tsallis non-extensive entropy. Here, we discuss the role of its entropic index in determining the thresholds. When this index is spanning the interval (0,1), for some images, the values of thresholds can have large leaps. In this manner, we observe abrupt transitions in the appearance of corresponding bi-level or multi-level images. These gray-level image transitions are analogous to order or texture transitions observed in physical systems, transitions which are driven by the temperature or by other physical quantities.


title: Leader-follower Tracking Control with Guaranteed Consensus Performance for Interconnected Systems with Linear Dynamic Uncertain Coupling

abstract: This paper considers the leader-follower tracking control problem for linear interconnected systems with undirected topology and linear dynamic coupling. Interactions between the systems are treated as linear dynamic uncertainty and are described in terms of integral quadratic constraints (IQCs). A consensus-type tracking control protocol is proposed for each system based on its state relative its neighbors. In addition a selected set of subsystems uses for control their relative states with respect to the leader. Two methods are proposed for the design of this control protocol. One method uses a coordinate transformation to recast the protocol design problem as a decentralized robust control problem for an auxiliary interconnected large scale system. Another method is direct, it does not employ coordinate transformation; it also allows for more general linear uncertain interactions. Using these methods, sufficient conditions are obtained which guarantee that the system tracks the leader. These conditions guarantee a suboptimal bound on the system consensus and tracking performance. The proposed methods are compared using a simulation example, and their effectiveness is discussed. Also, algorithms are proposed for computing suboptimal controllers.


title: Combining probability distributions: Extending the logarithmic pooling approach

abstract: Combining distributions is an important issue in decision theory and Bayesian inference. Logarithmic pooling is a popular method to aggregate expert opinions by using a set of weights that reflect the reliability of each information source. However, the resulting pooled distribution depends heavily on set of weights given to each opinion/prior and thus careful consideration must be given to the choice of weights. In this paper we review and extend the statistical theory of logarithmic pooling, focusing on the assignment of the weights using a hierarchical prior distribution. We explore several statistical applications, such as the estimation of survival probabilities, meta-analysis and Bayesian melding of deterministic models of population growth and epidemics. We show that it is possible learn the weights from data, although identifiability issues may arise for some configurations of priors and data. Furthermore, we show how the hierarchical approach leads to posterior distributions that are able to accommodate prior-data conflict in complex models.


title: Largest Eigenvalue of the Laplacian Matrix

abstract: We study the eigenspace of the Laplacian matrix of a simple graph corresponding to the largest eigenvalue, subsequently arriving at the theory of modular decomposition of T. Gallai.


title: Anomalous charge pumping in a one-dimensional optical superlattice

abstract: We model atomic motion in a sliding superlattice potential to explore topological "charge pumping" and to find optimal parameters for experimental observation of this phenomenon. We analytically study the band-structure, finding how the Wannier states evolve as two sinusoidal lattices are moved relative to one-another, and relate this evolution to the center of mass motion of an atomic cloud. We pay particular attention to counterintuitive or anomalous regimes, such as when the atomic motion is opposite to that of the lattice.


title: Integer points and their orthogonal lattices

abstract: Linnik proved in the late 1950's the equidistribution of integer points on large spheres under a congruence condition. The congruence condition was lifted in 1988 by Duke (building on a break-through by Iwaniec) using completely different techniques. We conjecture that this equidistribution result also extends to the pairs consisting of a vector on the sphere and the shape of the lattice in its orthogonal complement. We use a joining result for higher rank diagonalizable actions to obtain this conjecture under an additional congruence condition.


title: Grassmannian Codes as Lifts of Matrix Codes Derived as Images of Linear Block Codes over Finite Fields

abstract: Let $p$ be a prime such that $p \equiv 2$ or $3$ mod $5$. Linear block codes over the non-commutative matrix ring of $2 \times 2$ matrices over the prime field $GF(p)$ endowed with the Bachoc weight are derived as isometric images of linear block codes over the Galois field $GF(p^2)$ endowed with the Hamming metric. When seen as rank metric codes, this family of matrix codes satisfies the Singleton bound and thus are maximum rank distance codes, which are then lifted to form a special class of subspace codes, the Grassmannian codes, that meet the anticode bound. These so-called anticode-optimal Grassmannian codes are associated in some way with complete graphs. New examples of these maximum rank distance codes and anticode-optimal Grassmannian codes are given.


title: Parameter Estimation of Ground Moving Targets Based on SKT-DLVT Processing

abstract: It is well known that the motion of a ground moving target may induce the range cell migration, spectrum spread and velocity ambiguity during the imaging time, which makes the image smeared. To eliminate the influence of these factors on image focusing, a novel method for parameter estimation of ground moving targets, known as SKT-DLVT, is proposed in this paper. In this method, the segmental keystone transform (SKT) is used to correct the range walk of targets simultaneously, and a new transform, namely, Doppler Lv's transform (LVT) is applied on the azimuth signal to estimate the parameters. Theoretical analysis confirms that no interpolation is needed for the proposed method and the targets can be well focused within limited searching range of the ambiguity number. The proposed method is capable of obtaining the accurate parameter estimates efficiently in the low signal-to-noise ratio (SNR) scenario with low computational burden and memory cost, making it suitable to be applied in memory-limited and real-time processing systems. The effectiveness of the proposed method is demonstrated by both simulated and real data.


title: Singular perturbation by bending for an adhesive obstacle problem

abstract: A free boundary problem arising from materials science is studied in one-dimensional case. The problem studied here is an obstacle problem for the non-convex energy consisting of a bending energy, tension and an adhesion energy. If the bending energy, which is a higher order term, is deleted then "edge" singularities of the solutions (surfaces) may occur at the free boundary as Alt-Caffarelli type variational problems. The main result of this paper is to give a singular limit of the energy utilizing the notion of $\Gamma$-convergence, when the bending energy can be regarded as a perturbation. This singular limit energy only depends on the state of surfaces at the free boundary as seen in singular perturbations for phase transition models.


title: Index of refraction, Rayleigh scattering length, and Sellmeier coefficients in solid and liquid argon and xenon

abstract: Large liquid argon detectors have become widely used in low rate experiments, including dark matter and neutrino research. However, the optical properties of liquid argon are not well understood at the large scales relevant for current and near-future detectors.The index of refraction of liquid argon at the scin- tillation wavelength has not been measured, and current Rayleigh scattering length calculations disagree with measurements. Furthermore, the Rayleigh scattering length and index of refraction of solid argon and solid xenon at their scintillation wavelengths have not been previously measured or calculated. We introduce a new calculation using existing data in liquid and solid argon and xenon to extrapolate the optical properties at the scintillation wavelengths using the Sellmeier dispersion relationship.


title: A phenomenological analysis of azimuthal asymmetries in unpolarized semi-inclusive deep inelastic scattering

abstract: We present a phenomenological analysis of the cos-phi and cos-2phi asymmetries in unpolarized semi-inclusive deep inelastic scattering, based on the recent multidimensional data released by the COMPASS and HERMES Collaborations. In the TMD framework, valid at relatively low transverse momenta, these asymmetries arise from intrinsic transverse momentum and transverse spin effects, and from their correlations. The role of the Cahn and Boer-Mulders effects in both azimuthal moments is explored up to order 1/Q. As the kinematics of the present experiments is dominated by the low-Q^2 region, higher-twist contributions turn out to be important, affecting the results of our fits.


title: Simple loops on 2-bridge spheres in Heckoid orbifolds for the trivial knot

abstract: In this paper, we give a necessary and sufficient condition for an essential simple loop on a $2$-bridge sphere in an even Heckoid orbifold for the trivial knot to be null-homotopic, peripheral or torsion in the orbifold. We also give a necessary and sufficient condition for two essential simple loops on a $2$-bridge sphere in an even Heckoid orbifold for the trivial knot to be homotopic in the orbifold.


title: Algebraic and geometric aspects of rational $\Gamma$-inner functions

abstract: The set \[ \Gamma {\stackrel{\rm def}{=}} \{(z+w,zw):|z|\leq 1,|w|\leq 1\} \subset {\mathbb{C}}^2 \] has intriguing complex-geometric properties; it has a 3-parameter group of automorphisms, its distinguished boundary is a ruled surface homeomorphic to the M\"obius band and it has a special subvariety which is the only complex geodesic that is invariant under all automorphisms. We exploit this geometry to develop an explicit and detailed structure theory for the rational maps from the unit disc to $\Gamma$ that map the boundary of the disc to the distinguished boundary of $\Gamma$.


title: Nonconforming finite element method applied to the driven cavity problem

abstract: A cheapest stable nonconforming finite element method is presented for solving the incompressible flow in a square cavity without smoothing the corner singularities. The stable cheapest nonconforming finite element pair based on $P_1\times \widetilde{{\mathcal P}_{c}^h}$ on rectangular meshes \cite{stab-cheapest} is employed with a minimal modification for the discontinuous Dirichlet data on the top boundary, where $\widetilde{{\mathcal P}_{c}^h}$ is the finite element space of piecewise constant pressures with the globally one-dimensional checker-board pattern subspace eliminated. The proposed Stokes elements have least number of degrees of freedom compared to those of known stable Stokes elements. Three accuracy indications for our elements are analyzed and numerically verified. Also, various numerous computational results obtained by using our proposed element show excellent accuracy.


title: G\"artner-Ellis condition for squared asymptotically stationary Gaussian processes

abstract: The G\"artner-Ellis condition for the square of an asymptotically stationary Gaussian process is established. The same limit holds for the conditional distri-bution given any fixed initial point, which entails weak multiplicative ergodicity. The limit is shown to be the Laplace transform of a convolution of Gamma distributions with Poisson compound of exponentials. A proof based on Wiener-Hopf factorization induces a probabilistic interpretation of the limit in terms of a regression problem.


title: Canonical traces and directly finite Leavitt path algebras

abstract: Motivated by the study of traces on graph $C^*$-algebras, we consider traces (additive, central maps) on Leavitt path algebras, the algebraic counterparts of graph $C^*$-algebras. In particular, we consider traces which vanish on nonzero graded components of a Leavitt path algebra and refer to them as {\em canonical} since they are uniquely determined by their values on the vertices. A desirable property of a $\mathbb{C}$-valued trace on a $C^*$-algebra is that the trace of an element of the positive cone is nonnegative. We adapt this property to traces on a Leavitt path algebra $L_K(E)$ with values in any involutive ring. We refer to traces with this property as positive. If a positive trace is injective on positive elements, we say that it is faithful. We characterize when a canonical, $K$-linear trace is positive and when it is faithful in terms of its values on the vertices. As a consequence, we obtain a bijective correspondence between the set of faithful, gauge invariant, $\mathbb{C}$-valued (algebra) traces on $L_{\mathbb{C}}(E)$ of a countable graph $E$ and the set of faithful, semifinite, lower semicontinuous, gauge invariant (operator theory) traces on the corresponding graph $C^*$-algebra $C^*(E)$. With the direct finite condition (i.e $xy=1$ implies $yx=1$) for unital rings adapted to rings with local units, we characterize directly finite Leavitt path algebras as exactly those having the underlying graphs in which no cycle has an exit. Our proof involves consideration of "local" Cohn-Leavitt subalgebras of finite subgraphs. Lastly, we show that, while related, the class of locally noetherian, the class of directly finite, and the class of Leavitt path algebras which admit a faithful trace are different in general.


title: Exploring Hierarchies in Online Social Networks

abstract: Social hierarchy (i.e., pyramid structure of societies) is a fundamental concept in sociology and social network analysis. The importance of social hierarchy in a social network is that the topological structure of the social hierarchy is essential in both shaping the nature of social interactions between individuals and unfolding the structure of the social networks. The social hierarchy found in a social network can be utilized to improve the accuracy of link prediction, provide better query results, rank web pages, and study information flow and spread in complex networks. In this paper, we model a social network as a directed graph G, and consider the social hierarchy as DAG (directed acyclic graph) of G, denoted as GD. By DAG, all the vertices in G can be partitioned into different levels, the vertices at the same level represent a disjoint group in the social hierarchy, and all the edges in DAG follow one direction. The main issue we study in this paper is how to find DAG GD in G. The approach we take is to find GD by removing all possible cycles from G such that G = U(G) + GD where U(G) is a maximum Eulerian subgraph which contains all possible cycles. We give the reasons for doing so, investigate the properties of GD found, and discuss the applications. In addition, we develop a novel two-phase algorithm, called Greedy-&-Refine, which greedily computes an Eulerian subgraph and then refines this greedy solution to find the maximum Eulerian subgraph. We give a bound between the greedy solution and the optimal. The quality of our greedy approach is high. We conduct comprehensive experimental studies over 14 real-world datasets. The results show that our algorithms are at least two orders of magnitude faster than the baseline algorithm.


title: A Row-parallel 8$\times$8 2-D DCT Architecture Using Algebraic Integer Based Exact Computation

abstract: An algebraic integer (AI) based time-multiplexed row-parallel architecture and two final-reconstruction step (FRS) algorithms are proposed for the implementation of bivariate AI-encoded 2-D discrete cosine transform (DCT). The architecture directly realizes an error-free 2-D DCT without using FRSs between row-column transforms, leading to an 8$\times$8 2-D DCT which is entirely free of quantization errors in AI basis. As a result, the user-selectable accuracy for each of the coefficients in the FRS facilitates each of the 64 coefficients to have its precision set independently of others, avoiding the leakage of quantization noise between channels as is the case for published DCT designs. The proposed FRS uses two approaches based on (i) optimized Dempster-Macleod multipliers and (ii) expansion factor scaling. This architecture enables low-noise high-dynamic range applications in digital video processing that requires full control of the finite-precision computation of the 2-D DCT. The proposed architectures and FRS techniques are experimentally verified and validated using hardware implementations that are physically realized and verified on FPGA chip. Six designs, for 4- and 8-bit input word sizes, using the two proposed FRS schemes, have been designed, simulated, physically implemented and measured. The maximum clock rate and block-rate achieved among 8-bit input designs are 307.787 MHz and 38.47 MHz, respectively, implying a pixel rate of 8$\times$307.787$\approx$2.462 GHz if eventually embedded in a real-time video-processing system. The equivalent frame rate is about 1187.35 Hz for the image size of 1920$\times$1080. All implementations are functional on a Xilinx Virtex-6 XC6VLX240T FPGA device.


title: Energy dissipation in single-domain ferromagnetic nanoparticles: Dynamical approach

abstract: We study, both analytically and numerically, the phenomenon of energy dissipation in single-domain ferromagnetic nanoparticles driven by an alternating magnetic field. Our interest is focused on the power loss resulting from the Landau-Lifshitz-Gilbert equation, which describes the precessional motion of the nanoparticle magnetic moment. We determine the power loss as a function of the field amplitude and frequency and analyze its dependence on different regimes of forced precession induced by circularly and linearly polarized magnetic fields. The conditions to maximize the nanoparticle heating are also analyzed.


title: Theoretical Study of Phosphorene Tunneling Field Effect Transistors

abstract: In this work, device performances of tunneling field effect transistors (TFETs) based on phosphorene are explored via self-consistent atomistic quantum transport simulations. Phosphorene is an ultra-thin two-dimensional (2-D) material with a direct band gap suitable for TFETs applications. Our simulation shows that phosphorene TFETs exhibit subthreshold slope (SS) below 60 mV/dec and a wide range of on-current depending on the transport direction due to highly anisotropic band structures of phosphorene. By benchmarking with monolayer MoTe2 TFETs, we predict that phosphorene TFETs oriented in the small effective mass direction can yield much larger on-current at the same on-current/off-current ratio than monolayer MoTe2 TFETs. It is also observed that a gate underlap structure is required for scaling down phosphorene TFETs in the small effective mass direction to suppress the source-to-drain direct tunneling leakage current.


title: Discretising the velocity distribution for directional dark matter experiments

abstract: Dark matter (DM) direct detection experiments which are directionally-sensitive may be the only method of probing the full velocity distribution function (VDF) of the Galactic DM halo. We present an angular basis for the DM VDF which can be used to parametrise the distribution in order to mitigate astrophysical uncertainties in future directional experiments and extract information about the DM halo. This basis consists of discretising the VDF in a series of angular bins, with the VDF being only a function of the DM speed $v$ within each bin. In contrast to other methods, such as spherical harmonic expansions, the use of this basis allows us to guarantee that the resulting VDF is everywhere positive and therefore physical. We present a recipe for calculating the event rates corresponding to the discrete VDF for an arbitrary number of angular bins $N$ and investigate the discretisation error which is introduced in this way. For smooth, Standard Halo Model-like distribution functions, only $N=3$ angular bins are required to achieve an accuracy of around $10-30\%$ in the number of events in each bin. Shortly after confirmation of the DM origin of the signal with around 50 events, this accuracy should be sufficient to allow the discretised velocity distribution to be employed reliably. For more extreme VDFs (such as streams), the discretisation error is typically much larger, but can be improved with increasing $N$. This method paves the way towards an astrophysics-independent analysis framework for the directional detection of dark matter.


title: On quantification of systemic redundancy in reliable systems

abstract: In this paper, we consider the problem of quantifying systemic redundancy in reliable systems having multiple controllers with overlapping functionality. In particular, we consider a multi-channel system with multi-controller configurations -- where controllers are required to respond optimally, in the sense of best-response correspondence, a {\it reliable-by-design} requirement, to non-faulty controllers so as to ensure or maintain some system properties. Here we introduce a mathematical framework, based on the notion of relative entropy of probability measures associated with steady-state solutions of Fokker-Planck equations for a family of stochastically perturbed multi-channel systems, that provides useful information towards a systemic assessment of redundancy in the system.


title: Width Hierarchy for k-OBDD of Small Width

abstract: In this paper was explored well known model k-OBDD. There are proven width based hierarchy of classes of boolean functions which computed by k-OBDD. The proof of hierarchy is based on sufficient condition of Boolean function's non representation as k-OBDD and complexity properties of Boolean function SAF. This function is modification of known Pointer Jumping (PJ) and Indirect Storage Access (ISA) functions.


title: Remnant quantum resources of collapsed macroscopic quantum superpositions

abstract: We consider the collapse of a macroscopic quantum superposition occurring due to the measurement which optimally distinguishes its branches. Given a macroscopic superposition of N spin-1/2 particles, we use such a Helstrom measurement to construct the local unitary operator which maximizes the usefulness of the superposition for Heisenberg-limited phase estimation (i.e., with quantum Cram\'{e}r Rao bound proportional to 1/N). In contrast, the collapsed state is not useful as a probe for phase estimation below the standard quantum limit. For the case N=2, we compute the entanglement entropy of the collapsed state and show that it is reduced below that of the initial superposition when the superposition is macroscopic. We consider the remnant quantum resources of collapsed hierarchical Schr\"{o}dinger cat states to show that collapsed macroscopic superpositions can still be useful for ultraprecise quantum metrology.


title: Bridge numbers of knots in the page of an open book

abstract: Given any closed, connected, orientable $3$--manifold and integers $g\geq g(M), D > 0$, we show the existence of knots in $M$ whose genus $g$ bridge number is greater than $D$. These knots lie in a page of an open book decomposition of $M$, and the proof proceeds by examining the action of the map induced by the monodromy on the arc and curve complex of a page. A corollary is that there are Berge knots of arbitrarily large genus one bridge number.


title: The Possibility of Kelvin-Helmholtz Instability in Solar Spicules

abstract: Transversal oscillations of spicules axes may be related to the propagation of magnetohydrodynamic waves along them. These waves may become unstable and the instability can be of the Kelvin-Helmholtz type. We use the dispersion relation of kink mode derived from linearized magnetohydrodynamic equations. The input parameters of the derived dispersion equation, namely, spicules and their ambient medium densities ratios as well as their corresponding magnetic fields ratios, are considered to be within the range $0-1$. By solving the dispersion equation numerically, we show that for higher densities and lower magnetic fields ratios within the range mentioned, the KHI onset in type II spicules conditions is possible. This possibility decreases with an increase in Alfven velocity inside spicules. A rough criterion for appearing of Kelvin-Helmholtz instability is obtained. We also drive a more reliable and exact criterion for KHI onset of kink waves.


title: From the Hartree dynamics to the Vlasov equation

abstract: We consider the evolution of quasi-free states describing $N$ fermions in the mean field limit, as governed by the nonlinear Hartree equation. In the limit of large $N$, we study the convergence towards the classical Vlasov equation. For a class of regular interaction potentials, we establish precise bounds on the rate of convergence.


title: A two-dimensional continued fractions algorithm with Lagrange and Dirichlet properties

abstract: A Lagrange Theorem in dimension 2 is proved, for a particular two-dimensional algorithm, with a very natural geometrical definition. Dirichlet-type properties for the convergence of the algorithm are also proved. These properties procced from a geometrical quality of the algorithm. Some refercences are given to the works of various authors, in the domain of multidimensional continued fractions algorithms.


title: Sketch-based Shape Retrieval using Pyramid-of-Parts

abstract: We present a multi-scale approach to sketch-based shape retrieval. It is based on a novel multi-scale shape descriptor called Pyramidof- Parts, which encodes the features and spatial relationship of the semantic parts of query sketches. The same descriptor can also be used to represent 2D projected views of 3D shapes, allowing effective matching of query sketches with 3D shapes across multiple scales. Experimental results show that the proposed method outperforms the state-of-the-art method, whether the sketch segmentation information is obtained manually or automatically by considering each stroke as a semantic part.


title: Stability and instability of the Einstein-Lichnerowicz constraint system

abstract: We investigate the relevance of the conformal method by investigating stability issues for the Einstein-Lichnerowicz conformal constraint system in a nonlinear scalar-field setting. We prove the stability of the system with respect to arbitrary perturbations of generic focusing physics data on closed locally conformally flat manifolds, in any dimension. We also show that our stability result is sharp by constructing explicit instability examples when its assumptions are not satisfied. Our results apply to a more general class of constraint-like systems.


title: Barkhausen noise in the Random Field Ising Magnet Nd$_2$Fe$_{14}$B

abstract: With sintered needles aligned and a magnetic field applied transverse to its easy axis, the rare-earth ferromagnet Nd$_2$Fe$_{14}$B becomes a room-temperature realization of the Random Field Ising Model. The transverse field tunes the pinning potential of the magnetic domains in a continuous fashion. We study the magnetic domain reversal and avalanche dynamics between liquid helium and room temperatures at a series of transverse fields using a Barkhausen noise technique. The avalanche size and energy distributions follow power-law behavior with a cutoff dependent on the pinning strength dialed in by the transverse field, consistent with theoretical predictions for Barkhausen avalanches in disordered materials. A scaling analysis reveals two regimes of behavior: one at low temperature and high transverse field, where the dynamics are governed by the randomness, and the second at high temperature and low transverse field where thermal fluctuations dominate the dynamics.


title: Dynamics of the Cosmological Apparent Horizon: Surface Gravity & Temperature

abstract: In the context of thermodynamics applied to our cosmological apparent horizon, we explicit in greater details our previous work which established the Friedmann Equations from projection of Hayward's Unified First Law. In particular, we show that the dynamical Hayward-Kodama surface gravity is perfectly well-defined and is suitable for this derivation. We then relate this surface gravity to a physical notion of temperature, and show this has constant, positive sign for any kind of past-inner trapping horizons. Hopefully this will clarify the choice of temperature in a dynamical Friedmann-Lema\^itre-Roberston-Walker spacetime.


title: Impacts of Bad Data on the PMU based Line Outage Detection

abstract: Power systems become more prone to cyber-attacks due to the high integration of information technologies. In this paper, we demonstrate that the outages of some lines can be masked by injecting false data into a set of measurements. The success of the topology attack can be guaranteed by making that: 1)the injected false data obeys KCL and KVL to avoid being detected by the bad data detection program in the state estimation; 2)the residual is increased such that the line outage cannot be detected by PMU data. A quadratic programming problem is set up to determine the optimal attack vector that can maximize the residual of the outaged line. The IEEE 39-bus system is used to demonstrate the masking scheme.


title: Are Discoveries Spurious? Distributions of Maximum Spurious Correlations and Their Applications

abstract: Over the last two decades, many exciting variable selection methods have been developed for finding a small group of covariates that are associated with the response from a large pool. Can the discoveries from these data mining approaches be spurious due to high dimensionality and limited sample size? Can our fundamental assumptions about the exogeneity of the covariates needed for such variable selection be validated with the data? To answer these questions, we need to derive the distributions of the maximum spurious correlations given a certain number of predictors, namely, the distribution of the correlation of a response variable $Y$ with the best $s$ linear combinations of $p$ covariates $\mathbf{X}$, even when $\mathbf{X}$ and $Y$ are independent. When the covariance matrix of $\mathbf{X}$ possesses the restricted eigenvalue property, we derive such distributions for both a finite $s$ and a diverging $s$, using Gaussian approximation and empirical process techniques. However, such a distribution depends on the unknown covariance matrix of $\mathbf{X}$. Hence, we use the multiplier bootstrap procedure to approximate the unknown distributions and establish the consistency of such a simple bootstrap approach. The results are further extended to the situation where the residuals are from regularized fits. Our approach is then used to construct the upper confidence limit for the maximum spurious correlation and to test the exogeneity of the covariates. The former provides a baseline for guarding against false discoveries and the latter tests whether our fundamental assumptions for high-dimensional model selection are statistically valid. Our techniques and results are illustrated with both numerical examples and real data analysis.


title: Sharp thresholds for Gibbs-non-Gibbs transition in the fuzzy Potts model with a Kac-type interaction

abstract: We investigate the Gibbs properties of the fuzzy Potts model on the d-dimensional torus with Kac interaction. We use a variational approach for profiles inspired by that of Fernandez, den Hollander and Mart{\i}nez for their study of the Gibbs-non-Gibbs transitions of a dynamical Kac-Ising model on the torus. As our main result, we show that the mean-field thresholds dividing Gibbsian from non-Gibbsian behavior are sharp in the fuzzy Kac-Potts model with class size unequal two. On the way to this result we prove a large deviation principle for color profiles with diluted total mass densities and use monotocity arguments.


title: Maximum Likelihood de novo reconstruction of viral populations using paired end sequencing data

abstract: We present MLEHaplo, a maximum likelihood de novo assembly algorithm for reconstructing viral haplotypes in a virus population from paired-end next generation sequencing (NGS) data. Using the pairing information of reads in our proposed Viral Path Reconstruction Algorithm (ViPRA), we generate a small subset of paths from a De Bruijn graph of reads that serve as candidate paths for true viral haplotypes. Our proposed method MLEHaplo then generates a maximum likelihood estimate of the viral population using the paths reconstructed by ViPRA. We evaluate and compare MLEHaplo on simulated datasets of 1200 base pairs at different sequence coverage, on HCV strains with sequencing errors, and on a lab mixture of five HIV-1 strains. MLEHaplo reconstructs full length viral haplotypes having a 100% sequence identity to the true viral haplotypes in most of the small genome simulated viral populations at 250x sequencing coverage. While reference based methods either under-estimate or over-estimate the viral haplotypes, MLEHaplo limits the over-estimation to 3 times the size of true viral haplotypes, reconstructs the full phylogeny in the HCV to greater than 99% sequencing identity and captures more sequencing variation for the HIV-1 strains dataset compared to their known consensus sequences.


title: Scheduling of unit-length jobs with cubic incompatibility graphs on three uniform machines

abstract: In the paper we consider the problem of scheduling $n$ identical jobs on 3 uniform machines with speeds $s_1, s_2,$ and $s_3$ to minimize the schedule length. We assume that jobs are subjected to some kind of mutual exclusion constraints, modeled by a cubic incompatibility graph. We show that if the graph is 2-chromatic then the problem can be solved in $O(n^2)$ time. If the graph is 3-chromatic, the problem becomes NP-hard even if $s_1>s_2=s_3$. However, in this case there exists a $4/3$-approximation algorithm running in $O(n^3)$ time. Moreover, this algorithm solves the problem almost surely to optimality if $3s_1/4 \leq s_2 = s_3$.


title: PhaseLift is robust to a constant fraction of arbitrary errors

abstract: Consider the task of recovering an unknown $n$-vector from phaseless linear measurements. This task is the phase retrieval problem. Through the technique of lifting, this nonconvex problem may be convexified into a semidefinite rank-one matrix recovery problem, known as PhaseLift. Under a linear number of exact Gaussian measurements, PhaseLift recovers the unknown vector exactly with high probability. Under noisy measurements, the solution to a variant of PhaseLift has error proportional to the $\ell_1$ norm of the noise. In the present paper, we study the robustness of this variant of PhaseLift to a case with noise and gross, arbitrary corruptions. We prove that PhaseLift can tolerate a small, fixed fraction of gross errors, even in the highly underdetermined regime where there are only $O(n)$ measurements. The lifted phase retrieval problem can be viewed as a rank-one robust Principal Component Analysis (PCA) problem under generic rank-one measurements. From this perspective, the proposed convex program is simpler that the semidefinite version of the sparse-plus-low-rank formulation standard in the robust PCA literature. Specifically, the rank penalization through a trace term is unnecessary, and the resulting optimization program has no parameters that need to be chosen. The present work also achieves the information theoretically optimal scaling of $O(n)$ measurements without the additional logarithmic factors that appear in existing general robust PCA results.


title: Classification of catalytic branching processes and structure of the criticality set

abstract: We study a catalytic branching process (CBP) with any finite set of catalysts. This model describes a system of particles where the movement is governed by a Markov chain with arbitrary finite or countable state space and the branching may only occur at the points of catalysis. The results obtained generalize and strengthen those known in cases of CBP with a single catalyst and of branching random walk on d-dimensional integer lattice with a finite number of sources of particles generation. We propose to classify CBP with N catalysts as supercritical, critical or subcritical according to the value of the Perron root of a specified NxN matrix. Such classification agrees with the moment analysis performed here for local and total particles numbers. By introducing the criticality set C we also consider the influence of catalysts parameters on the process behavior. The proof is based on construction of auxiliary multi-type Bellman-Harris processes with the help of hitting times under taboo and on application of multidimensional renewal theorems. Keywords and phrases: catalytic branching process, classification, hitting times under taboo, moment analysis, multi-type Bellman-Harris process.


title: Bayesian Inference of Arrival Rate and Substitution Behavior from Sales Transaction Data with Stockouts

abstract: When an item goes out of stock, sales transaction data no longer reflect the original customer demand, since some customers leave with no purchase while others substitute alternative products for the one that was out of stock. Here we develop a Bayesian hierarchical model for inferring the underlying customer arrival rate and choice model from sales transaction data and the corresponding stock levels. The model uses a nonhomogeneous Poisson process to allow the arrival rate to vary throughout the day, and allows for a variety of choice models. Model parameters are inferred using a stochastic gradient MCMC algorithm that can scale to large transaction databases. We fit the model to data from a local bakery and show that it is able to make accurate out-of-sample predictions, and to provide actionable insight into lost cookie sales.


title: Optimal cyclic codes with generalized Niho type zeroes and the weight distribution

abstract: In this paper we extend the works \cite{gegeng2,XLZD} further in two directions and compute the weight distribution of these cyclic codes under more relaxed conditions. It is interesting to note that many cyclic codes in the family are optimal and have only a few non-zero weights. Besides using similar ideas from \cite{gegeng2,XLZD}, we carry out some subtle manipulation of certain exponential sums.


title: Using and Designing Platforms for In Vivo Education Experiments

abstract: In contrast to typical laboratory experiments, the everyday use of online educational resources by large populations and the prevalence of software infrastructure for A/B testing leads us to consider how platforms can embed in vivo experiments that do not merely support research, but ensure practical improvements to their educational components. Examples are presented of randomized experimental comparisons conducted by subsets of the authors in three widely used online educational platforms Khan Academy, edX, and ASSISTments. We suggest design principles for platform technology to support randomized experiments that lead to practical improvements enabling Iterative Improvement and Collaborative Work and explain the benefit of their implementation by WPI co-authors in the ASSISTments platform.


title: Stable Leader Election in Population Protocols Requires Linear Time

abstract: A population protocol *stably elects a leader* if, for all $n$, starting from an initial configuration with $n$ agents each in an identical state, with probability 1 it reaches a configuration $\mathbf{y}$ that is correct (exactly one agent is in a special leader state $\ell$) and stable (every configuration reachable from $\mathbf{y}$ also has a single agent in state $\ell$). We show that any population protocol that stably elects a leader requires $\Omega(n)$ expected "parallel time" --- $\Omega(n^2)$ expected total pairwise interactions --- to reach such a stable configuration. Our result also informs the understanding of the time complexity of chemical self-organization by showing an essential difficulty in generating exact quantities of molecular species quickly.


title: Supporting Instructors in Collaborating with Researchers using MOOClets

abstract: Most education and workplace learning takes place in classroom contexts far removed from laboratories or field sites with special arrangements for scientific research. But digital online resources provide a novel opportunity for large scale efforts to bridge the real world and laboratory settings which support data collection and randomized A/B experiments comparing different versions of content or interactions [2]. However, there are substantial technological and practical barriers in aligning instructors and researchers to use learning technologies like blended lessons/exercises & MOOCs as both a service for students and a realistic context to conduct research. This paper explains how the concept of a MOOClet can facilitate research-practitioner collaborations. MOOClets [3] are defined as modular components of a digital resource that can be implemented in technology to: (1) allow modification to create multiple versions, (2) allow experimental comparison and personalization of different versions, (3) reliably specify what data are collected. We suggest a framework in which instructors specify what kinds of changes to lessons, exercises, and emails they would be willing to adopt, and what data they will collect and make available. Researchers can then: (1) specify or design experiments that compare the effects of different versions on quantifiable outcomes. (2) Explore algorithms for maximizing particular outcomes by choosing alternative versions of a MOOClet based on the input variables available. We present a prototype survey tool for instructors intended to facilitate practitioner researcher matches and successful collaborations.


title: Asymptotic Justification of Bandlimited Interpolation of Graph signals for Semi-Supervised Learning

abstract: Graph-based methods play an important role in unsupervised and semi-supervised learning tasks by taking into account the underlying geometry of the data set. In this paper, we consider a statistical setting for semi-supervised learning and provide a formal justification of the recently introduced framework of bandlimited interpolation of graph signals. Our analysis leads to the interpretation that, given enough labeled data, this method is very closely related to a constrained low density separation problem as the number of data points tends to infinity. We demonstrate the practical utility of our results through simple experiments.


title: Fourier Monte Carlo Renormalization Group Approach to Crystalline Membranes

abstract: The computation of the critical exponent eta characterizing the universal elastic behavior of crystalline membranes in the flat phase continues to represent challenges to theorists as well as computer simulators that manifest themselves in a considerable spread of numerical results for eta published in the literature. We present new insight to this problem that results from combining Wilson's momentum shell renormalization group method with the power of modern computer simulations based on the Fourier Monte Carlo algorithm. After discussing the ideas and difficulties underlying this combined scheme, we present a calculation of the renormalization group flow of the effective 2d Young modulus for momentum shells of different thickness. Extrapolation to infinite shell thickness allows to produce results in reasonable agreement with those obtained by functional renormalization group or by Fourier Monte Carlo simulations in combination with finite size scaling. Moreover, our new method allows for the first time to obtain a decent estimate for the value of the Wegner exponent omega that determines the leading correction to scaling, which in turn allows to refine our numerical estimate for eta previously obtained from precise finite size scaling data.