title: Growing supermassive black holes in the late stages of galaxy mergers are heavily obscured

abstract: Mergers of galaxies are thought to cause significant gas inflows to the inner parsecs, which can activate rapid accretion onto supermassive black holes (SMBHs), giving rise to Active Galactic Nuclei (AGN). During a significant fraction of this process, SMBHs are predicted to be enshrouded by gas and dust. Studying 52 galactic nuclei in infrared-selected local Luminous and Ultra-luminous infrared galaxies in different merger stages in the hard X-ray band, where radiation is less affected by absorption, we find that the amount of material around SMBHs increases during the last phases of the merger. We find that the fraction of Compton-thick (CT, $N_{\rm\,H}\geq 10^{24}\rm\,cm^{-2}$) AGN in late merger galaxies is higher ($f_{\rm\,CT}=65^{+12}_{-13}\%$) than in local hard X-ray selected AGN ($f_{\rm\,CT}=27\pm 4\%$), and that obscuration reaches its maximum when the nuclei of the two merging galaxies are at a projected distance of $D_{12}\simeq0.4-10.8$ kiloparsecs ($f_{\rm\,CT}=77_{-17}^{+13}\%$). We also find that all AGN of our sample in late merger galaxies have $N_{\rm\,H}> 10^{23}\rm\,cm^{-2}$, which implies that the obscuring material covers $95^{+4}_{-8}\%$ of the X-ray source. These observations show that the material is most effectively funnelled from the galactic scale to the inner tens of parsecs during the late stages of galaxy mergers, and that the close environment of SMBHs in advanced mergers is richer in gas and dust with respect to that of SMBHs in isolated galaxies, and cannot be explained by the classical AGN unification model in which the torus is responsible for the obscuration.


title: Hot Molecular Gas in the Circumnuclear Disk

abstract: We present an analysis of archival ISO observations of pure-rotational lines of H2 in three pointings in the central 3 parsecs of the Galaxy: toward the Southwest region and Northeast region of the Galactic center Circumnuclear Disk, and toward the supermassive black hole Sgr A*. We detect pure rotational lines from 0-0 S(0) to S(13), as well as a number of rovibrationally excited transitions. From the pure rotational lines, we are able to describe the molecular gas with three discrete temperature components: a `hot' component between 500-600 K, a `hotter' component at 1250-1350 K, and a `hottest' component at > 2600 K. Toward Sgr A*, likely due to a combination of poorer baselines and weaker emission, we only detect a single hot component, at 1100 K. The observed excitation is consistent with heating via C-shocks. We also fit a continuous temperature distribution to the the S(1) through S(7) lines by assuming a power-law distribution of temperatures. We measure power law indices of n = 3.22 for the Northeast region and n = 2.83 for the Southwest region, with a smaller index indicating a higher fraction of warm gas. These indices are lower than those measured for other galaxies or other Galactic center clouds, which are measured to have n= 4-6. If we extrapolate this temperature distribution down to a cutoff temperature of 50 K, then the total molecular gas mass that we would measure for the Southwest (Northeast) region is 32 (140) % of the total molecular gas mass inferred from the dust emission, and 26 (125) % of the total molecular gas mass inferred from the CO emission for these regions. Ultimately, we find that the disagreement in the amount of mass recovered for these two regions means that this method cannot yet be verified to yield a reliable and independent estimate of the mass in the local region of the Circumnuclear Disk.


title: Scaling Theory of Entanglement at the Many-Body Localization Transition

abstract: We study the universal properties of eigenstate entanglement entropy across the transition between many-body localized (MBL) and thermal phases. We develop an improved real space renormalization group approach that enables numerical simulation of large system sizes and systematic extrapolation to the infinite system size limit. For systems smaller than the correlation length, the average entanglement follows a sub-thermal volume law, whose coefficient is a universal scaling function. Furthermore, the full distribution of entanglement follows a universal scaling form, and exhibits a bimodal structure that produces universal subleading power-law corrections to the leading volume-law. For systems larger than the correlation length, the short interval entanglement exhibits a discontinuous jump across the transition from fully thermal volume-law on the thermal side, to a pure area-law on the MBL side.


title: Curtailing the Dark Side in Non-Standard Neutrino Interactions

abstract: In presence of non-standard neutrino interactions the neutrino flavor evolution equation is affected by a degeneracy which leads to the so-called LMA-Dark solution. It requires a solar mixing angle in the second octant and implies an ambiguity in the neutrino mass ordering. Non-oscillation experiments are required to break this degeneracy. We perform a combined analysis of data from oscillation experiments with the neutrino scattering experiments CHARM and NuTeV. We find that the degeneracy can be lifted if the non-standard neutrino interactions take place with down quarks, but it remains for up quarks. However, CHARM and NuTeV constraints apply only if the new interactions take place through mediators not much lighter than the electroweak scale. For light mediators we consider the possibility to resolve the degeneracy by using data from future coherent neutrino-nucleus scattering experiments. We find that, for an experiment using a stopped-pion neutrino source, the LMA-Dark degeneracy will either be resolved, or the presence of new interactions in the neutrino sector will be established with high significance.


title: Resolved magnetic structures in the disk-halo interface of NGC 628

abstract: Magnetic fields are essential to fully understand the interstellar medium (ISM) and its role in the disk-halo interface of galaxies is still poorly understood. Star formation is known to expel hot gas vertically into the halo and these outflows have important consequences for mean-field dynamo theory in that they can be efficient in removing magnetic helicity. We perform new observations of the nearby face-on spiral galaxy NGC 628 with the Karl G. Jansky Very Large Array (JVLA) at S-band and the Effelsberg 100-m telescope at frequencies of 2.6 GHz and 8.35 GHz. We obtain some of the most sensitive radio continuum images in both total and linearly polarised intensity of any external galaxy observed so far in addition to high-quality images of Faraday depth and polarisation angle from which we obtained evidence for drivers of magnetic turbulence in the disk-halo connection. Such drivers include a superbubble detected via a significant Faraday depth gradient coinciding with a HI hole. We observe an azimuthal periodic pattern in Faraday depth with a pattern wavelength of 3.7$\pm$ 0.1 kpc, indicating Parker instabilities. The lack of a significant anti-correlation between Faraday depth and magnetic pitch angle indicates that these loops are vertical in nature with little helical twisting, unlike in IC 342. We find that the magnetic pitch angle is systematically larger than the morphological pitch angle of the polarisation arms which gives evidence for the action of a large-scale dynamo where the regular magnetic field is not coupled to the gas flow and obtains a significant radial component. We additionally discover a lone region of ordered magnetic field to the north of the galaxy with a high degree of polarisation and a small pitch angle, a feature that has not been observed in any other galaxy so far and is possibly caused by an asymmetric HI hole.


title: Revisiting the dilatation operator of the Wilson-Fisher fixed point

abstract: We revisit the order $\varepsilon$ dilatation operator of the Wilson-Fisher fixed point obtained by Kehrein, Pismak, and Wegner in light of recent results in conformal field theory. Our approach is algebraic and based only on symmetry principles. The starting point of our analysis is that the first correction to the dilatation operator is a conformal invariant, which implies that its form is fixed up to an infinite set of coefficients associated with the scaling dimensions of higher-spin currents. These coefficients can be fixed using well-known perturbative results, however, they were recently re-obtained using CFT arguments without relying on perturbation theory. Our analysis then implies that all order-$\varepsilon$ scaling dimensions of the Wilson-Fisher fixed point can be fixed by symmetry.


title: Equivalence of restricted Boltzmann machines and tensor network states

abstract: The restricted Boltzmann machine (RBM) is one of the fundamental building blocks of deep learning. RBM finds wide applications in dimensional reduction, feature extraction, and recommender systems via modeling the probability distributions of a variety of input data including natural images, speech signals, and customer ratings, etc. We build a bridge between RBM and tensor network states (TNS) widely used in quantum many-body physics research. We devise efficient algorithms to translate an RBM into the commonly used TNS. Conversely, we give sufficient and necessary conditions to determine whether a TNS can be transformed into an RBM of given architectures. Revealing these general and constructive connections can cross-fertilize both deep learning and quantum many-body physics. Notably, by exploiting the entanglement entropy bound of TNS, we can rigorously quantify the expressive power of RBM on complex data sets. Insights into TNS and its entanglement capacity can guide the design of more powerful deep learning architectures. On the other hand, RBM can represent quantum many-body states with fewer parameters compared to TNS, which may allow more efficient classical simulations.


title: Gaussian-based coupled-cluster theory for the ground state and band structure of solids

abstract: We present the results of Gaussian-based ground-state and excited-state equation-of-motion coupled-cluster theory with single and double excitations for three-dimensional solids. We focus on diamond and silicon, which are paradigmatic covalent semiconductors. In addition to ground-state properties (the lattice constant, bulk modulus, and cohesive energy), we compute the quasiparticle band structure and band gap. We sample the Brillouin zone with up to 64 k-points using norm-conserving pseudopotentials and polarized double- and triple-zeta basis sets, leading to canonical coupled-cluster calculations with as many as 256 electrons in 2,176 orbitals.


title: How to break the density-anisotropy degeneracy in spherical stellar systems

