title: The turbulent life of dust grains in the supernova-driven, multi-phase interstellar medium

abstract: Dust grains are an important component of the interstellar medium (ISM) of galaxies. We present the first direct measurement of the residence times of interstellar dust in the different ISM phases, and of the transition rates between these phases, in realistic hydrodynamical simulations of the multi-phase ISM. Our simulations include a time-dependent chemical network that follows the abundances of H^+, H, H_2, C^+ and CO and take into account self-shielding by gas and dust using a tree-based radiation transfer method. Supernova explosions are injected either at random locations, at density peaks, or as a mixture of the two. For each simulation, we investigate how matter circulates between the ISM phases and find more sizeable transitions than considered in simple mass exchange schemes in the literature. The derived residence times in the ISM phases are characterised by broad distributions, in particular for the molecular, warm and hot medium. The most realistic simulations with random and mixed driving have median residence times in the molecular, cold, warm and hot phase around 17, 7, 44 and 1 Myr, respectively. The transition rates measured in the random driving run are in good agreement with observations of Ti gas-phase depletion in the warm and cold phases in a simple depletion model, although the depletion in the molecular phase is under-predicted. ISM phase definitions based on chemical abundance rather than temperature cuts are physically more meaningful, but lead to significantly different transition rates and residence times because there is no direct correspondence between the two definitions.

categories: astro-ph.GA astro-ph.SR physics.comp-ph physics.flu-dyn


title: The Ly$\alpha$ emission from high-$z$ galaxies hosting strong Damped Ly$\alpha$ systems

abstract: We study the average Ly$\alpha$ emission associated with high-$z$ strong (log $N$(H I) $\ge$ 21) damped Ly$\alpha$ systems (DLAs). We report Ly$\alpha$ luminosities ($L_{\rm Ly\alpha}$) for the full as well as various sub-samples based on $N$(H I), $z$, $(r-i)$ colours of QSOs and rest equivalent width of Si II$\lambda$1526 line (i.e., $W_{1526}$). For the full sample, we find $L_{\rm Ly\alpha}$$< 10^{41} (3\sigma)\ \rm erg\ s^{-1}$ with a $2.8\sigma$ level detection of Ly$\alpha$ emission in the red part of the DLA trough. The $L_{\rm Ly\alpha}$ is found to be higher for systems with higher $W_{1526}$ with its peak, detected at $\geq 3\sigma$, redshifted by about 300-400 $\rm km\ s^{-1}$ with respect to the systemic absorption redshift, as seen in Lyman Break Galaxies (LBGs) and Ly$\alpha$ emitters. A clear signature of a double-hump Ly$\alpha$ profile is seen when we consider $W_{1526} \ge 0.4$ \AA\ and $(r-i) < 0.05$. Based on the known correlation between metallicity and $W_{1526}$, we interpret our results in terms of star formation rate (SFR) being higher in high metallicity (mass) galaxies with high velocity fields that facilitates easy Ly$\alpha$ escape. The measured Ly$\alpha$ surface brightness requires local ionizing radiation that is 4 to 10 times stronger than the metagalactic UV background at these redshifts. The relationship between the SFR and surface mass density of atomic gas seen in DLAs is similar to that of local dwarf and metal poor galaxies. We show that the low luminosity galaxies will contribute appreciably to the stacked spectrum if the size-luminosity relation seen for H I at low-$z$ is also present at high-$z$. Alternatively, large Ly$\alpha$ halos seen around LBGs could also explain our measurements.


title: Is the dark matter particle its own antiparticle?

abstract: We propose a test based on direct detection data that allows to determine if the dark matter particle is different from its antiparticle. The test requires the precise measurement of the dark matter spin-independent direct detection cross sections off $\mathrm{three}$ different nuclei, and consists of interpreting such signals in terms of self-conjugate (particle $=$ antiparticle) dark matter to see if such interpretation is consistent. If it is not, the dark matter must be different from its antiparticle. We illustrate this procedure for two sets of target nuclei, $\mathrm{\{Xe, Ar, Si\}}$ and $\mathrm{\{Xe, Ar, Ge\}}$, identifying the regions of the parameter space where it is particularly feasible. Our results indicate that future signals in direct detection experiments, if sufficiently accurate, might be used to establish that the dark matter particle is not its own antiparticle --a major step towards the determination of the fundamental nature of the dark matter.


title: Exotic colored scalars at the LHC

abstract: We study the phenomenology of exotic color-triplet scalar particles $X$ with charge $|Q|=2/3, 4/3,5/3,7/3,8/3$ and $10/3$. If $X$ is a non-singlet of $SU(2)_W$ representation, mass splitting within the multiplet allows for cascade decays of the members into the lightest state. We study examples where the lightest state, in turn, decays into a three-body $W^\pm jj$ final state, and show that in such case the entire multiplet is compatible with existing direct collider searches and indirect precision tests down to $m_X\sim250$~GeV. However, bound states $S$, made of $XX^\dag$ pairs at $m_S\approx2m_X$, form under rather generic conditions and their decay to diphoton can be the first discovery channel of the model. Furthermore, for $SU(2)_W$-non-singlets, the mode $S\to W^+W^-$ may be observable and the width of $S\to\gamma\gamma$ and $S\to jj$ may appear large as a consequence of mass splittings within the $X$-multiplet. As an example we study in detail the case of an $SU(2)_W$ quartet, finding that $m_X\simeq450$~GeV is allowed by all current searches.


title: Constraining axionlike particles using the distance-duality relation

abstract: One of the fundamental results used in observational cosmology is the distance duality relation (DDR), which relates the luminosity distance, $D_L$, with angular diameter distance, $D_A$, at a given redshift $z$. We employ the observed limits of this relation to constrain the coupling of axion like particles (ALPs) with photons. With our detailed $3D$ ALP-photon mixing simulation in standard $\Lambda$CDM universe and latest DDR limits observed in Holand \& Barros (2016) we limit the coupling constant $g_{\phi}\le 6 \times 10^{-13} GeV^{-1}\left(\frac{nG}{\langle B \rangle_{\rm Mpc}}\right)$ for ALPs of mass $\le 10^{-15}$ eV. The DDR observations can provide very stringent constraint on ALPs mixing in future. Also any deviation in DDR can be conventionally explained as photons decaying to axions or vice-versa.


title: A Novel Averaging Technique for Discrete Entropy-Stable Dissipation Operators for Ideal MHD

abstract: Entropy stable schemes can be constructed with a specific choice of the numerical flux function. First, an entropy conserving flux is constructed. Secondly, an entropy stable dissipation term is added to this flux to guarantee dissipation of the discrete entropy. Present works in the field of entropy stable numerical schemes are concerned with thorough derivations of entropy conservative fluxes for ideal MHD. However, as we show in this work, if the dissipation operator is not constructed in a very specific way, it cannot lead to a generally stable numerical scheme. The two main findings presented in this paper are that the entropy conserving flux of Ismail & Roe can easily break down for certain initial conditions commonly found in astrophysical simulations, and that special care must be taken in the derivation of a discrete dissipation matrix for an entropy stable numerical scheme to be robust. We present a convenient novel averaging procedure to evaluate the entropy Jacobians of the ideal MHD and the compressible Euler equations that yields a discretization with favorable robustness properties.


title: Assessment of the Information Content of the Power Spectrum and Bispectrum

abstract: The covariance matrix of the matter and halo power spectrum and bispectrum are studied. Using a large suite of simulations, we find that the non-Gaussianity in the covariance is significant already at mildly nonlinear scales. We compute the leading disconnected non-Gaussian correction to the matter bispectrum covariance using perturbation theory, and find that the corrections result in good agreement in the mildly nonlinear regime. The shot noise contribution to the halo power spectrum and bispectrum covariance is computed using the Poisson model, and the model yields decent agreement with simulation results. However, when the shot noise is estimated from the individual realization, which is usually done in reality, we find that the halo covariance is substantially reduced and gets close to the Gaussian covariance. This is because most of the non-Gaussianity in the covariance arises from the fluctuations in the Poisson shot noise. We use the measured non-Gaussian covariance to access the information content of the power spectrum and bispectrum. The signal-to-noise ratio, S/N, of the matter and halo power spectrum levels off in the mildly nonlinear regime, $k \sim 0.1 - 0.2 \,\mathrm{Mpc} h^{-1}$. In the nonlinear regime the S/N of the matter and halo bispectrum increases but much slower than the Gaussian results suggest. We find that both the S/N for power spectrum and bispectrum are overestimated by the Gaussian covariances, but the problem being much more serious for the bispectrum. Because the bispectrum is affected strongly by nonlinearity and shot noise, inclusion of the bispectrum only adds modest amount of S/N compared to that of the power spectrum.


title: Parameter spaces for algebraic equivalence

abstract: A cycle is algebraically trivial if it can be exhibited as the difference of two fibers in a family of cycles parameterized by a smooth scheme. Over an algebraically closed field, it is a result of Weil that it suffices to consider families of cycles parameterized by curves, or by abelian varieties. In this paper, we extend these results to arbitrary base fields. The strengthening of these results turns out to be a key step in our work elsewhere extending Murre's results on algebraic representatives for varieties over algebraically closed fields to arbitrary perfect fields.


title: A Call for New Physics : The Muon Anomalous Magnetic Moment and Lepton Flavor Violation

abstract: We review how the muon anomalous magnetic moment ($g-2$) and the quest for lepton flavor violation are intimately correlated. Indeed the decay $\mu \to e \gamma$ is induced by the same amplitude for different choices of in- and outgoing leptons. In this work, we try to address some intriguing questions such as: Which hierarchy in the charged lepton sector one should have in order to reconcile possible signals coming simultaneously from $g-2$ and lepton flavor violation? What can we learn if the $g-2$ anomaly is confirmed by the upcoming flagship experiments at FERMILAB and J-PARC, and no signal is seen in the decay $\mu \rightarrow e\gamma$ in the foreseeable future? On the other hand, if the $\mu \rightarrow e\gamma$ decay is seen in the upcoming years, do we need to necessarily observe a signal also in $g-2$?. In this attempt, we generally study the correlation between these observables in a detailed analysis of simplified models. We derive master integrals and fully analytical and exact expressions for both phenomena, and adress other flavor violating signals. We investigate under which conditions the observations can be made compatible and discuss their implications. Lastly, we discuss in this context several extensions of the SM, such as the Minimal Supersymmetric Standard Model, Left-Right symmetric model, $B-L$ model, scotogenic model, two Higgs doublet model, Zee-Babu model, 331 model, and $L_{\mu}-L_{\tau}$, dark photon, seesaw models type~I, II and III, and also address the interplay with $\mu \to eee$ decay and $\mu-e$ conversion.


