title: Spatial distribution of FIR rotationally excited CH+ and OH emission lines in the Orion Bar PDR

abstract: The abundance of CH+ and OH and excitation are predicted to be enhanced by the presence of vibrationally excited H2 or hot gas (~500-1000 K) in PDRs with high incident FUV radiation field. The excitation may also originate in dense gas (>10^5 cm-3) followed by nonreactive collisions. Previous observations suggest that the CH+ and OH correlate with dense and warm gas, and formation pumping contributes to CH+ excitation. We examine the spatial distribution of the CH+ and OH emission in the Orion Bar to establish their physical origin and main formation and excitation mechanisms. We present spatially sampled maps of the CH+ J=3-2 transition at 119.8 {\mu}m and the OH {\Lambda}-doublet at 84 {\mu}m in the Orion Bar over an area of 110"x110" with Herschel (PACS). We compare the spatial distribution of these molecules with those of their chemical precursors, C+, O and H2, and tracers of warm and dense gas. We assess the spatial variation of CH+ J=2-1 velocity-resolved line profile observed with Herschel (HIFI). The OH and CH+ lines correlate well with the high-J CO emission and delineate the warm and dense molecular region. While similar, the differences in the CH+ and OH morphologies indicate that CH+ formation and excitation are related to the observed vibrationally excited H2. This indicates that formation pumping contributes to the excitation of CH+. Interestingly, the peak of the rotationally excited OH 84 {\mu}m emission coincides with a bright young object, proplyd 244-440, which shows that OH can be an excellent tracer of UV-irradiated dense gas. The spatial distribution of CH+ and OH revealed in our maps is consistent with previous modeling studies. Both formation pumping and nonreactive collisions in a UV-irradiated dense gas are important CH+ J=3-2 excitation processes. The excitation of the OH {\Lambda}-doublet at 84 {\mu}m is mainly sensitive to the temperature and density.


title: How to Compute Modulo Prime-Power Sums

abstract: The problem of computing modulo prime-power sums is investigated in distributed source coding as well as computation over Multiple-Access Channel (MAC). We build upon group codes and present a new class of codes called Quasi Group Codes (QGC). A QGC is a subset of a group code. These codes are not closed under the group addition. We investigate some properties of QGC's, and provide a packing and a covering bound. Next, we use these bounds to derived achievable rates for distributed source coding as well as computation over MAC. We show that strict improvements over the previously known schemes can be obtained using QGC's.


title: Inversion formulas and range characterizations for the attenuated geodesic ray transform

abstract: We present two range characterizations for the attenuated geodesic X-ray transform defined on pairs of functions and one-forms on simple surfaces. Such characterizations are based on first isolating the range over sums of functions and one-forms, then separating each sub-range in two ways, first by implicit conditions, second by deriving new inversion formulas for sums of functions and one-forms.


title: Direct and Central Products of Localities

abstract: We develop a theory of direct and central products of partial groups and localities.


title: On topological approach to local theory of surfaces in Calabi-Yau threefolds

abstract: We study the web of dualities relating various enumerative invariants, notably Gromov-Witten invariants and invariants that arise in topological gauge theory. In particular, we study Donaldson-Thomas gauge theory and its reductions to D=4 and D=2 which are relevant to the local theory of surfaces in Calabi-Yau threefolds.


title: Implementation of a discrete Immersed Boundary Method in OpenFOAM

abstract: In this paper, the Immersed Boundary Method (IBM) proposed by Pinelli is implemented for finite volume approximations of incompressible Navier-Stokes equations solutions in the open source toolbox OpenFOAM version 2.2. Solid obstacles are described using a discrete forcing approach for boundary conditions. Unlike traditional approaches encompassing the presence of a solid body using conformal meshes and imposing no-slip boundary conditions on the boundary faces of the mesh, the solid body is here represented on the Eulerian Cartesian mesh through an ad-hoc body force evaluated on a set of Lagrangian markers. The markers can move across the Eulerian mesh, hence allowing for a straightforward analysis of motion or deformation of the body. The IBM method is described and implemented in PisoFOAM, whose Pressure-Implicit Split-Operator (PISO) solver is modified accordingly. The presence of the solid body and the divergence-free of the fluid velocity are imposed simultaneously by sub-iterating between IBM and the pressure correction step. This scheme allows for the use of fast optimized Poisson solvers while granting excellent accuracy with respect to the previously mentioned constraints. Various 2D and 3D well-documented test cases of flows around fixed or moving circular cylinders are simulated and carefully validated against existing data from the literature. The capability of the new solver is discussed in terms of accuracy and numerical performances.


title: Rare event simulation via importance sampling for linear SPDE's

abstract: The goal of this paper is to develop provably efficient importance sampling Monte Carlo methods for the estimation of rare events within the class of linear stochastic partial differential equations (SPDEs). We find that if a spectral gap of appropriate size exists, then one can identify a lower dimensional manifold where the rare event takes place. This allows one to build importance sampling changes of measures that perform provably well even pre-asymptotically (i.e. for small but non-zero size of the noise) without degrading in performance due to infinite dimensionality or due to long simulation time horizons. Simulation studies supplement and illustrate the theoretical results.


title: Non-interacting fermions at finite temperature in a $d$-dimensional trap: universal correlations

abstract: We study a system of $N$ non-interacting spin-less fermions trapped in a confining potential, in arbitrary dimensions $d$ and arbitrary temperature $T$. The presence of the trap introduces an edge where the average density of fermions vanishes. Far from the edge, near the center of the trap (the so called "bulk regime"), physical properties of the fermions have traditionally been understood using the Local Density Approximation. However, this approximation drastically fails near the edge where the density vanishes. In this paper we show that, even near the edge, novel universal properties emerge, independently of the details of the confining potential. We show that for large $N$, these fermions in a confining trap, in arbitrary dimensions and at finite temperature, form a determinantal point process. As a result, any $n$-point correlation function can be expressed as an $n \times n$ determinant whose entry is called the kernel. Near the edge, we derive the large $N$ scaling form of the kernels. In $d=1$ and $T=0$, this reduces to the so called Airy kernel, that appears in the Gaussian Unitary Ensemble (GUE) of random matrix theory. In $d=1$ and $T>0$ we show a remarkable connection between our kernel and the one appearing in the $1+1$-dimensional Kardar-Parisi-Zhang equation at finite time. Consequently our result provides a finite $T$ generalization of the Tracy-Widom distribution, that describes the fluctuations of the rightmost fermion at $T=0$. In $d>1$ and $T \geq 0$, while the connection to GUE no longer holds, the process is still determinantal whose analysis provides a new class of kernels, generalizing the $1d$ Airy kernel at $T=0$ obtained in random matrix theory. Some of our finite temperature results should be testable in present-day cold atom experiments, most notably our detailed predictions for the temperature dependence of the fluctuations near the edge.

categories: cond-mat.stat-mech cond-mat.dis-nn cond-mat.quant-gas math-ph math.MP math.PR


title: Explaining the $R_{K}$ and $R_{D^{(*)}}$ anomalies with vector leptoquarks

abstract: Recently the $B$ factories BaBar and Belle as well as the LHCb experiment have reported several anomalies in the semileptonic $B$ meson decays such as $R_{K}$ and $R_{D^{(*)}}$ etc. We investigate these deviations by considering the vector leptoquarks relevant for both $b \to s l^+ l^-$ and $b \to c l \bar \nu_l$ transitions. The leptoquark parameter space is constrained by using the experimentally measured branching ratios of $B_s \to l^+ l^-$, $\bar B \to X_s l^+ l^- (\nu \bar \nu)$ and $B_u^+ \to l^+ \nu_l$ processes. Using the constrained leptoquark couplings, we compute the branching ratios, forward-backward asymmetries, $\tau$ and $D^*$ polarization parameters in the $\bar B \to D^{(*)} l \bar \nu_l$ processes. We find that the vector leptoquarks can explain both $R_{D^{(*)}}$ and $R_K$ anomalies simultaneously. Furthermore, we study the rare leptonic $B_{u, c}^* \to l \bar\nu$ decay processes in this model.


title: On the energy landscape of the mixed even $p$-spin model

abstract: We investigate the energy landscape of the mixed even $p$-spin model with Ising spin configurations. We show that for any given energy level between zero and the maximal energy, with overwhelming probability there exist exponentially many distinct spin configurations such that their energies stay near this energy level. Furthermore, their magnetizations and overlaps are concentrated around some fixed constants. In particular, at the level of maximal energy, we prove that the Hamiltonian exhibits exponentially many orthogonal peaks. This improves the results of Chatterjee and Ding-Eldan-Zhai, where the former established a logarithmic size of the number of the orthogonal peaks, while the latter proved a polynomial size. Our second main result obtains disorder chaos at zero temperature and at any external field. As a byproduct, this implies that the fluctuation of the maximal energy is superconcentrated when the external field vanishes and obeys a Gaussian limit law when the external field is present.


title: Eta-Eta' mixing and the derivative of the topological susceptibility at zero momentum transfer

abstract: The couplings of the isosinglet axial-vector currents to the Eta and Eta' mesons are evaluated in a stable, model independent way by use of polynomial kernels in dispersion integrals. The corrections to the Gell-Mann-Oakes-Renner relation in the isoscalar channel are deduced. The derivative of the topological susceptibility at the origin is calculated taking into account instantons and instanton screening.


title: Analyzing signal attenuation in PFG anomalous diffusion via a non-Gaussian phase distribution approximation approach by fractional derivatives

abstract: Anomalous diffusion exists widely in polymer and biological systems. Pulsed field gradient (PFG) techniques have been increasingly used to study anomalous diffusion in NMR and MRI. However, the interpretation of PFG anomalous diffusion is complicated. Moreover, there is not an exact signal attenuation expression based on fractional derivatives for PFG anomalous diffusion, which includes the finite gradient pulse width effect. In this paper, a new method, a Mainardi-Luchko-Pagnini (MLP) phase distribution approximation, is proposed to describe PFG fractional diffusion. MLP phase distribution is a non-Gaussian phase distribution. From the fractional diffusion equation based on fractional derivatives in both real space and phase space, the obtained probability distribution function is a MLP distribution. The MLP distribution leads to a Mittag-Leffler function based PFG signal attenuation rather than the exponential or stretched exponential attenuation that is obtained from a Gaussian phase distribution (GPD) under a short gradient pulse approximation. The MLP phase distribution approximation is employed to get a complete signal attenuation expression E{\alpha}(-Dfb*{\alpha},\b{eta}) that includes the finite gradient pulse width effect for all three types of PFG fractional diffusion. The results obtained in this study are in agreement with the results from the literature. These results provide a new, convenient approximation formalism to interpret complex PFG fractional diffusion experiments.