abstract: We present a new non-parametric Jeans code, GravSphere, that recovers the density $\rho(r)$ and velocity anisotropy $\beta(r)$ of spherical stellar systems, assuming only that they are in a steady-state. Using a large suite of mock data, we confirm that with only line-of-sight velocity data, GravSphere provides a good estimate of the density at the projected stellar half mass radius, $\rho(R_{1/2})$, but is not able to measure $\rho(r)$ or $\beta(r)$, even with 10,000 tracer stars. We then test three popular methods for breaking this $\rho-\beta$ degeneracy: using multiple populations with different $R_{1/2}$; using higher order `Virial Shape Parameters' (VSPs); and including proper motion data. We find that two populations provide an excellent recovery of $\rho(r)$ in-between their respective $R_{1/2}$. However, even with a total of $\sim 7,000$ tracers, we are not able to well-constrain $\beta(r)$ for either population. By contrast, using 1000 tracers with higher order VSPs we are able to measure $\rho(r)$ over the range $0.5 < r/R_{1/2} < 2$ and broadly constrain $\beta(r)$. Including proper motion data for all stars gives an even better performance, with $\rho$ and $\beta$ well-measured over the range $0.25 < r/R_{1/2} < 4$. Finally, we test GravSphere on a triaxial mock galaxy that has axis ratios typical of a merger remnant, $[1:0.8:0.6]$. In this case, GravSphere can become slightly biased. However, we find that when this occurs the data are poorly fit, allowing us to detect when such departures from spherical symmetry become problematic.


title: High-order schemes for the Euler equations in cylindrical/spherical coordinates

abstract: We consider implementations of high-order finite difference Weighted Essentially Non-Oscillatory (WENO) schemes for the Euler equations in cylindrical and spherical coordinate systems with radial dependence only. The main concern of this work lies in ensuring both high-order accuracy and conservation. Three different spatial discretizations are assessed: one that is shown to be high-order accurate but not conservative, one conservative but not high-order accurate, and a new approach that is both high-order accurate and conservative. For cylindrical and spherical coordinates, we present convergence results for the advection equation and the Euler equations with an acoustics problem; we then use the Sod shock tube and the Sedov point-blast problems in cylindrical coordinates to verify our analysis and implementations.

categories: math.NA math-ph math.MP physics.comp-ph physics.flu-dyn


title: The SLUGGS Survey: a catalog of over 4000 globular cluster radial velocities in 27 nearby early-type galaxies

abstract: Here we present positions and radial velocities for over 4000 globular clusters (GCs) in 27 nearby early-type galaxies from the SLUGGS survey. The SLUGGS survey is designed to be representative of elliptical and lenticular galaxies in the stellar mass range 10 $<$ log M$_{\ast}$/M$_{\odot}$ $<$ 11.7. The data have been obtained over many years, mostly using the very stable multi-object spectrograph DEIMOS on the Keck II 10m telescope. Radial velocities are measured using the calcium triplet lines with a velocity accuracy of $\pm$ 10-15 km/s. We use phase space diagrams (i.e. velocity--position diagrams) to identify contaminants such as foreground stars and background galaxies, and to show that the contribution of GCs from neighboring galaxies is generally insignificant. Likely ultra-compact dwarfs are tabulated separately. We find that the mean velocity of the GC system is close to that of the host galaxy systemic velocity, indicating that the GC system is in overall dynamical equilibrium within the galaxy potential. We also find that the GC system velocity dispersion scales with host galaxy stellar mass in a similar manner to the Faber-Jackson relation for the stellar velocity dispersion. Publication of these GC radial velocity catalogs should enable further studies in many areas, such as GC system substructure, kinematics, and host galaxy mass measurements.


title: Evaluating (and Improving) Estimates of the Solar Radial Magnetic Field Component from Line-of-Sight Magnetograms

abstract: Although for many solar physics problems the desirable or meaningful boundary is the radial component of the magnetic field $B_{\rm r}$, the most readily available measurement is the component of the magnetic field along the line-of-sight to the observer, $B_{\rm los}$. As this component is only equal to the radial component where the viewing angle is exactly zero, some approximation is required to estimate $B_{\rm r}$ at all other observed locations. In this study, a common approximation known as the "$\mu$-correction", which assumes all photospheric field to be radial, is compared to a method which invokes computing a potential field that matches the observed $B_{\rm los}$, from which the potential field radial component, $B_{\rm r}^{\rm pot}$ is recovered. We demonstrate that in regions that are truly dominated by radially-oriented field at the resolution of the data employed, the $\mu$-correction performs acceptably if not better than the potential-field approach. However, it is also shown that for any solar structure which includes horizontal fields, i.e. active regions, the potential-field method better recovers both the strength of the radial field and the location of magnetic neutral line.


title: A characterization of the Radon-Nikodym property for vector valued measures

abstract: If $\mu_1,\mu_2,\dots$ are positive measures on a measurable space $(X,\Sigma)$ and $v_1,v_2, \dots$ are elements of a Banach space ${\mathbb E}$ such that $\sum_{n=1}^\infty \|v_n\| \mu_n(X) < \infty$, then $\omega (S)= \sum_{n=1}^\infty v_n \mu_n(S)$ defines a vector measure of bounded variation on $(X,\Sigma)$. We show ${\mathbb E}$ has the Radon-Nikodym property if and only if every ${\mathbb E}$-valued measure of bounded variation on $(X,\Sigma)$ is of this form. As an application of this result we show that under natural conditions an operator defined on positive measures, has a unique extension to an operator defined on ${\mathbb E}$-valued measures for any Banach space ${\mathbb E}$ that has the Radon-Nikodym property.


title: Pseudogaps in strongly interacting Fermi gases

abstract: A central challenge in modern condensed matter physics is developing the tools for understanding nontrivial yet unordered states of matter. One important idea to emerge in this context is that of a "pseudogap": the fact that under appropriate circumstances the normal state displays a suppression of the single particle spectral density near the Fermi level, reminiscent of the gaps seen in ordered states of matter. While these concepts arose in a solid state context, it is now being explored in cold gases. This article reviews the current experimental and theoretical understanding of the normal state of strongly interacting Fermi gases, with particular focus on the phenomonology which is traditionally associated with the pseudogap.


title: A non-Archimedean Ohsawa-Takegoshi extension theorem

abstract: We prove an Ohsawa-Takegoshi-type extension theorem on the Berkovich closed unit disc over a complete non-Archimedean field. As an application, we establish a non-Archimedean analogue of Demailly's regularization theorem for quasisubharmonic functions on the Berkovich unit disc.


title: Optical patterning of magnetic domains and defects in ferromagnetic liquid crystal colloids

abstract: A promising approach in designing composite materials with unusual physical behavior combines solid nanostructures and orientationally ordered soft matter at the mesoscale. Such composites not only inherit properties of their constituents but also can exhibit emergent behavior, such as ferromagnetic ordering of colloidal metal nanoparticles forming mesoscopic magnetization domains when dispersed in a nematic liquid crystal. Here we demonstrate the optical patterning of domain structures and topological defects in such ferromagnetic liquid crystal colloids, which allows for altering their response to magnetic fields. Our findings reveal the nature of the defects in this soft matter system which is different as compared to non-polar nematics and ferromagnets alike.


title: Entropy dissipation of Fokker-Planck equations on graphs

abstract: We study the nonlinear Fokker-Planck equation on graphs, which is the gradient flow in the space of probability measures supported on the nodes with respect to the discrete Wasserstein metric. The energy functional driving the gradient flow consists of a Boltzmann entropy, a linear potential and a quadratic interaction energy. We show that the solution converges to the Gibbs measures exponentially fast with a rate that can be given analytically. The continuous analog of this asymptotic rate is related to the Yano's formula.


title: CSI 2264: Simultaneous optical and X-ray variability in pre-Main Sequence stars. I: Time resolved X-ray spectral analysis during optical dips and accretion bursts in stars with disks

abstract: Pre-main sequence stars are variable sources. In stars with disks, this variability is related to the morphology of the inner circumstellar region (<0.1 AU) and that of the photosphere and corona, all impossible to be spatially resolved with present day techniques. This has been the main motivation for the Coordinated Synoptic Investigation of NGC 2264. In this paper, we focus on the stars with disks. We analyze the X-ray spectral properties extracted during optical bursts and dips in order to unveil the nature of these phenomena. We analyze simultaneous CoRoT and Chandra/ACIS-I observations to search for coherent optical and X-ray flux variability in stars with disks. Then, stars are analyzed in two different samples. In stars with variable extinction, we look for a simultaneous increase of optical extinction and X-ray absorption during the optical dips; in stars with accretion bursts, we search for soft X-ray emission and increasing X-ray absorption during the bursts. Results. We find evidence for coherent optical and X-ray flux variability among the stars with variable extinction. In 9/24 stars with optical dips, we observe a simultaneous increase of X-ray absorption and optical extinction. In seven dips, it is possible to calculate the NH/AV ratio in order to infer the composition of the obscuring material. In 5/20 stars with optical accretion bursts, we observe increasing soft X-ray emission during the bursts that we associate to the emission of accreting gas. It is not surprising that these properties are not observed in all the stars with dips and bursts, since favorable geometric configurations are required. The observed variable absorption during the dips is mainly due to dust-free material in accretion streams. In stars with accretion bursts, we observe on average a larger soft X-ray spectral component not observed in non accreting stars.


title: Novel mid-infrared dispersive wave generation in gas-filled PCF by transient ionization-driven changes in dispersion

abstract: Gas-filled hollow-core photonic crystal fibre (PCF) is being used to generate ever wider supercontinuum spectra, in particular via dispersive wave (DW) emission in the deep and vacuum ultraviolet, with a multitude of applications. DWs are the result of the resonant transfer of energy from a self-compressed soliton, a process which relies crucially on phase matching. It was recently predicted that, in the strong-field regime, the additional transient anomalous dispersion introduced by gas ionization would allow phase-matched DW generation in the mid-infrared (MIR)-something that is forbidden in the absence of free electrons. Here we report for the first time the experimental observation of such MIR DWs, embedded in a 4.7-octave-wide supercontinuum that uniquely reaches simultaneously to the vacuum ultraviolet, with up to 1.7 W of total average power.


title: Quantum Entanglement in Neural Network States

abstract: Machine learning, one of today's most rapidly growing interdisciplinary fields, promises an unprecedented perspective for solving intricate quantum many-body problems. Understanding the physical aspects of the representative artificial neural-network states is recently becoming highly desirable in the applications of machine learning techniques to quantum many-body physics. Here, we study the quantum entanglement properties of neural-network states, with a focus on the restricted-Boltzmann-machine (RBM) architecture. We prove that the entanglement of all short-range RBM states satisfies an area law for arbitrary dimensions and bipartition geometry. For long-range RBM states we show by using an exact construction that such states could exhibit volume-law entanglement, implying a notable capability of RBM in representing efficiently quantum states with massive entanglement. We further examine generic RBM states with random weight parameters. We find that their averaged entanglement entropy obeys volume-law scaling and meantime strongly deviates from the Page-entropy of the completely random pure states. We show that their entanglement spectrum has no universal part associated with random matrix theory and bears a Poisson-type level statistics. Using reinforcement learning, we demonstrate that RBM is capable of finding the ground state (with power-law entanglement) of a model Hamiltonian with long-range interaction. In addition, we show, through a concrete example of the one-dimensional symmetry-protected topological cluster states, that the RBM representation may also be used as a tool to analytically compute the entanglement spectrum. Our results uncover the unparalleled power of artificial neural networks in representing quantum many-body states, which paves a novel way to bridge computer science based machine learning techniques to outstanding quantum condensed matter physics problems.


title: Galactic Dark Matter Halos and Globular Cluster Populations. III: Extension to Extreme Environments

abstract: The total mass M_GCS in the globular cluster (GC) system of a galaxy is empirically a near-constant fraction of the total mass M_h = M_bary + M_dark of the galaxy, across a range of 10^5 in galaxy mass. This trend is radically unlike the strongly nonlinear behavior of total stellar mass M_star versus M_h. We discuss extensions of this trend to two more extreme situations: (a) entire clusters of galaxies, and (b) the Ultra-Diffuse Galaxies (UDGs) recently discovered in Coma and elsewhere. Our calibration of the ratio \eta_M = M_GCS / M_h from normal galaxies, accounting for new revisions in the adopted mass-to-light ratio for GCs, now gives \eta_M = 2.9 \times 10^{-5} as the mean absolute mass fraction. We find that the same ratio appears valid for galaxy clusters and UDGs. Estimates of \eta_M in the four clusters we examine tend to be slightly higher than for individual galaxies, butmore data and better constraints on the mean GC mass in such systems are needed to determine if this difference is significant. We use the constancy of \eta_M to estimate total masses for several individual cases; for example, the total mass of the Milky Way is calculated to be M_h = 1.1 \times 10^{12} M_sun. Physical explanations for the uniformity of \eta_M are still descriptive, but point to a picture in which massive, dense star clusters in their formation stages were relatively immune to the feedback that more strongly influenced lower-density regions where most stars form.


title: Beyond Whittle: Nonparametric correction of a parametric likelihood with a focus on Bayesian time series analysis

abstract: The Whittle likelihood is widely used for Bayesian nonparametric estimation of the spectral density of stationary time series. However, the loss of efficiency for non-Gaussian time series can be substantial. On the other hand, parametric methods are more powerful if the model is well-specified, but may fail entirely otherwise. Therefore, we suggest a nonparametric correction of a parametric likelihood taking advantage of the efficiency of parametric models while mitigating sensitivities through a nonparametric amendment. Using a Bernstein-Dirichlet prior for the nonparametric spectral correction, we show posterior consistency and illustrate the performance of our procedure in a simulation study and with LIGO gravitational wave data.


title: Atomic Data Revisions for Transitions Relevant to Observations of Interstellar, Circumgalactic, and Intergalactic Matter

abstract: Measurements of element abundances in galaxies from astrophysical spectroscopy depend sensitively on the atomic data used. With the goal of making the latest atomic data accessible to the community, we present a compilation of selected atomic data for resonant absorption lines at wavelengths longward of 911.753 {\AA} (the \ion{H}{1} Lyman limit), for key heavy elements (heavier than atomic number 5) of astrophysical interest. In particular, we focus on the transitions of those ions that have been observed in the Milky Way interstellar medium (ISM), the circumgalactic medium (CGM) of the Milky Way and/or other galaxies, and the intergalactic medium (IGM). We provide wavelengths, oscillator strengths, associated accuracy grades, and references to the oscillator strength determinations. We also attempt to compare and assess the recent oscillator strength determinations. For about 22\% of the lines that have updated oscillator strength values, the differences between the former values and the updated ones are $\gtrsim$~0.1 dex. Our compilation will be a useful resource for absorption line studies of the ISM, as well as studies of the CGM and IGM traced by sight lines to quasars and gamma-ray bursts. Studies (including those enabled by future generations of extremely large telescopes) of absorption by galaxies against the light of background galaxies will also benefit from our compilation.


title: On a conjecture of Demailly and new bounds on Waldschmidt constants in ${\mathbb P}^N$

abstract: In the present note we prove a conjecture of Demailly for finite sets of sufficiently many very general points in projective spaces. This gives a lower bound on Waldschmidt constants of such sets. Waldschmidt constants are asymptotic invariants of subschemes receiving recently considerable amount of attention.


title: Optical manipulation of shape-morphing elastomeric liquid crystal microparticles doped with gold nanocrystals

abstract: We demonstrate facile optical manipulation of shape of birefringent colloidal microparticles made from liquid crystal elastomers. Using soft lithography and polymerization, we fabricate elastomeric microcylinders with weakly undulating director oriented on average along their long axes. These particles are infiltrated with gold nanospheres acting as heat transducers that allow for an efficient localized transfer of heat from a focused infrared laser beam to a submicrometer region within a microparticle. Photothermal control of ordering in the liquid crystal elastomer using scanned beams allows for a robust control of colloidal particles, enabling both reversible and irreversible changes of shape. Possible applications include optomechanics, microfluidics, and reconfigurable colloidal composites with shape-dependent self-assembly.


title: Selection of quasi-stationary states in the Navier-Stokes equation on the torus

abstract: The two dimensional incompressible Navier-Stokes equation on $D_\delta := [0, 2\pi\delta] \times [0, 2\pi]$ with $\delta \approx 1$, periodic boundary conditions, and viscosity $0 < \nu \ll 1$ is considered. Bars and dipoles, two explicitly given quasi-stationary states of the system, evolve on the time scale $\mathcal{O}(e^{-\nu t})$ and have been shown to play a key role in its long-time evolution. Of particular interest is the role that $\delta$ plays in selecting which of these two states is observed. Recent numerical studies suggest that, after a transient period of rapid decay of the high Fourier modes, the bar state will be selected if $\delta \neq 1$, while the dipole will be selected if $\delta = 1$. Our results support this claim and seek to mathematically formalize it. We consider the system in Fourier space, project it onto a center manifold consisting of the lowest eight Fourier modes, and use this as a model to study the selection of bars and dipoles. It is shown for this ODE model that the value of $\delta$ controls the behavior of the asymptotic ratio of the low modes, thus determining the likelihood of observing a bar state or dipole after an initial transient period. Moreover, in our model, for all $\delta \approx 1$, there is an initial time period in which the high modes decay at the rapid rate $\mathcal{O}(e^{-t/\nu})$, while the low modes evolve at the slower $\mathcal{O}(e^{-\nu t})$ rate. The results for the ODE model are proven using energy estimates and invariant manifolds and further supported by formal asymptotic expansions and numerics.


title: Synthesizing Normalized Faces from Facial Identity Features

abstract: We present a method for synthesizing a frontal, neutral-expression image of a person's face given an input face photograph. This is achieved by learning to generate facial landmarks and textures from features extracted from a facial-recognition network. Unlike previous approaches, our encoding feature vector is largely invariant to lighting, pose, and facial expression. Exploiting this invariance, we train our decoder network using only frontal, neutral-expression photographs. Since these photographs are well aligned, we can decompose them into a sparse set of landmark points and aligned texture maps. The decoder then predicts landmarks and textures independently and combines them using a differentiable image warping operation. The resulting images can be used for a number of applications, such as analyzing facial attributes, exposure and white balance adjustment, or creating a 3-D avatar.


title: Eight new luminous z > 6 quasars selected via SED model fitting of VISTA, WISE and Dark Energy Survey Year 1 Observations

abstract: We present the discovery and spectroscopic confirmation with the ESO NTT and Gemini South telescopes of eight new 6.0 < z < 6.5 quasars with z$_{AB}$ < 21.0. These quasars were photometrically selected without any star-galaxy morphological criteria from 1533 deg$^{2}$ using SED model fitting to photometric data from the Dark Energy Survey (g, r, i, z, Y), the VISTA Hemisphere Survey (J, H, K) and the Wide-Field Infrared Survey Explorer (W1, W2). The photometric data was fitted with a grid of quasar model SEDs with redshift dependent Lyman-{\alpha} forest absorption and a range of intrinsic reddening as well as a series of low mass cool star models. Candidates were ranked using on a SED-model based $\chi^{2}$-statistic, which is extendable to other future imaging surveys (e.g. LSST, Euclid). Our spectral confirmation success rate is 100% without the need for follow-up photometric observations as used in other studies of this type. Combined with automatic removal of the main types of non-astrophysical contaminants the method allows large data sets to be processed without human intervention and without being over run by spurious false candidates. We also present a robust parametric redshift estimating technique that gives comparable accuracy to MgII and CO based redshift estimators. We find two z $\sim$ 6.2 quasars with HII near zone sizes < 3 proper Mpc which could indicate that these quasars may be young with ages < 10$^6$ - 10$^7$ years or lie in over dense regions of the IGM. The z = 6.5 quasar VDESJ0224-4711 has J$_{AB}$ = 19.75 is the second most luminous quasar known with z > 6.5.


title: Continuous diffusion model for concentration dependence of nitroxide EPR parameters in normal and supercooled water

abstract: We measured electron paramagnetic resonance (EPR) spectra of 14N- and 15N-labeled perdeuterated TEMPONE radicals in normal and supercooled water at various radical concentrations. By fitting the EPR spectra to spectral shape functions based on the modified Bloch equations, we obtained concentration dependences of EPR parameters of radicals at each measured temperature. From concentration dependences of the EPR parameters quantifying spin dephasing, coherence transfer, and hyperfine splitting, we determined linear concentration coefficients, whose values depend on the relative motion of radicals due to modulation of the Heisenberg spin exchange (HSE) and dipole-dipole (DD) interactions between them. We applied the continuous diffusion model for relative motion of radicals and we evaluated the diffusion coefficients of radicals from the concentration coefficients using the standard relations and the relations derived from kinetic equations for the spin evolution of interacting radical pair. It was found that the latter equations lead to the better agreement between the diffusion coefficients calculated from different concentration coefficients. The calculated diffusion coefficients of 14N- and 15N-labeled radicals show similar values, which is an expected result that supports the presented method. Upon lowering the temperature into the supercooled state, the calculated diffusion coefficients decrease slower than is predicted by the Stokes-Einstein relation and slower than the rotational diffusion coefficient. Similar effects were detected in NMR studies of the rotational and translational motion of water molecules in supercooled water.


title: Conservation laws and symmetries of a generalized Kawahara equation

abstract: The generalized Kawahara equation $u_t=a(t) u_{xxxxx} +b(t)u_{xxx} +c(t)f(u) u_x$ appears in many physical applications. A complete classification of low-order conservation laws and point symmetries is obtained for this equation, which includes as a special case the usual Kawahara equation $u_t = \alpha u u_x+\beta u^2u_x +\gamma u_{xxx}+\mu u_{xxxxx}$. A general connection between conservation laws and symmetries for the generalized Kawahara equation is derived through the Hamiltonian structure of this equation and its relationship to Noether's theorem using a potential formulation.


title: Some Connections Between Cycles and Permutations that Fix a Set and Touchard Polynomials and Covers of Multisets

abstract: We present a new proof of a fundamental result concerning cycles of random permutations which gives some intuition for the connection between Touchard polynomials and the Poisson distribution. We also introduce a rather novel permutation statistic and study its distribution. This quantity, indexed by $m$, is the number of sets of size $m$ fixed by the permutation. This leads to a new and simpler derivation of the exponential generating function for the number of covers of certain multisets.


title: On integrable field theories as dihedral affine Gaudin models

abstract: We introduce the notion of a classical dihedral affine Gaudin model, associated with an untwisted affine Kac-Moody algebra $\widetilde{\mathfrak{g}}$ equipped with an action of the dihedral group $D_{2T}$, $T \geq 1$ through (anti-)linear automorphisms. We show that a very broad family of classical integrable field theories can be recast as examples of such classical dihedral affine Gaudin models. Among these are the principal chiral model on an arbitrary real Lie group $G_0$ and the $\mathbb{Z}_T$-graded coset $\sigma$-model on any coset of $G_0$ defined in terms of an order $T$ automorphism of its complexification. Most of the multi-parameter integrable deformations of these $\sigma$-models recently constructed in the literature provide further examples. The common feature shared by all these integrable field theories, which makes it possible to reformulate them as classical dihedral affine Gaudin models, is the fact that they are non-ultralocal. In particular, we also obtain affine Toda field theory in its lesser-known non-ultralocal formulation as another example of this construction. We propose that the interpretation of a given classical non-ultralocal integrable field theory as a classical dihedral affine Gaudin model provides a natural setting within which to address its quantisation. At the same time, it may also furnish a general framework for understanding the massive ODE/IM correspondence since the known examples of integrable field theories for which such a correspondence has been formulated can all be viewed as dihedral affine Gaudin models.


title: Measurement of the antineutrino to neutrino charged-current interaction cross section ratio in MINERvA

abstract: We present measurements of the neutrino and antineutrino total charged-current cross sections on carbon and their ratio using the MINERvA scintillator-tracker. The measurements span the energy range 2-22 GeV and were performed using forward and reversed horn focusing modes of the Fermilab low-energy NuMI beam to obtain large neutrino and antineutrino samples. The flux is obtained using a sub-sample of charged-current events at low hadronic energy transfer along with precise higher energy external neutrino cross section data overlapping with our energy range between 12-22 GeV. We also report on the antineutrino-neutrino cross section ratio, Rcc, which does not rely on external normalization information. Our ratio measurement, obtained within the same experiment using the same technique, benefits from the cancellation of common sample systematic uncertainties and reaches a precision of 5% at low energy. Our results for the antineutrino-nucleus scattering cross section and for Rcc are the most precise to date in the energy range $E_{\nu} <$ 6GeV.


title: Mixed Effects Models are Sometimes Terrible

abstract: Mixed-effects models have emerged as the gold standard of statistical analysis in different sub-fields of linguistics (Baayen, Davidson & Bates, 2008; Johnson, 2009; Barr, et al, 2013; Gries, 2015). One problematic feature of these models is their failure to converge under maximal (or even near-maximal) random effects structures. The lack of convergence is relatively unaddressed in linguistics and when it is addressed has resulted in statistical practices (e.g. Jaeger, 2009; Gries, 2015; Bates, et al, 2015b) that are premised on the idea that non-convergence is an indication that a random effects structure is over-specified (or not parsimonious), the parsimonious convergence hypothesis (PCH). We test the PCH by running simulations in lme4 under two sets of assumptions for both a linear dependent variable and a binary dependent variable in order to assess the rate of non-convergence for both types of mixed effects models when a known maximal effect structure is used to generate the data (i.e. when non-convergence cannot be explained by random effects with zero variance). Under the PCH, lack of convergence is treated as evidence against a more maximal random effects structure, but that result is not upheld with our simulations. We provide an alternative model, fully specified Bayesian models implemented in rstan (Stan Development Team, 2016; Carpenter, et al, in press) that removed the convergence problems almost entirely in simulations of the same conditions. These results indicate that when there is known non-zero variance for all slopes and intercepts, under realistic distributions of data and with moderate to severe imbalance, mixed effects models in lme4 have moderate to high non-convergence rates which can cause linguistic researchers to wrongfully exclude random effect terms.


title: A proper motions study of the globular cluster NGC 3201

abstract: With a high value of heliocentric radial velocity, a retrograde orbit, and being suspected to have an extragalactic origin, NGC 3201 is an interesting globular cluster for kinematical studies. Our purpose is to calculate the relative proper motions (PMs) and membership probability for the stars in the wide region of globular cluster NGC 3201. Proper motion based membership probabilities are used to isolate the cluster sample from the field stars. The membership catalogue will help address the question of chemical inhomogeneity in the cluster. Archive CCD data taken with a wide-field imager (WFI) mounted on the ESO 2.2m telescope are reduced using the high-precision astrometric software developed by Anderson et al. for the WFI images. The epoch gap between the two observational runs is $\sim$14.3 years. To standardize the $BVI$ photometry, Stetson's secondary standard stars are used. The CCD data with an epoch gap of $\sim$14.3 years enables us to decontaminate the cluster stars from field stars efficiently. The median precision of PMs is better than $\sim$0.8 mas~yr$^{-1}$ for stars having $V<$18 mag that increases up to $\sim$1.5 mas~yr$^{-1}$ for stars with $18<V<20$ mag. Kinematic membership probabilities are calculated using proper motions for stars brighter than $V\sim$20 mag. An electronic catalogue of positions, relative PMs, $BVI$ magnitudes and membership probabilities in $\sim$19.7$\times$17 arcmin$^2$ region of NGC 3201 is presented. We use our membership catalogue to identify probable cluster members among the known variables and $X$-ray sources in the direction of NGC 3201.


title: A Necessary and Suffcient Condition for Uniqueness of the Trivial Solution in Semilinear Parabolic Equations

abstract: In their 1968 paper Fujita and Watanabe considered the issue of uniqueness of the trivial solution of semilinear parabolic equations with respect to the class of bounded, non-negative solutions. In particular they showed that if the underlying ODE has non-unique solutions (as characterised via an Osgood-type condition) {\em and} the nonlinearity $f$ satisfies a concavity condition, then the parabolic PDE also inherits the non-uniqueness property. This concavity assumption has remained in place either implicitly or explicitly in all subsequent work in the literature relating to this and other, similar, non-uniqueness phenomena in parabolic equations. In this paper we provide an elementary proof of non-uniqueness for the PDE without any such concavity assumption on $f$. An important consequence of our result is that uniqueness of the trivial solution of the PDE is equivalent to uniqueness of the trivial solution of the corresponding ODE, which in turn is known to be equivalent to an Osgood-type integral condition on $f$.


title: Dynamical coupling between protein conformational fluctuation and hydration water: Heterogeneous dynamics of biological water

abstract: We investigate dynamical coupling between water and amino acid side-chain residues in solvation dynamics by selecting residues often used as natural probes, namely tryptophan, tyrosine and histidine, located at different positions on protein surface and having various degrees of solvent exposure. Such differently placed residues are found to exhibit different timescales of relaxation. The total solvation response, as measured by the probe is decomposed in terms of its interactions with (i) protein core, (ii) side-chain atoms and (iii) water molecules. Significant anti cross-correlations among these contributions are observed as a result of side-chain assisted energy flow between protein core and hydration layer, which is important for the proper functionality of a protein. It is also observed that there are rotationally faster as well as slower water molecules than that of bulk solvent, which are considered to be responsible for the multitude of timescales that are observed in solvation dynamics. We also establish that slow solvation derives a significant contribution from protein side-chain fluctuations. When the motion of the protein side-chains is forcefully quenched, solvation either becomes faster or slower depending on the location of the probe.


title: Towards Principled Methods for Training Generative Adversarial Networks

abstract: The goal of this paper is not to introduce a single algorithm or method, but to make theoretical steps towards fully understanding the training dynamics of generative adversarial networks. In order to substantiate our theoretical analysis, we perform targeted experiments to verify our assumptions, illustrate our claims, and quantify the phenomena. This paper is divided into three sections. The first section introduces the problem at hand. The second section is dedicated to studying and proving rigorously the problems including instability and saturation that arize when training generative adversarial networks. The third section examines a practical and theoretically grounded direction towards solving these problems, while introducing new tools to study them.


title: Natural Encoding of Information through Interacting Impulses

abstract: How to artificially encode observer in universal information coding structure like DNA ? It requires naturally creating information Bits and natural encoding triplet code enables recognizing other encoded information. These Bits become standard units of different information languages in modern communications. Fundamental interactions build structure of Universe. Numerous multilevel inter-species interactions selforganize biosystems. Human interactions unify these and many others. Physical reality is only interactions identified or not yet. Each interaction is elementary yes-no action of impulse which models a natural Bit. Natural interactive process, transferring Bits, models information process. Information is universal physical substance a phenomenon of interaction which not only originates information but transfers it sequentially. Mutually interacting processes enable creating new elements like chemical chain reactions. The elements enclosing components of reaction memorize the interactive yes-no result similar to encoding. Energy quantity and quality of specific interaction determine sequence of transferring information, its encoding, and limit the code length. The introduced formalism of natural emergence information and its encoding also shows advantage over non-natural encoding. The impulse sequential natural encoding merges memory with the time of memorizing information and compensates the cost by running time intervals of encoding. Information process binds the encoding impulse reversible microprocesses in multiple impulses macroprocess of information irreversible dynamics establishing interactive integrated information dynamics. The encoding process integrates geometrical triplet coding structure rotating double helix of sequencing cells Bits, which commands cognition, intelligence including conscience. The results validate computer simulation, and experiments.


title: Shear banding in metallic glasses described by alignments of Eshelby quadrupoles

abstract: Plastic deformation of metallic glasses performed well below the glass transition temperature leads to the formation of shear bands as a result of shear localization. It is believed that shear banding originates from individual stress concentrators having quadrupolar symmetry. To elucidate the underlying mechanisms of shear band formation, microstructural investigations were carried out on sheared zones using transmission electron microscopy. Here we show evidence of a characteristic signature present in shear bands manifested in the form of sinusoidal density variations. We present an analytical solution for the observed post-deformation state derived from continuum mechanics using an alignment of quadrupolar stress field perturbations for the plastic events. Since we observe qualitatively similar features for three different types of metallic glasses that span the entire range of characteristic properties of metallic glasses, we conclude that the reported deformation behavior is generic for all metallic glasses, and thus has far-reaching consequences for the deformation behavior of amorphous solids in general.


title: Limits on the Lorentz Invariance Violation from UHECR astrophysics

abstract: In this paper, Lorentz Invariance Violation (LIV) is introduced in the calculations of photon propagation in the Universe. LIV is considered in the photon sector and the mean free path of the $\gamma \gamma \rightarrow e^{+} e^{-}$ interaction is calculated. The corresponding photon horizon including LIV effects is used to predict major changes in the propagation of photons with energy above $10^{18}$ eV. The flux of GZK photons on Earth considering LIV is calculated for several source models of ultra-high energy cosmic ray (UHECR). The predicted flux of GZK gamma-rays is compared to the new upper limits on the photon flux obtained by the Pierre Auger Observatory in order to impose upper limits on the LIV coefficients of order $n =$ 0, 1 and 2. The limits on the LIV coefficients derived here are more realistic than previous works and in some cases more restrictive. The analysis resulted in LIV upper limits in the photon sector of $\delta_{\gamma,0}^{limit} \sim -10^{-20}$, $\delta_{\gamma,1}^{limit} \sim -10^{-38} \; \mathrm{eV^{-1}}$ and $\delta_{\gamma,2}^{limit} \sim -10^{-56} \; \mathrm{eV^{-2}}$ in the astrophysical scenario which best describes UHECR data.


title: Cosmic Ray Antiprotons at High Energies

abstract: Cosmic ray antiprotons provide a powerful tool to probe dark matter annihilations in our galaxy. The sensitivity of this important channel is, however, diluted by sizable uncertainties in the secondary antiproton background. In this work, we improve the calculation of secondary antiproton production with a particular focus on the high energy regime. We employ the most recent collider data and identify a substantial increase of antiproton cross sections with energy. This increase is driven by the violation of Feynman scaling as well as by an enhanced strange hyperon production. The updated antiproton production cross sections are made publicly available for independent use in cosmic ray studies. In addition, we provide the correlation matrix of cross section uncertainties for the AMS-02 experiment. At high energies, the new cross sections improve the compatibility of the AMS-02 data with a pure secondary origin of antiprotons in cosmic rays.


title: \'Etale homotopy groups of algebraic groups and homogeneous spaces

abstract: We show the vanishing of the second homotopy group of the \'etale homotopy type of a smooth connected algebraic group over a separably closed field, completed away from the characteristic. This is an algebraic analogue of a classical theorem of Elie Cartan. Based on this result, we establish an explicit formula for the similarly completed second homotopy group of a homogeneous space.


title: Toward the design principles of molecular machines

abstract: Living things avoid equilibrium using molecular machines. Such microscopic soft-matter objects encounter relatively large friction and fluctuations. We discuss design principles for effective molecular machine operation in this unfamiliar context.


title: 3D Morphology Prediction of Progressive Spinal Deformities from Probabilistic Modeling of Discriminant Manifolds

abstract: We introduce a novel approach for predicting the progression of adolescent idiopathic scoliosis from 3D spine models reconstructed from biplanar X-ray images. Recent progress in machine learning have allowed to improve classification and prognosis rates, but lack a probabilistic framework to measure uncertainty in the data. We propose a discriminative probabilistic manifold embedding where locally linear mappings transform data points from high-dimensional space to corresponding low-dimensional coordinates. A discriminant adjacency matrix is constructed to maximize the separation between progressive and non-progressive groups of patients diagnosed with scoliosis, while minimizing the distance in latent variables belonging to the same class. To predict the evolution of deformation, a baseline reconstruction is projected onto the manifold, from which a spatiotemporal regression model is built from parallel transport curves inferred from neighboring exemplars. Rate of progression is modulated from the spine flexibility and curve magnitude of the 3D spine deformation. The method was tested on 745 reconstructions from 133 subjects using longitudinal 3D reconstructions of the spine, with results demonstrating the discriminatory framework can identify between progressive and non-progressive of scoliotic patients with a classification rate of 81% and prediction differences of 2.1$^{o}$ in main curve angulation, outperforming other manifold learning methods. Our method achieved a higher prediction accuracy and improved the modeling of spatiotemporal morphological changes in highly deformed spines compared to other learning methods.


title: Positive feedback in coordination games: stochastic evolutionary dynamics and the logit choice rule

abstract: We study the problem of stochastic stability for evolutionary dynamics under the logit choice rule. We consider general classes of coordination games, symmetric or asymmetric, with an arbitrary number of strategies, which satisfies the marginal bandwagon property (i.e., there is positive feedback to coordinate). Our main result is that the most likely evolutionary escape paths from a status quo convention consist of a series of identical mistakes. As an application of our result, we show that the Nash bargaining solution arises as the long run convention for the evolutionary Nash demand game under the usual logit choice rule. We also obtain a new bargaining solution if the logit choice rule is combined with intentional idiosyncratic plays. The new bargaining solution is more egalitarian than the Nash bargaining solution, demonstrating that intentionality implies equality under the logit choice model.


title: Optimal control of linear systems with limited control actions: threshold-based event-triggered control

abstract: We consider a finite-horizon linear-quadratic optimal control problem where only a limited number of control messages are allowed for sending from the controller to the actuator. To restrict the number of control actions computed and transmitted by the controller, we employ a threshold-based event-triggering mechanism that decides whether or not a control message needs to be calculated and delivered. Due to the nature of threshold-based event-triggering algorithms, finding the optimal control sequence requires minimizing a quadratic cost function over a non-convex domain. In this paper, we firstly provide an exact solution to the non-convex problem mentioned above by solving an exponential number of quadratic programs. To reduce computational complexity, we, then, propose two efficient heuristic algorithms based on greedy search and the Alternating Direction Method of Multipliers (ADMM) method. Later, we consider a receding horizon control strategy for linear systems controlled by event-triggered controllers, and we also provide a complete stability analysis of receding horizon control that uses finite horizon optimization in the proposed class. Numerical examples testify to the viability of the presented design technique.


title: Formality for g-manifolds

abstract: To any $\mathfrak{g}$-manifold $M$ are associated two dglas $\operatorname{tot}\big(\Lambda^{\bullet} \mathfrak{g}^\vee \otimes_{\Bbbk} T_{\operatorname{poly}}^{\bullet} \big)$ and $\operatorname{tot} \big(\Lambda^{\bullet} \mathfrak{g}^\vee\otimes_{\Bbbk} D_{\operatorname{poly}}^{\bullet} \big)$, whose cohomologies $H_{\operatorname{CE}}(\mathfrak{g}, T_{\operatorname{poly}}^{\bullet} \xrightarrow{0} T_{\operatorname{poly}}^{\bullet+1})$ and $H_{\operatorname{CE}}(\mathfrak{g}, D_{\operatorname{poly}}^{\bullet} \xrightarrow{0} D_{\operatorname{poly}}^{\bullet+1})$ are Gerstenhaber algebras. We establish a formality theorem for $\mathfrak{g}$-manifolds: there exists an $L_\infty$ quasi-isomorphism $\Phi: \operatorname{tot}\big(\Lambda^{\bullet} \mathfrak{g}^\vee \otimes_{\Bbbk} T_{\operatorname{poly}}^{\bullet} \big) \to \operatorname{tot} \big(\Lambda^{\bullet} \mathfrak{g}^\vee\otimes_{\Bbbk} D_{\operatorname{poly}}^{\bullet} \big)$ whose first `Taylor coefficient' (1) is equal to the Hochschild-Kostant-Rosenberg map twisted by the square root of the Todd cocycle of the $\mathfrak{g}$-manifold $M$ and (2) induces an isomorphism of Gerstenhaber algebras on the level of cohomology. Consequently, the Hochschild-Kostant-Rosenberg map twisted by the square root of the Todd class of the $\mathfrak{g}$-manifold $M$ is an isomorphism of Gerstenhaber algebras from $H_{\operatorname{CE}}(\mathfrak{g}, T_{\operatorname{poly}}^{\bullet} \xrightarrow{0} T_{\operatorname{poly}}^{\bullet+1})$ to $H_{\operatorname{CE}}(\mathfrak{g}, D_{\operatorname{poly}}^{\bullet} \xrightarrow{0} D_{\operatorname{poly}}^{\bullet+1})$.