title: A composite tale of two right-handed neutrinos

abstract: In this work, we develop a model for Higgs-like composites based on two generations of right handed neutrinos which condense. We analyze the Spontaneous Symmetry Breaking of the theory with two explicit breakings, setting the different scales of the model and obtaining massive bosons as a result. Finally, we calculate the gravitational wave imprint left by the phase transition associated to the symmetry breaking of a generic potential dictated by the symmetries of the composites.


title: Large bulk photovoltaic effect and spontaneous polarization of single-layer monochalcogenides

abstract: We use a first-principles density functional theory approach to calculate the shift current and linear absorption of uniformly illuminated single-layer Ge and Sn monochalcogenides. We predict strong absorption in the visible spectrum and a large effective three-dimensional shift current ($\sim$100 $\mu$A/V$^2$), larger than has been previously observed in other polar systems. Moreover, we show that the integral of the shift-current tensor is correlated to the large spontaneous effective three-dimensional electric polarization ($\sim$1.9 C/m$^2$). Our calculations indicate that the shift current will be largest in the visible spectrum, suggesting that these monochalcogenides may be promising for polar optoelectronic devices. A Rice-Mele tight-binding model is used to rationalize the shift-current response for these systems, and its dependence on polarization, in general terms with implications for other polar materials


title: Scaling Laws of Passive-Scalar Diffusion in the Interstellar Medium

abstract: Passive scalar mixing (metals, molecules, etc.) in the turbulent interstellar medium (ISM) is critical for abundance patterns of stars and clusters, galaxy and star formation, and cooling from the circumgalactic medium. However, the fundamental scaling laws remain poorly understood in the highly supersonic, magnetized, shearing regime relevant for the ISM. We therefore study the full scaling laws governing passive-scalar transport in idealized simulations of supersonic turbulence. Using simple phenomenological arguments for the variation of diffusivity with scale based on Richardson diffusion, we propose a simple fractional diffusion equation to describe the turbulent advection of an initial passive scalar distribution. These predictions agree well with the measurements from simulations, and vary with turbulent Mach number in the expected manner, remaining valid even in the presence of a large-scale shear flow (e.g. rotation in a galactic disk). The evolution of the scalar distribution is not the same as obtained using simple, constant "effective diffusivity" as in Smagorinsky models, because the scale-dependence of turbulent transport means an initially Gaussian distribution quickly develops highly non-Gaussian tails. We also emphasize that these are mean scalings that only apply to ensemble behaviors (assuming many different, random scalar injection sites): individual Lagrangian "patches" remain coherent (poorly-mixed) and simply advect for a large number of turbulent flow-crossing times.

categories: astro-ph.GA math.NA physics.comp-ph physics.flu-dyn


title: Di-boson signatures as Standard Candles for Partial Compositeness

abstract: Composite Higgs Models are often constructed including fermionic top partners with a mass around the TeV scale, with the top partners playing the role of stabilizing the Higgs potential and enforcing partial compositeness for the top quark. A class of models of this kind can be formulated in terms of fermionic strongly coupled gauge theories. A common feature they all share is the presence of specific additional scalar resonances, namely two neutral singlets and a colored octet, described by a simple effective Lagrangian. We study the phenomenology of these scalars, both in a model independent and model dependent way, including the bounds from all the available searches in the relevant channels with di-boson and di-top final states. We develop a generic framework which can be used to constrain any model containing pseudo-scalar singlets or octets. Using it, we find that such signatures provide strong bounds on the compositeness scale complementary to the traditional EWPT and Higgs couplings deviations. In many cases a relatively light scalar can be on the verge of discovery as a first sign of new physics.


title: Defects of liquid crystals with variable degree of orientation

abstract: We prove that the zero set of a minimizer of the modified Ericksen energy for nematic liquid crystals with variable degree of orientation consists locally of isolated points and a finite union of H\"older continuous curves with only finitely many crossings.


title: Nonlinear molecular deformations give rise to stress stiffening, yielding and non-uniform stress propagation in actin networks

abstract: We use optical tweezers microrheology and fluorescence microscopy to apply nonlinear microscale strains to entangled and crosslinked actin networks, and measure the resulting stress and actin filament deformations. We couple nonlinear stress response and subsequent relaxation to the velocity profiles of individual fluorescent-labeled actin segments at varying times throughout the strain and varying distances from the strain path to determine the underlying molecular dynamics that give rise to the debated nonlinear response and stress propagation of crosslinked and entangled actin networks at the microscale. We show that initial stress stiffening arises from acceleration of strained filaments due to molecular extension along the strain, while softening and yielding is coupled to filament deceleration, halting and recoil. We demonstrate a surprising non-monotonic dependence of velocity profiles on crosslinker concentration. Namely, networks with no crosslinks or substantial crosslinks both exhibit fast initial filament velocities and reduced molecular recoil while intermediate crosslinker concentrations display reduced velocities and increased recoil. We show that these collective results are due to a balance of network elasticity and transient crosslinker unbinding and rebinding. We further show that elasticity and relaxation timescales display an exponential dependence on crosslinker concentration that reveal that crosslinks dominate entanglement dynamics when the length between crosslinkers becomes smaller than the length between entanglements. In accord with recent simulations, our relaxation dynamics demonstrate that post-strain stress can be long-lived in crosslinked networks by distributing stress to a small fraction of highly-strained connected filaments that span the network and sustain the load, thereby allowing the rest of the network to recoil and relax.


title: Efficient numerical solution to vacuum decay with many fields

abstract: Finding numerical solutions describing bubble nucleation is notoriously difficult in more than one field space dimension. Traditional shooting methods fail because of the extreme non-linearity of field evolution over a macroscopic distance as a function of initial conditions. Minimization methods tend to become either slow or imprecise for larger numbers of fields due to their dependence on the high dimensionality of discretized function spaces. We present a new method for finding solutions which is both very efficient and able to cope with the non-linearities. Our method directly integrates the equations of motion except at a small number of junction points, so we do not need to introduce a discrete domain for our functions. The method, based on multiple shooting, typically finds solutions involving three fields in under a minute, and can find solutions for eight fields in about an hour. We include a numerical package for Mathematica which implements the method described here.


title: Structure and Physical Conditions in the Huygens Region of the Orion Nebula

abstract: HST images, MUSE maps of emission-lines, and an atlas of high velocity resolution emission-line spectra have been used to establish for the firrst time correlations of the electron temperature, electron density, radial velocity, turbulence, and orientation within the main ionization front of the nebula. From the study of the combined properties of multiple features, it is established that variations in the radial velocity are primarily caused by the photo-evaporating ionization front being viewed at different angles. There is a progressive increase of the electron temperature and density with decreasing distance from the dominant ionizing star Theta1 Ori C. The product of these characteristics (NexTe) is the most relevant parameter in modeling a blister-type nebula like the Huygens Region, where this quantity should vary with the surface brightness in Halpha. Several lines of evidence indicate that small-scale structure and turbulence exists down to the level of our resolution of a few arcseconds. Although photo-evaporative ow must contribute at some level to the well-known non-thermal broadening of the emission lines, comparison of quantitative predictions with the observed optical line widths indicate that it is not the major additive broad- ening component. Derivation of Te values for H+ from radio+optical and optical-only ionized hydro- gen emission showed that this temperature is close to that derived from [Nii] and that the transition from the well-known at extinction curve that applies in the Huygens Region to a more normal steep extinction curve occurs immediately outside of the Bright Bar feature of the nebula.


title: H-theorem and Thermodynamics for generalized entropies that depend only on the probability

abstract: We consider a previously proposed non-extensive statistical mechanics in which the entropy depends only on the probability, this was obtained from a f(\beta) distribution and its corresponding Boltzmann factor. We show that the first term correcting the usual entropy also arises from several f(\beta) distributions, we also construct the corresponding H-function and demonstrate that a generalized H-theorem is fulfilled. Furthermore, expressing this H-function as function of the simplest Maxwellian state we find, up to a first approximation some modified thermodynamic quantities for an ideal gas. In order to gain some insight about the behavior of the proposed generalized statistics, we present some simulation results for the case of a square-well and Lennard-Jones potentials, showing that an effective repulsive interaction is obtained with the new formalism.


title: Standard Model Background of the Cosmological Collider

abstract: The inflationary universe can be viewed as a "Cosmological Collider" with energy of Hubble scale, producing very massive particles and recording their characteristic signals in primordial non-Gaussianities. To utilize this collider to explore any new physics at very high scales, it is a prerequisite to understand the background signals from the particle physics Standard Model. In this paper we describe the Standard Model background of the Cosmological Collider.


title: When do Correlations Increase with Firing Rate?

abstract: A central question in neuroscience is to understand how noisy firing patterns are used to transmit information. Because neural spiking is noisy, spiking patterns are often quantified via pairwise correlations, or the probability that two cells will spike coincidentally, above and beyond their baseline firing rate. One observation frequently made in experiments, is that correlations can increase systematically with firing rate. Theoretical studies have determined that stimulus-dependent correlations that increase with firing rate can have beneficial effects on information coding; however, we still have an incomplete understanding of what circuit mechanisms do, or do not, produce this correlation-firing rate relationship. Here, we study the relationship between pairwise correlations and firing rates in recurrently coupled excitatory-inhibitory spiking networks with conductance-based synapses. We found that with stronger excitatory coupling, a positive relationship emerges between pairwise correlations and firing rates. To explain these findings, we used linear response theory to predict the full correlation matrix and to decompose correlations in terms of graph motifs. We then used this decomposition to explain why covariation of correlations with firing rate -- a relationship previously explained in feedforward networks driven by correlated input -- emerges in some recurrent networks but not in others. Furthermore, when correlations covary with firing rate, this relationship is reflected in low-rank structure in the correlation matrix.