title: Finite LTL Synthesis is EXPTIME-complete

abstract: LTL synthesis -- the construction of a function to satisfy a logical specification formulated in Linear Temporal Logic -- is a 2EXPTIME-complete problem with relevant applications in controller synthesis and a myriad of artificial intelligence applications. In this research note we consider De Giacomo and Vardi's variant of the synthesis problem for LTL formulas interpreted over finite rather than infinite traces. Rather surprisingly, given the existing claims on complexity, we establish that LTL synthesis is EXPTIME-complete for the finite interpretation, and not 2EXPTIME-complete as previously reported. Our result coincides nicely with the planning perspective where non-deterministic planning with full observability is EXPTIME-complete and partial observability increases the complexity to 2EXPTIME-complete; a recent related result for LTL synthesis shows that in the finite case with partial observability, the problem is 2EXPTIME-complete.


title: QUBIC Technical Design Report

abstract: QUBIC is an instrument aiming at measuring the B mode polarisation anisotropies at medium scales angular scales (30-200 multipoles). The search for the primordial CMB B-mode polarization signal is challenging, because of many difficulties: smallness of the expected signal, instrumental systematics that could possibly induce polarization leakage from the large E signal into B, brighter than anticipated polarized foregrounds (dust) reducing to zero the initial hope of finding sky regions clean enough to have a direct primordial B-modes observation. The QUBIC instrument is designed to address all aspects of this challenge with a novel kind of instrument, a Bolometric Interferometer, combining the background-limited sensitivity of Transition-Edge-Sensors and the control of systematics allowed by the observation of interference fringe patterns, while operating at two frequencies to disentangle polarized foregrounds from primordial B mode polarization. Its characteristics are described in details in this Technological Design Report.


title: Consistency of Social Sensing Signatures Across Major US Cities

abstract: Previous studies have shown that Twitter users have biases to tweet from certain locations, locational bias, and during certain hours, temporal bias. We used three years of geolocated Twitter Data to quantify these biases and test our central hypothesis that Twitter users biases are consistent across US cities. Our results suggest that temporal and locational bias of Twitter users are inconsistent between three US metropolitan cities. We derive conclusions about the role of the complexity of the underlying data producing process on its consistency and argue for the potential research avenue for Geospatial Data Science to test and quantify these inconsistencies in the class of organically evolved Big Data.


title: The cross-over from tunnelling to multiphoton ionization of atoms

abstract: We present a theory illuminating the cross-over from strong-field tunnelling ionization to weak-field multiphoton ionization in the interaction of a classical laser field with a hydrogen atom. A simple formula is derived in which the ionization amplitude appears as a product of two separate amplitudes. The first describes the initial polarization of the atom by virtual multiphoton absorption and the second the subsequent tunnelling out of the polarized atom. Tunnelling directly from the ground state and multiphoton absorption without tunnelling appear naturally as the limits of the theory.


title: A Perspective on Deep Imaging

abstract: The combination of tomographic imaging and deep learning, or machine learning in general, promises to empower not only image analysis but also image reconstruction. The latter aspect is considered in this perspective article with an emphasis on medical imaging to develop a new generation of image reconstruction theories and techniques. This direction might lead to intelligent utilization of domain knowledge from big data, innovative approaches for image reconstruction, and superior performance in clinical and preclinical applications. To realize the full impact of machine learning on medical imaging, major challenges must be addressed.


title: Super-resolution community detection for layer-aggregated multilayer networks

abstract: Applied network science often involves preprocessing network data before applying a network-analysis method, and there is typically a theoretical disconnect between these steps. For example, it is common to aggregate time-varying network data into windows prior to analysis, and the tradeoffs of this preprocessing are not well understood. Focusing on the problem of detecting small communities in multilayer networks, we study the effects of layer aggregation by developing random-matrix theory for modularity matrices associated with layer-aggregated networks with $N$ nodes and $L$ layers, which are drawn from an ensemble of Erd\H{o}s-R\'enyi networks. We study phase transitions in which eigenvectors localize onto communities (allowing their detection) and which occur for a given community provided its size surpasses a detectability limit $K^*$. When layers are aggregated via a summation, we obtain $K^*\varpropto \mathcal{O}(\sqrt{NL}/T)$, where $T$ is the number of layers across which the community persists. Interestingly, if $T$ is allowed to vary with $L$ then summation-based layer aggregation enhances small-community detection even if the community persists across a vanishing fraction of layers, provided that $T/L$ decays more slowly than $ \mathcal{O}(L^{-1/2})$. Moreover, we find that thresholding the summation can in some cases cause $K^*$ to decay exponentially, decreasing by orders of magnitude in a phenomenon we call super-resolution community detection. That is, layer aggregation with thresholding is a nonlinear data filter enabling detection of communities that are otherwise too small to detect. Importantly, different thresholds generally enhance the detectability of communities having different properties, illustrating that community detection can be obscured if one analyzes network data using a single threshold.


title: Dynamical History Of The Local Group In $\Lambda$CDM $\rm{II}$ $-$ Including External Perturbers In 3D

abstract: We attempt to fit the observed radial velocities (RVs) of ${\,{\sim}\,}30$ Local Group (LG) galaxies using a 3D dynamical model of it and its immediate environment within the context of the standard cosmological paradigm, $\Lambda$CDM. This extends and confirms the basic results of our previous axisymmetric investigation of the LG (MNRAS, 459, 2237). We find that there remains a tendency for observed RVs to exceed those predicted by our best-fitting model. The typical mismatch is slightly higher than in our 2D model, with a root mean square value of ${\,{\sim}\,}50$ km/s. Our main finding is that including the 3D distribution of massive perturbing dark matter halos is unlikely to help greatly with the high velocity galaxy problem. Nonetheless, the 2D and 3D results differ in several other ways such as which galaxies' RVs are most problematic and the preferred values of parameters common to both models. The anomalously high RVs of several LG dwarfs may be better explained if the Milky Way (MW) and Andromeda (M31) were once moving much faster than in our models. This would allow LG dwarfs to gain very high RVs via gravitational slingshot encounters with a massive fast-moving galaxy. Such a scenario is possible in some modified gravity theories, especially those which require the MW and M31 to have previously undergone a close flyby. In a $\Lambda$CDM context, however, this scenario is not feasible as the resulting dynamical friction would cause a rapid merger.


title: An improved current potential method for fast computation of stellarator coil shapes

abstract: Several fast methods for computing stellarator coil shapes are compared, including the classical NESCOIL procedure [Merkel, Nucl. Fusion 27, 867 (1987)], its generalization using truncated singular value decomposition, and a Tikhonov regularization approach we call REGCOIL in which the squared current density is included in the objective function. Considering W7-X and NCSX geometries, and for any desired level of regularization, we find the REGCOIL approach simultaneously achieves lower surface-averaged and maximum values of both current density (on the coil winding surface) and normal magnetic field (on the desired plasma surface). This approach therefore can simultaneously improve the free-boundary reconstruction of the target plasma shape while substantially increasing the minimum distances between coils, preventing collisions between coils while improving access for ports and maintenance. The REGCOIL method also allows finer control over the level of regularization, it preserves convexity to ensure the local optimum found is the global optimum, and it eliminates two pathologies of NESCOIL: the resulting coil shapes become independent of the arbitrary choice of angles used to parameterize the coil surface, and the resulting coil shapes converge rather than diverge as Fourier resolution is increased. We therefore contend that REGCOIL should be used instead of NESCOIL for applications in which a fast and robust method for coil calculation is needed, such as when targeting coil complexity in fixed-boundary plasma optimization, or for scoping new stellarator geometries.


title: Threaded-Field-Lines Model for the Low Solar Corona Powered by the Alfven Wave Turbulence

abstract: We present an updated global model of the solar corona, including the transition region. We simulate the realistic tree-dimensional (3D) magnetic field using the data from the photospheric magnetic field measurements and assume the magnetohydrodynamic (MHD) Alfv\'en wave turbulence and its non-linear dissipation to be the only source for heating the coronal plasma and driving the solar wind. In closed field regions the dissipation efficiency in a balanced turbulence is enhanced. In the coronal holes we account for a reflection of the outward propagating waves, which is accompanied by generation of weaker counter-propagating waves. The non-linear cascade rate degrades in strongly imbalanced turbulence, thus resulting in colder coronal holes. The distinctive feature of the presented model is the description of the low corona as almost-steady-state low-beta plasma motion and heat flux transfer along the magnetic field lines. We trace the magnetic field lines through each grid point of the lower boundary of the global corona model, chosen at some heliocentric distance, $R=R_{b}\sim1.1\ R_\odot$ well above the transition region. One can readily solve the plasma parameters along the magnetic field line from 1D equations for the plasma motion and heat transport together with the Alfv\'en wave propagation, which adequately describe physics within the heliocentric distances range, $R_{\odot}<R<R_{b}$, in the low solar corona. By interfacing this threaded-field-lines model with the full MHD global corona model at $r=R_{b}$, we find the global solution and achieve a faster-than-real-time performance of the model on $\sim200$ cores.


title: Light scattering study of the "pseudo-layer" compression elastic constant in a twist-bend nematic liquid crystal

abstract: The nematic twist-bend (TB) phase, exhibited by certain achiral thermotropic liquid crystalline (LC) dimers, features a nanometer-scale, heliconical rotation of the average molecular long axis (director) with equally probable left- and right-handed domains. On meso to macroscopic scales, the TB phase may be considered as a stack of equivalent slabs or "pseudo-layers", each one helical pitch in thickness. The long wavelength fluctuation modes should then be analogous to those of a smectic-A phase, and in particular the hydrodynamic mode combining "layer" compression and bending ought to be characterized by an effective layer compression elastic constant $B_{eff}$ and average director splay constant $K_1^{eff}$. The magnitude of $K_1^{eff}$ is expected to be similar to the splay constant of an ordinary nematic LC, but due to the absence of a true mass density wave, $B_{eff}$ could differ substantially from the typical value of $\sim 10^6$ Pa in a conventional smectic-A. Here we report the results of a dynamic light scattering study, which confirms the "pseudo-layer" structure of the TB phase with $B_{eff}$ in the range $\sim 10^3-10^4$ Pa. We show additionally that the temperature dependence of $B_{eff}$ at the TB to nematic transition is accurately described by a coarse-grained free energy density, which is based on a Landau-deGennes expansion in terms of a heli-polar order parameter that characterizes the TB state and is linearly coupled to bend distortion of the director.