title: Tree Structured Synthesis of Gaussian Trees

abstract: A new synthesis scheme is proposed to effectively generate a random vector with prescribed joint density that induces a (latent) Gaussian tree structure. The quality of synthesis is measured by total variation distance between the synthesized and desired statistics. The proposed layered and successive encoding scheme relies on the learned structure of tree to use minimal number of common random variables to synthesize the desired density. We characterize the achievable rate region for the rate tuples of multi-layer latent Gaussian tree, through which the number of bits needed to simulate such Gaussian joint density are determined. The random sources used in our algorithm are the latent variables at the top layer of tree, the additive independent Gaussian noises, and the Bernoulli sign inputs that capture the ambiguity of correlation signs between the variables.


title: Phase diagram and neutron spin resonance of superconducting NaFe$_{1-x}$Cu$_x$As

abstract: We use transport and neutron scattering to study the electronic phase diagram and spin excitations of NaFe$_{1-x}$Cu$_x$As single crystals. Similar to Co- and Ni-doped NaFeAs, a bulk superconducting phase appears near $x\approx2\%$ with the suppression of stripe-type magnetic order in NaFeAs. Upon further increasing Cu concentration the system becomes insulating, culminating in an antiferromagnetically ordered insulating phase near $x\approx 50\%$. Using transport measurements, we demonstrate that the resistivity in NaFe$_{1-x}$Cu$_x$As exhibits non-Fermi-liquid behavior near $x\approx1.8\%$. Our inelastic neutron scattering experiments reveal a single neutron spin resonance mode exhibiting weak dispersion along $c$-axis in NaFe$_{0.98}$Cu$_{0.02}$As. The resonance is high in energy relative to the superconducting transition temperature $T_{\rm c}$ but weak in intensity, likely resulting from impurity effects. These results are similar to other iron pnictides superconductors despite the superconducting phase in NaFe$_{1-x}$Cu$_x$As is continuously connected to an antiferromagnetically ordered insulating phase near $x\approx 50\%$ with significant electronic correlations. Therefore, electron correlations is an important ingredient of superconductivity in NaFe$_{1-x}$Cu$_x$As and other iron pnictides.