title: Euclidean distance matrix completion and point configurations from the minimal spanning tree

abstract: The paper introduces a special case of the Euclidean distance matrix completion problem (edmcp) of interest in statistical data analysis where only the minimal spanning tree distances are given and the matrix completion must preserve the minimal spanning tree. Two solutions are proposed, one an adaptation of a more general method based on a dissimilarity parameterized formulation, the other an entirely novel method which constructs the point configuration directly through a guided random search. These methods as well as three standard edcmp methods are described and compared experimentally on real and synthetic data. It is found that the constructive method given by the guided random search algorithm clearly outperforms all others considered here. Notably, standard methods including the adaptation force peculiar, and generally unwanted, geometric structure on the point configurations their completions produce.


title: Completeness of the Nearby White Dwarf Sample: Let Us Count the Ways

abstract: We have recently extended our ongoing survey of the local white dwarf population, effectively doubling the sample volume. Based upon the latest distance determinations, Holberg et al. (2016) estimated the present 20 pc and 25 pc samples were about 86 and 68 percent complete, respectively. Here we examine how the completeness of the 25 pc sample depends upon other observables such as apparent magnitude, proper motion, photometric color index, etc. The results may provide additional clues to why Sirius-Like systems are underrepresented in the extended 25 pc sample and how additional nearby single white dwarf stars may be found.


title: An Approach to Speed-up the Word Sense Disambiguation Procedure through Sense Filtering

abstract: In this paper, we are going to focus on speed up of the Word Sense Disambiguation procedure by filtering the relevant senses of an ambiguous word through Part-of-Speech Tagging. First, this proposed approach performs the Part-of-Speech Tagging operation before the disambiguation procedure using Bigram approximation. As a result, the exact Part-of-Speech of the ambiguous word at a particular text instance is derived. In the next stage, only those dictionary definitions (glosses) are retrieved from an online dictionary, which are associated with that particular Part-of-Speech to disambiguate the exact sense of the ambiguous word. In the training phase, we have used Brown Corpus for Part-of-Speech Tagging and WordNet as an online dictionary. The proposed approach reduces the execution time upto half (approximately) of the normal execution time for a text, containing around 200 sentences. Not only that, we have found several instances, where the correct sense of an ambiguous word is found for using the Part-of-Speech Tagging before the Disambiguation procedure.


title: Iterative Refinement for Machine Translation

abstract: Existing machine translation decoding algorithms generate translations in a strictly monotonic fashion and never revisit previous decisions. As a result, earlier mistakes cannot be corrected at a later stage. In this paper, we present a translation scheme that starts from an initial guess and then makes iterative improvements that may revisit previous decisions. We parameterize our model as a convolutional neural network that predicts discrete substitutions to an existing translation based on an attention mechanism over both the source sentence as well as the current translation output. By making less than one modification per sentence, we improve the output of a phrase-based translation system by up to 0.4 BLEU on WMT15 German-English translation.


title: Combinatorial Multi-Armed Bandit with General Reward Functions

abstract: In this paper, we study the stochastic combinatorial multi-armed bandit (CMAB) framework that allows a general nonlinear reward function, whose expected value may not depend only on the means of the input random variables but possibly on the entire distributions of these variables. Our framework enables a much larger class of reward functions such as the $\max()$ function and nonlinear utility functions. Existing techniques relying on accurate estimations of the means of random variables, such as the upper confidence bound (UCB) technique, do not work directly on these functions. We propose a new algorithm called stochastically dominant confidence bound (SDCB), which estimates the distributions of underlying random variables and their stochastically dominant confidence bounds. We prove that SDCB can achieve $O(\log{T})$ distribution-dependent regret and $\tilde{O}(\sqrt{T})$ distribution-independent regret, where $T$ is the time horizon. We apply our results to the $K$-MAX problem and expected utility maximization problems. In particular, for $K$-MAX, we provide the first polynomial-time approximation scheme (PTAS) for its offline problem, and give the first $\tilde{O}(\sqrt T)$ bound on the $(1-\epsilon)$-approximation regret of its online problem, for any $\epsilon>0$.


title: Production of exotic charmonium in $\gamma \gamma$ interactions at hadronic colliders

abstract: In this paper we investigate the Exotic Charmonium (EC) production in $\gamma \gamma$ interactions present in proton-proton, proton-nucleus and nucleus-nucleus collisions at the CERN Large Hadron Collider (LHC) energies as well as for the proposed energies of the Future Circular Collider (FCC). Our results demonstrate that the experimental study of these processes is feasible and can be used to constrain the theoretical decay widths and shed some light on the configuration of the considered multiquark states.


title: Mathematical analysis of inertial waves in rectangular basins with one sloping boundary

abstract: Here we consider the problem of small oscillations of a rotating inviscid incompressible fluid. From a mathematical point of view, new exact solutions to the two-dimensional Poincar\'e-Sobolev equation in a class of domains including trapezoid are found in an explicit form and their main properties are described. These solutions correspond to the absolutely continuous spectrum of a linear operator that is associated with this system of equations. For specialists in Astrophysics and Geophysics the existence of these solutions signifies the existence of some previously unknown type of inertial waves corresponding to the continuous spectrum of inertial oscillations. A fundamental distinction between monochromatic inertial waves and waves of the new type is shown: usual characteristics (frequency, amplitude, wave vector, dispersion relation, direction of energy propagation, and so on) are not applicable to the last. Main properties of these waves are described. In particular it is proved that they are progressive. Main features of their energy transfer are described. The existence of such inertial waves enables us to explain in a new way a lot of experimental data that were obtained in Geophysics in the past two decades and to predict the occurrence of such oscillations in natural waters.


title: Morphological and Biochemical Adaptive Changes Associated With A Short-period Starvation of Adult Male Japanese Quail (Coturnix japonica)

abstract: Background The morphological and biochemical impact of a short-period of starvation on Japanese quail was investigated. Materials and methods Ten adult male Japanese quail were divided into two groups; control fed and starved. The control-fed group was offered food and water ad libitum and the starved group was subjected to a short-period of food deprivation. After 2.5 days, the serum was obtained and different parameters including the total protein, AST, ALT, triglyceride, HDL, LDL, creatinine and urea were assessed. Gastrointestinal tract, stomach and liver were excised and their masses were estimated. Paraffin and resin embedded sections from the proventriculus, gizzard, liver, duodenum, kidney and pancreas were examined with a light microscopy. Results Significant decreases in the masses of body, gastrointestinal tract, stomach and liver of the starved group were recorded. The liver and duodenum were the most affected organs. The liver showed depletion of glycogen, vacuolation, hyperemia and cellular infiltrations. Duodenal villi showed degenerative changes in lamina epithelialis and cellular infiltrations in the lamina propria. Biochemical analysis revealed a decreased level of total protein, AST and ALT, increased cholesterol, triglycerides and LDL, and unchanged HDL, urea and creatinine by starvation. Conclusion The current study described in details the effect of short time starvation on quail organs. Time-point adaptive responses of male quail to starvation and refeeding on quail organs will be investigated in future studies.


title: Optimal Bitter Coil Solenoid

abstract: Bitter coil is an electromagnet used for the generation of exceptionally strong magnetic fields. The upper bound of magnet flux density is restricted by several factors. One principal restriction is the high stresses due to Lorentz forces in the coil. The Lorentz forces generate the distributed body force, which acts as the pressure of magnetic field. The common radial thickness profile of the Bitter coil is constant. In this paper the possibility of optimization by means of non-constant radial thickness profile of the Bitter coil is studied. The close form expression for optimal thickness profile is obtained. Both designs are compared and the considerable improvement of magnetic flux density is demonstrated. Moreover, the optimal design improves the shape of cooling channels. Namely, the highest cross-section of cooling channel is at the most thermally loaded inner surface of the coil.


title: Lumbar degenerative spondylolisthesis epidemiology: a systemic review with a focus on gender-specific and age-specific prevalence

abstract: The epidemiology of lumbar degenerative spondylolisthesis (DS) remains controversial. We performed a systemic review with the aim to have a better understanding of DS's prevalence in general population. The results showed the prevalence of DS is very gender specific and age specific. Both women and men have few DS before 50 years old, after 50 years old both women and men start to develop DS, with women having a faster developing rate than men. For elderly Chinese (>=65 yrs, mean age: 72.5 yrs), large population based studies (MsOS(Hong Kong) and MrOS (Hong Kong), females n=2000 and males n=2000) showed DS prevalence was for 25.0% for women and 19.1% for men, and the prevalence F:M (women:men) ratio was 1.3:1. The published data (MsOS(USA) and MrOS(USA) studies) seem to show elderly Caucasian American has a higher DS prevalence, being approximately 60-70% higher than elderly Chinese; however the prevalence F:M ratio was similar to elderly Chinese population. Patient data showed female patients more often received treatment than men; and preliminary data show the ratio of numbers of female patients received treatment compared with men did not differ between Northeast Asians (Chinese, Japanese, and Korean) and European and American Caucasians, being around 2:1 in elderly population. The existing data also suggest that menopause may be a contributing factor for the accelerated DS development in post-menopausal


title: BP Piscium: its flaring disk imaged with SPHERE/ZIMPOL

abstract: Whether BP Piscium (BP Psc) is either a pre-main sequence T Tauri star at d ~ 80 pc, or a post-main sequence G giant at d ~ 300 pc is still not clear. As a first-ascent giant, it is the first to be observed with a molecular and dust disk. Alternatively, BP Psc would be among the nearest T Tauri stars with a protoplanetary disk (PPD). We investigate whether the disk geometry resembles typical PPDs, by comparing polarimetric images with radiative transfer models. Our VLT/SPHERE/ZIMPOL observations allow us to perform Polarimetric Differential Imaging; Reference Star Differential Imaging; and Richardson-Lucy deconvolution. We present the first visible light polarization and intensity images of the disk of BP Psc. Our deconvolution confirms the disk shape as detected before, mainly showing the southern side of the disk. In polarized intensity the disk is imaged at larger detail and also shows the northern side, giving it the typical shape of high inclination flared disks. We explain the observed disk features by retrieving the large-scale geometry with MCMax radiative transfer modeling, which yields a strongly flared model, atypical for disks of T Tauri stars.