title: Layer-dependent properties of SnS2 and SnSe2 novel two-dimensional materials

abstract: The layer dependent structural, electronic and vibrational properties of SnS2 and SnSe2 are investigated using first-principles density functional theory (DFT). The in-plane lattice constants, interlayer distances and binding energies are found to be layer-independent. Bulk SnS2 and SnSe2 are both indirect band gap semiconductors with Eg = 2.18 eV and 1.07 eV, respectively. Few-layer and monolayer 2D systems also possess an indirect band gap, which is increased to 2.41 eV and 1.69 eV for single layers of SnS2 and SnSe2. The effective mass theory of 2D excitons, which takes into account the combined effect of the anisotropy, non-local 2D screening and layer-dependent 3D screening, predicts strong excitonic effects. The binding energy of indirect excitons in monolayer samples, Ex~0.9 eV, is substantially reduced to Ex = 0.14 eV in bulk SnS2 and Ex = 0.09 eV in bulk SnSe2. The layer-dependent Raman spectra display a strong decrease of intensities of the Raman active A1g mode upon decreasing the number of layers down to a monolayer, by a factor of 7 in the case of SnS2 and a factor of 20 in the case of SnSe2 which can be used to identify number of layers in a 2D sample.


title: Warped Convolutions: Efficient Invariance to Spatial Transformations

abstract: Convolutional Neural Networks (CNNs) are extremely efficient, since they exploit the inherent translation-invariance of natural images. However, translation is just one of a myriad of useful spatial transformations. Can the same efficiency be attained when considering other spatial invariances? Such generalized convolutions have been considered in the past, but at a high computational cost. We present a construction that is simple and exact, yet has the same computational complexity that standard convolutions enjoy. It consists of a constant image warp followed by a simple convolution, which are standard blocks in deep learning toolboxes. With a carefully crafted warp, the resulting architecture can be made equivariant to a wide range of two-parameter spatial transformations. We show encouraging results in realistic scenarios, including the estimation of vehicle poses in the Google Earth dataset (rotation and scale), and face poses in Annotated Facial Landmarks in the Wild (3D rotations under perspective).


title: Probabilistic Population Projections for Countries with Generalized HIV/AIDS Epidemics

abstract: The United Nations (UN) issued official probabilistic population projections for all countries to 2100 in July 2015. This was done by simulating future levels of total fertility and life expectancy from Bayesian hierarchical models, and combining the results using a standard cohort-component projection method. The 40 countries with generalized HIV/AIDS epidemics were treated differently from others, in that the projections used the highly multistate Spectrum/EPP model, a complex 15-compartment model that was designed for short-term projections of quantities relevant to policy for the epidemic. Here we propose a simpler approach that is more compatible with the existing UN probabilistic projection methodology for other countries. Changes in life expectancy are projected probabilistically using a simple time series regression model on current life expectancy, HIV prevalence and ART coverage. These are then converted to age- and sex-specific mortality rates using a new family of model life tables designed for countries with HIV/AIDS epidemics that reproduces the characteristic hump in middle adult mortality. These are then input to the standard cohort-component method, as for other countries. The method performed well in an out-of-sample cross-validation experiment. It gives similar population projections to Spectrum/EPP in the short run, while being simpler and avoiding multistate modeling.


title: Electrohydrodynamic Quincke rotation of a prolate ellipsoid

abstract: We experimentally study the occurrence of spontaneous spinning (Quincke rotation) of an ellipsoid in a uniform DC electric field. For an ellipsoid suspended in an unbounded fluid, we find two stable states characterized by the orientation of the ellipsoid long axis relative to the applied electric field : spinless (parallel) and spinning (perpendicular). The phase diagram of ellipsoid behavior as a function of field strength and aspect ratio is in close agreement with the theory of Cebers et al. Phys. Rev .E 63:016301 (2000). We also investigated the dynamics of the ellipsoidal Quincke rotor resting on a planar surface with normal perpendicular to the field direction. We find novel behaviors, such as swinging (long axis oscillating around the applied field direction) and tumbling, due to the confinement.


title: Behavior of the maximum likelihood in quantum state tomography

abstract: Quantum state tomography on a d-dimensional system demands resources that grow rapidly with d. They may be reduced by using model selection to tailor the number of parameters in the model (i.e., the size of the density matrix). Most model selection methods typically rely on a test statistic and a null theory that describes its behavior when two models are equally good. Here, we consider the loglikelihood ratio. Because of the positivity constraint $\rho \geq 0$, quantum state space does not generally satisfy local asymptotic normality, meaning the classical null theory for the loglikelihood ratio (the Wilks theorem) should not be used. Thus, understanding and quantifying how positivity affects the null behavior of this test statistic is necessary for its use in model selection for state tomography. We define a new generalization of local asymptotic normality, metric-projected local asymptotic normality, show that quantum state space satisfies it, and derive a replacement for the Wilks theorem. In addition to enabling reliable model selection, our results shed more light on the qualitative effects of the positivity constraint on state tomography.


title: Dispersive Stiffness of Dzyaloshinskii Domain Walls

abstract: It is well documented that subjecting perpendicular magnetic films which exhibit the interfacial Dzyaloshinskii-Moriya interaction (DMI) to an in-plane magnetic field results in a domain wall (DW) energy, $\sigma$, that is highly anisotropic with respect to the orientation of the DW in the film plane, $\Theta$. We demonstrate that this anisotropy has a profound impact on the elastic response of the DW as characterized by the surface stiffness, $\tilde{\sigma}(\Theta) = \sigma(\Theta) + \sigma"(\Theta)$, and evaluate its dependence on the length scale of deformation. The influence of the stiffness on DW mobility in the creep regime is assessed, with analytic and numerical calculations showing trends in $\tilde{\sigma}$ that better represent experimental measurements of domain wall velocity in magnetic thin films compared to $\sigma$ alone. Our treatment provides experimental support for theoretical models of the mobility of anisotropic elastic manifolds and makes progress toward a more complete understanding of magnetic domain wall creep.


title: 3D Face Reconstruction by Learning from Synthetic Data

abstract: Fast and robust three-dimensional reconstruction of facial geometric structure from a single image is a challenging task with numerous applications. Here, we introduce a learning-based approach for reconstructing a three-dimensional face from a single image. Recent face recovery methods rely on accurate localization of key characteristic points. In contrast, the proposed approach is based on a Convolutional-Neural-Network (CNN) which extracts the face geometry directly from its image. Although such deep architectures outperform other models in complex computer vision problems, training them properly requires a large dataset of annotated examples. In the case of three-dimensional faces, currently, there are no large volume data sets, while acquiring such big-data is a tedious task. As an alternative, we propose to generate random, yet nearly photo-realistic, facial images for which the geometric form is known. The suggested model successfully recovers facial shapes from real images, even for faces with extreme expressions and under various lighting conditions.


title: Relativistic Monte Carlo

abstract: Hamiltonian Monte Carlo (HMC) is a popular Markov chain Monte Carlo (MCMC) algorithm that generates proposals for a Metropolis-Hastings algorithm by simulating the dynamics of a Hamiltonian system. However, HMC is sensitive to large time discretizations and performs poorly if there is a mismatch between the spatial geometry of the target distribution and the scales of the momentum distribution. In particular the mass matrix of HMC is hard to tune well. In order to alleviate these problems we propose relativistic Hamiltonian Monte Carlo, a version of HMC based on relativistic dynamics that introduce a maximum velocity on particles. We also derive stochastic gradient versions of the algorithm and show that the resulting algorithms bear interesting relationships to gradient clipping, RMSprop, Adagrad and Adam, popular optimisation methods in deep learning. Based on this, we develop relativistic stochastic gradient descent by taking the zero-temperature limit of relativistic stochastic gradient Hamiltonian Monte Carlo. In experiments we show that the relativistic algorithms perform better than classical Newtonian variants and Adam.


title: Accurate, Empirical Radii and Masses of Planets and their Host Stars with Gaia Parallaxes

abstract: We present empirical measurements of the radii of 116 stars that host transiting planets. These radii are determined using only direct observables-the bolometric flux at Earth, the effective temperature, and the parallax provided by the Gaia first data release-and thus are virtually model independent, extinction being the only free parameter. We also determine each star's mass using our newly determined radius and the stellar density, itself a virtually model independent quantity from previously published transit analyses. These stellar radii and masses are in turn used to redetermine the transiting planet radii and masses, again using only direct observables. The median uncertainties on the stellar radii and masses are ~8% and ~30%, respectively, and the resulting uncertainties on the planet radii and masses are ~9% and ~22%, respectively. These accuracies are generally larger than previously published model-dependent precisions of ~5% and ~6% on the planet radii and masses, respectively, but the newly determined values are purely empirical. We additionally report radii for 242 stars hosting radial-velocity (non-transiting) planets, with median achieved accuracy of ~2%. Using our empirical stellar masses we verify that the majority of putative "retired A stars" in the sample are indeed more massive than ~1.2 Msun. Most importantly, the bolometric fluxes and angular radii reported here for a total of 498 planet host stars-with median accuracies of 1.7% and 1.8%, respectively-serve as a fundamental dataset to permit the re-determination of transiting planet radii and masses with the Gaia second data release to ~3% and ~5% accuracy, better than currently published precisions, and determined in an entirely empirical fashion.


title: Nonlinear dynamics and chaos in an optomechanical beam