title: Nowhere differentiable functions with respect to the position

abstract: Let $\Omega$ be a bounded domain in $\mathbb{C}$ such that $\partial \Omega$ does not contain isolated points. Let $R(\Omega)$ be the space of uniform limits on $\overline{\Omega}$ of rational functions with poles off $\overline{\Omega}$, endowed with the supremum norm. We prove that either generically all functions $f$ in $R(\Omega)$ satisfy % $$ \limsup_{\substack{z \to z_0 z \in \partial \Omega}} \Big| \frac{f(z) - f(z_0)}{z - z_0} \Big| = + \infty $$ for every $z_0 \in \partial \Omega$ or no such function in $R(\Omega)$ meets this requirement. In the first case, the generic function $f \in R(\Omega)$ is nowhere differentiable on $\partial \Omega$ with respect to the position. We give specific examples where each case of the previous dichotomy holds. We also extend the previous result to unbounded domains.


title: Are Continuum Predictions of Clustering Chaotic?

abstract: Gas-solid multiphase flows are prone to develop an instability known as clustering. Two-fluid models, which treat the particulate phase as a continuum, are known to reproduce the qualitative features of this instability, producing highly-dynamic, spatiotemporal patterns. However, it is unknown whether such simulations are truly aperiodic or a type of complex periodic behavior. By showing that the system possesses a sensitive dependence on initial conditions and a positive largest Lyapunov exponent, $\lambda_1 \approx 1/\tau$, we provide a tentative answer: continuum predictions of clustering are chaotic. We further demonstrate that the chaotic behavior is dimensionally dependent, a conclusion which unifies previous results and strongly suggests that the chaotic behavior is not a result of the fundamental kinematic instability, but of the secondary (inherently multidimensional) instability.