title: Degrees of regular sequences with a symmetric group action

abstract: We consider ideals in a polynomial ring that are generated by regular sequences of homogeneous polynomials and are stable under the action of the symmetric group permuting the variables. In previous work, we determined the possible isomorphism types for these ideals. Following up on that work, we now analyze the possible degrees of the elements in such regular sequences. For each case of our classification, we provide some criteria guaranteeing the existence of regular sequences in certain degrees.


title: The Broadband Spectral Variability of Holmberg IX X-1

abstract: We present results from four new broadband X-ray observations of the extreme ultraluminous X-ray source Holmberg IX X-1 ($L_{\rm{X}} > 10^{40}$ erg/s), performed by $Suzaku$ and $NuSTAR$ in coordination. Combined with the archival data, we now have broadband observations of this remarkable source from six separate epochs. Two of these new observations probe lower fluxes than seen previously, allowing us to extend our knowledge of the broadband spectral variability exhibited. The spectra are well fit by two thermal blackbody components, which dominate the emission below 10 keV, as well as a steep ($\Gamma \sim 3.5$) powerlaw tail which dominates above $\sim$15 keV. Remarkably, while the 0.3-10.0 keV flux varies by a factor of $\sim$3 between all these epochs, the 15-40 keV flux varies by only $\sim$20%. Although the spectral variability is strongest in the $\sim$1-10 keV band, both of the thermal components are required to vary when all epochs are considered. We also re-visit the search for iron absorption features, leveraging the high-energy $NuSTAR$ data to improve our sensitivity to extreme velocity outflows in light of the ultra-fast outflow recently detected in NGC 1313 X-1. Iron absorption from a similar outflow along our line of sight can be ruled out in this case. We discuss these results in the context of super-Eddington accretion models that invoke a funnel-like geometry for the inner flow, and propose a scenario in which we have an almost face-on view of a funnel that expands to larger radii with increasing flux, resulting in an increasing degree of geometrical collimation for the emission from intermediate temperature regions.


title: Toric surfaces over an arbitrary field

abstract: We study toric varieties over an arbitrary field with an emphasis on toric surfaces in the Merkurjev-Panin motivic category of "K-motives". We explore the decomposition of certain toric varieties as K-motives into products of central simple algebras, the geometric and topological information encoded in these central simple algebras, and the relationship between the decomposition of the K-motives and the semiorthogonal decomposition of the derived categories. We obtain the information mentioned above for toric surfaces by explicitly classifying all minimal smooth projective toric surfaces using toric geometry.


title: The fixation probability and time for a doubly beneficial mutant

abstract: For a highly beneficial mutant $A$ entering a randomly reproducing population of constant size, we study the situation when a second beneficial mutant $B$ arises before $A$ has fixed. If the selection coefficient of $B$ is greater than the selection coefficient of $A$, and if $A$ and $B$ can recombine at some rate $\rho$, there is a chance that the double beneficial mutant $AB$ forms and eventually fixes. We give a convergence result for the fixation probability of $AB$ and its fixation time for large selection coefficients.


title: Solving Multi-Objective Optimization via Adaptive Stochastic Search with Domination Measure

abstract: For general multi-objective optimization problems, we propose a novel performance metric called domination measure to measure the quality of a solution, which can be intuitively interpreted as the size of the portion of the solution space that dominates that solution. As a result, we reformulate the original multi-objective problem into a stochastic single-objective one and propose a model-based approach to solve it. We show that an ideal version algorithm of the proposed approach converges to a set representation of the global optima of the reformulated problem. We also investigate the numerical performance of an implementable version algorithm by comparing it with numerous existing multi-objective optimization methods on popular benchmark test functions. The numerical results show that the proposed approach is effective in generating a finite and uniformly spread approximation of the Pareto optimal set of the original multi-objective problem, and is competitive to the tested existing methods.


title: A single electron transistor with quantum rings

abstract: We Have developed the concept of a new kind of single-electron transistor in which the transport of the electron through a quantum wire is controlled by charged quantum rings. Using a 2D harmonic potential as the transverse constraint, we numerically investigated the transport of the electron through the wire. We have shown that in the low energy limit, for a suitable configuration of the rings, called the quadrupole configuration, we are able to adjust the conductance of the wire and therefore control the switching process.


title: Measurement of the cross section and longitudinal double-spin asymmetry for di-jet production in polarized $pp$ collisions at $\sqrt{s}$ = 200 GeV

abstract: We report the first measurement of the longitudinal double-spin asymmetry $A_{LL}$ for mid-rapidity di-jet production in polarized $pp$ collisions at a center-of-mass energy of $\sqrt{s} = 200$ GeV. The di-jet cross section was measured and is shown to be consistent with next-to-leading order (NLO) perturbative QCD predictions. $A_{LL}$ results are presented for two distinct topologies, defined by the jet pseudorapidities, and are compared to predictions from several recent NLO global analyses. The measured asymmetries, the first such correlation measurements, support those analyses that find positive gluon polarization at the level of roughly 0.2 over the region of Bjorken-$x > 0.05$.


title: Generating invariant rings of quivers in arbitrary characteristic

abstract: It is well known that the ring of polynomial invariants of a reductive group is finitely generated. However, it is difficult to give strong upper bounds on the degrees of the generators, especially over fields of positive characteristic. In this paper, we make use of the theory of good filtrations along with recent results on the null cone to provide polynomial bounds for matrix semi-invariants in arbitrary characteristic, and consequently for matrix invariants. Our results generalize to invariants and semi-invariants of quivers.


title: An Efficient Optimal Algorithm for the Successive Minima Problem

abstract: In many applications including integer-forcing linear multiple-input and multiple-output (MIMO) receiver design, one needs to solve a successive minima problem (SMP) on an $n$-dimensional lattice to get an optimal integer coefficient matrix $\A^\star\in \mathbb{Z}^{n\times n}$. In this paper, we first propose an efficient optimal SMP algorithm with an $\bigO(n^2)$ memory complexity. The main idea behind the new algorithm is it first initializes with a suitable suboptimal solution, which is then updated, via a novel algorithm with only $\bigO(n^2)$ flops in each updating, until $\A^{\star}$ is obtained. Different from existing algorithms which find $\A^\star$ column by column through using a sphere decoding search strategy $n$ times, the new algorithm uses a search strategy once only. We then rigorously prove the optimality of the proposed algorithm. Furthermore, we theoretically analyze its complexity. In particular, we not only show that the new algorithm is $\Omega(n)$ times faster than the most efficient existing algorithm with polynomial memory complexity, but also assert that it is even more efficient than the most efficient existing algorithm with exponential memory complexity. Finally, numerical simulations are presented to illustrate the optimality and efficiency of our novel SMP algorithm.


title: Functional Asynchronous Networks: Factorization of Dynamics and Function

abstract: In this note we describe the theory of functional asynchronous networks and one of the main results, the Modularization of Dynamics Theorem, which for a large class of functional asynchronous networks gives a factorization of dynamics in terms of constituent subnetworks. For these networks we can give a complete description of the network function in terms of the function of the events comprising the network and thereby answer a question originally raised by Alon in the context of biological networks.


title: Proposing Plausible Answers for Open-ended Visual Question Answering

abstract: Answering open-ended questions is an essential capability for any intelligent agent. One of the most interesting recent open-ended question answering challenges is Visual Question Answering (VQA) which attempts to evaluate a system's visual understanding through its answers to natural language questions about images. There exist many approaches to VQA, the majority of which do not exhibit deeper semantic understanding of the candidate answers they produce. We study the importance of generating plausible answers to a given question by introducing the novel task of `Answer Proposal': for a given open-ended question, a system should generate a ranked list of candidate answers informed by the semantics of the question. We experiment with various models including a neural generative model as well as a semantic graph matching one. We provide both intrinsic and extrinsic evaluations for the task of Answer Proposal, showing that our best model learns to propose plausible answers with a high recall and performs competitively with some other solutions to VQA.


title: Non-reciprocal quantum interactions and devices via autonomous feed-forward

abstract: In a recent work [A. Metelmann and A. A. Clerk, Phys. Rev. X 5, 021025 (2015)], a general reservoir-engineering approach for generating non-reciprocal quantum interactions and devices was described. We show here how in many cases this general recipe can be viewed as an example of autonomous feed-forward: the full dissipative evolution is identical to the unconditional evolution in a setup where an observer performs an ideal quantum measurement of one system, and then uses the results to drive a second system. We also extend the application of this approach to non-reciprocal quantum amplifiers, showing the added functionality possible when using two engineered reservoirs. In particular, we demonstrate how to construct an ideal phase-preserving cavity-based amplifier which is full non-reciprocal, quantum-limited and free of any fundamental gain-bandwidth constraint.


title: Double stage nematic bond-ordering above double stripe magnetism: application to BaTi$_2$Sb$_2$O

abstract: Spin-driven nematicity, or the breaking of the point-group symmetry of the lattice without long-range magnetic order, is clearly quite important in iron-based superconductors. From a symmetry point of view, nematic order can be described as a coherent locking of spin fluctuations in two interpenetrating N\'eel sublattices with ensuing nearest-neighbor bond order and an absence of static magnetism. Here, we argue that the low-temperature state of the recently discovered superconductor BaTi$_2$Sb$_2$O is a strong candidate for a more exotic form of spin-driven nematic order, in which fluctuations occurring in four N\'eel sublattices promote both nearest- and next-nearest neighbor bond order. We develop a low-energy field theory of this state and show that it can have, as a function of temperature, up to two separate bond-order phase transitions -- namely, one that breaks rotation symmetry and one that breaks reflection and translation symmetries of the lattice. The resulting state has an orthorhombic lattice distortion, an intra-unit-cell charge density wave, and no long-range magnetic order, all consistent with reported measurements of the low-temperature phase of BaTi$_2$Sb$_2$O. We then use density functional theory calculations to extract exchange parameters to confirm that the model is applicable to BaTi$_2$Sb$_2$O.