abstract: Optical non-linearities, such as thermo-optic effects and free-carrier-dispersion, are often considered as undesired effects in silicon-based resonators and, more specifically, optomechanical (OM) cavities, affecting the relative detuning between an optical resonance and the excitation laser. However, the interplay between such mechanisms could also enable unexpected physical phenomena to be used in new applications. In the present work, we exploit those non-linearities and their intercoupling with the mechanical degrees of freedom of a silicon OM nanobeam to unveil a rich set of fundamentally different complex dynamics. By smoothly changing the parameters of the excitation laser, namely its power and wavelength, we demonstrate accurate control for activating bi-dimensional and tetra-dimensional limit-cycles, a period doubling route and chaos. In addition, by scanning the laser parameters in opposite senses we demonstrate bistability and hysteresis between bi-dimensional and tetra-dimensional limit-cycles, between different coherent mechanical states and between tetra-dimensional limit-cycles and chaos. As a result of implementing the Rosenstein algorithm to the experimental time series we have extracted a Largest Lyapunov Exponent of LLE=1.3x10^5 s-1, which is between one and two orders of magnitude lower than the typical oscillating frequencies of the system. Most of the experimental features can be well reproduced with a model of first-order non-linear differential equations coupled through the number of intracavity photons. Our findings open several routes towards exploiting silicon-based OM photonic crystals for systems with memory and, more specifically, for chaos based applications such as secure information transfer or sensing.


title: When is $a^{n} + 1$ the sum of two squares?

abstract: Using Fermat's two squares theorem and properties of cyclotomic polynomials, we prove assertions about when numbers of the form $a^{n}+1$ can be expressed as the sum of two integer squares. We prove that $a^n + 1$ is the sum of two squares for all $n \in \mathbb{N}$ if and only if $a$ is a perfect square. We also prove that for $a\equiv 0,1,2\pmod{4},$ if $a^{n} + 1$ is the sum of two squares, then $a^{\delta} + 1$ is the sum of two squares for all $\delta | n, \ \delta>1$. Using Aurifeuillian factorization, we show that if $a$ is a prime and $a\equiv 1 \pmod{4}$, then there are either zero or infinitely many odd $n$ such that $a^n+1$ is the sum of two squares. When $a\equiv 3\pmod{4},$ we define $m$ to be the least positive integer such that $\frac{a+1}{m}$ is the sum of two squares, and prove that if $a^n+1$ is the sum of two squares for any odd integer $n,$ then $m | n$, and both $a^m+1$ and $\frac{n}{m}$ are sums of two squares.


title: Learning Robust Features for Gait Recognition by Maximum Margin Criterion

abstract: In the field of gait recognition from motion capture data, designing human-interpretable gait features is a common practice of many fellow researchers. To refrain from ad-hoc schemes and to find maximally discriminative features we may need to explore beyond the limits of human interpretability. This paper contributes to the state-of-the-art with a machine learning approach for extracting robust gait features directly from raw joint coordinates. The features are learned by a modification of Linear Discriminant Analysis with Maximum Margin Criterion so that the identities are maximally separated and, in combination with an appropriate classifier, used for gait recognition. Experiments on the CMU MoCap database show that this method outperforms eight other relevant methods in terms of the distribution of biometric templates in respective feature spaces expressed in four class separability coefficients. Additional experiments indicate that this method is a leading concept for rank-based classifier systems.


title: Subalgebras of Solomon's descent algebra based on alternating runs

abstract: The number of alternating runs is a natural permutation statistic. We show it can be used to define some commutative subalgebras of the symmetric group algebra, and more precisely of the descent algebra. The Eulerian peak algebras naturally appear as subalgebras of our run algebras. We also calculate the orthogonal idempotents for run algebras in terms of noncommutative symmetric functions.


title: Cardinal kinematics: I. Rotation fields of the APOGEE Survey

abstract: Correlations between stellar chemistry and kinematics have long been used to gain insight into the evolution of the Milky Way Galaxy. Orbital angular momentum is a key physical parameter and it is often estimated from three-dimensional space motions. We here demonstrate the lower uncertainties that can be achieved in the estimation of one component of velocity through selection of stars in key directions and use of line-of-sight velocity alone (i.e. without incorporation of proper motion data). In this first paper we apply our technique to stars observed in the direction of Galactic rotation in the APOGEE survey. We first derive the distribution of azimuthal velocities, $v_\phi$, then from these and observed radial coordinates, estimate the stellar guiding centre radii, $R_g$, within $6.9\leq R \leq 10$ kpc with uncertainties smaller than (or of the order of) 1kpc. We show that there is no simple way to select a clean stellar sample based on low errors on proper motions and distances to obtain high-quality 3D velocities and hence one should pay particular attention when trying to identify kinematically peculiar stars based on velocities derived using the proper motions. Using our azimuthal velocity estimations, we investigate the joint distribution of elemental abundances and rotational kinematics free from the blurring effects of epicyclic motions, and we derive the $\partial v_\phi / \partial R_g$ trends for the thin and thick discs as a function of radius. Our analysis provides further evidence for radial migration within the thin disc and hints against radial migration playing a significant role in the evolution of the thick disc.


title: First Gaia Local Group Dynamics: Magellanic Clouds Proper Motion and Rotation

abstract: We use the Gaia data release 1 (DR1) to study the proper motion (PM) fields of the Large and Small Magellanic Clouds (LMC, SMC). This uses the Tycho-Gaia Astrometric Solution (TGAS) PMs for 29 Hipparcos stars in the LMC and 8 in the SMC. The LMC PM in the West and North directions is inferred to be $(\mu_W,\mu_N) = (-1.872 \pm 0.045, 0.224 \pm 0.054)$ mas/yr, and the SMC PM $(\mu_W,\mu_N) = (-0.874 \pm 0.066, -1.229 \pm 0.047)$ mas/yr. These results have similar accuracy and agree to within the uncertainties with existing Hubble Space Telescope (HST) PM measurements. Since TGAS uses different methods with different systematics, this provides an external validation of both data sets and their underlying approaches. Residual DR1 systematics may affect the TGAS results, but the HST agreement implies this must be below the random errors. Also in agreement with prior HST studies, the TGAS LMC PM field clearly shows the clockwise rotation of the disk, even though it takes the LMC disk in excess of $10^8$ years to complete one revolution. The implied rotation curve amplitude for young LMC stars is consistent with that inferred from line-of-sight (LOS) velocity measurements. Comparison of the PM and LOS rotation curves implies a kinematic LMC distance modulus $m-M = 18.54 \pm 0.39$, consistent but not yet competitive with photometric methods. These first results from Gaia on the topic of Local Group dynamics provide an indication of how its future data releases will revolutionize this field.


title: Quantum Speed Limit Bounds in an Open Quantum Evolution

abstract: Quantum mechanics dictates bounds for the minimal evolution time between predetermined initial and final states. Several of these Quantum Speed Limit (QSL) bounds were derived for non-unitary dynamics using different approaches. Here, we perform a systematic analysis of the most common QSL bounds in the damped Jaynes-Cummings model, covering the Markovian and non-Markovian regime. We show that only one of the analysed bounds cleaves to the essence of the QSL theory outlined in the pioneer works of Mandelstam \& Tamm and Margolus \& Levitin in the context of unitary evolutions. We also show that all of QSL bounds analysed reflect the fact that in our model non-Markovian effects speed up the quantum evolution. However, it is not possible to infer the Markovian or non-Markovian behaviour of the dynamics only analysing the QSL bounds.


title: Soft Black Hole Absorption Rates as Conservation Laws

abstract: The absorption rate of low-energy, or soft, electromagnetic radiation by spherically symmetric black holes in arbitrary dimensions is shown to be fixed by conservation of energy and large gauge transformations. We interpret this result as the explicit realization of the Hawking-Perry-Strominger Ward identity for large gauge transformations in the background of a non-evaporating black hole. Along the way we rederive and extend previous analytic results regarding the absorption rate for the minimal scalar and the photon.


title: Star Formation Quenching Timescale of Central Galaxies in a Hierarchical Universe

abstract: Central galaxies make up the majority of the galaxy population, including the majority of the quiescent population at $\mathcal{M}_* > 10^{10}\mathrm{M}_\odot$. Thus, the mechanism(s) responsible for quenching central galaxies plays a crucial role in galaxy evolution as whole. We combine a high resolution cosmological $N$-body simulation with observed evolutionary trends of the "star formation main sequence," quiescent fraction, and stellar mass function at $z < 1$ to construct a model that statistically tracks the star formation histories and quenching of central galaxies. Comparing this model to the distribution of central galaxy star formation rates in a group catalog of the SDSS Data Release 7, we constrain the timescales over which physical processes cease star formation in central galaxies. Over the stellar mass range $10^{9.5}$ to $10^{11} \mathrm{M}_\odot$ we infer quenching e-folding times that span $1.5$ to $0.5\; \mathrm{Gyr}$ with more massive central galaxies quenching faster. For $\mathcal{M}_* = 10^{10.5}\mathrm{M}_\odot$, this implies a total migration time of $\sim 4~\mathrm{Gyrs}$ from the star formation main sequence to quiescence. Compared to satellites, central galaxies take $\sim 2~\mathrm{Gyrs}$ longer to quench their star formation, suggesting that different mechanisms are responsible for quenching centrals versus satellites. Finally, the central galaxy quenching timescale we infer provides key constraints for proposed star formation quenching mechanisms. Our timescale is generally consistent with gas depletion timescales predicted by quenching through strangulation. However, the exact physical mechanism(s) responsible for this still remain unclear.


title: Stellar clusterings around "Isolated" Massive YSOs in the LMC

abstract: Observations suggest that there is a significant fraction of O-stars in the field of the Milky Way that appear to have formed in isolation or in low mass clusters ($<$100 $M_\odot$). The existence of these high-mass stars that apparently formed in the field challenges the generally accepted paradigm, which requires star formation to occur in clustered environments. In order to understand the physical conditions for the formation of these stars, it is necessary to observe isolated high-mass stars while they are still forming. With the $Hubble$ $Space$ $Telescope$, we observe the seven most isolated massive ($>$8 $M_\odot$) young stellar objects (MYSOs) in the Large Magellanic Cloud (LMC). The observations show that while these MYSOs are remote from other MYSOs, OB associations, and even from known giant molecular clouds, they are actually not isolated at all. Imaging reveals $\sim$100 to several hundred pre--main-sequence (PMS) stars in the vicinity of each MYSO. These previously undetected PMS stars form prominent compact clusters around the MYSOs, and in most cases they are also distributed sparsely across the observed regions. Contrary to what previous high-mass field star studies show, these observations suggest that high-mass stars may not be able to form in clusters with masses less than 100 $M_\odot$. If these MYSOs are indeed the best candidates for isolated high-mass star formation, then the lack of isolation is at odds with random sampling of the IMF. Moreover, while isolated MYSOs may not exist, we find evidence that isolated clusters containing O-stars can exist, which in itself is rare.