title: The Motivic Cofiber of $\tau$

abstract: Consider the Tate twist $\tau \in H^{0,1}(S^{0,0})$ in the mod 2 cohomology of the motivic sphere. After 2-completion, the motivic Adams spectral sequence realizes this element as a map $\tau \colon S^{0,-1} \to S^{0,0}$, with cofiber $C\tau$. We show that this motivic 2-cell complex can be endowed with a unique $E_{\infty}$ ring structure. Moreover, this promotes the known isomorphism $\pi_{\ast,\ast} C\tau \cong \mathrm{Ext}^{\ast,\ast}_{BP_{\ast}BP}(BP_{\ast},BP_{\ast})$ to an isomorphism of rings which also preserves higher products. We then consider the closed symmetric monoidal category $({ }_{C\tau}\textbf{Mod}, - \wedge_{C\tau} -)$ which lives in the kernel of Betti realization. Given a motivic spectrum $X$, the $C\tau$-induced spectrum $X \wedge C\tau$ is usually better behaved and easier to understand than $X$ itself. We specifically illustrate this concept in the examples of the mod 2 Eilenberg-Maclane spectrum $H\mathbb{F}_2$, the mod 2 Moore spectrum $S^{0,0}/2$ and the connective hermitian $K$-theory spectrum $kq$.