title: Fractional Schrodinger-Poisson Equations With General Nonlinearities

abstract: In this paper we investigate the existence of positive solutions and ground state solution for a class of fractional Schr\"odinger-Poisson equations in $\mathbb R^3$ with general nonlinearities.


title: Degree of W-operator and Noncrossing Partition

abstract: Goulden and Jackson first introduced the cut-and-join operator. The cut-and-join is widely used in studying the Hurwitz number and many other topological recursion problems. Mironov, Morosov and Natanzon give a more general construction and call it W-operator W([n]). As a special case, the cut-and-join operator is W([2]). In this paper, we study the structure of W([n]). We prove that W([n]) can be written as the sum of n! terms and each term corresponds uniquely to a permutation in Sn. We also define the degree of each term. We prove that there is a correspondence between the terms of W([n]) with highest degree and the noncrossing partitions.


title: M@TE - Monitoring at TeV Energies

abstract: Blazars are extremely variable objects emitting radiation across the electromagnetic spectrum and showing variability on time scales from minutes to years. For the understanding of the emission mechanisms, simultaneous multi-wavelength observations are crucial. Various models for flares predict simultaneous flux increases in the X-ray and gamma-ray band or more complex variability patterns, depending on the dominant process responsible for the gamma-ray emission. Monitoring at TeV energies is providing important information to distinguish between different models. To study duty cycle and variability time scales of an object, an unbiased data sample is essential, and good sensitivity and continuous monitoring are needed to resolve smaller time scales. A dedicated long-term monitoring program at TeV energies has been started by the FACT project. Its success clearly illustrated that the usage of silicon based photo sensors (SIPMs) is ideal for long-term monitoring. They provide not only an excellent and stable detector performance, but also allow for observations during bright ambient light minimizing observational gaps and increasing the instrument's duty cycle. The observation time in a single longitude is limited to 6 hours. To study typical variability time scales of few hours to one day, the ultimate goal is 24/7 monitoring with a network of small telescopes around the globe (DWARF project). The installation of an Imaging Air Cherenkov Telescope is planned in San Pedro Martir, Mexico. For the M@TE (Monitoring at TeV energies) telescope, a mount from a previous experiment is being refurbished to be equipped with a camera using the new generation of SiPMs. In the presentation, the status of the M@TE project will be reported outlining the scientific potential, including the possibility to extend monitoring campaigns to 12 hours by coordinated observations together with FACT.


title: Chimera states in a network-organized public goods game with destructive agents

abstract: We found that a network-organized metapopulation of cooperators, defectors and destructive agents playing the public goods game with mutations, can collectively reach global synchronization or chimera states. Global synchronization is accompanied by a collective periodic burst of cooperation, whereas chimera states reflect the tendency of the networked metapopulation to be fragmented in clusters of synchronous and incoherent bursts of cooperation. Numerical simulations have shown that the system's dynamics alternates between these two steady states through a first order transition. Depending on the parameters determining the dynamical and topological properties, chimera states with different numbers of coherent and incoherent clusters are observed. Our results present the first systematic study of chimera states and their characterization in the context of evolutionary game theory. This provides a valuable insight into the details of their occurrence, extending the relevance of such states to natural and social systems.


title: First Study of Combined Blazar Light Curves with FACT and HAWC

abstract: For studying variable sources like blazars, it is crucial to achieve unbiased monitoring, either with dedicated telescopes in pointing mode or survey instruments. At TeV energies, the High Altitude Water Cherenkov (HAWC) observatory monitors approximately two thirds of the sky every day. It uses the water Cherenkov technique, which provides an excellent duty cycle independent of weather and season. The First G-APD Cherenkov Telescope (FACT) monitors a small sample of sources with better sensitivity, using the imaging air Cherenkov technique. Thanks to its camera with silicon-based photosensors, FACT features an excellent detector performance and stability and extends its observations to times with strong moonlight, increasing the duty cycle compared to other imaging air Cherenkov telescopes. As FACT and HAWC have overlapping energy ranges, a joint study can exploit the longer daily coverage given that the observatories' locations are offset by 5.3 hours. Furthermore, the better sensitivity of FACT adds a finer resolution of features on hour-long time scales, while the continuous duty cycle of HAWC ensures evenly sampled long-term coverage. Thus, the two instruments complement each other to provide a more complete picture of blazar variability. In this presentation, the first joint study of light curves from the two instruments will be shown, correlating long-term measurements with daily sampling between air and water Cherenkov telescopes. The presented results focus on the study of the variability of the bright blazars Mrk 421 and Mrk 501 during the last two years featuring various flaring activities.


title: The determinism and boundedness of self-assembling structures

abstract: Self-assembly processes are widespread in nature, and lie at the heart of many biological and physical phenomena. The characteristics of self-assembly building blocks determine the structures that they form. Two crucial properties are the determinism and boundedness of the self-assembly. The former tells us whether the same set of building blocks always generates the same structure, and the latter whether it grows indefinitely. These properties are highly relevant in the context of protein structures, as the difference between deterministic protein self-assembly and nondeterministic protein aggregation is central to a number of diseases. Here we introduce a graph theoretical approach that can determine the determinism and boundedness for several geometries and dimensionalities of self-assembly more accurately and quickly than conventional methods. We apply this methodology to a previously studied lattice self-assembly model and discuss generalizations to a wide range of other self-assembling systems


title: Solvable normal subgroups of 2-knot groups

abstract: If $X$ is an orientable, strongly minimal $PD_4$-complex and $\pi_1(X)$ has one end then it has no nontrivial locally-finite normal subgroup. Hence if $\pi$ is a 2-knot group then (a) if $\pi$ is virtually solvable then either $\pi$ has two ends or $\pi\cong\Phi$, with presentation $\langle{a,t}|ta=a^2t\rangle$, or $\pi$ is torsion-free and polycyclic of Hirsch length 4; (b) either $\pi$ has two ends, or $\pi$ has one end and the centre $\zeta\pi$ is torsion-free, or $\pi$ has infinitely many ends and $\zeta\pi$ is finite; and (c) the Hirsch-Plotkin radical $\sqrt\pi$ is nilpotent.


title: Selective addressing of solid-state spins at the nanoscale via magnetic resonance frequency encoding

abstract: The nitrogen-vacancy (NV) centre in diamond is a leading platform for nanoscale sensing and imaging, as well as quantum information processing in the solid state. To date, individual control of two NV electronic spins at the nanoscale has been demonstrated. However, a key challenge is to scale up such control to arrays of NV spins. Here we apply nanoscale magnetic resonance frequency encoding to realize site-selective addressing and coherent control of a four-site array of NV spins. Sites in the array are separated by 100 nm, with each site containing multiple NVs separated by ~15 nm. Microcoils fabricated on the diamond chip provide electrically tuneable magnetic-field gradients ~0.1 G/nm. Tailored application of gradient fields and resonant microwaves allow site-selective NV spin manipulation and sensing applications, including Rabi oscillations, imaging, and nuclear magnetic resonance (NMR) spectroscopy with nanoscale resolution. Microcoil-based magnetic resonance of solid-state spins provides a practical platform for quantum-assisted sensing, quantum information processing, and the study of nanoscale spin networks.


title: Inverse Power Flow Problem

abstract: This paper formulates an inverse power flow problem which is to infer a nodal admittance matrix (hence the network structure of a power system) from voltage and current phasors measured at a number of buses. We show that the admittance matrix can be uniquely identified from a sequence of measurements corresponding to different steady states when every node in the system is equipped with a measurement device, and a Kron-reduced admittance matrix can be determined even if some nodes in the system are not monitored (hidden nodes). Furthermore, we propose effective algorithms based on graph theory to uncover the actual admittance matrix of radial systems with hidden nodes. We provide theoretical guarantees for the recovered admittance matrix and demonstrate that the actual admittance matrix can be fully recovered even from the Kron-reduced admittance matrix under some mild assumptions. Simulations on standard test systems confirm that these algorithms are capable of providing accurate estimates of the admittance matrix from noisy sensor data.


title: Sampling hyperparameters in hierarchical models: improving on Gibbs for high-dimensional latent fields and large data sets

abstract: We consider posterior sampling in the very common Bayesian hierarchical model in which observed data depends on high-dimensional latent variables that, in turn, depend on relatively few hyperparameters. When the full conditional over the latent variables has a known form, the marginal posterior distribution over hyperparameters is accessible and can be sampled using a Markov chain Monte Carlo (MCMC) method on a low-dimensional parameter space. This may improve computational efficiency over standard Gibbs sampling since computation is not over the high-dimensional space of latent variables and correlations between hyperparameters and latent variables become irrelevant. When the marginal posterior over hyperparameters depends on a fixed-dimensional sufficient statistic, precomputation of the sufficient statistic renders the cost of the low-dimensional MCMC independent of data size. Then, when the hyperparameters are the primary variables of interest, inference may be performed in big-data settings at modest cost. Moreover, since the form of the full conditional for the latent variables does not depend on the form of the hyperprior distribution, the method imposes no restriction on the hyperprior, unlike Gibbs sampling that typically requires conjugate distributions. We demonstrate these efficiency gains in four computed examples.


title: Novelty Learning via Collaborative Proximity Filtering

abstract: The vast majority of recommender systems model preferences as static or slowly changing due to observable user experience. However, spontaneous changes in user preferences are ubiquitous in many domains like media consumption and key factors that drive changes in preferences are not directly observable. These latent sources of preference change pose new challenges. When systems do not track and adapt to users' tastes, users lose confidence and trust, increasing the risk of user churn. We meet these challenges by developing a model of novelty preferences that learns and tracks latent user tastes. We combine three innovations: a new measure of item similarity based on patterns of consumption co-occurrence; model for {\em spontaneous} changes in preferences; and a learning agent that tracks each user's dynamic preferences and learns individualized policies for variety. The resulting framework adaptively provides users with novelty tailored to their preferences for change per se.


title: Characteristic Polynomials of Symmetric Matrices over the Univariate Polynomial Ring