title: Extending Semi-numeric Reionisation Models to the First Stars and Galaxies

abstract: Semi-numeric methods have made it possible to efficiently model the epoch of reionisation (EoR). While most implementations involve a reduction to a simple three-parameter model, we introduce a new mass-dependent ionising efficiency parameter that folds in physical parameters that are constrained by the latest numerical simulations. This new parameterization enables the effective modeling of a broad range of host halo masses containing ionising sources, extending from the smallest Population III host halos with $M \sim 10^6 M_\odot$, which are often ignored, to the rarest cosmic peaks with $M \sim 10^{12} M_\odot$ during EoR. We compare the resulting ionising histories with a typical three-parameter model and also compare with the latest constraints from the Planck mission. Our model results in a optical depth due to Thomson scattering, $\tau_{\mathrm{e}}$ = 0.057, that is consistent with Planck. The largest difference in our model is shown in the resulting bubble size distributions which peak at lower characteristic sizes and are broadened. We also consider the uncertainties of the various physical parameters and comparing the resulting ionising histories broadly disfavors a small contribution from galaxies. As the smallest haloes cease a meaningful contribution to the ionising photon budget after $z = 10$, implying they play a role in determining the start of EoR and little else.


title: Spherical Harmonic Analyses of Intensity Mapping Power Spectra

abstract: Intensity mapping is a promising technique for surveying the large scale structure of our Universe from $z=0$ to $z \sim 150$, using the brightness temperature field of spectral lines to directly observe previously unexplored portions of out cosmic timeline. Examples of targeted lines include the $21\,\textrm{cm}$ hyperfine transition of neutral hydrogen, rotational lines of carbon monoxide, and fine structure lines of singly ionized carbon. Recent efforts have focused on detections of the power spectrum of spatial fluctuations, but have been hindered by systematics such as foreground contamination. This has motivated the decomposition of data into Fourier modes perpendicular and parallel to the line-of-sight, which has been shown to be a particularly powerful way to diagnose systematics. However, such a method is well-defined only in the limit of a narrow-field, flat-sky approximation. This limits the sensitivity of intensity mapping experiments, as it means that wide surveys must be separately analyzed as a patchwork of smaller fields. In this paper, we develop a framework for analyzing intensity mapping data in a spherical Fourier-Bessel basis, which incorporates curved sky effects without difficulty. We use our framework to generalize a number of techniques in intensity mapping data analysis from the flat sky to the curved sky. These include visibility-based estimators for the power spectrum, treatments of interloper lines, and the "foreground wedge" signature of spectrally smooth foregrounds.


title: The Most Massive Galaxies and Black Holes Allowed by $\Lambda$CDM

abstract: Given a galaxy's stellar mass, its host halo mass has a lower limit from the cosmic baryon fraction and known baryonic physics. At z>4, galaxy stellar mass functions place lower limits on halo number densities that approach expected $\Lambda$CDM halo mass functions. High-redshift galaxy stellar mass functions can thus place interesting limits on number densities of massive haloes, which are otherwise very difficult to measure. Although halo mass functions at z<8 are consistent with observed galaxy stellar masses if galaxy baryonic conversion efficiencies increase with redshift, JWST and WFIRST will more than double the redshift range over which useful constraints are available. We calculate maximum galaxy stellar masses as a function of redshift given expected halo number densities from $\Lambda$CDM. We apply similar arguments to black holes. If their virial mass estimates are accurate, number density constraints alone suggest that the quasars SDSS J1044-0125 and SDSS J010013.02+280225.8 likely have black hole mass -- stellar mass ratios higher than the median z=0 relation, confirming the expectation from Lauer bias. Finally, we present a public code to evaluate the probability of an apparently $\Lambda$CDM-inconsistent high-mass halo being detected given the combined effects of multiple surveys and observational errors.


title: The Long Term Evolution of ASASSN-14li

abstract: We present late-time optical spectroscopy taken with the Large Binocular Telescope's Multi-Object Double Spectrograph, late-time SWIFT UVOT and XRT observations, as well as improved ASAS-SN pre-discovery limits on the nearby (d=90.3 Mpc, z=0.0206) tidal disruption event (TDE) ASASSN-14li. The late-time optical spectra show H$\alpha$ emission well in excess of that seen in the SDSS host galaxy spectrum, indicating that the processes powering the luminous flares associated with TDEs can operate for several hundreds of days. The SWIFT observations reveal the presence of lingering apparently thermal UV (T$_{\rm UV} \sim 3.5\times10^4$ K) and X-ray (T$_{\rm X} \sim 7\times10^5$ K) emission. The characteristic temperatures evolve by, at most, a factor of $\sim2$ over the 600 day follow-up campaign. The X-ray, UV, and H$\alpha$ luminosities evolve roughly in tandem and at a rate that is consistent with a power-law decay at late times. This behavior is in stark contrast with the majority of optically discovered TDEs, which are X-ray faint and evolve on shorter timescales. Finally we address how the unique properties of ASASSN-14li can be used to probe the relationship between the TDE rate and host galaxy properties.


title: Geometric constraints on the space of N=2 SCFTs III: enhanced Coulomb branches and central charges

abstract: This is the third in a series of three papers on the systematic analysis of rank 1 four dimensional $\mathcal{N}=2$ SCFTs. In the first two papers we developed and carried out a strategy for classifying and constructing physical planar rank-1 Coulomb branch geometries of $\mathcal{N}=2$ SCFTs. Here we describe general features of the Higgs and mixed branch geometries of the moduli space of these SCFTs, and use this, along with their Coulomb branch geometry, to compute their conformal and flavor central charges. We conclude with a summary of the state of the art for rank-1 $\mathcal{N}=2$ SCFTs.


title: Are Cosmological Gas Accretion Streams Multiphase and Turbulent?

abstract: Simulations of cosmological filamentary accretion reveal flows ("streams") of warm gas, ~$10^4$ K, which are efficient in bringing gas into galaxies. We present a phenomenological scenario where gas in such flows -- if it is shocked as it enters the halo as we assume -- become biphasic and, as a result, turbulent. We consider a collimated stream of warm gas that flows into a halo from an over dense filament of the cosmic web. The post-shock streaming gas expands because it has a higher pressure than the ambient halo gas, and fragments as it cools. The fragmented stream forms a two phase medium: a warm cloudy phase embedded in hot post-shock gas. We argue that the hot phase sustains the accretion shock. A fraction of the initial kinetic energy of the infalling gas is converted into turbulence among and within the warm clouds. The thermodynamic evolution of the post-shock gas is largely determined by the relative timescales of several processes -- the cooling, the expansion of the post-shock gas, the amount of turbulence in the clouds, and the halo dynamics. We expect the gas to become multiphase when the cooling and dynamical times are of the same order-of-magnitude. In this framework, we show that this occurs in the important mass range of ~$10^{11}$ to $10^{13}$ M$_\odot$ , where the bulk of stars have formed in galaxies. Gas accreting along cosmic web filaments may eventually lose coherence and mix with the ambient halo gas. Through both the phase separation and "disruption" of the stream, the accretion efficiency onto a galaxy in a halo dynamical time is lowered. De-collimating flows make the direct interaction between galaxy feedback and accretion streams more likely, thereby further reducing the overall accretion efficiency. Moderating the gas accretion efficiency through these mechanisms may help to alleviate a number of significant challenges in theoretical galaxy formation.


title: Vortices and Vermas

abstract: In three-dimensional gauge theories, monopole operators create and destroy vortices. We explore this idea in the context of 3d N=4 gauge theories in the presence of an Omega-background. In this case, monopole operators generate a non-commutative algebra that quantizes the Coulomb-branch chiral ring. The monopole operators act naturally on a Hilbert space, which is realized concretely as the equivariant cohomology of a moduli space of vortices. The action furnishes the space with the structure of a Verma module for the Coulomb-branch algebra. This leads to a new mathematical definition of the Coulomb-branch algebra itself, related to that of Braverman-Finkelberg-Nakajima. By introducing additional boundary conditions, we find a construction of vortex partition functions of 2d N=(2,2) theories as overlaps of coherent states (Whittaker vectors) for Coulomb-branch algebras, generalizing work of Braverman-Feigin-Finkelberg-Rybnikov on a finite version of the AGT correspondence. In the case of 3d linear quiver gauge theories, we use brane constructions to exhibit vortex moduli spaces as handsaw quiver varieties, and realize monopole operators as interfaces between handsaw-quiver quantum mechanics, generalizing work of Nakajima.


title: Star formation triggered by galaxy interactions in modified gravity

abstract: Together with interstellar turbulence, gravitation is one key player in star formation. It acts both at galactic scales in the assembly of gas into dense clouds, and inside those structures for their collapse and the formation of pre-stellar cores. To understand to what extent the large scale dynamics govern the star formation activity of galaxies, we present hydrodynamical simulations in which we generalise the behaviour of gravity to make it differ from Newtonian dynamics in the low acceleration regime. We focus on the extreme cases of interacting galaxies, and compare the evolution of galaxy pairs in the dark matter paradigm to that in the Milgromian Dynamics (MOND) framework. Following up on the seminal work by Tiret & Combes, this paper documents the first simulations of galaxy encounters in MOND with a detailed Eulerian hydrodynamical treatment of baryonic physics, including star formation and stellar feedback. We show that similar morphologies of the interacting systems can be produced by both the dark matter and MOND formalisms, but require a much slower orbital velocity in the MOND case. Furthermore, we find that the star formation activity and history are significantly more extended in space and time in MOND interactions, in particular in the tidal debris. Such differences could be used as observational diagnostics and make interacting galaxies prime objects in the study of the nature of gravitation at galactic scales.


title: Unbounded memory advantage in stochastic simulation using quantum mechanics

abstract: Simulating the stochastic evolution of real quantities on a digital computer requires a trade-off between the precision to which these quantities are approximated, and the memory required to store them. The statistical accuracy of the simulation is thus generally limited by the internal memory available to the simulator. Here, using tools from computational mechanics, we show that quantum processors with a fixed finite memory can simulate stochastic processes of real variables to arbitrarily high precision. This demonstrates a provable, unbounded memory advantage that a quantum simulator can exhibit over its best possible classical counterpart.


title: The Razor's Edge of Collapse: The Transition Point from Lognormal to Powerlaw in Molecular Cloud PDFs