title: Radiation from transmission lines PART I: free space transmission lines

abstract: This work derives exact expressions for the radiation from two conductors non isolated TEM transmission lines of any cross section in free space. We cover the cases of infinite, semi-infinite and finite transmission lines and show that while an infinite transmission line does not radiate, there is a smooth transition between the radiation from a finite to a semi-infinite transmission line. Our analysis is in the frequency domain and we consider transmission lines carrying any combination of forward and backward waves. The analytic results are validated by successful comparison with ANSYS commercial software simulation results, and successful comparisons with other published results.


title: Backward Raman Amplification in the Long-wavelength Infrared

abstract: The wealth of work in backward Raman amplification in plasma has focused on the extreme intensity limit, however backward Raman amplification may also provide an effective and practical mechanism for generating intense, broad bandwidth, long-wavelength infrared radiation (LWIR). An electromagnetic simulation coupled with a relativistic cold fluid plasma model is used to demonstrate the generation of picosecond pulses at a wavelength of 10 microns with terawatt powers through backward Raman amplification. The effects of collisional damping, Landau damping, pump depletion, and wave breaking are examined, as well as the resulting design considerations for a LWIR Raman amplifier.


title: A New Family of Asymmetric Distributions for Modeling Light-Tailed and Right-Skewed Data

abstract: A new three-parameter cumulative distribution function defined on $(\alpha,\infty)$, for some $\alpha\geq0$, with asymmetric probability density function and showing exponential decays at its both tails, is introduced. The new distribution is near to familiar distributions like the gamma and log-normal distributions, but this new one shows own elements and thus does not generalize neither of these distributions. Hence, the new distribution constitutes a new alternative to fit values showing light-tailed behaviors. Further, this new distribution shows great flexibility to fit the bulk of data by tuning some parameters. We refer to this new distribution as the generalized exponential log-squared distribution (GEL-S). Statistical properties of the GEL-S distribution are discussed. The maximum likelihood method is proposed for estimating the model parameters, but incorporating adaptations in computational procedures due to difficulties in the manipulation of the parameters. The perfomance of the new distribution is studied using simulations. Applications to real data sets coming from different domains are showed.


title: $K\pi$ $I=1/2$ $S$-wave from $\eta_c$ decay data at BaBar and classic Meson-Meson scattering from LASS

abstract: A recent analysis of data on the two photon production of the $\eta_c$ and its decay to $K(K\pi)$ has determined the $K\pi$ $S$-wave amplitude in a "model-independent" way assuming primarily that the additional kaon is a spectator in this decay. The purpose of this paper is to fit these results, together with classic $K\pi$ production data from LASS, within a formalism that implements unitarity for the di-meson interaction. This fixes the $I=1/2$ $K\pi\to K\pi$ $S$-wave amplitude up to 2.4 GeV. This resolves the Barrelet ambiguity in the original LASS analysis, and constrains the amount of inelasticity in $K\pi$ scattering, highlighting that this becomes significant beyond 1.8 GeV. This result needs to be checked by experimental information on the many inelastic channels, in particular $K\eta^\prime$ and $K\pi\pi\pi$. Our analysis provides a single representation for the $K\pi$ $S$-wave from threshold, controlled by Chiral Perturbation Theory, through the broad $\kappa$, $K_0^*(1430)$ and $K_0^*(1950)$ resonances. There is no arbitrary sum of Breit-Wigner forms and random backgrounds for real $K\pi$ masses. Rather the form provides a representation that can be translated to other processes with $K\pi$ interactions with their own coupling functions, while automatically maintaining consistency with the chiral dynamics near threshold, with the LASS data and the new results on $\eta_c$ decay.


title: A Unifying Perspective: Solitary Traveling Waves As Discrete Breathers And Energy Criteria For Their Stability

abstract: In this work, we provide two complementary perspectives for the (spectral) stability of solitary traveling waves in Hamiltonian nonlinear dynamical lattices, of which the Fermi-Pasta-Ulam and the Toda lattice are prototypical examples. One is as an eigenvalue problem for a stationary solution in a co-traveling frame, while the other is as a periodic orbit modulo shifts. We connect the eigenvalues of the former with the Floquet multipliers of the latter and based on this formulation derive an energy-based spectral stability criterion. It states that a sufficient (but not necessary) condition for a change in the wave stability occurs when the functional dependence of the energy (Hamiltonian) $H$ of the model on the wave velocity $c$ changes its monotonicity. Moreover, near the critical velocity where the change of stability occurs, we provide explicit leading-order computation of the unstable eigenvalues, based on the second derivative of the Hamiltonian $H"(c_0)$ evaluated at the critical velocity $c_0$. We corroborate this conclusion with a series of analytically and numerically tractable examples and discuss its parallels with a recent energy-based criterion for the stability of discrete breathers.


title: Fractional Parts and their Relations to the Values of the Riemann Zeta Function

abstract: A well-known result, due to Dirichlet and later generalized by de la Vallee-Poussin, expresses a relationship between the sum of fractional parts and the Euler-Mascheroni constant. In this paper, we prove an asymptotic relationship between the summation of the products of fractional parts with powers of integers on one hand, and the values of the Riemann zeta function, on the other hand. Dirichlet's classical result falls as a particular case of this more general theorem.


title: Probing Radiative Neutrino Mass Models At The LHC Via Trilepton Events