abstract: Viewing a bivariate polynomial f in R[x,t] as a family of univariate polynomials in t parametrized by real numbers x, we call f real rooted if this family consists of monic polynomials with only real roots. If f is the characteristic polynomial of a symmetric matrix with entries in R[x], it is obviously real rooted. In this article the converse is established, namely that every real rooted bivariate polynomial is the characteristic polynomial of a symmetric matrix over the univariate real polynomial ring. As a byproduct we present a purely algebraic proof of the Helton-Vinnikov Theorem which solved the 60 year old Lax conjecture on the existence of definite determinantal representation of ternary hyperbolic forms.


title: Observation of vortex-antivortex pairing in decaying 2D turbulence of a superfluid gas

abstract: In a two-dimensional (2D) classical fluid, a large-scale flow structure emerges out of turbulence, which is known as the inverse energy cascade where energy flows from small to large length scales. An interesting question is whether this phenomenon can occur in a superfluid, which is inviscid and irrotational by nature. Atomic Bose-Einstein condensates (BECs) of highly oblate geometry provide an experimental venue for studying 2D superfluid turbulence, but their full investigation has been hindered due to a lack of the circulation sign information of individual quantum vortices in a turbulent sample. Here, we demonstrate a vortex sign detection method by using Bragg scattering, and we investigate decaying turbulence in a highly oblate BEC at low temperatures, with our lowest being $\sim 0.5 T_c$, where $T_c$ is the superfluid critical temperature. We observe that weak spatial pairing between vortices and antivortices develops in the turbulent BEC, which corresponds to the vortex-dipole gas regime predicted for high dissipation. Our results provide a direct quantitative marker for the survey of various 2D turbulence regimes in the BEC system.


title: Three phase classification of an uninterrupted traffic flow: a $k$-means clustering study

abstract: We investigate the speed time series of the vehicles recorded by a camera at a section of a highway in the city of Isfahan, Iran. Using $k$-means clustering algorithm, we find that the natural number of clustering for this set of data is $3$. This is in agreement with the three-phase theory of uninterrupted traffic flows. According to this theory, the three traffic phases are categorized as {\em free flow} (F), {\em synchronized} (S) and {\em wide moving jam} (J). We obtain the transition speeds and densities at F-S and also S-J transitions. We also apply the Shannon entropy analysis on the speed time series over finite windows, which equips us to monitor in areal time the instant state of a traffic flow.


title: Decoherence Can Relax Cosmic Acceleration: an Example

abstract: We investigate back reaction in de Sitter space in an approach where only states that are observationally accessible are included in the density matrix. Using the Bunch-Davies vacuum as the initial condition we find for a conformal scalar field and a cosmological constant that tracing over the unobservable states beyond the cosmological horizon leads to a thermal spectrum of particles and that such a configuration is unstable under semi-classical back reaction. It is concluded that this prescription results in an instability of de Sitter space with a gradually increasing horizon size.


title: Additive unit structure of endomorphism rings and invariance of modules

abstract: We use the type theory for rings of operators due to Kaplansky to describe the structure of modules that are invariant under automorphisms of their injective envelopes. Also, we highlight the importance of Boolean rings in the study of such modules. As a consequence of this approach, we are able to further the study initiated by Dickson and Fuller regarding when a module invariant under automorphisms of its injective envelope is invariant under any endomorphism of it. In particular, we find conditions for several classes of noetherian rings which ensure that modules invariant under automorphisms of their injective envelopes are quasi-injective. In the case of a commutative noetherian ring, we show that any automorphism-invariant module is quasi-injective. We also provide multiple examples that show that our conditions are the best possible, in the sense that if we relax them further then there exist automorphism-invariant modules which are not quasi-injective. We finish this paper by dualizing our results to the automorphism-coinvariant case.


title: Localized $L^p$-estimates for eigenfunctions: II

abstract: If $(M,g)$ is a compact Riemannian manifold of dimension $n\ge 2$ we give necessary and sufficient conditions for improved $L^p(M)$-norms of eigenfunctions for all $2<p\ne p_c=\tfrac{2(n+1)}{n-1}$, the critical exponent. Since improved $L^{p_c}(M)$ bounds imply improvement all other exponents, these conditions are necessary for improved bounds for the critical space. We also show that improved $L^{p_c}(M)$ bounds are valid if these conditions are met and if the half-wave operators, $U(t)$, have no caustics when $t\ne 0$. The problem of finding a necessary and sufficient condition for $L^{p_c}(M)$ improvement remains an interesting open problem.


title: A Generalized Correlation Index for Quantifying Signal Morphological Similarity

abstract: In biomedical applications, the similarity between a signal measured from an injured subject and a reference signal measured from a normal subject can be used to quantify the injury severity. This paper proposes a generalization of the adaptive signed correlation index (ASCI) to account for specific signal features of interest and extend the trichotomization of conventional ASCI to an arbitrary number of levels. In the context of spinal cord injury assessment, a computational example is presented to illustrate the enhanced resolution of the proposed measure and its ability to offer a more refined measure of the level of injury.


title: Non-uniform sampling schemes of the Brillouin zone for many-electron perturbation-theory calculations in reduced dimensionality

abstract: First principles calculations based on many-electron perturbation theory methods, such as the \textit{ab initio} GW and GW plus Bethe-Salpeter equation (GW-BSE) approach, are reliable ways to predict quasiparticle and optical properties of materials, respectively. However, these methods involve more care in treating the electron-electron interaction and are considerably more computationally demanding when applied to systems with reduced dimensionality, since the electronic confinement leads a slower convergence of sums over the Brillouin zone due to a much more complicated screening environment that manifests in the "head" and "neck" elements of the dielectric matrix. Here, we present two new schemes to sample the Brillouin zone for GW and GW-BSE calculations: the non-uniform neck subsampling method and the clustered sampling interpolation method, which can respectively be used for a family of single-particle problems, such as GW calculations, and for problems involving the scattering of two-particle states, such as when solving the BSE. We tested these methods on several few-layer semiconductors and graphene and show that they perform a much more efficient sampling of the Brillouin zone and yield two to three orders of magnitude reduction in the computer time. These two methods can be readily incorporated into several \textit{ab initio} packages that compute electronic and optical properties through the GW and GW-BSE approaches.


title: Gas and stellar spiral arms and their offsets in the grand-design spiral galaxy M51