abstract: We derive an analytic expression for the transitional column density value ($s_t$) between the lognormal and power-law form of the probability distribution function (PDF) in star-forming molecular clouds. Our expression for $s_t$ depends on the mean column density, the variance of the lognormal portion of the PDF, and the slope of the power-law portion of the PDF. We show that $s_t$ can be related to physical quantities such as the sonic Mach number of the flow and the power-law index for a self-gravitating isothermal sphere. This implies that the transition point between the lognormal and power-law density/column density PDF represents the critical density where turbulent and thermal pressure balance, the so-called "post-shock density." We test our analytic prediction for the transition column density using dust PDF observations reported in the literature as well as numerical MHD simulations of self-gravitating supersonic turbulence with the Enzo code. We find excellent agreement between the analytic $s_t$ and the measured values from the numerical simulations and observations (to within 1.5 A$_V$). We discuss the utility of our expression for determining the properties of the PDF from unresolved low density material in dust observations, for estimating the post-shock density, and for determining the HI-H$_2$ transition in clouds.


title: On-sky single-mode fiber coupling measurements at the Large Binocular Telescope

abstract: The demonstration of efficient single-mode fiber (SMF) coupling is a key requirement for the development of a compact, ultra-precise radial velocity (RV) spectrograph. iLocater is a next generation instrument for the Large Binocular Telescope (LBT) that uses adaptive optics (AO) to inject starlight into a SMF. In preparation for commissioning iLocater, a prototype SMF injection system was installed and tested at the LBT in the Y-band (0.970-1.065 $\mu$m). This system was designed to verify the capability of the LBT AO system as well as characterize on-sky SMF coupling efficiencies. SMF coupling was measured on stars with variable airmasses, apparent magnitudes, and seeing conditions for six half-nights using the Large Binocular Telescope Interferometer. We present the overall optical and mechanical performance of the SMF injection system, including details of the installation and alignment procedure. A particular emphasis is placed on analyzing the instrument's performance as a function of telescope elevation to inform the final design of the fiber injection system for iLocater.


title: The iLocater cryostat: design and thermal control strategy for precision radial velocity measurements

abstract: The current generation of precision radial velocity (RV) spectrographs are seeing-limited instruments. In order to achieve high spectral resolution on 8m class telescopes, these spectrographs require large optics and in turn, large instrument volumes. Achieving milli-Kelvin thermal stability for these systems is challenging but is vital in order to obtain a single measurement RV precision of better than 1m/s. This precision is crucial to study Earth-like exoplanets within the habitable zone. iLocater is a next generation RV instrument being developed for the Large Binocular Telescope. Unlike seeing-limited RV instruments, iLocater uses adaptive optics (AO) to inject a diffraction-limited beam into single-mode fibers. These fibers illuminate the instrument spectrograph, facilitating a diffraction-limited design and a small instrument volume compared to present-day instruments. This enables intrinsic instrument stability and facilitates precision thermal control. We present the current design of the iLocater cryostat which houses the instrument spectrograph and the strategy for its thermal control. The spectrograph is situated within a pair of radiation shields mounted inside an MLI lined vacuum chamber. The outer radiation shield is actively controlled to maintain instrument stability at the sub-mK level and minimize effects of thermal changes from the external environment. An inner shield passively dampens any residual temperature fluctuations and is radiatively coupled to the optical board. To provide intrinsic stability, the optical board and optic mounts will be made from Invar and cooled to 58K to benefit from a zero coefficient of thermal expansion (CTE) value at this temperature. Combined, the small footprint of the instrument spectrograph, the use of Invar, and precision thermal control will allow long-term sub-milliKelvin stability to facilitate precision RV measurements.


title: iLocater: A Diffraction-limited Doppler Spectrometer for the Large Binocular Telescope

abstract: We are developing a stable and precise spectrograph for the Large Binocular Telescope (LBT) named "iLocater." The instrument comprises three principal components: a cross-dispersed echelle spectrograph that operates in the YJ-bands (0.97-1.30 microns), a fiber-injection acquisition camera system, and a wavelength calibration unit. iLocater will deliver high spectral resolution (R~150,000-240,000) measurements that permit novel studies of stellar and substellar objects in the solar neighborhood including extrasolar planets. Unlike previous planet-finding instruments, which are seeing-limited, iLocater operates at the diffraction limit and uses single mode fibers to eliminate the effects of modal noise entirely. By receiving starlight from two 8.4m diameter telescopes that each use "extreme" adaptive optics (AO), iLocater shows promise to overcome the limitations that prevent existing instruments from generating sub-meter-per-second radial velocity (RV) precision. Although optimized for the characterization of low-mass planets using the Doppler technique, iLocater will also advance areas of research that involve crowded fields, line-blanketing, and weak absorption lines.


title: Three-dimensional distribution of ejecta in Supernova 1987A at 10 000 days

abstract: Due to its proximity, SN 1987A offers a unique opportunity to directly observe the geometry of a stellar explosion as it unfolds. Here we present spectral and imaging observations of SN 1987A obtained ~10,000 days after the explosion with HST/STIS and VLT/SINFONI at optical and near-infrared wavelengths. These observations allow us to produce the most detailed 3D map of H-alpha to date, the first 3D maps for [Ca II] \lambda \lambda 7292, 7324, [O I] \lambda \lambda 6300, 6364 and Mg II \lambda \lambda 9218, 9244, as well as new maps for [Si I]+[Fe II] 1.644 \mu m and He I 2.058 \mu m. A comparison with previous observations shows that the [Si I]+[Fe II] flux and morphology have not changed significantly during the past ten years, providing evidence that it is powered by 44Ti. The time-evolution of H-alpha shows that it is predominantly powered by X-rays from the ring, in agreement with previous findings. All lines that have sufficient signal show a similar large-scale 3D structure, with a north-south asymmetry that resembles a broken dipole. This structure correlates with early observations of asymmetries, showing that there is a global asymmetry that extends from the inner core to the outer envelope. On smaller scales, the two brightest lines, H-alpha and [Si I]+[Fe II] 1.644 \mu m, show substructures at the level of ~ 200 - 1000 km/s and clear differences in their 3D geometries. We discuss these results in the context of explosion models and the properties of dust in the ejecta.


title: A convergent $\frac{1}{N}$ expansion for GUE

abstract: We show that the asymptotic $1/N$ expansion for the averages of linear statistics of the GUE is convergent when the test function is an entire function of order two and finite type. This allows to fully recover the mean eigenvalue density function for finite $N$ from the coefficients of the expansion thus providing a resummation procedure. As an intermediate result we compute the bilateral Laplace transform of the GUE reproducing kernel in the half-sum variable, generalizing a formula of Haagerup and Thorbj{\o}rnsen.


title: Fabrication of fully automated setup for high temperature thermal conductivity measurement

abstract: In this work, we report the fabrication of experimental setup for high temperature thermal conductivity ($\kappa$) measurement. It can characterize samples with various dimensions and shapes. Steady state based axial heat flow technique is used for $\kappa$ measurement. Heat loss is measured using parallel thermal conductance technique. Simple design, lightweight and small size sample holder is developed by using a thin heater and limited components. Low heat loss value is achieved by using very low thermal conductive insulator block with small cross sectional area. Power delivered to the heater is measured accurately by using 4-wire technique and for this, heater is developed with 4-wire. This setup is validated by using $Bi_{0.36}Sb_{1.45}Te_3$, polycrystalline bismuth, gadolinium, and alumina samples. The data obtained for these samples are found to be in good agreement with the reported data. The maximum deviation of 6 \% in the value $\kappa$ are observed. This maximum deviation is observed with the gadolinium sample. We also report the thermal conductivity of polycrystalline tellurium from 320 K to 550 K and the non monotonous behavior of $\kappa$ with temperature is observed.


title: No hydrogen exosphere detected around the super-Earth HD97658 b

abstract: The exoplanet HD97658b provides a rare opportunity to probe the atmospheric composition and evolution of moderately irradiated super-Earths. It transits a bright K star at a moderate orbital distance of 0.08 au. Its low density is compatible with a massive steam envelope that could photodissociate at high altitudes and become observable as escaping hydrogen. Our analysis of 3 transits with HST/STIS at Ly-alpha reveals no such signature, suggesting that the thermosphere is not hydrodynamically expanding and is subjected to a low escape of neutral hydrogen (<10^8 g/s at 3 sigma). Using HST Ly-alpha and Chandra & XMM-Newton observations at different epochs, we find that HD97658 is a weak and soft X-ray source with signs of chromospheric variability in the Ly-alpha line core. We determine an average reference for the intrinsic Ly-alpha line and XUV spectrum of the star, and show that HD97658 b is in mild conditions of irradiation compared to other known evaporating exoplanets with an XUV irradiation about 3 times lower than the evaporating warm Neptune GJ436 b. This could be why the thermosphere of HD97658b is not expanding: the low XUV irradiation prevents an efficient photodissociation of any putative steam envelope. Alternatively, it could be linked to a low hydrogen content or inefficient conversion of the stellar energy input. The HD97658 system provides clues for understanding the stability of low-mass planet atmospheres. Our study of HD97658 b can be seen as a control experiment of our methodology, confirming that it does not bias detections of atmospheric escape and underlining its strength and reliability. Our results show that stellar activity can be efficiently discriminated from absorption signatures by a transiting exospheric cloud. They also highlight the potential of observing the upper atmosphere of small transiting planets to probe their physical and chemical properties


title: An Adaptive Psychoacoustic Model for Automatic Speech Recognition

abstract: Compared with automatic speech recognition (ASR), the human auditory system is more adept at handling noise-adverse situations, including environmental noise and channel distortion. To mimic this adeptness, auditory models have been widely incorporated in ASR systems to improve their robustness. This paper proposes a novel auditory model which incorporates psychoacoustics and otoacoustic emissions (OAEs) into ASR. In particular, we successfully implement the frequency-dependent property of psychoacoustic models and effectively improve resulting system performance. We also present a novel double-transform spectrum-analysis technique, which can qualitatively predict ASR performance for different noise types. Detailed theoretical analysis is provided to show the effectiveness of the proposed algorithm. Experiments are carried out on the AURORA2 database and show that the word recognition rate using our proposed feature extraction method is significantly increased over the baseline. Given models trained with clean speech, our proposed method achieves up to 85.39% word recognition accuracy on noisy data.


title: Synchronous Lagrangian variational principles in General Relativity