abstract: Trilepton event represents one of the probes of the new physics at high energy colliders. In this talk, we consider the search for processes with final states $\ell_{\alpha}^{\pm }\ell_{\beta}^{\pm}\ell_{\gamma}^{\mp}$ + $\slashed{E}_{T}$ where ${\alpha}$, ${\beta}$, ${\gamma}$= $e,\mu,\tau$, via the production of singlet charged scalar $S^{\pm}$ which arise in a class of radiative neutrino mass models. We discuss the opposite sign same flavor leptons signal, as well as the background free channel in view to get a significant excess at $\sqrt{s}$= 8 TeV and $\sqrt{s}$ = 14 TeV at the hadron collider LHC.


title: Interpolating sequences in spaces with the complete Pick property

abstract: We characterize interpolating sequences for multiplier algebras of spaces with the complete Pick property. Specifically, we show that a sequence is interpolating if and only if it is separated and generates a Carleson measure. This generalizes results of Carleson for the Hardy space and of Bishop, Marshall and Sundberg for the Dirichlet space. Furthermore, we investigate interpolating sequences for pairs of Hilbert function spaces.


title: Simplicial Homotopy Theory, Link Homology and Khovanov Homology

abstract: The purpose of this note is to point out that simplicial methods and the well-known Dold-Kan construction in simplicial homotopy theory can be fruitfully applied to convert link homology theories into homotopy theories. Dold and Kan prove that there is a functor from the category of chain complexes over a commutative ring with unit to the category of simplicial objects over that ring such that chain homotopic maps go to homotopic maps in the simplicial category. Furthermore, this is an equivalence of categories. In this way, given a link homology theory, we construct a mapping taking link diagrams to a category of simplicial objects such that up to looping or delooping, link diagrams related by Reidemeister moves will give rise to homotopy equivalent simplicial objects, and the homotopy groups of these objects will be equal to the link homology groups of the original link homology theory. The construction is independent of the particular link homology theory. A simplifying point in producing a homotopy simplicial object in relation to a chain complex occurs when the chain complex is itself derived (via face maps) from a simplicial object that satisfies the Kan extension condition. Under these circumstances one can use that simplicial object rather than apply the Dold-Kan functor to the chain complex. We will give examples of this situation in regard to Khovanov homology. We will investigate detailed working out of this correspondence in separate papers. The purpose of this note is to announce the basic relationships for using simplicial methods in this domain. Thus we do more than just quote the Dold-Kan Theorem. We give a review of simplicial theory and we point to specific constructions, particularly in relation to Khovanov homology, that can be used to make simplicial homotopy types directly.


title: Cauchy's residue sore thumb

abstract: Cauchy's method from two centuries ago for computing integrals along the real axis by passing into the complex plane is not rigorous by present-day standards. Yet when properly formulated, his original approach is simpler than modern presentations of the residue calculus.


title: The connectivity of graphs of graphs with self-loops and a given degree sequence