abstract: Theoretical studies on the response of interstellar gas to a gravitational potential disc with a quasi-stationary spiral arm pattern suggest that the gas experiences a sudden compression due to standing shock waves at spiral arms. This mechanism, called a galactic shock wave, predicts that gas spiral arms move from downstream to upstream of stellar arms with increasing radius inside a corotation radius. In order to investigate if this mechanism is at work in the grand-design spiral galaxy M51, we have measured azimuthal offsets between the peaks of stellar mass and gas mass distributions in its two spiral arms. The stellar mass distribution is created by the spatially resolved spectral energy distribution fitting to optical and near infrared images, while the gas mass distribution is obtained by high-resolution CO and HI data. For the inner region (r < 150"), we find that one arm is consistent with the galactic shock while the other is not. For the outer region, results are less certain due to the narrower range of offset values, the weakness of stellar arms, and the smaller number of successful offset measurements. The results suggest that the nature of two inner spiral arms are different, which is likely due to an interaction with the companion galaxy.


title: Theoretical Study Of The Effects Of Magnetic Field Geometry On The High-Energy Emission Of Blazars

abstract: The knowledge of the structure of the magnetic field inside a blazar jet, as deduced from polarization observations at radio to optical wavelengths, is closely related to the formation and propagation of relativistic jets that result from accretion onto supermassive black holes. However, a largely unexplored aspect of the theoretical understanding of radiation transfer physics in blazar jets has been the magnetic field geometry as revealed by the polarized emission and the connection between the variability in polarization and flux across the spectrum. Here, we explore the effects of various magnetic geometries that can exist inside a blazar jet: parallel, oblique, toroidal, and tangled. We investigate the effects of changing the orientation of the magnetic field, according to the above-mentioned geometries, on the resulting high-energy spectral energy distributions (SEDs) and spectral variability patterns (SVPs) of a typical blazar. We use the MUlti-ZOne Radiation Feedback (MUZORF) model of Joshi et al. (2014) to carry out this study and to relate the geometry of the field to the observed SEDs at X-ray and gamma-ray energies. One of the goals of the study is to understand the relationship between synchrotron and inverse Compton peaks in blazar SEDs and the reason for the appearance of gamma-ray "orphan flares" observed in some blazars. This can be associated with the directionality of the magnetic field, which creates a difference in the radiation field as seen by an observer versus that seen by the electrons in the emission region.


title: Minimal Skew energy of oriented bicyclic graphs with a given diameter

abstract: Let $S(G^{\sigma})$ be the skew-adjacency matrix of the oriented graph $G^{\sigma}$, which is obtained from a simple undirected graph $G$ by assigning an orientation $\sigma$ to each of its edges. The skew energy of an oriented graph $G^{\sigma}$ is defined as the sum of absolute values of all eigenvalues of $S(G^{\sigma})$. For any positive integer $d$ with $3\leq d\leq n-3$, we determine the graph with minimal skew energy among all oriented bicyclic graphs that contain no vertex disjoint odd cycle of lengths $s$ and $l$ with $s+l\equiv 2(mod 4)$ on $n$ vertices with a given diameter $d$.


title: The role of phases in detecting three-qubit entanglement

abstract: We propose separability criteria for three-qubit states in terms of diagonal and anti-diagonal entries to detect entanglement with positive partial transposes. We report that the phases, that is, the angular parts of anti-diagonal entries, play a crucial role in determining whether a given three-qubit state is separable or entangled, and they must obey even an identity for separability in some cases. These criteria are strong enough to detect PPT (positive partial transpose) entanglement with nonzero volume. In several cases when all the entries are zero except for diagonal and anti-diagonal entries, we characterize separability using phases. These include the cases when anti-diagonal entries of such states share a common magnitude, and when ranks are less than or equal to six. We also compute the lengths of rank six cases, and find three-qubit separable states with lengths $8$ whose maximum ranks of partial transposes are $7$.


title: The Computational Complexity of Ball Permutations

abstract: Inspired by connections to two dimensional quantum theory, we define several models of computation based on permuting distinguishable particles (which we call balls), and characterize their computational complexity. In the quantum setting, we find that the computational power of this model depends on the initial input states. More precisely, with a standard basis input state, we show how to approximate the amplitudes of this model within additive error using the model DQC1 (the class of problems solvable with one clean qubit), providing evidence that the model in this case is weaker than universal quantum computing. However, for specific choices of input states, the model is shown to be universal for BQP in an encoded sense. We use representation theory of the symmetric group to partially classify the computational complexity of this model for arbitrary input states. Interestingly, we find some input states which yield a model intermediate between DQC1 and BQP. Furthermore, we consider a restricted version of this model based on an integrable scattering problem in 1+1 dimensions. We show it is universal under postselection, if we allow intermediate destructive measurements and specific input states. Therefore, the existence of any classical procedure to sample from the output distribution of this model within multiplicative error implies collapse of polynomial hierarchy to its third level. Finally, we define a classical version of this model in which one can probabilistically permute balls. We find this yields a complexity class which is intermediate between L and BPP. Moreover, we find a nondeterministic version of this model is NP-complete.


title: Photoproduction of Upsilon states in ultraperipheral collisions at the CERN Large Hadron Collider within the color dipole approach

abstract: The exclusive photoproduction of upsilon states $\Upsilon(1S)$ and its radially excited states $\Upsilon(2S,3S)$ is investigated in the context of ultra-peripheral collisions at the LHC energies. Predictions are presented for their production in proton-proton, proton-nucleus and nucleus-nucleus collision at the energies available at the LHC run 2. The rapidity and transverse momentum distributions are shown, and the robustness of the model is tested against the experimental results considering $\psi(1S,2S)$ and $\Upsilon(1S)$ states. The theoretical framework considered in the analysis is the light-cone color dipole formalism, which includes consistently parton saturation effects and nuclear shadowing corrections.


title: KMS states on the operator algebras of reducible higher-rank graphs

abstract: We study the equilibrium or KMS states of the Toeplitz C*-algebra of a finite higher-rank graph which is reducible. The Toeplitz algebra carries a gauge action of a higher-dimensional torus, and a dynamics arises by choosing an embedding of the real numbers in the torus. Here we use an embedding which leads to a dynamics which has previously been identified as "preferred", and we scale the dynamics so that 1 is a critical inverse temperature. As with 1-graphs, we study the strongly connected components of the vertices of the graph. The behaviour of the KMS states depends on both the graphical relationships between the components and the relative size of the spectral radii of the vertex matrices of the components. We test our theorems on graphs with two connected components. We find that our techniques give a complete analysis of the KMS states with inverse temperatures down to a second critical temperature beta_c<1.


title: Electric dipole moments of charged leptons at one loop in the presence of massive neutrinos

abstract: Violation of $CP$ invariance is a quite relevant phenomenon that is found in the Standard Model (SM), though in small amounts. This has been an incentive to look for high-energy descriptions in which $CP$ violation is increased, thus enhancing effects that are suppressed in the SM, such as the electric dipole moments (EDMs) of elementary particles. In the present investigation, we point out that charged currents in which axial couplings are different from vector couplings are able to produce one-loop contributions to EDMs of charged leptons if neutrinos are massive and if these currents violate $CP$. We develop our discussion around charged currents involving heavy neutrinos and a $W'$ gauge boson coupling to SM charged leptons. Using the most stringent bound on the electron EDM, provided by the ACME Collaboration, we determine that the upper bound on the difference between axial and vector currents lies within $\sim10^{-10}$ and $\sim10^{-7}$ for heavy-neutrino masses between $0.5\,{\rm TeV}$ and $6\,{\rm TeV}$ and if the $W'$ mass is within $0.45\,{\rm TeV}-7\,{\rm TeV}$. This possibility is analyzed altogether with the anomalous magnetic moments of charged leptons, among which we estimate, for the $\tau$ lepton, an anomalous magnetic moment contribution between $\sim10^{-8}$ and $\sim10^{-10}$ for neutrino masses ranging from $0.5\,{\rm TeV}$ to $6\,{\rm TeV}$ and a $W'$ mass between $0.45\,{\rm TeV}$ and $7\,{\rm TeV}$. The general charged currents are also used to calculate the branching ratio for $\mu\to e\gamma$, which gets suppressed if the set of masses of heavy neutrinos is quasidegenerate. In a scenario of nondegenerate neutrino masses, we find that regions of neutrino and $W'$ masses in which the contributions to this flavor changing branching ratio are lower than the current upper bound exist. We show that such regions can be widened if the $W'$ gauge boson mass is larger.


title: Parametric Auto-Excitation of Magnetic Droplet Soliton Perimeter Modes

abstract: Recent experiments performed in current-driven nano-contacts with strong perpendicular anisotropy have shown that spin-transfer torque can drive self-localized spin waves [1, 2] that above a certain threshold intensity can condense into a highly nonlinear magnetodynamic and nano-sized state known as a magnetic droplet soliton [3]. Here we demonstrate analytically, numerically, and experimentally that at sufficiently large driving currents, and for a spin polarization that is tilted away from the film normal, the circular droplet soliton can become unstable to periodic excitations of its perimeter. We furthermore show that these perimeter excitation modes (PEMs) are parametrically excited when the fundamental droplet soliton precession frequency is close to twice the frequency of one or more of the PEMs. As a consequence, for increasing applied fields, progressively higher PEMs can be excited. Quantitative agreement with experiment confirms this picture.


title: Chromospheric polarimetry through multi-line observations of the 850 nm spectral region

abstract: Future solar missions and ground-based telescopes aim to understand the magnetism of the solar chromosphere. We performed a supporting study in Quintero Noda et al. (2016) focused on the infrared Ca II 8542 A line and we concluded that is one of the best candidates because it is sensitive to a large range of atmospheric heights, from the photosphere to the middle chromosphere. However, we believe that it is worth to try improving the results produced by this line observing additional spectral lines. In that regard, we examined the neighbour solar spectrum looking for spectral lines that could increase the sensitivity to the atmospheric parameters. Interestingly, we discovered several photospheric lines that greatly improve the photospheric sensitivity to the magnetic field vector. Moreover, they are located close to a second chromospheric line that also belongs to the Ca II infrared triplet, i.e. the Ca II 8498 A line, and enhances the sensitivity to the atmospheric parameters at chromospheric layers. We conclude that the lines in the vicinity of the Ca II 8542 A line not only increase its sensitivity to the atmospheric parameters at all layers, but also they constitute an excellent spectral window for chromospheric polarimetry.


title: Enhancing quantum correlations and entanglement in an optomechanical system via cross-Kerr nonlinearity

abstract: In this work, we theoretically study the quantum correlations present in an optomechanical system by invoking an additional cross-Kerr coupling between the optical and mechanical mode. Under experimentally achievable conditions, we first show that a significant enhancement of the steady-state entanglement could be achieved at a considerably lower driving power, which is also extremely robust with respect to system parameters and environmental temperature. Then, we employ Gaussian quantum discord as a witness of the genuine quantumness of the correlation present in the system and discuss its dependence on the cross-Kerr nonlinearity.


title: Achiral symmetry breaking and positive Gaussian modulus lead to scalloped colloidal membranes

abstract: In the presence of a non-adsorbing polymer, monodisperse rod-like particles assemble into colloidal membranes, which are one rod-length thick liquid-like monolayers of aligned rods. Unlike 3D edgeless bilayer vesicles, colloidal monolayer membranes form open structures with an exposed edge, thus presenting an opportunity to study physics of thin elastic sheets. Membranes assembled from single-component chiral rods form flat disks with uniform edge twist. In comparison, membranes comprised of mixture of rods with opposite chiralities can have the edge twist of either handedness. In this limit disk-shaped membranes become unstable, instead forming structures with scalloped edges, where two adjacent lobes with opposite handedness are separated by a cusp-shaped point defect. Such membranes adopt a 3D configuration, with cusp defects alternatively located above and below the membrane plane. In the achiral regime the cusp defects have repulsive interactions, but away from this limit we measure effective long-ranged attractive binding. A phenomenological model shows that the increase in the edge energy of scalloped membranes is compensated by concomitant decrease in the deformation energy due to Gaussian curvature associated with scalloped edges, demonstrating that colloidal membranes have positive Gaussian modulus. A simple excluded volume argument predicts the sign and magnitude of the Gaussian curvature modulus that is in agreement with experimental measurements. Our results provide insight into how the interplay between membrane elasticity, geometrical frustration and achiral symmetry breaking can be used to fold colloidal membranes into 3D shapes.


title: Partitions of Equiangular Tight Frames

abstract: We present a new efficient algorithm to construct partitions of a special class of equiangular tight frames (ETFs) that satisfy the operator norm bound established by a theorem of Marcus, Spielman, and Srivastava (MSS), which they proved as a corollary yields a positive solution to the Kadison-Singer problem. In particular, we prove that certain diagonal partitions of complex ETFs generated by recursive skew-symmetric conference matrices yield a refinement of the MSS bound. Moreover, we prove that all partitions of ETFs whose largest subset has cardinality three or less also satisfy the MSS bound.


title: Immune Therapeutic Strategies Using Optimal Controls with $L^1$ and $L^2$ Type Objectives

abstract: Therapeutic strategies to correct an excessive immune response to pathogenic infection is investigated as an optimal control problem. The control problem is formulated around a four dimensional mathematical model describing the inflammatory response to a pathogenic insult with two therapeutic control inputs which have either a direct pro- or anti-inflammatory effect in the given system. We use Pontryagin's maximum principle and discuss necessary optimality conditions. We consider both an $L^1$ type objective functional as well as an $L^2$ type objective. For the former, the presence of singular control will be addressed. For each case, numerical simulations using a nonlinear programming optimization solver to acquire different drug treatment strategies are presented and discussed. The results provide insight for possible treatment strategies and the methods could be a relevant tool for future practice to assist in better prediction of clinical outcomes and subsequently better treatment for patients.


title: Single Pass PCA of Matrix Products

abstract: In this paper we present a new algorithm for computing a low rank approximation of the product $A^TB$ by taking only a single pass of the two matrices $A$ and $B$. The straightforward way to do this is to (a) first sketch $A$ and $B$ individually, and then (b) find the top components using PCA on the sketch. Our algorithm in contrast retains additional summary information about $A,B$ (e.g. row and column norms etc.) and uses this additional information to obtain an improved approximation from the sketches. Our main analytical result establishes a comparable spectral norm guarantee to existing two-pass methods; in addition we also provide results from an Apache Spark implementation that shows better computational and statistical performance on real-world and synthetic evaluation datasets.


title: Multispectral Kerr frequency comb initiated by Faraday ripples

abstract: In a uniform microresonator, the generation of a broadband Kerr frequency comb is triggered by Turing patterns. Here, we study a distinctly different route to initiate the Kerr frequency comb by Faraday ripples. Momentum conservation is ensured by azimuthal modulation of the cavity dispersion. With a good agreement with the theoretical analysis, we demonstrate a multispectral Kerr frequency comb covering telecommunication O, C, L, and 2 {\mu}m bands. Comb coherence and absence of a subcomb offset are confirmed by cw heterodyne beat note and amplitude noise spectra measurements. The device can be used for achieving broadband optical frequency synthesizer and high-capacity coherent telecommunication.


title: Public Computation & Boundary Play

abstract: In this paper, we introduce 'public computation' as a genre of learning environments that can be used to radically broaden public participation in authentic, computation-enabled STEM disciplinary practices. Our paradigmatic approach utilizes open source software designed for professional scientists, engineers and digital artists, and situates them in an undiluted form, alongside live and archived expert support, in a public space. We present a case study of DigiPlay, a prototypical public computation space we designed at the University of Calgary, where users can interact directly with scientific simulations as well as the underlying open source code using an array of massive multi- touch screens. We argue that in such a space, public interactions with the code can be thought of as boundary work and play, through which public participation becomes legitimate scientific act, as the public engages in scientific creation through truly open-ended explorations with the code.


title: Analysis of One-Bit Quantized Precoding for the Multiuser Massive MIMO Downlink

abstract: We present a mathematical analysis of linear precoders for downlink massive MIMO multiuser systems that employ one-bit digital-to-analog converters at the basestation in order to reduce complexity and mitigate power usage. The analysis is based on the Bussgang theorem, and applies generally to any linear precoding scheme. We examine in detail the special case of the quantized zero-forcing (ZF) precoder, and derive a simple asymptotic expression for the resulting symbol error rate at each terminal. Our analysis illustrates that the performance of the quantized ZF precoder depends primarily on the ratio of the number of antennas to the number of users, and our simulations show that it can outperform the much more complicated maximum likelihood encoder for low-to-moderate signal to noise ratios, where massive MIMO systems are presumed to operate. We also use the Bussgang theorem to derive a new linear precoder optimized for the case of one-bit quantization, and illustrate its improved performance.


title: Semistochastic Heat-bath Configuration Interaction method: selected configuration interaction with semistochastic perturbation theory

abstract: We extend the recently proposed heat-bath configuration interaction (HCI) method [Holmes, Tubman, Umrigar, J. Chem. Theory Comput. 12, 3674 (2016)], by introducing a semistochastic algorithm for performing multireference Epstein-Nesbet perturbation theory, in order to completely eliminate the severe memory bottleneck of the original method. The proposed algorithm has several attractive features. First, there is no sign problem that plagues several quantum Monte Carlo methods. Second, instead of using Metropolis-Hastings sampling, we use the Alias method to directly sample determinants from the reference wavefunction, thus avoiding correlations between consecutive samples. Third, in addition to removing the memory bottleneck, semistochastic HCI (SHCI) is faster than the deterministic variant for many systems if a stochastic error of 0.1 mHa is acceptable. Fourth, within the SHCI algorithm one can trade memory for a modest increase in computer time. Fifth, the perturbative calculation is embarrassingly parallel. The SHCI algorithm extends the range of applicability of the original algorithm, allowing us to calculate the correlation energy of very large active spaces. We demonstrate this by performing calculations on several first row dimers including F2 with an active space of (14e, 108o), Mn-Salen cluster with an active space of (28e, 22o), and Cr2 dimer with up to a quadruple-zeta basis set with an active space of (12e, 190o). For these systems we were able to obtain better than 1 mHa accuracy with a wall time of merely 55 seconds, 37 seconds, and 56 minutes on 1, 1, and 4 nodes, respectively.


title: Spin transport and dynamics in all-oxide perovskite La$_{2/3}$Sr$_{1/3}$MnO$_3$/SrRuO$_3$ bilayers probed by ferromagnetic resonance

abstract: Thin films of perovskite oxides offer the possibility of combining emerging concepts of strongly correlated electron phenomena and spin current in magnetic devices. However, spin transport and magnetization dynamics in these complex oxide materials are not well understood. Here, we experimentally quantify spin transport parameters and magnetization damping in epitaxial perovskite ferromagnet/paramagnet bilayers of La$_{2/3}$Sr$_{1/3}$MnO$_3$/SrRuO$_3$ (LSMO/SRO) by broadband ferromagnetic resonance spectroscopy. From the SRO thickness dependence of Gilbert damping, we estimate a short spin diffusion length of $\lesssim$1 nm in SRO and an interfacial spin-mixing conductance comparable to other ferromagnet/paramagnetic-metal bilayers. Moreover, we find that anisotropic non-Gilbert damping due to two-magnon scattering also increases with the addition of SRO. Our results demonstrate LSMO/SRO as a spin-source/spin-sink system that may be a foundation for examining spin-current transport in various perovskite heterostructures.


title: Unveiling the Multi-fractal Structure of Complex Networks

abstract: The fractal nature of graphs has traditionally been investigated by using the nodes of networks as the basic units. Here, instead, we propose to concentrate on the graph edges, and introduce a practical and computationally not demanding method for revealing changes in the fractal behavior of networks, and particularly for allowing distinction between mono-fractal, quasi mono-fractal, and multi-fractal structures. We show that degree homogeneity plays a crucial role in determining the fractal nature of the underlying network, and report on six different protein-protein interaction networks along with their corresponding random networks. Our analysis allows to identify varying levels of complexity in the species.


title: On torsion-free nilpotent loops

abstract: We show that a torsion-free nilpotent loop (that is, a loop nilpotent with respect to the dimension filtration) has a torsion-free nilpotent left multiplication group of, at most, the same class. We also prove that a free loop is residually torsion-free nilpotent and that the same holds for any free commutative loop. Although this last result is much stronger than the usual residual nilpotency of the free loop proved by Higman, it is proved, essentially, by the same method.


title: Stochastic Gradient MCMC with Stale Gradients

abstract: Stochastic gradient MCMC (SG-MCMC) has played an important role in large-scale Bayesian learning, with well-developed theoretical convergence properties. In such applications of SG-MCMC, it is becoming increasingly popular to employ distributed systems, where stochastic gradients are computed based on some outdated parameters, yielding what are termed stale gradients. While stale gradients could be directly used in SG-MCMC, their impact on convergence properties has not been well studied. In this paper we develop theory to show that while the bias and MSE of an SG-MCMC algorithm depend on the staleness of stochastic gradients, its estimation variance (relative to the expected estimate, based on a prescribed number of samples) is independent of it. In a simple Bayesian distributed system with SG-MCMC, where stale gradients are computed asynchronously by a set of workers, our theory indicates a linear speedup on the decrease of estimation variance w.r.t. the number of workers. Experiments on synthetic data and deep neural networks validate our theory, demonstrating the effectiveness and scalability of SG-MCMC with stale gradients.


title: On the Convergence of Stochastic Gradient MCMC Algorithms with High-Order Integrators

abstract: Recent advances in Bayesian learning with large-scale data have witnessed emergence of stochastic gradient MCMC algorithms (SG-MCMC), such as stochastic gradient Langevin dynamics (SGLD), stochastic gradient Hamiltonian MCMC (SGHMC), and the stochastic gradient thermostat. While finite-time convergence properties of the SGLD with a 1st-order Euler integrator have recently been studied, corresponding theory for general SG-MCMCs has not been explored. In this paper we consider general SG-MCMCs with high-order integrators, and develop theory to analyze finite-time convergence properties and their asymptotic invariant measures. Our theoretical results show faster convergence rates and more accurate invariant measures for SG-MCMCs with higher-order integrators. For example, with the proposed efficient 2nd-order symmetric splitting integrator, the {\em mean square error} (MSE) of the posterior average for the SGHMC achieves an optimal convergence rate of $L^{-4/5}$ at $L$ iterations, compared to $L^{-2/3}$ for the SGHMC and SGLD with 1st-order Euler integrators. Furthermore, convergence results of decreasing-step-size SG-MCMCs are also developed, with the same convergence rates as their fixed-step-size counterparts for a specific decreasing sequence. Experiments on both synthetic and real datasets verify our theory, and show advantages of the proposed method in two large-scale real applications.


title: Short-term prediction of localized cloud motion using ground-based sky imagers

abstract: Fine-scale short-term cloud motion prediction is needed for several applications, including solar energy generation and satellite communications. In tropical regions such as Singapore, clouds are mostly formed by convection; they are very localized, and evolve quickly. We capture hemispherical images of the sky at regular intervals of time using ground-based cameras. They provide a high resolution and localized cloud images. We use two successive frames to compute optical flow and predict the future location of clouds. We achieve good prediction accuracy for a lead time of up to 5 minutes.


title: Detecting Rainfall Onset Using Sky Images

abstract: Ground-based sky cameras (popularly known as Whole Sky Imagers) are increasingly used now-a-days for continuous monitoring of the atmosphere. These imagers have higher temporal and spatial resolutions compared to conventional satellite images. In this paper, we use ground-based sky cameras to detect the onset of rainfall. These images contain additional information about cloud coverage and movement and are therefore useful for accurate rainfall nowcast. We validate our results using rain gauge measurement recordings and achieve an accuracy of 89% for correct detection of rainfall onset.