abstract: The problem of formulating synchronous variational principles in the context of General Relativity is discussed. Based on the analogy with classical relativistic particle dynamics, the existence of variational principles is pointed out in relativistic classical field theory which are either asynchronous or synchronous. The historical Einstein-Hilbert and Palatini variational formulations are found to belong to the first category. Nevertheless, it is shown that an alternative route exists which permits one to cast these principles in terms of equivalent synchronous Lagrangian variational formulations. The advantage is twofold. First, synchronous approaches allow one to overcome the lack of gauge symmetry of the asynchronous principles. Second, the property of manifest covariance of the theory is also restored at all levels, including the symbolic Euler-Lagrange equations, with the variational Lagrangian density being now identified with a $4-$scalar. As an application, a joint synchronous variational principle holding both for the non-vacuum Einstein and Maxwell equations is displayed, with the matter source being described by means of a Vlasov kinetic treatment.


title: Inhibitory loop robustly induces anticipated synchronization in neuronal microcircuits

abstract: We investigate the synchronization properties between two excitatory coupled neurons in the presence of an inhibitory loop mediated by an interneuron. Dynamical inhibition together with noise independently applied to each neuron provide phase diversity in the dynamics of the neuronal motif. We show that the interplay between the coupling strengths and external noise controls the phase relations between the neurons in a counter-intuitive way. For a master-slave configuration (unidirectional coupling) we find that the slave can anticipate the master, on average, if the slave is subject to the inhibitory feedback. In this non-usual regime, called anticipated synchronization (AS), the phase of the post-synaptic neuron is advanced with respect to that of the pre-synaptic neuron. We also show that the AS regime survives even in the presence of unbalanced bidirectional excitatory coupling. Moreover, for the symmetric mutually coupled situation, the neuron that is subject to the inhibitory loop leads in phase.


title: On a Local Version of the Bak-Sneppen Model

abstract: A major difficulty in studying the Bak-Sneppen model is in effectively comparing it with well-understood models. This stems from the use of two geometries: complete graph geometry to locate the global fitness minimizer, and graph geometry to replace the species in the neighborhood of the minimizer. We present a variant in which only the graph geometry is used. This allows to obtain the stationary distribution through random walk dynamics. We use this to show that for constant-degree graphs, the stationary fitness distribution converges to an IID law as the number of vertices tends to infinity. We also discuss exponential ergodicity through coupling, and avalanches for the model.


title: Scaling of tripartite entanglement at impurity quantum phase transitions

abstract: The emergence of a diverging length scale in many-body systems at a quantum phase transition implies that total entanglement has to reach its maximum there. In order to fully characterize this, one has to consider multipartite entanglement as, for instance, bipartite entanglement between individual particles fails to signal this effect. However, quantification of multipartite entanglement is very hard and detecting it may not be possible due to the lack of accessibility to all individual particles. For these reasons it will be more sensible to partition the system into relevant subsystems, each containing few to many spins, and study entanglement between those constituents as a coarse-grain picture of multipartite entanglement between individual particles. In impurity systems, famously exemplified by two-impurity and two-channel Kondo models, it is natural to divide the system into three parts, namely, impurities and the left and right bulks. By exploiting two tripartite entanglement measures, based on negativity, we show that at impurity quantum phase transitions the tripartite entanglement diverges and shows scaling behavior. While the critical exponents are different for each tripartite entanglement measure they both provide very similar critical exponents for the two-impurity and the two-channel Kondo models suggesting that they belong to the same universality class.


title: Manifest Covariant Hamiltonian Theory of General Relativity

abstract: The problem of formulating a manifest covariant Hamiltonian theory of General Relativity in the presence of source fields is addressed, by extending the so-called "DeDonder-Weyl" formalism to the treatment of classical fields in curved space-time. The theory is based on a synchronous variational principle for the Einstein equation, formulated in terms of superabundant variables. The technique permits one to determine the continuum covariant Hamiltonian structure associated with the Einstein equation. The corresponding continuum Poisson bracket representation is also determined. The theory relies on first-principles, in the sense that the conclusions are reached in the framework of a non-perturbative covariant approach, which allows one to preserve both the 4-scalar nature of Lagrangian and Hamiltonian densities as well as the gauge invariance property of the theory.


title: The Close AGN Reference Survey (CARS) - Mrk 1018 returns to the shadows after 30 years as a Seyfert 1