abstract: `Double edge swaps' transform one graph into another while preserving the graph's degree sequence, and have thus been used in a number of popular Markov chain Monte Carlo (MCMC) sampling techniques. However, while double edge-swaps can transform, for any fixed degree sequence, any two graphs inside the classes of simple graphs, multigraphs, and pseudographs, this is not true for graphs which allow self-loops but not multiedges (loopy graphs). Indeed, we exactly characterize the degree sequences where double edge swaps cannot reach every valid loopy graph and develop an efficient algorithm to determine such degree sequences. The same classification scheme to characterize degree sequences can be used to prove that, for all degree sequences, loopy graphs are connected by a combination of double and triple edge swaps. Thus, we contribute the first MCMC sampler that uniformly samples loopy graphs with any given sequence.

categories: math.CO cs.SI physics.data-an physics.soc-ph


title: Efficient and Adaptive Linear Regression in Semi-Supervised Settings

abstract: We consider the linear regression problem under semi-supervised settings wherein the available data typically consists of: (i) a small or moderate sized 'labeled' data, and (ii) a much larger sized 'unlabeled' data. Such data arises naturally from settings where the outcome, unlike the covariates, is expensive to obtain, a frequent scenario in modern studies involving large databases like electronic medical records (EMR). Supervised estimators like the ordinary least squares (OLS) estimator utilize only the labeled data. It is often of interest to investigate if and when the unlabeled data can be exploited to improve estimation of the regression parameter in the adopted linear model. In this paper, we propose a class of 'Efficient and Adaptive Semi-Supervised Estimators' (EASE) to improve estimation efficiency. The EASE are two-step estimators adaptive to model mis-specification, leading to improved (optimal in some cases) efficiency under model mis-specification, and equal (optimal) efficiency under a linear model. This adaptive property, often unaddressed in the existing literature, is crucial for advocating 'safe' use of the unlabeled data. The construction of EASE primarily involves a flexible 'semi-non-parametric' imputation, including a smoothing step that works well even when the number of covariates is not small; and a follow up 'refitting' step along with a cross-validation (CV) strategy both of which have useful practical as well as theoretical implications towards addressing two important issues: under-smoothing and over-fitting. We establish asymptotic results including consistency, asymptotic normality and the adaptive properties of EASE. We also provide influence function expansions and a 'double' CV strategy for inference. The results are further validated through extensive simulations, followed by application to an EMR study on auto-immunity.


title: Blind Deconvolution with Additional Autocorrelations via Convex Programs

abstract: In this work we characterize all ambiguities of the linear (aperiodic) one-dimensional convolution on two fixed finite-dimensional complex vector spaces. It will be shown that the convolution ambiguities can be mapped one-to-one to factorization ambiguities in the $z-$domain, which are generated by swapping the zeros of the input signals. We use this polynomial description to show a deterministic version of a recently introduced masked Fourier phase retrieval design. In the noise-free case a (convex) semi-definite program can be used to recover exactly the input signals if they share no common factors (zeros). Then, we reformulate the problem as deterministic blind deconvolution with prior knowledge of the autocorrelations. Numerically simulations show that our approach is also robust against additive noise.


title: Data-pattern tomography of entangled states

abstract: We discuss the data-pattern tomography for reconstruction of entangled states of light. We show that for a moderate number of probe coherent states it is possible to achieve high accuracy of representation not only for single-mode states but also for two-mode entangled states. We analyze the stability of these representations to the noise and demonstrate the conservation of the purity and entanglement. Simulating the probe and signal measurements, we show that systematic error inherent for representation of realistic signal response with finite sets of responses from probe states still allows one to infer reliably the signal states preserving entanglement.


title: The Effect of Protoplanetary Disk Cooling Times on the Formation of Gas Giant Planets by Gravitational Instability

abstract: Observational evidence exists for the formation of gas giant planets on wide orbits around young stars by disk gravitational instability, but the roles of disk instability and core accretion for forming gas giants on shorter period orbits are less clear. The controversy extends to population synthesis models of exoplanet demographics and to hydrodynamical models of the fragmentation process. The latter refers largely to the handling of radiative transfer in three dimensional (3D) hydrodynamical models, which controls heating and cooling processes in gravitationally unstable disks, and hence dense clump formation. A suite of models using the $\beta$ cooling approximation is presented here. The initial disks have masses of 0.091 $M_\odot$ and extend from 4 to 20 AU around a 1 $M_\odot$ protostar. The initial minimum Toomre $Q_i$ values range from 1.3 to 2.7, while $\beta$ ranges from 1 to 100. We show that the choice of $Q_i$ is equal in importance to the $\beta$ value assumed: high $Q_i$ disks can be stable for small $\beta$, when the initial disk temperature is taken as a lower bound, while low $Q_i$ disks can fragment for high $\beta$. These results imply that the evolution of disks toward low $Q_i$ must be taken into account in assessing disk fragmentation possibilities, at least in the inner disk, i.e., inside about 20 AU. The models suggest that if low $Q_i$ disks can form, there should be an as yet largely undetected population of gas giants orbiting G dwarfs between about 6 AU and 16 AU.


title: Grid cells show field-to-field variability and this explains the aperiodic response of inhibitory interneurons

abstract: Research on network mechanisms and coding properties of grid cells assume that the firing rate of a grid cell in each of its fields is the same. Furthermore, proposed network models predict spatial regularities in the firing of inhibitory interneurons that are inconsistent with experimental data. In this paper, by analyzing the response of grid cells recorded from rats during free navigation, we first show that there are strong variations in the mean firing rate of the fields of individual grid cells and thus show that the data is inconsistent with the theoretical models that predict similar peak magnitudes. We then build a two population excitatory-inhibitory network model in which sparse spatially selective input to the excitatory cells, presumed to arise from e.g. salient external stimuli, hippocampus or a combination of both, leads to the variability in the firing field amplitudes of grid cells. We show that, when combined with appropriate connectivity between the excitatory and inhibitory neurons, the variability in the firing field amplitudes of grid cells results in inhibitory neurons that do not exhibit regular spatial firing, consistent with experimental data. Finally, we show that even if the spatial positions of the fields are maintained, variations in the firing rates of the fields of grid cells are enough to cause remapping of hippocampal cells.


title: Bose-Einstein Condensation and Symmetry Breaking of a Complex Charged Scalar Field

abstract: In this work the Klein-Gordon (KG) equation for a complex scalar field with U(1) symmetry endowed in a mexican-hat scalar field potential with thermal and electromagnetic contributions is written as a Gross-Pitaevskii (GP)-like equation. This equation is interpreted as a charged generalization of the GP equation at finite temperatures found in previous works. Its hydrodynamical representation is obtained and the corresponding thermodynamical properties are derived and related to measurable quantities. The condensation temperature in the non-relativistic regime associated with the aforementioned system within the semiclassical approximation is calculated. Also, a generalized equation for the conservation of energy for a charged bosonic gas is found when electromagnetic fields are introduced, and it is studied how under certain circumstances its breaking of symmetry can give some insight on the phase transition of the system not just into the condensed phase but also on other related systems.


title: Unknowable Manipulators: Social Network Curator Algorithms

abstract: For a social networking service to acquire and retain users, it must find ways to keep them engaged. By accurately gauging their preferences, it is able to serve them with the subset of available content that maximises revenue for the site. Without the constraints of an appropriate regulatory framework, we argue that a sufficiently sophisticated curator algorithm tasked with performing this process may choose to explore curation strategies that are detrimental to users. In particular, we suggest that such an algorithm is capable of learning to manipulate its users, for several qualitative reasons: 1. Access to vast quantities of user data combined with ongoing breakthroughs in the field of machine learning are leading to powerful but uninterpretable strategies for decision making at scale. 2. The availability of an effective feedback mechanism for assessing the short and long term user responses to curation strategies. 3. Techniques from reinforcement learning have allowed machines to learn automated and highly successful strategies at an abstract level, often resulting in non-intuitive yet nonetheless highly appropriate action selection. In this work, we consider the form that these strategies for user manipulation might take and scrutinise the role that regulation should play in the design of such systems.


title: Three-band Hubbard model for Na$_2$IrO$_3$: Topological insulator, zigzag antiferromagnet, and Kitaev-Heisenberg material

abstract: Na$_2$IrO$_3$ was one of the first materials proposed to feature the Kane-Mele type topological insulator phase. Contemporaneously it was claimed that the very same material is in a Mott insulating phase which is described by the Kitaev-Heisenberg (KH) model. First experiments indeed revealed Mott insulating behavior in conjunction with antiferromagnetic long-range order. Further refined experiments established antiferromagnetic order of zigzag type which is not captured by the KH model. Since then several extensions and modifications of the KH model were proposed in order to describe the experimental findings. Here we suggest that adding charge fluctuations to the KH model represents an alternative explanation of zigzag antiferromagnetism. Moreover, a phenomenological three-band Hubbard model unifies all the pieces of the puzzle: topological insulator physics for weak and KH model for strong electron-electron interactions as well as a zigzag antiferromagnet at intermediate interaction strength.


title: Grain Boundary Resistance in Copper Interconnects from an Atomistic Model to a Neural Network

abstract: Orientation effects on the resistivity of copper grain boundaries are studied systematically with two different atomistic tight binding methods. A methodology is developed to model the resistivity of grain boundaries using the Embedded Atom Model, tight binding methods and non-equilibrum Green's functions (NEGF). The methodology is validated against first principles calculations for small, ultra-thin body grain boundaries (<5nm) with 6.4% deviation in the resistivity. A statistical ensemble of 600 large, random structures with grains is studied. For structures with three grains, it is found that the distribution of resistivities is close to normal. Finally, a compact model for grain boundary resistivity is constructed based on a neural network.


title: The JCMT Gould Belt Survey: A First Look at IC 5146

abstract: We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.


title: Scattering of a composite quasiparticle by an impurity on a lattice

abstract: We study scattering of a composite quasiparticle, which possesses a degree of freedom corresponding to relative separation between two bound excitations, by a delta-like impurity potential on a one-dimensional discrete lattice. Firstly, we show that, due to specific properties of their dispersion, lattice excitations bind to impurities with both negative and positive potentials. We demonstrate that the finite size of the composite excitation leads to formation of multiple excitation-impurity bound states. The number and the degree of localization of these bound states depend on the signs and relative magnitudes of the impurity potential and the binding strength of two quasiparticles. We also report the existence of excitation-impurity bound states whose energies are located in the continuum band. Secondly, we study a change in the entanglement between the centre of mass and relative coordinate degrees of freedom of a biexciton wave packet during single impurity scattering and decoherence caused by it. For a composite quasiparticle on a lattice, the entanglement between its relative and centre of mass coordinate degrees of freedom arises naturally due to inseparability of the two-particle Hamiltonian. One of the main focuses of our study is to investigate how this inseparability affects the creation of the biexciton-impurity bound states and the entanglement dynamics.

categories: quant-ph cond-mat.other gr-qc math-ph math.MP


title: Asynchronous Parallel Algorithms for Nonconvex Big-Data Optimization. Part II: Complexity and Numerical Results

abstract: We present complexity and numerical results for a new asynchronous parallel algorithmic method for the minimization of the sum of a smooth nonconvex function and a convex nonsmooth regularizer, subject to both convex and nonconvex constraints. The proposed method hinges on successive convex approximation techniques and a novel probabilistic model that captures key elements of modern computational architectures and asynchronous implementations in a more faithful way than state-of-the-art models. In the companion paper we provided a detailed description on the probabilistic model and gave convergence results for a diminishing stepsize version of our method. Here, we provide theoretical complexity results for a fixed stepsize version of the method and report extensive numerical comparisons on both convex and nonconvex problems demonstrating the efficiency of our approach.


title: Multi-channel Sensing And Resource Allocation in Energy Constrained Cognitive Radio Networks

abstract: We consider a cognitive radio network in a multi-channel licensed environment. Secondary user transmits in a channel if the channel is sensed to be vacant. This results in a tradeoff between sensing time and transmission time. When secondary users are energy constrained, energy available for transmission is less if more energy is used in sensing. This gives rise to an energy tradeoff. For multiple primary channels, secondary users must decide appropriate sensing time and transmission power in each channel to maximize average aggregate-bit throughput in each frame duration while ensuring quality-of-service of primary users. Considering time and energy as limited resources, we formulate this problem as a resource allocation problem. Initially a single secondary user scenario is considered and solution is obtained using decomposition and alternating optimization techniques. Later we extend the analysis for the case of multiple secondary users. Simulation results are presented to study effect of channel occupancy, fading and energy availability on performance of proposed method.


title: AdS Poisson homogeneous spaces and Drinfel'd doubles

abstract: The correspondence between Poisson homogeneous spaces over a Poisson-Lie group $G$ and Lagrangian Lie subalgebras of the classical double $D({\mathfrak g})$ is revisited and explored in detail for the case in which ${\mathfrak g}=D(\mathfrak a)$ is a classical double itself. We apply these results to give an explicit description of some coisotropic 2d Poisson homogeneous spaces over the group $\mathrm{SL}(2,R)\cong\mathrm{SO}(2,1)$, namely 2d anti de Sitter space, 2d hyperbolic space and the lightcone in 3d Minkowski space. We show how each of these spaces is obtained as a quotient with respect to a Poisson-subgroup for one of the three inequivalent Lie bialgebra structures on ${sl}(2,R)$ and as a coisotropic one for the others. We then construct families of coisotropic Poisson homogeneous structures for 3d anti de Sitter space $\mathrm{AdS}_3$ and show that the ones that are quotients by a Poisson subgroup are determined by a three-parameter family of classical $r$-matrices for ${so}(2,2)$, while the non Poisson-subgroup cases are much more numerous. In particular, we present the two Poisson homogeneous structures on $\mathrm{AdS}_3$ that arise from two Drinfel'd double structures on $\mathrm{SO}(2,2)$. The first one realises $\mathrm{AdS}_3$ as a quotient of $\mathrm{SO}(2,2)$ by the Poisson-subgroup $\mathrm{SL}(2,R)$, while the second one, the non-commutative spacetime of the twisted $\kappa$-AdS deformation, realises $\mathrm{AdS}_3$ as a coisotropic Poisson homogeneous space.


title: The Landau-Lifshitz equation, the NLS, and the magnetic rogue wave as a by-product of two colliding regular "positons"

abstract: In this article we present a new method for construction of exact solutions of the Landau-Lifshitz-Gilbert equation (LLG) for ferromagnetic nanowires. The method is based on the established relationship between the LLG and the nonlinear Schr\"odinger equation (NLS), and is aimed at resolving an old problem: how to produce multiple-rogue wave solutions of NLS using just the Darboux-type transformations. The solutions of this type - known as P-breathers - have been proven to exist by Dubard and Matveev, but their technique heavily relied on using the solutions of yet another nonlinear equation, Kadomtsev-Petviashvili I equation (KP-I), and its relationship with NLS. We have shown that in fact one doesn't have to use KP-I but can instead reach the same results just with NLS solutions, but only if they are dressed via the binary Darboux transformation. In particular, our approach allows to construct all the Dubard-Matveev P-breathers. Furthermore, the new method can lead to some completely new, previously unknown solutions. One particular solution that we have constructed describes two positon-like waves, colliding with each other and in the process producing a new, short-lived rogue wave. We called this unusual solution (rogue wave begotten after the impact of two solitons) the "impacton".


title: Superfluid-Mott-insulator phase transition of light in a two-mode cavity array with ultrastrong coupling

abstract: In this paper we construct a new type of cavity array, in each cavity of which multiple two-level atoms interact with two independent photon modes. This system can be totally governed by a two-mode Dicke-lattice model, which includes all of the counter-rotating terms and therefore works well in the ultrastrong coupling regime achieved in recent experiments. Attributed to its special atom-photon coupling scheme, this model supports a global conserved excitation and a continuous $U(1)$ symmetry, rather than the discrete $Z_{2}$ symmetry in the standard Dicke-lattice model. This distinct change of symmetry via adding an extra photon mode strongly impacts the nature of photon localization/delocalization behavior. Specifically, the atom-photon interaction features stable Mott-lobe structures of photons and a second-order superfluid-Mott-insulator phase transition, which share similarities with the Jaynes-Cummings-lattice and Bose-Hubbard models. More interestingly, the Mott-lobe structures predicted here depend crucially on the atom number of each site. We also show that our model can be mapped into a continuous $XX$ spin model. Finally, we propose a scheme to implement the introduced cavity array in circuit quantum electrodynamics. This work broadens our understanding of strongly-correlated photons.


title: Feedforward Architectures Driven by Inhibitory Interactions

abstract: Directed information transmission is paramount for many social, physical, and biological systems. For neural systems, scientists have studied this problem under the paradigm of feedforward networks for decades. In most models of feedforward networks, activity is exclusively driven by excitatory neurons and the wiring patterns between them, while inhibitory neurons play only a stabilizing role for the network dynamics. Motivated by recent experimental discoveries of hippocampal circuitry, cortical circuitry, and the diversity of inhibitory neurons throughout the brain, here we illustrate that one can construct such networks even if the connectivity between the excitatory units in the system remains random. This is achieved by endowing inhibitory nodes with a more active role in the network. Our findings demonstrate that apparent feedforward activity can be caused by a much broader network-architectural basis than often assumed.