abstract: We report the discovery that the known `changing look' AGN Mrk 1018 has changed spectral type for a second time. New VLT-MUSE data taken in 2015 as part of the Close AGN Reference Survey (CARS) shows that the AGN has returned to its original Seyfert 1.9 classification. The CARS sample is selected to contain only bright type 1 AGN, but Mrk 1018's broad emission lines and continuum, typical of type 1 AGN, have almost entirely disappeared. We use spectral fitting of the MUSE spectrum and previously available spectra to determine the drop in broad line flux and the Balmer decrement. We find that the broad line flux has decreased by a factor of 4.75+/- 0.5 in H{\alpha} since an SDSS spectrum was taken in 2000. The Balmer decrement has not changed significantly implying no enhanced reddening with time, but the remaining broad lines are more asymmetric than those present in the type 1 phase. We posit that the change is due to an intrinsic drop in flux from the accretion disk rather than variable extinction or a tidal disruption event.


title: Ornstein-Uhlenbeck approximation of one-step processes: a differential equation approach

abstract: The steady state of the Fokker-Planck equation corresponding to a density dependent one-step process is approximated by a suitable normal distribution. Starting from the master equations of the process, written in terms of the time dependent probabilities, $p_k(t)$ of the states $k=0,1,\ldots , N$, their continuous (in space) version, the Fokker-Planck equation is formulated. This PDE approximation enables us to create analytic approximation formulas for the steady state distribution. These formulas are derived based on heuristic reasoning and then their accuracy is proved to be of order $1/N^{\beta}$ with some power $\beta <1$.


title: The Close AGN Reference Survey (CARS) - What is causing Mrk1018's return to the shadows after 30 years?

abstract: We recently discovered that the active galactic nucleus (AGN) of Mrk 1018 has changed optical type again after 30 years as a type 1 AGN. Here we combine Chandra, NuStar, Swift, Hubble Space Telescope and ground-based observations to explore the cause of this change. The 2-10keV flux declines by a factor of ~8 between 2010 and 2016. We show with our X-ray observation that this is not caused by varying neutral hydrogen absorption along the line-of-sight up to the Compton-thick level. The optical-UV spectral energy distributions are well fit with a standard geometrically thin optically thick accretion disc model that seems to obey the expected $L\sim T^4$ relation. It confirms that a decline in accretion disc luminosity is the primary origin for the type change. We detect a new narrow-line absorber in Lya blue-shifted by ~700km/s with respect to the systemic velocity of the galaxy. This new Lya absorber could be evidence for the onset of an outflow or a companion black hole with associated gas that could be related to the accretion rate change. However, the low column density of the absorber means that it is not the direct cause for Mrk 1018's changing-look nature.


title: Hamiltonian approach to GR - Part 1: covariant theory of classical gravity

abstract: A challenging issue in General Relativity concerns the determination of the manifestly-covariant continuum Hamiltonian structure underlying the Einstein field equations and the related formulation of the corresponding covariant Hamilton-Jacobi theory. The task is achieved by adopting a synchronous variational principle requiring distinction between the prescribed deterministic metric tensor $\widehat{g}(r)\equiv \left\{ \widehat{g}_{\mu \nu }(r)\right\} $ solution of the Einstein field equations which determines the geometry of the background space-time and suitable variational fields $x\equiv \left\{ g,\pi \right\} $ obeying an appropriate set of continuum Hamilton equations, referred to here as GR-Hamilton equations$.$ It is shown that a prerequisite for reaching such a goal is that of casting the same equations in evolutionary form by means of a Lagrangian parametrization for a suitably-reduced canonical state. As a result, the corresponding Hamilton-Jacobi theory is established in manifestly-covariant form. Physical implications of the theory are discussed. These include the investigation of the structural stability of the GR-Hamilton equations with respect to vacuum solutions of the Einstein equations, assuming that wave-like perturbations are governed by the canonical evolution equations.


title: Connecting the quantum and classical worlds

abstract: By considering (non-relativistic) quantum mechanics as it is done in practice in particular in condensed-matter physics, it is argued that a deterministic, unitary time evolution within a chosen Hilbert space always has a limited scope, leaving a lot of room for embedding the quantum-classical transition into our current theories without recurring to difficult-to-accept interpretations of quantum mechanics. Nonunitary projections to initial and final states, the breaking of time-reversal symmetry, a change of Hilbert space, and the introduction of classical concepts such as external potentials or localized atomic nuclei are widespread in quantum mechanical calculations. Furthermore, quantum systems require classical environments that enable the symmetry breaking that is necessary for creating the atomic configurations of molecules and crystals. This paper argues that such classical environments are provided by finite-temperature macroscopic systems in which the range of quantum correlations and entanglement is limited. This leads to classical behavior on larger scales, and to collapse-like events in all dynamical processes that become coupled to the thermalized degrees of freedom.


title: Hamiltonian approach to GR - Part 2: covariant theory of quantum gravity

abstract: A non-perturbative quantum field theory of General Relativity is presented which leads to a new realization of the theory of Covariant Quantum-Gravity (CQG-theory). The treatment is founded on the recently-identified Hamiltonian structure associated with the classical space-time, i.e., the corresponding manifestly-covariant Hamilton equations and the related Hamilton-Jacobi theory. The quantum Hamiltonian operator and the CQG-wave equation for the corresponding CQG-state and wave-function are realized in $% 4-$scalar form. The new quantum wave equation is shown to be equivalent to a set of quantum hydrodynamic equations which warrant the consistency with the classical GR Hamilton-Jacobi equation in the semiclassical limit. A perturbative approximation scheme is developed, which permits the adoption of the harmonic oscillator approximation for the treatment of the Hamiltonian potential. As an application of the theory, the stationary vacuum CQG-wave equation is studied, yielding a stationary equation for the CQG-state in terms of the $4-$scalar invariant-energy eigenvalue associated with the corresponding approximate quantum Hamiltonian operator. The conditions for the existence of a discrete invariant-energy spectrum are pointed out. This yields a possible estimate for the graviton mass together with a new interpretation about the quantum origin of the cosmological constant.


title: Extended Harmonic Map Equations and the Chaotic Soliton Solutions

abstract: In this paper, the theory of harmonic maps is extended. The soliton or traveling wave solutions of Euler's equations of the extended harmonic maps are studied. In certain cases, the chaotic behaviors of these partial equations can be found for the particular case of the metrics and the potential functions of the extended harmonic equations.


title: Soft-Collinear Supersymmetry

abstract: Soft-Collinear Effective Theory (SCET) is a framework for modeling the infrared structure of theories whose long distance behavior is dominated by soft and collinear divergences. This paper demonstrates that SCET can be made compatible with supersymmetry (SUSY). Explicitly, the effective Lagrangian for $\mathcal{N} = 1$ SUSY Yang-Mills is constructed and shown to be a complete description for the infrared of this model. For contrast, we also construct the effective Lagrangian for chiral SUSY theories with Yukawa couplings, specifically the single flavor Wess-Zumino model. Only a subset of the infrared divergences are reproduced by the Lagrangian -- to account for the complete low energy description requires the inclusion of local operators. SCET is formulated by expanding fields along a light-like direction and then subsequently integrating out degrees-of-freedom that are away from the light-cone. Defining the theory with respect to a specific frame obfuscates Lorentz invariance -- given that SUSY is a space-time symmetry, this presents a possible obstruction. The cleanest language with which to expose the congruence between SUSY and SCET requires exploring two novel formalisms: collinear fermions as two-component Weyl spinors, and SCET in light-cone gauge. By expressing SUSY Yang-Mills in "collinear superspace", a slice of superspace derived by integrating out half the fermionic coordinates, the light-cone gauge SUSY SCET theory can be written in terms of superfields. As a byproduct, bootstrapping up to the full theory yields the first algorithmic approach for determining the SUSY Yang-Mills on-shell superspace action. This work paves the way toward discovering the effective theory for the collinear limit of $\mathcal{N} = 4$ SUSY Yang-Mills.


title: Note on the Bondal-Orlov Functors for Toric DM Stacks

abstract: We calculate explicit formulas for the general equivariant Bondal-Orlov functors on the localized K-theory groups for a crepant birational transformation of toric DM stacks. We recall some facts that the Bondal-Orlov functors give equivalences on the bounded derived categories. Applying twice of these functors we get the Seidel-Thomas spherical twists for the derived category.


title: One-loop gravity divergences in D>4 cannot all be removed

abstract: Unlike the one-loop QG divergences in D=4, which can all be transformed away, those in arbitrary higher (even) dimension cannot.


title: $L^p$ Expander Graphs

abstract: We discuss how graph expansion is related to the behavior of $L^{p}$-functions on the covering tree. We show that the non-trivial eigenvalues of the adjacency operator on aa $(q+1)$-regular graph are bounded by $q^{1/p}+q^{(p-1)/p}$ - the $L^{p}$-norm of the operator on the covering tree - if and only if properly averaged lifts of functions from the graph to the tree lie in $L^{p+\epsilon}$ for every $\epsilon>0$. We generalize the result to operators on edges and to bipartite graphs. The work is based on a combinatorial interpretation of representation-theoretic ideas.


title: A Panchromatic View of Relativistic Jets in Narrow-Line Seyfert 1 Galaxies

abstract: The discovery by the Large Area Telescope on board Fermi of variable gamma-ray emission from radio-loud narrow-line Seyfert 1 (NLSy1) galaxies revealed the presence of a possible third class of Active Galactic Nuclei (AGN) with relativistic jets in addition to blazars and radio galaxies. Considering that NLSy1 are usually hosted in spiral galaxies, this finding poses intriguing questions about the nature of these objects and the formation of relativistic jets. We report on a systematic investigation of the gamma-ray properties of a sample of radio-loud NLSy1, including the detection of new objects, using 7 years of Fermi-LAT data with the new Pass 8 event-level analysis. In addition we discuss the radio-to-very-high-energy properties of the gamma-ray emitting NLSy1, their host galaxy, and black hole mass in the context of the blazar scenario and the unification of relativistic jets at different scales.


title: Pulse reflection in a random waveguide with a turning point

abstract: We present an analysis of wave propagation and reflection in an acoustic waveguide with random sound soft boundary and a turning point. The waveguide has slowly bending axis and variable cross section. The variation consists of a slow and monotone change of the width of the waveguide and small and rapid fluctuations of the boundary, on the scale of the wavelength. These fluctuations are modeled as random. The turning point is many wavelengths away from the source, which emits a pulse that propagates toward the turning point, where it is reflected. To focus attention on this reflection, we assume that the waveguide supports a single propagating mode from the source to the turning point, beyond which all the waves are evanescent. We consider a scaling regime where scattering at the random boundary has a significant effect on the reflected pulse. In this regime scattering from the random boundary away from the turning point is negligible, while scattering from the random boundary around the turning point results in a strong, deterministic pulse deformation. The reflected pulse shape is not the same as the emitted one. It is damped, due to scattering at the boundary, and is deformed by dispersion in the waveguide. The reflected pulse also carries a random phase.


title: Bayesian Reinforcement Learning: A Survey

abstract: Bayesian methods for machine learning have been widely investigated, yielding principled methods for incorporating prior information into inference algorithms. In this survey, we provide an in-depth review of the role of Bayesian methods for the reinforcement learning (RL) paradigm. The major incentives for incorporating Bayesian reasoning in RL are: 1) it provides an elegant approach to action-selection (exploration/exploitation) as a function of the uncertainty in learning; and 2) it provides a machinery to incorporate prior knowledge into the algorithms. We first discuss models and methods for Bayesian inference in the simple single-step Bandit model. We then review the extensive recent literature on Bayesian methods for model-based RL, where prior information can be expressed on the parameters of the Markov model. We also present Bayesian methods for model-free RL, where priors are expressed over the value function or policy class. The objective of the paper is to provide a comprehensive survey on Bayesian RL algorithms and their theoretical and empirical properties.


title: Robust recovery-type a posteriori error estimators for streamline upwind/Petrov Galerkin discretizations for singularly perturbed problems

abstract: In this paper, we investigate adaptive streamline upwind/Petrov Galerkin (SUPG) methods for singularly perturbed convection-diffusion-reaction equations in a new dual norm presented in [Du and Zhang, J. Sci. Comput. (2015)]. The flux is recovered by either local averaging in conforming $H({\rm div})$ spaces or weighted global $L^2$ projection onto conforming $H({\rm div})$ spaces. We further introduce a recovery stabilization procedure, and develop completely robust a posteriori error estimators with respect to the singular perturbation parameter $\varepsilon$. Numerical experiments are reported to support the theoretical results and to show that the estimated errors depend on the degrees of freedom uniformly in $\varepsilon$.


title: Local approximation of arbitrary functions by solutions of nonlocal equations

abstract: We show that any function can be locally approximated by solutions of prescribed linear equations of nonlocal type. In particular, we show that every function is locally $s$-caloric, up to a small error. The case of non-elliptic and non-parabolic operators is taken into account as well.


title: Unrestricted State Complexity of Binary Operations on Regular and Ideal Languages

abstract: We study the state complexity of binary operations on regular languages over different alphabets. It is known that if $L'_m$ and $L_n$ are languages of state complexities $m$ and $n$, respectively, and restricted to the same alphabet, the state complexity of any binary boolean operation on $L'_m$ and $L_n$ is $mn$, and that of product (concatenation) is $m 2^n - 2^{n-1}$. In contrast to this, we show that if $L'_m$ and $L_n$ are over different alphabets, the state complexity of union and symmetric difference is $(m+1)(n+1)$, that of difference is $mn+m$, that of intersection is $mn$, and that of product is $m2^n+2^{n-1}$. We also study unrestricted complexity of binary operations in the classes of regular right, left, and two-sided ideals, and derive tight upper bounds. The bounds for product of the unrestricted cases (with the bounds for the restricted cases in parentheses) are as follows: right ideals $m+2^{n-2}+2^{n-1}$ ($m+2^{n-2}$); left ideals $mn+m+n$ ($m+n-1$); two-sided ideals $m+2n$ ($m+n-1$). The state complexities of boolean operations on all three types of ideals are the same as those of arbitrary regular languages, whereas that is not the case if the alphabets of the arguments are the same. Finally, we update the known results about most complex regular, right-ideal, left-ideal, and two-sided-ideal languages to include the unrestricted cases.


title: $\downarrow$-posets

abstract: We investigate a certain class of posets arising from semilattice actions. Let $S$ be a semilattice with identity. Let $S$ act on a set $C$. For $c,d\in C$ put $c\leq d$ iff there is some $s\in S$ with $ds=c$. Then $(C,\leq)$ is a poset. Let's call the posets that arise in this way $\downarrow$-posets. We give a reasonable second order characterization of $\downarrow$-posets and show that there is no first order characterization.


title: Long-time behavior for crystal dislocation dynamics

abstract: We describe the asymptotic states for the solutions of a nonlocal equation of evolutionary type, which have the physical meaning of the atom dislocation function in a periodic crystal. More precisely, we can describe accurately the "smoothing effect" on the dislocation function occurring slightly after a "particle collision" (roughly speaking, two opposite transitions layers average out) and, in this way, we can trap the atom dislocation function between a superposition of transition layers which, as time flows, approaches either a constant function or a single heteroclinic (depending on the algebraic properties of the orientations of the initial transition layers). The results are endowed of explicit and quantitative estimates and, as a byproduct, we show that the ODE systems of particles that governs the evolution of the transition layers does not admit stationary solutions (i.e., roughly speaking, transition layers always move).


title: Acoustic quasi-lossless transmission in arbitrary pathway of network

abstract: Acoustic metamaterials have exhibited extraordinary possibilities to manipulate the propagation of the sound wave. Up to now, it is still a challenge to control the propagation of the sound wave in an arbitrary pathway of a network. Here, we design a symmetry breaking cross-shape metamaterial comprised of Helmholtz resonant cells and a square column. The square column is eccentrically arranged. The sound wave can quasi-lossless transmit through channels along the eccentric direction with compressed spaces, which breaks through the general transmission phenomenon. This exotic propagation characteristic is verified by the band structure and the mode of the metamaterial. Two acoustic networks, including a 2x2 network and an 8x8 network, show that the sound wave can quasi-lossless propagate along various arbitrary pathways, such as the Great Wall shape, the stair step shape and the serpentine shape, by reconfiguring eccentric directions. This ability opens up a new venue to route the sound wave and exhibits promising applications, from the acoustic communication to the energy transmission.