title: Solvable Sachdev-Ye-Kitaev models in higher dimensions: from diffusion to many-body localization

abstract: Many aspects of many-body localization (MBL) transitions remain elusive so far. Here, we propose a higher-dimensional generalization of the Sachdev-Ye-Kitaev (SYK) model and show that it exhibits a MBL transition. The model on a bipartite lattice has $N$ Majorana fermions with SYK interactions on each site of the $A$ sublattice and $M$ free Majorana fermions on each site the of $B$ sublattice, where $N$ and $M$ are large and finite. For $r$$\equiv$$M/N\!<\!r_c$=1, it describes a diffusive metal exhibiting maximal chaos. Remarkably, its diffusive constant $D$ vanishes [$D$$\propto$$ (r_c-r)^{1/2}$] as $r$$\rightarrow$$r_c$, implying a dynamical transition to a MBL phase. It is further supported by numerical calculations of level statistics which changes from Wigner-Dyson ($r$$<$$r_c$) to Poisson ($r$$>$$r_c$) distributions. Note that no subdiffusive phase intervenes between diffusive and MBL phases. Moreover, the critical exponent $\nu$=0, violating the Harris criterion. Our higher-dimensional SYK model may provide a promising arena to explore exotic MBL transitions.


title: \"Uber-Gravity and the Cosmological Constant Problem

abstract: Recently, the idea of taking ensemble average over gravity models has been introduced. Based on this idea, we study the ensemble average over (effectively) all the gravity models (constructed from Ricci scalar) dubbing the name \"uber-gravity which is a {\it{fixed point}} in the model space. The \"uber-gravity has interesting universal properties, independent from the choice of basis: $i)$ it mimics Einstein-Hilbert gravity for high-curvature regime, $ii)$ it predicts stronger gravitational force for an intermediate-curvature regime, $iii)$ surprisingly, for low-curvature regime, i.e. $R<R_0$ where $R$ is Ricci scalar and $R_0$ is a given scale, the Lagrangian vanishes automatically and $iiii)$ there is a sharp transition between low- and intermediate-curvature regimes at $R=R_0$. We show that the \"uber-gravity response is robust to all values of vacuum energy, $\rho_{vac}$ when there is no other matter. So as a toy model, \"uber-gravity, gives a way to think about the hierarchy problems e.g. the cosmological constant problem. Due to the transition at $R=R_0$ there is a chance for \"uber-gravity to bypass Weinberg's no-go theorem. The cosmology of this model is also promising because of its non-trivial predictions for small curvature scales in comparison to $\Lambda$CDM model.


title: Galaxy Zoo: star-formation versus spiral arm number

abstract: Spiral arms are common features in low-redshift disc galaxies, and are prominent sites of star-formation and dust obscuration. However, spiral structure can take many forms: from galaxies displaying two strong `grand design' arms, to those with many `flocculent' arms. We investigate how these different arm types are related to a galaxy's star-formation and gas properties by making use of visual spiral arm number measurements from Galaxy Zoo 2. We combine UV and mid-IR photometry from GALEX and WISE to measure the rates and relative fractions of obscured and unobscured star formation in a sample of low-redshift SDSS spirals. Total star formation rate has little dependence on spiral arm multiplicity, but two-armed spirals convert their gas to stars more efficiently. We find significant differences in the fraction of obscured star-formation: an additional $\sim 10$ per cent of star-formation in two-armed galaxies is identified via mid-IR dust emission, compared to that in many-armed galaxies. The latter are also significantly offset below the IRX-$\beta$ relation for low-redshift star-forming galaxies. We present several explanations for these differences versus arm number: variations in the spatial distribution, sizes or clearing timescales of star-forming regions (i.e., molecular clouds), or contrasting recent star-formation histories.


title: Independence by Random Scaling

abstract: We give conditions under which a scalar random variable T can be coupled to a random scaling factor $\xi$ such that T and $\xi$T are rendered stochastically independent. A similar result is obtained for random measures. One consequence is a generalization of a result by Pitman and Yor on the Poisson-Dirichlet distribution to its negative parameter range. Another application are diffusion excursions straddling an exponential random time.


title: Quasiparticle interference and resonant states in normal and superconducting line nodal semimetals

abstract: We study impurity scattering in the normal and d-wave superconducting states of line nodal semimetals and show that, due to additional scattering phase space available for impurities on the surface, the quasiparticle interference pattern acquires an extended character instead of a discrete collection of delta function peaks. Moreover, using the T-matrix formalism, we demonstrate that the conventional behavior of a scalar impurity in a d-wave superconductor breaks down on the surface of a line nodal semimetal in the quasi flat band limit.


title: Transiting Planets with LSST III: Detection Rate per Year of Operation

abstract: The Large Synoptic Survey Telescope (LSST) will generate light curves for approximately 1 billion stars. Our previous work has demonstrated that, by the end of the LSST 10 year mission, large numbers of transiting exoplanetary systems could be recovered using the LSST "deep drilling" cadence. Here we extend our previous work to examine how the recoverability of transiting planets over a range of orbital periods and radii evolves per year of LSST operation. As specific example systems we consider hot Jupiters orbiting solar-type stars and hot Neptunes orbiting K-Dwarfs at distances from Earth of several kpc, as well as super-Earths orbiting nearby low-mass M-dwarfs. The detection of transiting planets increases steadily with the accumulation of data over time, generally becoming large (greater than 10 percent) after 4 - 6 years of operation. However, we also find that short-period (less than 2 day) hot Jupiters orbiting G-dwarfs and hot Neptunes orbiting K-dwarfs can already be discovered within the first 1 - 2 years of LSST operation.


title: Spin-dependent $\mu \to e$ conversion

abstract: The experimental sensitivity to $\mu \to e$ conversion on nuclei is expected to improve by four orders of magnitude in coming years. We consider the impact of $\mu \to e$ flavour-changing tensor and axial-vector four-fermion operators which couple to the spin of nucleons. Such operators, which have not previously been considered, contribute to $\mu \to e$ conversion in three ways: in nuclei with spin they mediate a spin-dependent transition; in all nuclei they contribute to the coherent ($A^2$-enhanced) spin-independent conversion via finite recoil effects and via loop mixing with dipole, scalar, and vector operators. We estimate the spin-dependent rate in Aluminium (the target of the upcoming COMET and Mu2e experiments), show that the loop effects give the greatest sensitivity to tensor and axial-vector operators involving first-generation quarks, and discuss the complementarity of the spin-dependent and independent contributions to $\mu \to e$ conversion


title: Curved spacetimes with local $\kappa$-Poincar\'e dispersion relation

abstract: We use our previously developed identification of dispersion relations with Hamilton functions on phase space to locally implement the $\kappa$-Poincar\'e dispersion relation in the momentum spaces at each point of a generic curved spacetime. We use this general construction to build the most general Hamiltonian compatible with spherical symmetry and the Plank-scale-deformed one such that in the local frame it reproduces the $\kappa$-Poincar\'e dispersion relation. Specializing to Planck-scale-deformed Schwarzschild geometry, we find that the photon sphere around a black hole becomes a thick shell since photons of different energy will orbit the black hole on circular orbits at different altitudes. We also compute the redshift of a photon between different observers at rest, finding that there is a Planck-scale correction to the usual redshift only if the observers detecting the photon have different masses.


title: Cheshire: An Online Algorithm for Activity Maximization in Social Networks

abstract: User engagement in social networks depends critically on the number of online actions their users take in the network. Can we design an algorithm that finds when to incentivize users to take actions to maximize the overall activity in a social network? In this paper, we model the number of online actions over time using multidimensional Hawkes processes, derive an alternate representation of these processes based on stochastic differential equations (SDEs) with jumps and, exploiting this alternate representation, address the above question from the perspective of stochastic optimal control of SDEs with jumps. We find that the optimal level of incentivized actions depends linearly on the current level of overall actions. Moreover, the coefficients of this linear relationship can be found by solving a matrix Riccati differential equation, which can be solved efficiently, and a first order differential equation, which has a closed form solution. As a result, we are able to design an efficient online algorithm, Cheshire, to sample the optimal times of the users' incentivized actions. Experiments on both synthetic and real data gathered from Twitter show that our algorithm is able to consistently maximize the number of online actions more effectively than the state of the art.


title: The metallicity dependence of WR winds

abstract: Wolf-Rayet (WR) stars are the most advanced stage in the evolution of the most massive stars. The strong feedback provided by these objects and their subsequent supernova (SN) explosions are decisive for a variety of astrophysical topics such as the cosmic matter cycle. Consequently, understanding the properties of WR stars and their evolution is indispensable. A crucial but still not well known quantity determining the evolution of WR stars is their mass-loss rate. Since the mass loss is predicted to increase with metallicity, the feedback provided by these objects and their spectral appearance are expected to be a function of the metal content of their host galaxy. This has severe implications for the role of massive stars in general and the exploration of low metallicity environments in particular. Hitherto, the metallicity dependence of WR star winds was not well studied. In this contribution, we review the results from our comprehensive spectral analyses of WR stars in environments of different metallicities, ranging from slightly super-solar to SMC-like metallicities. Based on these studies, we derived empirical relations for the dependence of the WN mass-loss rates on the metallicity and iron abundance, respectively.


title: Maximum Entropy Inferences on the Axion Mass in Models with Axion-Neutrino Interaction

abstract: In this work we use the Maximum Entropy Principle (MEP) to infer the mass of an axion which interacts to photons and neutrinos in an effective low energy theory. The Shannon entropy function to be maximized is suitably defined in terms of the axion branching ratios. We show that MEP strongly constrains the axion mass taking into account the current experimental bounds on the neutrinos masses. Assuming that the axion is massive enough to decay into all the three neutrinos and that MEP fixes all the free parameters of the model, the inferred axion mass is in the interval $0.1\ $eV$\ <m_{A}<0.2$ eV, which can be tested by forthcoming experiments such as IAXO. However, even in the case where MEP fixes just the axion mass and no other parameter, we found that $0.1$ eV $< m_A < 6.3$ eV in the DFSZ model with right-handed neutrinos. Moreover, a light axion, allowed to decay to photons and the lightest neutrino only, is determined by MEP as a viable dark matter candidate.

categories: hep-ph cond-mat.stat-mech hep-th physics.data-an


title: Accelerating gravitational wave parameter estimation with multi-band template interpolation

abstract: Parameter estimation on gravitational wave signals from compact binary coalescence (CBC) requires the evaluation of computationally intensive waveform models, typically the bottleneck in the analysis. This cost will increase further as low frequency sensitivity in later second and third generation detectors motivates the use of longer waveforms. We describe a method for accelerating parameter estimation by exploiting the chirping behaviour of the signals to sample the waveform sparsely for portions where the full frequency resolution is not required. We demonstrate that the method can reproduce the original results with a waveform mismatch of $\leq 5\times 10^{-7}$, but with a waveform generation cost up to $\sim 50$ times lower for computationally costly frequency-domain waveforms starting from below 8 Hz.


title: Questions on mod p representations of reductive p-adic groups

abstract: This is a list of questions raised by our joint work arXiv:1412.0737 and its sequels.


title: Distinguishing A Higgs-Like Dilaton Scenario With A Complete Bosonic Effective Field Theory Basis

abstract: A Higgs-like dilaton owns couplings that differ from those of the Standard Model Higgs and of a generic Composite Higgs. The complete bosonic basis for a Higgs-like dilaton is presented at the first subleading order. A comparison with the Standard Model, the Standard Model Effective Field Theory and the generic Lagrangian for the minimal $SO(5)/SO(4)$ Composite Higgs model is performed. Observables that can disentangle the different hypotheses are identified.


title: On the Expressive Power of Overlapping Architectures of Deep Learning

abstract: Expressive efficiency refers to the relation between two architectures A and B, whereby any function realized by B could be replicated by A, but there exists functions realized by A, which cannot be replicated by B unless its size grows significantly larger. For example, it is known that deep networks are exponentially efficient with respect to shallow networks, in the sense that a shallow network must grow exponentially large in order to approximate the functions represented by a deep network of polynomial size. In this work, we extend the study of expressive efficiency to the attribute of network connectivity and in particular to the effect of "overlaps" in the convolutional process, i.e., when the stride of the convolution is smaller than its filter size (receptive field). To theoretically analyze this aspect of network's design, we focus on a well-established surrogate for ConvNets called Convolutional Arithmetic Circuits (ConvACs), and then demonstrate empirically that our results hold for standard ConvNets as well. Specifically, our analysis shows that having overlapping local receptive fields, and more broadly denser connectivity, results in an exponential increase in the expressive capacity of neural networks. Moreover, while denser connectivity can increase the expressive capacity, we show that the most common types of modern architectures already exhibit exponential increase in expressivity, without relying on fully-connected layers.


title: Membranes by the Numbers

abstract: Many of the most important processes in cells take place on and across membranes. With the rise of an impressive array of powerful quantitative methods for characterizing these membranes, it is an opportune time to reflect on the structure and function of membranes from the point of view of biological numeracy. To that end, in this article, I review the quantitative parameters that characterize the mechanical, electrical and transport properties of membranes and carry out a number of corresponding order of magnitude estimates that help us understand the values of those parameters.


title: General order conditions for stochastic partitioned Runge-Kutta methods

abstract: In this paper stochastic partitioned Runge-Kutta (SPRK) methods are considered. A general order theory for SPRK methods based on stochastic B-series and multicolored, multishaped rooted trees is developed. The theory is applied to prove the order of some known methods, and it is shown how the number of order conditions can be reduced in some special cases, especially that the conditions for preserving quadratic invariants can be used as simplifying assumptions.


title: Heralded single photons based on spectral multiplexing and feed-forward control

abstract: We propose and experimentally demonstrate a novel approach to a heralded single photon source based on spectral multiplexing (SMUX) and feed-forward-based spectral manipulation of photons created by means of spontaneous parametric down-conversion in a periodically-poled LiNbO3 crystal. As a proof-of-principle, we show that our 3-mode SMUX increases the heralded single-photon rate compared to that of the individual modes without compromising the quality of the emitted single-photons. We project that by adding further modes, our approach can lead to a deterministic SPS.


title: Magnetically induced Ferroelectricity in Bi$_2$CuO$_4$

abstract: The tetragonal copper oxide Bi$_2$CuO$_4$ has an unusual crystal structure with a three-dimensional network of well separated CuO$_4$ plaquettes. This material was recently predicted to host electronic excitations with an unconventional spectrum and the spin structure of its magnetically ordered state appearing at T$_N$ $\sim$43 K remains controversial. Here we present the results of detailed studies of specific heat, magnetic and dielectric properties of Bi$_2$CuO$_4$ single crystals grown by the floating zone technique, combined with the polarized neutron scattering and high-resolution X-ray measurements. Our polarized neutron scattering data show Cu spins are parallel to the $ab$ plane. Below the onset of the long range antiferromagnetic ordering we observe an electric polarization induced by an applied magnetic field, which indicates inversion symmetry breaking by the ordered state of Cu spins. For the magnetic field applied perpendicular to the tetragonal axis, the spin-induced ferroelectricity is explained in terms of the linear magnetoelectric effect that occurs in a metastable magnetic state. A relatively small electric polarization induced by the field parallel to the tetragonal axis may indicate a more complex magnetic ordering in Bi$_2$CuO$_4$.


title: Unbiased clustering estimation in the presence of missing observations

abstract: In order to be efficient, spectroscopic galaxy redshift surveys do not obtain redshifts for all galaxies in the population targeted. The missing galaxies are often clustered, commonly leading to a lower proportion of successful observations in dense regions. One example is the close-pair issue for SDSS spectroscopic galaxy surveys, which have a deficit of pairs of observed galaxies with angular separation closer than the hardware limit on placing neighbouring fibers. Spatially clustered missing observations will exist in the next generations of surveys. Various schemes have previously been suggested to mitigate these effects, but none works for all situations. We argue that the solution is to link the missing galaxies to those observed with statistically equivalent clustering properties, and that the best way to do this is to rerun the targeting algorithm, varying the angular position of the observations. Provided that every pair has a non-zero probability of being observed in one realisation of the algorithm, then a pair-upweighting scheme linking targets to successful observations, can correct these issues. We present such a scheme, and demonstrate its validity using realisations of an idealised simple survey strategy.


title: Using angular pair upweighting to improve 3D clustering measurements

abstract: Three dimensional galaxy clustering measurements provide a wealth of cosmological information. However, obtaining spectra of galaxies is expensive, and surveys often only measure redshifts for a subsample of a target galaxy population. Provided that the spectroscopic data is representative, we argue that angular pair upweighting should be used in these situations to improve the 3D clustering measurements. We present a toy model showing mathematically how such a weighting can improve measurements, and provide a practical example of its application using mocks created for the Baryon Oscillation Spectroscopic Survey (BOSS). Our analysis of mocks suggests that, if an angular clustering measurement is available over twice the area covered spectroscopically, weighting gives a $\sim$10-20% reduction of the variance of the monopole correlation function on the BAO scale.


title: Performance Analysis for Time-of-Arrival Estimation with Oversampled Low-Complexity 1-bit A/D Conversion

abstract: Analog-to-digtial (A/D) conversion plays a crucial role when it comes to the design of energy-efficient and fast signal processing systems. As its complexity grows exponentially with the number of output bits, significant savings are possible when resorting to a minimum resolution of a single bit. However, then the nonlinear effect which is introduced by the A/D converter results in a pronounced performance loss, in particular for the case when the receiver is operated outside the low signal-to-noise ratio (SNR) regime. By trading the A/D resolution for a moderately faster sampling rate, we show that for time-of-arrival (TOA) estimation under any SNR level it is possible to obtain a low-complexity $1$-bit receive system which features a smaller performance degradation then the classical low SNR hard-limiting loss of $2/\pi$ ($-1.96$ dB). Key to this result is the employment of a lower bound for the Fisher information matrix which enables us to approximate the estimation performance for coarsely quantized receivers with correlated noise models in a pessimistic way.


title: Coherent Hypernucleus Production in Antiproton-Nucleus Annihilation Reactions as a Probe for $\kappa$ meson Exchange

abstract: The hypernucleus production reaction $\bar p + {}^AZ \to {}^A_\Lambda(Z-1) + \bar{\Lambda}$ in the beam momentum range 1.5-20 GeV/c is addressed theoretically as a coherent process. The calculations are based on a covariant $t$-channel meson exchange model for the elementary $\bar p p \to \bar\Lambda \Lambda$ annihilation amplitude with parameters fixed by comparison with empirical data. Besides pseudo-scalar $K$ and vector $K^*$ mesons we also account for correlated $\pi K$ contributions, modelled by the scalar $K^*_0(800)$ or $\kappa$ meson. Initial and final state nuclear interactions are taken into account in eikonal approximation. The bound baryon wave functions are obtained self-consistently in a covariant mean-field approach. It is shown that the hypernucleus production cross sections populating discrete states are dominated by the vector and scalar interaction channels. The pronounced sensitivity on the scalar $\kappa$ meson exchange contributions indicates that these reactions are well suited as a probe for correlated $\pi K$ exchange in the scalar $\kappa/K^*_0$ interaction channel.


title: SCUBA-2 follow-up of Herschel-SPIRE observed Planck overdensities

abstract: We present SCUBA-2 follow-up of 61 candidate high-redshift Planck sources. Of these, 10 are confirmed strong gravitational lenses and comprise some of the brightest such submm sources on the observed sky, while 51 are candidate proto-cluster fields undergoing massive starburst events. With the accompanying Herschel-SPIRE observations and assuming an empirical dust temperature prior of $34^{+13}_{-9}$ K, we provide photometric redshift and far-IR luminosity estimates for 172 SCUBA-2-selected sources within these Planck overdensity fields. The redshift distribution of the sources peak between a redshift of 2 and 4, with one third of the sources having $S_{500}$/$S_{350} > 1$. For the majority of the sources, we find far-IR luminosities of approximately $10^{13}\,\mathrm{L}_\odot$, corresponding to star-formation rates of around $1000$ M$_\odot \mathrm{yr}^{-1}$. For $S_{850}>8$ mJy sources, we show that there is up to an order of magnitude increase in star-formation rate density and an increase in uncorrected number counts of $6$ for $S_{850}>8$ mJy when compared to typical cosmological survey fields. The sources detected with SCUBA-2 account for only approximately $5$ per cent of the Planck flux at 353 GHz, and thus many more fainter sources are expected in these fields.


title: Robust Motion Planning employing Signal Temporal Logic

abstract: Motion planning classically concerns the problem of accomplishing a goal configuration while avoiding obstacles. However, the need for more sophisticated motion planning methodologies, taking temporal aspects into account, has emerged. To address this issue, temporal logics have recently been used to formulate such advanced specifications. This paper will consider Signal Temporal Logic in combination with Model Predictive Control. A robustness metric, called Discrete Average Space Robustness, is introduced and used to maximize the satisfaction of specifications which results in a natural robustness against noise. The comprised optimization problem is convex and formulated as a Linear Program.


title: Lorentzian Quantum Cosmology

abstract: We argue that the Lorentzian path integral is a better starting point for quantum cosmology than the Euclidean version. In particular, we revisit the mini-superspace calculation of the Feynman path integral for quantum gravity with a positive cosmological constant. Instead of rotating to Euclidean time, we deform the contour of integration over metrics into the complex plane, exploiting Picard-Lefschetz theory to transform the path integral from a conditionally convergent integral into an absolutely convergent one. We show that this procedure unambiguously determines which semiclassical saddle point solutions are relevant to the quantum mechanical amplitude. Imposing "no-boundary" initial conditions, i.e., restricting attention to regular, complex metrics with no initial boundary, we find that the dominant saddle contributes a semiclassical exponential factor which is precisely the {\it inverse} of the famous Hartle-Hawking result.


title: A lattice formulation of the F4 completion procedure

abstract: We write a procedure for constructing noncommutative Groebner bases. Reductions are done by particular linear projectors, called reduction operators. The operators enable us to use a lattice construction to reduce simultaneously each S-polynomial into a unique normal form. We write an implementation as well as an example to illustrate our procedure. Moreover, the lattice construction is done by Gaussian elimination, which relates our procedure to the F4 algorithm for constructing commutative Groebner bases.


title: Cross-screening in observational studies that test many hypotheses

abstract: We discuss observational studies that test many causal hypotheses, either hypotheses about many outcomes or many treatments. To be credible an observational study that tests many causal hypotheses must demonstrate that its conclusions are neither artifacts of multiple testing nor of small biases from nonrandom treatment assignment. In a sense that needs to be defined carefully, hidden within a sensitivity analysis for nonrandom assignment is an enormous correction for multiple testing: in the absence of bias, it is extremely improbable that multiple testing alone would create an association insensitive to moderate biases. We propose a new strategy called "cross-screening", different from but motivated by recent work of Bogomolov and Heller on replicability. Cross-screening splits the data in half at random, uses the first half to plan a study carried out on the second half, then uses the second half to plan a study carried out on the first half, and reports the more favorable conclusions of the two studies correcting using the Bonferroni inequality for having done two studies. If the two studies happen to concur, then they achieve Bogomolov-Heller replicability; however, importantly, replicability is not required for strong control of the family-wise error rate, and either study alone suffices for firm conclusions. In randomized studies with a few hypotheses, cross-split screening is not an attractive method when compared with conventional methods of multiplicity control, but it can become attractive when hundreds or thousands of hypotheses are subjected to sensitivity analyses in an observational study. We illustrate the technique by comparing 46 biomarkers in individuals who consume large quantities of fish versus little or no fish.


title: Exact RG Invariance and Symmetry Improved 2PI Effective Potential

abstract: The Symmetry Improved Two-Particle-Irreducible (SI2PI) formalism is a powerful tool to calculate the effective potential beyond perturbation theory, whereby infinite sets of selective loop-graph topologies can be resummed in a systematic and consistent manner. In this paper we study the Renormalization-Group (RG) properties of this formalism, by proving for the first time a number of new field-theoretic results. First, the RG runnings of all proper 2PI couplings are found to be UV finite, in the Hartree-Fock and sunset approximations of the 2PI effective action. Second, the SI2PI effective potential is exactly RG invariant, in contrast to what happens in the ordinary One-Particle-Irreducible (1PI) perturbation theory, where the effective potential is RG invariant only up to higher orders. Finally, we show how the effective potential of an O(2) theory evaluated in the SI2PI framework, appropriately RG improved, can reach a higher level of accuracy, even up to one order of magnitude, with respect to the corresponding one obtained in the 1PI formalism.


title: Non-Cohen-Macaulay canonical singularities

abstract: The purpose of this note is to show that in characteristic two there exists a smooth Fano variety for which the square of the anti-canonical bundle violates Kodaira vanishing. This is used in turn to construct non-Cohen-Macaulay canonical singularities in characteristic two.


title: Estimation and prediction in sparse and unbalanced tables

abstract: We consider the problem where we have a multi-way table of means, indexed by several factors, where each factor can have a large number of levels. The entry in each cell is the mean of some response, averaged over the observations falling into that cell. Some cells may be very sparsely populated, and in extreme cases, not populated at all. We might still like to estimate an expected response in such cells. We propose here a novel hierarchical ANOVA (HANOVA) representation for such data. Sparse cells will lean more on the lower-order interaction model for the data. These in turn could have components that are poorly represented in the data, in which case they rely on yet lower-order models. Our approach leads to a simple hierarchical algorithm, requiring repeated calculations of sub-table means of modified counts. The algorithm has shown superiority over the unshrinked methods in both simulations and real data sets.


title: There are no P-points in Silver extensions

abstract: We prove that after adding a Silver real no ultrafilter from the ground model can be extended to a P-point, and this remains to be the case in any further extension which has the Sacks property. We conclude that there are no P-points in the Silver model. In particular, it is possible to construct a model without P-points by iterating Borel partial orders. This answers a question of Michael Hru\v{s}\'ak. We also show that the same argument can be used for the side-by-side product of Silver forcing. This provides a model without P-points with the continuum arbitrary large, answering a question of Wolfgang Wohofsky.


title: Auto-context Convolutional Neural Network (Auto-Net) for Brain Extraction in Magnetic Resonance Imaging

abstract: Brain extraction or whole brain segmentation is an important first step in many of the neuroimage analysis pipelines. The accuracy and robustness of brain extraction, therefore, is crucial for the accuracy of the entire brain analysis process. With the aim of designing a learning-based, geometry-independent and registration-free brain extraction tool in this study, we present a technique based on an auto-context convolutional neural network (CNN), in which intrinsic local and global image features are learned through 2D patches of different window sizes. In this architecture three parallel 2D convolutional pathways for three different directions (axial, coronal, and sagittal) implicitly learn 3D image information without the need for computationally expensive 3D convolutions. Posterior probability maps generated by the network are used iteratively as context information along with the original image patches to learn the local shape and connectedness of the brain, to extract it from non-brain tissue. The brain extraction results we have obtained from our algorithm are superior to the recently reported results in the literature on two publicly available benchmark datasets, namely LPBA40 and OASIS, in which we obtained Dice overlap coefficients of 97.42% and 95.40%, respectively. Furthermore, we evaluated the performance of our algorithm in the challenging problem of extracting arbitrarily-oriented fetal brains in reconstructed fetal brain magnetic resonance imaging (MRI) datasets. In this application our algorithm performed much better than the other methods (Dice coefficient: 95.98%), where the other methods performed poorly due to the non-standard orientation and geometry of the fetal brain in MRI. Our CNN-based method can provide accurate, geometry-independent brain extraction in challenging applications.


title: Magnetically aligned dust and SiO maser polarization in the envelope of the red supergiant VY CMa

abstract: We use Atacama Large Millimeter/submillimeter Array Band 5 science verification observations of the red supergiant VY CMa to study the polarization of SiO thermal/masers lines and dust continuum at ~1.7 mm wavelength. We analyse both linear and circular polarization and derive the magnetic field strength and structure, assuming the polarization of the lines originates from the Zeeman effect, and that of the dust originates from aligned dust grains. We also discuss other effects that could give rise to the observed polarization. We detect, for the first time, significant polarization (~3%) of the circumstellar dust emission at millimeter wavelengths. The polarization is uniform with an electric vector position angle of $\sim8^\circ$. Varying levels of linear polarization are detected for the J=4-3 28SiO v=0, 1, 2, and 29SiO v=0, 1 lines, with the strongest polarization fraction of ~30% found for the 29SiO v=1 maser. The linear polarization vectors rotate with velocity, consistent with earlier observations. We also find significant (up to ~1%) circular polarization in several lines, consistent with previous measurements. We conclude that the detection is robust against calibration and regular instrumental errors, although we cannot yet fully rule out non-standard instrumental effects. Emission from magnetically aligned grains is the most likely origin of the observed continuum polarization. This implies that the dust is embedded in a magnetic field >13 mG. The maser line polarization traces the magnetic field structure. The magnetic field in the gas and dust is consistent with an approximately toroidal field configuration, but only higher angular resolution observations will be able to reveal more detailed field structure. If the circular polarization is due to Zeeman splitting, it indicates a magnetic field strength of ~1-3 Gauss, consistent with previous maser observations.


title: Gluon Tomography from Deeply Virtual Compton Scattering at Small-x

abstract: We present a full evaluation of the deeply virtual Compton scattering (DVCS) cross section in the dipole framework in the small-x region. The result features the $\cos \phi$ and $\cos 2\phi$ azimuthal angular correlations which have been missing in previous studies based on the dipole model. In particular, the $\cos 2\phi$ term is generated by the elliptic gluon Wigner distribution whose measurement at the planned electron-ion collider (EIC) provides an important information about the gluon tomography at small-x. We also show the consistency with the standard collinear factorization approach based on the quark and gluon generalized parton distributions (GPDs).


title: The Most Massive Heartbeat: An In-depth Analysis of {\iota} Orionis

abstract: {\iota} Ori is a well studied massive binary consisting of an O9 III + B1 III/IV star. Due to its high eccentricity (e = 0.764) and short orbital period (P orb = 29.13376 d), it has been considered to be a good candidate to show evidence of tidal effects; however, none have previously been identified. Using photometry from the BRITE-Constellation space photometry mission we have confirmed the existence of tidal distortions through the presence of a heartbeat signal at periastron. We combine spectroscopic and light curve analyses to measure the masses and radii of the components, revealing {\iota} Ori to be the most massive heartbeat system known to date. In addition, using a thorough frequency analysis, we also report the unprecedented discovery of multiple tidally induced oscillations in an O star. The amplitudes of the pulsations allow us to empirically estimate the tidal circularization rate, yielding an effective tidal quality factor Q $\approx 4 \times 10^{4}$ .


title: Measurement of the W boson helicity using top pair events at $\sqrt{s}$ = 8 TeV with the CMS detector

abstract: This document gives an overview over the recent results on helicity measurement of W boson originated from top pair events. The results are obtained using data collected by the CMS detector at a center-of-mass energy of 8 TeV. The helicity measurements are confronted with the most precise theoretical predictions of the standard model.


title: The Naming Game on the complete graph

abstract: We consider a model of language development, known as the naming game, in which agents invent, share and then select descriptive words for a single object, in such a way as to promote local consensus. When formulated on a finite and connected graph, a global consensus eventually emerges in which all agents use a common unique word. Previous numerical studies of the model on the complete graph with $n$ agents suggest that when no words initially exist, the time to consensus is of order $n^{1/2}$, assuming each agent speaks at a constant rate. We show rigorously that the time to consensus is at least $n^{1/2-o(1)}$, and that it is at most constant times $\log n$ when only two words remain. In order to do so we develop sample path estimates for quasi-left continuous semimartingales with bounded jumps.


title: The variational Laplace approach to approximate Bayesian inference

abstract: Variational approaches to approximate Bayesian inference provide very efficient means of performing parameter estimation and model selection. Among these, so-called variational-Laplace or VL schemes rely on Gaussian approximations to posterior densities on model parameters. In this note, we review the main variants of VL approaches, that follow from considering nonlinear models of continuous and/or categorical data. En passant, we also derive a few novel theoretical results that complete the portfolio of existing analyses of variational Bayesian approaches, including investigations of their asymptotic convergence. We also suggest practical ways of extending existing VL approaches to hierarchical generative models that include (e.g., precision) hyperparameters.


title: The Feasibility of Dynamically Granted Permissions: Aligning Mobile Privacy with User Preferences

abstract: Current smartphone operating systems regulate application permissions by prompting users on an ask-on-first-use basis. Prior research has shown that this method is ineffective because it fails to account for context: the circumstances under which an application first requests access to data may be vastly different than the circumstances under which it subsequently requests access. We performed a longitudinal 131-person field study to analyze the contextuality behind user privacy decisions to regulate access to sensitive resources. We built a classifier to make privacy decisions on the user's behalf by detecting when context has changed and, when necessary, inferring privacy preferences based on the user's past decisions and behavior. Our goal is to automatically grant appropriate resource requests without further user intervention, deny inappropriate requests, and only prompt the user when the system is uncertain of the user's preferences. We show that our approach can accurately predict users' privacy decisions 96.8% of the time, which is a four-fold reduction in error rate compared to current systems.


title: Optimized Cost per Click in Taobao Display Advertising

abstract: Taobao, as the largest online retail platform in the world, provides billions of online display advertising impressions for millions of advertisers every day. For commercial purposes, the advertisers bid for specific spots and target crowds to compete for business traffic. The platform chooses the most suitable ads to display in tens of milliseconds. Common pricing methods include cost per mille (CPM) and cost per click (CPC). Traditional advertising systems target certain traits of users and ad placements with fixed bids, essentially regarded as coarse-grained matching of bid and traffic quality. However, the fixed bids set by the advertisers competing for different quality requests cannot fully optimize the advertisers' key requirements. Moreover, the platform has to be responsible for the business revenue and user experience. Thus, we proposed a bid optimizing strategy called optimized cost per click (OCPC) which automatically adjusts the bid to achieve finer matching of bid and traffic quality of page view (PV) request granularity. Our approach optimizes advertisers' demands, platform business revenue and user experience and as a whole improves traffic allocation efficiency. We have validated our approach in Taobao display advertising system in production. The online A/B test shows our algorithm yields substantially better results than previous fixed bid manner.


title: The Description and Scaling Behavior for the Inner Region of the Boundary Layer for 2-D Wall-bounded Flows

abstract: A second derivative-based moment method is proposed for describing the thickness and shape of the region where viscous forces are dominant in turbulent boundary layer flows. Rather than the fixed location sublayer model presently employed, the new method defines thickness and shape parameters that are experimentally accessible without differentiation. It is shown theoretically that one of the new length parameters used as a scaling parameter is also a similarity parameter for the velocity profile. In fact, we show that this new length scale parameter removes one of the theoretical inconsistencies present in the traditional Prandtl Plus scaling's. Furthermore, the new length parameter and the Prandtl Plus scaling parameters perform identically when operating on experimental datasets if the Rotta similarity constraint (u_tau/u_e = constant) holds. This means that many of the past successes ascribed to the Prandtl Plus scaling also apply to the new parameter set but without one of the theoretical inconsistencies. Examples are offered to show how the new description method is useful in exploring the actual physics of the boundary layer.


title: Approach to spatial estimation of meteorological droughts in the Coello River basin, Colombia

abstract: The space-temporal evaluation to characterize meteorological droughts was based on data accumulated monthly precipitation between 1996-2005 from 20 meteorological stations distributed in the Coello River basin. Data precipitation was performed preprocessing with data consistency tests to correct and delete data over- or under estimated. To estimate missing precipitation data are compared three geostatistical interpolation methods derived from Kriging, associated with secondary variables such as the Ordinary Kriging, CoKrigin Ordinary associated with secondary variables of a Digital Elevation Model and data satellite TRMM. To select the statistical method the setting of each interpolation was compared with respect to three reference stations through three quality tests, which were Root Mean Square Error (RMSE), Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC). In this investigation two of three tests favor the Ordinary CoKriging using as a secondary variable Altitude (COK+DEM). With the interpolated series of precipitation were evaluated and characterized by drought Standardized Precipitation Index (SPI) at monthly and quarterly scale, calculating the parameters of severity, duration, intensity and frequency of droughts. By mapping are delimited the regions where occur the more negative values of SPI. In analyzing spacetemporal the months of January, February, July and August are the driest of the year. In 1997 the meteorological drought greatest damage occurs in the Coello River basin generally concentrated in the middle and lower part of the basin, with a maximum intensity of SPI -2,57.


title: Magnetic and Hyperfine interactions in HoFe1-xCrxO3 compounds

abstract: We report on the magnetic and Mossbauer properties of polycrystalline HoFe1-xCrxO3 compounds. Magnetization data reveals the continuous tailoring of magnetic transition due to weakening of Ho3+ Fe3+ and Fe3+ Fe3+ interactions in the entire temperature range by replacing the Fe3+ ions with Cr3+ ions. The observed decrease in Neel temperature (TN) and increase in spin re-orientation transition temperature (TSR) with the replacement of Fe3+ with Cr3+ is ascribed to the weakening of Fe(Cr) O Fe(Cr) antiferromagnetic exchange interaction. In addition, we also attribute such a change in TN to the enhancement of ferromagnetic interaction of adjacent Cr3+ moments through t e hybridization as a result of the structural distortion. The decrease in isomer shift (IS) suggests enhancement of the interaction between nuclear charge with the 3s electrons as a result of decrease in radial part of 3d wave function with Cr addition. In this paper we also discuss about the variation of quadrupole splitting (QS) and hyperfine fields (Hhf) with Cr addition in HoFe1-xCrxO3 compounds.


title: In-situ Polarized Neutron Reflectometry: Epitaxial Thin Film Growth of Fe on Cu(001) by DC Magnetron Sputtering

abstract: The step-wise growth of epitaxial Fe on Cu(001)/Si(001), investigated by in-situ polarized neutron reflectometry is presented. A sputter deposition system was integrated into the neutron reflectometer AMOR at the Swiss neutron spallation source SINQ, which enables the analysis of the microstructure and magnetic moments during all deposition steps of the Fe layer. We report on the progressive evolution of the accessible parameters describing the microstructure and the magnetic properties of the Fe film, which reproduce known features and extend our knowledge on the behavior of ultrathin iron films.


title: Structural and emission properties of Tb3+-doped nitrogen-rich silicon oxynitride films

abstract: Terbium doped silicon oxynitride host matrix is suitable for various applications such as light emitters compatible with CMOS technology or frequency converter systems for photovoltaic cells. In this study, amorphous Tb3+ ion doped nitrogen-rich silicon oxynitride (NRSON) thin films were fabricated using a reactive magnetron co-sputtering method, with various N2 flows and annealing conditions, in order to study their structural and emission properties. Rutherford backscattering (RBS) measurements and refractive index values confirmed the silicon oxynitride nature of the films. An electron microscopy analysis conducted for different annealing temperatures (T A) was also performed up to 1200 {\textdegree}C. Transmission electron microscopy (TEM) images revealed two different sublayers. The top layer showed porosities coming from a degassing of oxygen during deposition and annealing, while in the region close to the substrate, a multilayer-like structure of SiO2 and Si3N4 phases appeared, involving a spinodal decomposition. Upon a 1200 {\textdegree}C annealing treatment, a significant density of Tb clusters was detected, indicating a higher thermal threshold of rare earth (RE) clusterization in comparison to the silicon oxide matrix. With an opposite variation of the N2 flow during the deposition, the nitrogen excess parameter (Nex) estimated by RBS measurements was introduced to investigate the Fourier transform infrared (FTIR) spectrum behavior and emission properties. Different vibration modes of the Si--N and Si--O bonds have been carefully identified from the FTIR spectra characterizing such host matrices, especially the 'out-of-phase' stretching vibration mode of the Si--O bond. The highest Tb3+ photoluminescence (PL) intensity was obtained by optimizing the N incorporation and the annealing conditions. In addition, according to these conditions, the integrated PL intensity variation confirmed that the silicon nitride-based host matrix had a higher thermal threshold of rare earth clusterization than its silicon oxide counterpart. Analysis of time-resolved PL intensity versus T A showed the impact of Tb clustering on decay times, in agreement with the TEM observations. Finally, PL and PL excitation (PLE) experiments and comparison of the related spectra between undoped and Tb-doped samples were carried out to investigate the impact of the band tails on the excitation mechanism of Tb3+ ions.


title: Absence of magnetic long range order in Y$_{2}$CrSbO$_{7}$: bond-disorder induced magnetic frustration in a ferromagnetic pyrochlore

abstract: The consequences of nonmagnetic-ion dilution for the pyrochlore family Y$_{2}$($M_{1-x}N_{x}$)$_{2}$O$_{7}$ ($M$ = magnetic ion, $N$ = nonmagnetic ion) have been investigated. As a first step, we experimentally examine the magnetic properties of Y$_{2}$CrSbO$_{7}$ ($x$ = 0.5), in which the magnetic sites (Cr$^{3+}$) are percolative. Although the effective Cr-Cr spin exchange is ferromagnetic, as evidenced by a positive Curie-Weiss temperature, $\Theta_\mathrm{{CW}}$ = 20.1(6) K, our high-resolution neutron powder diffraction measurements detect no sign of magnetic long range order down to 2 K. In order to understand our observations, we performed numerical simulations to study the bond-disorder introduced by the ionic size mismatch between $M$ and $N$. Based on these simulations, bond-disorder ($x_{b}$ $\simeq$ 0.23) percolates well ahead of site-disorder ($x_{s}$ $\simeq$ 0.61). This model successfully reproduces the critical region (0.2 < $x$ < 0.25) for the N\'eel to spin glass phase transition in Zn(Cr$_{1-x}$Ga$_{x}$)$_{2}$O$_{4}$, where the Cr/Ga-sublattice forms the same corner-sharing tetrahedral network as the $M/N$-sublattice in Y$_{2}$($M_{1-x}N_{x}$)$_{2}$O$_{7}$, and the rapid drop in magnetically ordered moment in the N\'eel phase [Lee $et$ $al$, Phys. Rev. B 77, 014405 (2008)]. Our study stresses the nonnegligible role of bond-disorder on magnetic frustration, even in ferromagnets.


title: The Impact of Road Configuration on V2V-based Cooperative Localization

abstract: Cooperative localization with map matching has been shown to reduce Global Navigation Satellite System (GNSS) localization error from several meters to sub-meter level by fusing the GNSS measurements of four vehicles in our previous work. While further error reduction is expected to be achievable by increasing the number of vehicles, the quantitative relationship between the estimation error and the number of connected vehicles has neither been systematically investigated nor analytically proved. In this work, a theoretical study is presented that analytically proves the correlation between the localization error and the number of connected vehicles in two cases of practical interest. More specifically, it is shown that, under the assumption of small non-common error, the expected square error of the GNSS common error correction is inversely proportional to the number of vehicles, if the road directions obey a uniform distribution, or inversely proportional to logarithm of the number of vehicles, if the road directions obey a Bernoulli distribution. Numerical simulations are conducted to justify these analytic results. Moreover, the simulation results show that the aforementioned error decrement rates hold even when the assumption of small non-common error is violated.


title: Balanced flux formulations for multidimensional Evans function computations for viscous shocks

abstract: The Evans function is a powerful tool for the stability analysis of viscous shock profiles; zeros of this function carry stability information. In the one-dimensional case, it is typical to compute the Evans function using Goodman's integrated coordinates [G1]; this device facilitates the search for zeros of the Evans function by winding number arguments. Although integrated coordinates are not available in the multidimensional case, we show here that there is a choice of coordinates which gives similar advantages.


title: Guarantees for Greedy Maximization of Non-submodular Functions with Applications

abstract: We investigate the performance of the standard Greedy algorithm for cardinality constrained maximization of non-submodular nondecreasing set functions. While there are strong theoretical guarantees on the performance of Greedy for maximizing submodular functions, there are few guarantees for non-submodular ones. However, Greedy enjoys strong empirical performance for many important non-submodular functions, e.g., the Bayesian A-optimality objective in experimental design. We prove theoretical guarantees supporting the empirical performance. Our guarantees are characterized by a combination of the (generalized) curvature $\alpha$ and the submodularity ratio $\gamma$. In particular, we prove that Greedy enjoys a tight approximation guarantee of $\frac{1}{\alpha}(1- e^{-\gamma\alpha})$ for cardinality constrained maximization. In addition, we bound the submodularity ratio and curvature for several important real-world objectives, including the Bayesian A-optimality objective, the determinantal function of a square submatrix and certain linear programs with combinatorial constraints. We experimentally validate our theoretical findings for both synthetic and real-world applications.


title: Optical Nanoscopy of High-Tc Cuprate Nano-Constriction Devices Patterned by Helium Ion Beams

abstract: Helium-ion beams (HIB) focused to sub-nanometer scales have emerged as powerful tools for high-resolution imaging as well as nano-scale lithography, ion milling or deposition. Quantifying irradiation effects is essential for reliable device fabrication but most of the depth profiling information is provided by computer simulations rather than experiment. Here, we use atomic force microscopy (AFM) combined with scanning near-field optical microscopy (SNOM) to provide three-dimensional (3D) dielectric characterization of high-temperature superconductor devices fabricated by HIB. By imaging the infrared dielectric response we find that amorphization caused by the nominally 0.5 nm HIB extends throughout the entire 26.5 nm thickness of the cuprate film and by about 500 nm laterally. This unexpectedly widespread structural and electronic damage can be attributed to a Helium depth distribution substantially modified by internal device interfaces. Our study introduces AFM-SNOM as a quantitative nano-scale tomographic technique for non-invasive 3D characterization of irradiation damage in a wide variety of devices.


title: Revisiting stochastic off-policy action-value gradients

abstract: Off-policy stochastic actor-critic methods rely on approximating the stochastic policy gradient in order to derive an optimal policy. One may also derive the optimal policy by approximating the action-value gradient. The use of action-value gradients is desirable as policy improvement occurs along the direction of steepest ascent. This has been studied extensively within the context of natural gradient actor-critic algorithms and more recently within the context of deterministic policy gradients. In this paper we briefly discuss the off-policy stochastic counterpart to deterministic action-value gradients, as well as an incremental approach for following the policy gradient in lieu of the natural gradient.


title: Leveling the playing field for Visually Impaired using Transport Assistant

abstract: Visually impaired people face numerous challenges when it comes to transportation. Not only must they circumvent obstacles while navigating, but they also need access to essential information related to available public transport, up-to-date weather forecast, and convenient method for booking private taxis. In this paper we introduce Transport Assistant - a voice based assistive technology prototype, built with a goal of leveling the playing field for the visually impaired to solve these problems that they face in their day to day life. Being voice enabled makes it seamlessly integrate into the environment, and can be invoked by saying a hotword - hello assistant. The paper explores this research question, followed by investigating existing technologies, explains the methodology and design, then concludes by presenting the prototype and results.


title: Long-run dynamics of the U.S. patent classification system

abstract: Almost by definition, radical innovations create a need to revise existing classification systems. In this paper, we argue that classification system changes and patent reclassification are common and reveal interesting information about technological evolution. To support our argument, we present three sets of findings regarding classification volatility in the U.S. patent classification system. First, we study the evolution of the number of distinct classes. Reconstructed time series based on the current classification scheme are very different from historical data. This suggests that using the current classification to analyze the past produces a distorted view of the evolution of the system. Second, we study the relative sizes of classes. The size distribution is exponential so classes are of quite different sizes, but the largest classes are not necessarily the oldest. To explain this pattern with a simple stochastic growth model, we introduce the assumption that classes have a regular chance to be split. Third, we study reclassification. The share of patents that are in a different class now than they were at birth can be quite high. Reclassification mostly occurs across classes belonging to the same 1-digit NBER category, but not always. We also document that reclassified patents tend to be more cited than non-reclassified ones, even after controlling for grant year and class of origin.


title: Network Structure and Naive Sequential Learning

abstract: We study a sequential-learning model featuring a network of naive agents with Gaussian information structures. Agents apply a heuristic rule to aggregate predecessors' actions. They weigh these actions according the strengths of their social connections to different predecessors. We show this rule arises endogenously when agents wrongly believe others act solely on private information and thus neglect redundancies among observations. We provide a simple linear formula expressing agents' actions in terms of network paths and use this formula to characterize the set of networks where naive agents eventually learn correctly. This characterization implies that, on all networks where later agents observe more than one neighbor, there exist disproportionately influential early agents who can cause herding on incorrect actions. Going beyond existing social-learning results, we compute the probability of such mislearning exactly. This allows us to compare likelihoods of incorrect herding, and hence expected welfare losses, across network structures. The probability of mislearning increases when link densities are higher and when networks are more integrated. In partially segregated networks, divergent early signals can lead to persistent disagreement between groups.


title: Sub-Microsecond X-Ray Imaging Using Hole-Collecting Schottky type CdTe with Charge-Integrating Pixel Array Detectors

abstract: CdTe is increasingly being used as the x-ray sensing material in imaging pixel array detectors for x-rays, generally above 20 keV, where silicon sensors become unacceptably transparent. Unfortunately CdTe suffers from polarization, which can alter the response of the material over time and with accumulated dose. Most prior studies used long integration times or CdTe that was not of the hole-collecting Schottky type. We investigated the temporal response of hole-collecting Schottky type CdTe sensors on timescales ranging from tens of nanoseconds to several seconds. We found that the material shows signal persistence on the timescale of hundreds of milliseconds attributed to the detrapping of a shallow trap, and additional persistence on sub-microsecond timescales after polarization. The results show that this type of CdTe can be used for time resolved studies down to approximately 100 ns. However quantitative interpretation of the signal requires careful attention to bias voltages, polarization and exposure history.


title: Encoding qubits into oscillators with atomic ensembles and squeezed light

abstract: The Gottesman-Kitaev-Preskill (GKP) encoding of a qubit within an oscillator provides a number of advantages when used in a fault-tolerant architecture for quantum computing, most notably that Gaussian operations suffice to implement all single- and two-qubit Clifford gates. The main drawback of the encoding is that the logical states themselves are challenging to produce. Here we present a method for generating optical GKP-encoded qubits by coupling an atomic ensemble to a squeezed state of light. Particular outcomes of a subsequent spin measurement of the ensemble herald successful generation of the resource state in the optical mode. We analyze the method in terms of the resources required (total spin and amount of squeezing) and the probability of success. We propose a physical implementation using a Faraday-based quantum non-demolition interaction.


title: Robust Estimation of Scattering in Pulsar Timing Analysis

abstract: We present a robust approach to incorporating models for the time-variable broadening of the pulse profile due to scattering in the ionized interstellar medium into profile-domain pulsar timing analysis. We use this approach to simultaneously estimate temporal variations in both the dispersion measure (DM) and scattering, together with a model for the pulse profile that includes smooth evolution as a function of frequency, and the pulsar's timing model. We show that fixing the scattering timescales when forming time-of-arrival estimates, as has been suggested in the context of traditional pulsar timing analysis, can significantly underestimate the uncertainties in both DM, and the arrival time of the pulse, leading to bias in the timing parameters. We apply our method using a new, publicly available, GPU accelerated code, both to simulations, and observations of the millisecond pulsar PSR J1643$-$1224. This pulsar is known to exhibit significant scattering variability compared to typical millisecond pulsars, and we find including low-frequency ($< 1$ GHz) data without a model for these scattering variations leads to significant periodic structure in the DM, and also biases the astrometric parameters at the $4\sigma$ level, for example, changing proper motion in right ascension by $0.50 \pm 0.12$. If low frequency observations are to be included when significant scattering variations are present, we conclude it is necessary to not just model those variations, but also to sample the parameters that describe the variations simultaneously with all other parameters in the model, a task for which profile domain pulsar timing is ideally suited.


title: The energy-momentum tensor of electromagnetic fields in matter

abstract: We present a complete resolution of the Abraham-Minkowski controversy . This is done by considering several new aspects which invalidate previous discussions. We show that: 1)For polarized matter the center of mass theorem is no longer valid in its usual form. A contribution related to microscopic spin should be considered. 2)The electromagnetic dipolar energy density contributes to the inertia of matter and should be incorporated covariantly to the the energy-momentum tensor of matter. Then there is an electromagnetic component in matter's momentum density whose variation explains the results of the only experiment which supports Abraham's force. 3)Averaging the microscopic Lorentz's force results in the unambiguos expression for the force density exerted by the field. This force density is consistent with all the experimental evidence. 4)Momentum conservation determines the electromagnetic energy-momentum tensor. This tensor is different from Abraham's and Minkowski's tensors, but one recovers Minkowski's expression for the momentum density. The energy density is different from Poynting's expression but Poynting's vector remains the same. Our tensor is non-symmetric which allows the field to exert a distributed torque on matter. We use our results to discuss momentum and angular momentum exchange in various situations of physical interest. We find complete consistency of our equations in the description of the systems considered. We also show that several alternative expressions of the field energy-momentum tensor and force-density cannot be successfully used in all our examples. In particular we verify in two of these examples that the center of mass and spin introduced by us moves with constant velocity, but that the standard center of mass does not.


title: Resolved, expanding jets in the Galactic black hole candidate XTE J1908+094

abstract: Black hole X-ray binaries undergo occasional outbursts caused by changing inner accretion flows. Here we report high-angular resolution radio observations of the 2013 outburst of the black hole candidate X-ray binary system J1908+094, using data from the VLBA and EVN. We show that following a hard-to-soft state transition, we detect moving jet knots that appear asymmetric in morphology and brightness, and expand to become laterally resolved as they move away from the core, along an axis aligned approximately $-11$\degree\ east of north. We initially see only the southern component, whose evolution gives rise to a 15-mJy radio flare and generates the observed radio polarization. This fades and becomes resolved out after 4 days, after which a second component appears to the north, moving in the opposite direction. From the timing of the appearance of the knots relative to the X-ray state transition, a 90\degree\ swing of the inferred magnetic field orientation, the asymmetric appearance of the knots, their complex and evolving morphology, and their low speeds, we interpret the knots as working surfaces where the jets impact the surrounding medium. This would imply a substantially denser environment surrounding XTE J1908+094 than has been inferred to exist around the microquasar sources GRS 1915+105 and GRO J1655-40.


title: Classification and clustering for observations of event time data using non-homogeneous Poisson process models

abstract: Data of the form of event times arise in various applications. A simple model for such data is a non-homogeneous Poisson process (NHPP) which is specified by a rate function that depends on time. We consider the problem of having access to multiple independent observations of event time data, observed on a common interval, from which we wish to classify or cluster the observations according to their rate functions. Each rate function is unknown but assumed to belong to a finite number of rate functions each defining a distinct class. We model the rate functions using a spline basis expansion, the coefficients of which need to be estimated from data. The classification approach consists of using training data for which the class membership is known, to calculate maximum likelihood estimates of the coefficients for each group, then assigning test observations to a group by a maximum likelihood criterion. For clustering, by analogy to the Gaussian mixture model approach for Euclidean data, we consider mixtures of NHPP and use the expectation-maximisation algorithm to estimate the coefficients of the rate functions for the component models and group membership probabilities for each observation. The classification and clustering approaches perform well on both synthetic and real-world data sets. Code associated with this paper is available at https://github.com/duncan-barrack/NHPP .


title: Process convolution approaches for modeling interacting trajectories

abstract: Gaussian processes are a fundamental statistical tool used in a wide range of applications. In the spatio-temporal setting, several families of covariance functions exist to accommodate a wide variety of dependence structures arising in different applications. These parametric families can be restrictive and are insufficient in some situations. In contrast, process convolutions represent a flexible, interpretable approach to defining the covariance of a Gaussian process and have modest requirements to ensure validity. We introduce a generalization of the process convolution approach that employs multiple convolutions sequentially to form a "process convolution chain." In our proposed multi-stage framework, complex dependencies that arise from a combination of different interacting mechanisms are decomposed into a series of interpretable kernel smoothers. We demonstrate an application of process convolution chains to model killer whale movement, in which the paths taken by multiple individuals are not independent, but reflect dynamic social interactions within the population. Our proposed model for dependent movement provides inference for the latent dynamic social structure in the study population. Additionally, by leveraging the positive dependence among individual paths, we achieve a reduction in uncertainty for the estimated locations of the whales, compared to a model that treats paths as independent.


title: A fresh look at effect aliasing and interactions: some new wine in old bottles

abstract: Interactions and effect aliasing are among the fundamental concepts in experimental design. In this paper, some new insights and approaches are provided on these subjects. In the literature, the "de-aliasing" of aliased effects is deemed to be impossible. We argue that this "impossibility" can indeed be resolved by employing a new approach which consists of reparametrization of effects and exploitation of effect non-orthogonality. This approach is successfully applied to three classes of designs: regular and nonregular two-level fractional factorial designs, and three-level fractional factorial designs. For reparametrization, the notion of conditional main effects (cme's) is employed for two-level regular designs, while the linear-quadratic system is used for three-level designs. For nonregular two-level designs, reparametrization is not needed because the partial aliasing of their effects already induces non-orthogonality. The approach can be extended to general observational data by using a new bi-level variable selection technique based on the cme's. A historical recollection is given on how these ideas were discovered.


title: The Ohm-Rush content function II. Noetherian rings, valuation domains, and base change

abstract: The notion of an Ohm-Rush algebra, and its associated content map, has connections with prime characteristic algebra, polynomial extensions, and the Ananyan-Hochster proof of Stillman's conjecture. As further restrictions are placed (creating the increasingly more specialized notions of weak content, semicontent, content, and Gaussian algebras), the construction becomes more powerful. Here we settle the question in the affirmative over a Noetherian ring from our previous article of whether a faithfully flat weak content algebra is semicontent (and over an Artinian ring of whether such an algebra is content), though both questions remain open in general. We show that in content algebra maps over Pr\"ufer domains, heights are preserved and a dimension formula is satisfied. We show that an inclusion of nontrivial valuation domains is a content algebra if and only if the induced map on value groups is an isomorphism, and that such a map induces a homeomorphism on prime spectra. Examples are given throughout, including results that show the subtle role played by properties of transcendental field extensions.


title: Computational Search for Novel Hard Chromium-Based Materials

abstract: Nitrides, carbides and borides of transition metals are an attractive class of hard materials. Our recent preliminary explorations of the binary chemical compounds indicated that chromium-based materials are among the hardest transition metal compounds. Motivated by this, here we explore in detail the binary Cr-B, Cr-C and Cr-N systems using global optimization techniques. Calculated enthalpy of formation and hardness of predicted materials were used for Pareto optimization to define the hardest materials with lowest energy. Our calculations recover all numerous known stable compounds (except $Cr_{23}C_6$ with its large unit cell) and discover a novel stable phase $Pmn2_1$-$Cr_2C$. We resolve the structure of $Cr_2N$ and find it to be of anti-$CaCl_2$ type (space group Pnnm). Many of these phases possess remarkable hardness, but only $CrB_4$ is superhard (Vickers hardness 48 GPa). Among chromium compounds, borides generally possess highest hardnesses and greatest stability. Under pressure, we predict stabilization of a TMDC-like phase of $Cr_2N$, a WC-type phase of CrN, and a new compound $CrN_4$. Nitrogen-rich chromium nitride $CrN_4$ is a high energy-density material featuring polymeric nitrogen chains. In the presence of metal atoms polymerization of nitrogen takes place at much lower pressures: $CrN_4$ becomes stable at 15 GPa.


title: Evaluation of Machine Learning Methods to Predict Coronary Artery Disease Using Metabolomic Data

abstract: Metabolomic data can potentially enable accurate, non-invasive and low-cost prediction of coronary artery disease. Regression-based analytical approaches however might fail to fully account for interactions between metabolites, rely on a priori selected input features and thus might suffer from poorer accuracy. Supervised machine learning methods can potentially be used in order to fully exploit the dimensionality and richness of the data. In this paper, we systematically implement and evaluate a set of supervised learning methods (L1 regression, random forest classifier) and compare them to traditional regression-based approaches for disease prediction using metabolomic data.


title: Theory Morphisms in Church's Type Theory with Quotation and Evaluation

abstract: ${\rm CTT}_{\rm qe}$ is a version of Church's type theory with global quotation and evaluation operators that is engineered to reason about the interplay of syntax and semantics and to formalize syntax-based mathematical algorithms. ${\rm CTT}_{\rm uqe}$ is a variant of ${\rm CTT}_{\rm qe}$ that admits undefined expressions, partial functions, and multiple base types of individuals. It is better suited than ${\rm CTT}_{\rm qe}$ as a logic for building networks of theories connected by theory morphisms. This paper presents the syntax and semantics of ${\rm CTT}_{\rm uqe}$, defines a notion of a theory morphism from one ${\rm CTT}_{\rm uqe}$ theory to another, and gives two simple examples that illustrate the use of theory morphisms in ${\rm CTT}_{\rm uqe}$.


title: Computing in Memory with Spin-Transfer Torque Magnetic RAM

abstract: In-memory computing is a promising approach to addressing the processor-memory data transfer bottleneck in computing systems. We propose Spin-Transfer Torque Compute-in-Memory (STT-CiM), a design for in-memory computing with Spin-Transfer Torque Magnetic RAM (STT-MRAM). The unique properties of spintronic memory allow multiple wordlines within an array to be simultaneously enabled, opening up the possibility of directly sensing functions of the values stored in multiple rows using a single access. We propose modifications to STT-MRAM peripheral circuits that leverage this principle to perform logic, arithmetic, and complex vector operations. We address the challenge of reliable in-memory computing under process variations by extending ECC schemes to detect and correct errors that occur during CiM operations. We also address the question of how STT-CiM should be integrated within a general-purpose computing system. To this end, we propose architectural enhancements to processor instruction sets and on-chip buses that enable STT-CiM to be utilized as a scratchpad memory. Finally, we present data mapping techniques to increase the effectiveness of STT-CiM. We evaluate STT-CiM using a device-to-architecture modeling framework, and integrate cycle-accurate models of STT-CiM with a commercial processor and on-chip bus (Nios II and Avalon from Intel). Our system-level evaluation shows that STT-CiM provides system-level performance improvements of 3.93x on average (upto 10.4x), and concurrently reduces memory system energy by 3.83x on average (upto 12.4x).


title: Almost euclidean Isoperimetric Inequalities in spaces satisfying local Ricci curvature lower bounds

abstract: Motivated by Perelman's Pseudo Locality Theorem for the Ricci flow, we prove that if a Riemannian manifold has Ricci curvature bounded below in a metric ball which moreover has almost maximal volume, then in a smaller ball (in a quantified sense) it holds an almost-euclidean isoperimetric inequality. The result is actually established in the more general framework of non-smooth spaces satisfying local Ricci curvature lower bounds in a synthetic sense via optimal transportation.


title: Root separation for reducible monic polynomials of odd degree

abstract: We study root separation of reducible monic integer polynomials of odd degree. Let h(P) be the naive height, sep(P) the minimal distance between two distinct roots of an integer polynomial P(x) and sep(P)=h(P)^{-e(P)}. Let e_r*(d)=limsup_{deg(P)=d, h(P)-> +infty} e(P), where the limsup is taken over the reducible monic integer polynomials P(x) of degree d. We prove that e_r*(d) <= d-2. We also obtain a lower bound for e_r*(d) for d odd, which improves previously known lower bounds for e_r*(d) when d = 5, 7, 9.


title: Modeling quasi-ballistic transient thermal transport with spatially sinusoidal heating: a McKelvey-Shockley flux approach

abstract: Ballistic phonon effects, arising on length scales comparable to the mean-free-path, result in non-diffusive heat flow and alter the thermal properties of materials. Simple theoretical models that accurately capture non-diffusive transport physics are valuable for experimental analysis, technology design, and providing physical insight. In this work, we utilize and extend the McKelvey-Shockley (McK-S) flux method, a simple and accurate framework, to investigate ballistic effects in transient phonon transport submitted to a spatially sinusoidal heating profile, simulating a transient thermal grating. We begin by extending a previous McK-S formulation to include inelastic scattering, then obtain an analytical solution in the single phonon energy case (gray approximation), and after show how this approach can readily support a full phonon dispersion and mean-free-path distribution. The results agree with experimental data and compare very well to solutions of the phonon Boltzmann transport equation in the diffusive and weakly quasi-ballistic transport regimes. We discuss the role of ballistic and non-equilibrium physics, and show that inelastic scattering is key to retrieving the heat equation solution in the diffusive limit. Overall the McK-S flux method, which takes the form of a diffusion-like equation, proves to be a simple and accurate framework that is applicable from the ballistic to diffusive transport regime.


title: Supersymmetric D-term Twin Higgs

abstract: We propose a new type of supersymmetric Twin Higgs model where the SU(4) invariant quartic term is provided by a D-term potential of a new U(1) gauge symmetry. In the model the 125 GeV Higgs mass can be obtained for stop masses below 1 TeV, and a tuning required to obtain the correct electroweak scale can be as low as 20 %. A stop mass of about 2 TeV is also possible with tuning of order O(10) %.


title: Combining symmetry collective states with coupled cluster theory: Lessons from the Agassi model Hamiltonian

abstract: The failures of single-reference coupled cluster for strongly correlated many-body systems is flagged at the mean-field level by the spontaneous breaking of one or more physical symmetries of the Hamiltonian. Restoring the symmetry of the mean-field determinant by projection reveals that coupled cluster fails because it factorizes high-order excitation amplitudes incorrectly. However, symmetry-projected mean-field wave functions do not account sufficiently for dynamic (or weak) correlation. Here we pursue a merger of symmetry projection and coupled cluster theory, following previous work along these lines that utilized the simple Lipkin model system as a testbed [J. Chem. Phys. 146, 054110 (2017)]. We generalize the concept of a symmetry-projected mean-field wave function to the concept of a symmetry projected state, in which the factorization of high-order excitation amplitudes in terms of low-order ones is guided by symmetry projection and is not exponential, and combine them with coupled cluster theory in order to model the ground state of the Agassi Hamiltonian. This model has two separate channels of correlation and two separate physical symmetries which are broken under strong correlation. We show how the combination of symmetry collective states and coupled cluster is effective in obtaining correlation energies and order parameters of the Agassi model throughout its phase diagram.


title: Non-line-of-sight tracking of people at long range

abstract: A remote-sensing system that can determine the position of hidden objects has applications in many critical real-life scenarios, such as search and rescue missions and safe autonomous driving. Previous work has shown the ability to range and image objects hidden from the direct line of sight, employing advanced optical imaging technologies aimed at small objects at short range. In this work we demonstrate a long-range tracking system based on single laser illumination and single-pixel single-photon detection. This enables us to track one or more people hidden from view at a stand-off distance of over 50~m. These results pave the way towards next generation LiDAR systems that will reconstruct not only the direct-view scene but also the main elements hidden behind walls or corners.


title: Dynamic Wavelength-Tunable Photodetector Using Subwavelength Graphene Field-Effect Transistors

abstract: Dynamic wavelength tunability has long been the holy grail of photodetector technology. Because of its atomic thickness and unique properties, graphene opens up new paradigms to realize this concept, but so far this has been elusive experimentally. Here we employ detailed quantum transport modeling of photocurrent in graphene field-effect transistors (including realistic electromagnetic fields) to show that wavelength tunability is possible by dynamically changing the gate voltage. We reveal the phenomena that govern the behavior of this type of device and show significant departure from the simple expectations based on vertical transitions. We find strong focusing of the electromagnetic fields at the contact edges over the same length scale as the band-bending. Both of these spatially-varying potentials lead to an enhancement of non-vertical optical transitions, which dominate even the absence of phonon or impurity scattering. We also show that the vanishing density of states near the Dirac point leads to contact blocking and a gate-dependent modulation of the photocurrent.


title: Magnetic disorder in superconductors: Enhancement by mesoscopic fluctuations

abstract: We study the density of states (DOS) and the transition temperature $T_c$ in a dirty superconducting film with rare classical magnetic impurities of an arbitrary strength described by the Poissonian statistics. We take into account that the potential disorder is a source for mesoscopic fluctuations of the local DOS, and, consequently, for the effective strength of magnetic impurities. We find that these mesoscopic fluctuations result in a non-zero DOS for all energies in the region of the phase diagram where without this effect the DOS is zero within the standard mean-field theory. This mechanism can be more efficient in filling the mean-field superconducting gap than rare fluctuations of the potential disorder (instantons). Depending on the magnetic impurity strength, the suppression of $T_c$ by spin-flip scattering can be faster or slower than in the standard mean-field theory.


title: On the arithmetic of a family of degree-two K3 surfaces

abstract: Let $\mathbb{P}$ denote the weighted projective space with weights $(1,1,1,3)$ over the rationals, with coordinates $x,y,z,$ and $w$; let $\mathcal{X}$ be the generic element of the family of surfaces in $\mathbb{P}$ given by \begin{equation*} X\colon w^2=x^6+y^6+z^6+tx^2y^2z^2. \end{equation*} The surface $\mathcal{X}$ is a K3 surface over the function field $\mathbb{Q}(t)$. In this paper, we explicitly compute the geometric Picard lattice of $\mathcal{X}$, together with its Galois module structure, as well as derive more results on the arithmetic of $\mathcal{X}$ and other elements of the family $X$.


title: Electric dipoles vs. magnetic dipoles - for two molecules in a harmonic trap

abstract: We study energy levels of two heteronuclear molecules moving in a spherically symmetric harmonic trap. A role of electric dipole interactions is compared and contrasted with our earlier results (https://arxiv.org/abs/1512.00631) for two magnetic dipolar atoms. We stress importance of a rotational energy with its value which is very high compared to the energy of dipolar interaction. We show that dipolar forces do not play a significant role in the ground state of the system under typical experimental conditions. However, there exist excited states that exhibit anticrossings similar to the ones observed for magnetic dipoles.


title: Matter-wave physics with nanoparticles and biomolecules

abstract: These lecture notes emerged from a contribution to the "Les Houches Summer School, Session CVII - Current Trends in Atomic Physics, July 2016". It is meant to serve students as a guide to a selection of topics that are currently at the focus of molecular quantum optics and complements an earlier lecture on related topics [Arndt et al., 2014]. In this review we discuss recent advances in molecular quantum optics with large organic molecules. In particular, we present successful experiments of molecules of biological importance, such as neurotransmitters and map out the route towards interferometry of large bio-matter such as proteins. The interaction of internally complex molecules with various beam splitters is discussed for both far-field diffraction at a single nanomechanical grating and the Kapitza-Dirac Talbot-Lau near-field interferometer - addressing recent progress, experimental challenges, and prospects for quantum enhanced metrology of highly complex systems in the gas phase. A central part of this review deals with the preparation of slow beams of neutral biomolecules, ranging from a single amino acid to proteins or peptides in a micro-environment.


title: Demand Estimation and Chance-Constrained Fleet Management for Ride Hailing

abstract: In autonomous Mobility on Demand (MOD) systems, customers request rides from a fleet of shared vehicles that can be automatically positioned in response to customer demand. Recent approaches to MOD systems have focused on environments where customers can only request rides through an app or by waiting at a station. This paper develops MOD fleet management approaches for ride hailing, where customers may instead request rides simply by hailing a passing vehicle, an approach of particular importance for campus MOD systems. The challenge for ride hailing is that customer demand is not explicitly provided as it would be with an app, but rather customers are only served if a vehicle happens to be located at the arrival location. This work focuses on maximizing the number of served hailing customers in an MOD system by learning and utilizing customer demand. A Bayesian framework is used to define a novel customer demand model which incorporates observed pedestrian traffic to estimate customer arrival locations with a quantification of uncertainty. An exploration planner is proposed which routes MOD vehicles in order to reduce arrival rate uncertainty. A robust ride hailing fleet management planner is proposed which routes vehicles under the presence of uncertainty using a chance-constrained formulation. Simulation of a real-world MOD system on MIT's campus demonstrates the effectiveness of the planners. The customer demand model and exploration planner are demonstrated to reduce estimation error over time and the ride hailing planner is shown to improve the fraction of served customers in the system by 73\% over a baseline exploration approach.


title: Han's conjecture and Hochschild homology for null-square projective algebras

abstract: Let $\mathcal H$ be the class of algebras verifying Han's conjecture. In this paper we analyse two types of algebras with the aim of providing an inductive step towards the proof of this conjecture. Firstly we show that if an algebra $\Lambda$ is triangular with respect to a system of non necessarily primitive idempotents, and if the algebras at the idempotents belong to $\mathcal H$, then $\Lambda$ is in $\mathcal H$. Secondly we consider a $2\times 2$ matrix algebra, with two algebras on the diagonal, two projective bimodules in the corners, and zero corner products. They are not triangular with respect to the system of the two diagonal idempotents. However, the analogous result holds, namely if both algebras on the diagonal belong to $\mathcal H$, then the algebra itself is in $\mathcal H$.


title: Differences between Doppler velocities of ions and neutral atoms in a solar prominence

abstract: In astrophysical systems with partially ionized plasma the motion of ions is governed by the magnetic field while the neutral particles can only feel the magnetic field's Lorentz force indirectly through collisions with ions. The drift in the velocity between ionized and neutral species plays a key role in modifying important physical processes like magnetic reconnection, damping of magnetohydrodynamic waves, transport of angular momentum in plasma through the magnetic field, and heating. This paper investigates the differences between Doppler velocities of calcium ions and neutral hydrogen in a solar prominence to look for velocity differences between the neutral and ionized species. We simultaneously observed spectra of a prominence over an active region in H I 397 nm, H I 434 nm, Ca II 397 nm, and Ca II 854 nm using a high dispersion spectrograph of the Domeless Solar Telescope at Hida observatory, and compared the Doppler velocities, derived from the shift of the peak of the spectral lines presumably emitted from optically-thin plasma. There are instances when the difference in velocities between neutral atoms and ions is significant, e.g. 1433 events (~ 3 % of sets of compared profiles) with a difference in velocity between neutral hydrogen atoms and calcium ions greater than 3sigma of the measurement error. However, we also found significant differences between the Doppler velocities of two spectral lines emitted from the same species, and the probability density functions of velocity difference between the same species is not significantly different from those between neutral atoms and ions. We interpreted the difference of Doppler velocities as a result of motions of different components in the prominence along the line of sight, rather than the decoupling of neutral atoms from plasma.


title: Covering systems with restricted divisibility

abstract: We prove that every distinct covering system has a modulus divisible by either 2 or 3.


title: Global behaviour of radially symmetric solutions stable at infinity for gradient systems

abstract: This paper is concerned with radially symmetric solutions of systems of the form \[ u_t = -\nabla V(u) + \Delta_x u \] where space variable $x$ and and state-parameter $u$ are multidimensional, and the potential $V$ is coercive at infinity. For such systems, under generic assumptions on the potential, the asymptotic behaviour of solutions "stable at infinity", that is approaching a spatially homogeneous equilibrium when $|x|$ approaches $+\infty$, is investigated. It is proved that every such solutions approaches a stacked family of radially symmetric bistable fronts travelling to infinity. This behaviour is similar to the one of bistable solutions for gradient systems in one unbounded spatial dimension, described in a companion paper. It is expected (but unfortunately not proved at this stage) that behind these travelling fronts the solution again behaves as in the one-dimensional case (that is, the time derivative approaches zero and the solution approaches a pattern of stationary solutions).


title: Scalable Underapproximation for the Stochastic Reach-Avoid Problem for High-Dimensional LTI Systems using Fourier Transforms

abstract: We present a scalable underapproximation of the terminal hitting time stochastic reach-avoid probability at a given initial condition, for verification of high-dimensional stochastic LTI systems. While several approximation techniques have been proposed to alleviate the curse of dimensionality associated with dynamic programming, these techniques are limited and cannot handle larger, more realistic systems. We present a scalable method that uses Fourier transforms to compute an underapproximation of the reach-avoid probability for systems with disturbances with arbitrary probability densities. We characterize sufficient conditions for Borel-measurability of the value functions. We exploit fixed control sequences parameterized by the initial condition (an open-loop control policy) to generate the underapproximation. For Gaussian disturbances, the underapproximation can be obtained using existing efficient algorithms by solving a convex optimization problem. Our approach produces non-trivial lower bounds and is demonstrated on a chain of integrators with 40 states.


title: English Conversational Telephone Speech Recognition by Humans and Machines

abstract: One of the most difficult speech recognition tasks is accurate recognition of human to human communication. Advances in deep learning over the last few years have produced major speech recognition improvements on the representative Switchboard conversational corpus. Word error rates that just a few years ago were 14% have dropped to 8.0%, then 6.6% and most recently 5.8%, and are now believed to be within striking range of human performance. This then raises two issues - what IS human performance, and how far down can we still drive speech recognition error rates? A recent paper by Microsoft suggests that we have already achieved human performance. In trying to verify this statement, we performed an independent set of human performance measurements on two conversational tasks and found that human performance may be considerably better than what was earlier reported, giving the community a significantly harder goal to achieve. We also report on our own efforts in this area, presenting a set of acoustic and language modeling techniques that lowered the word error rate of our own English conversational telephone LVCSR system to the level of 5.5%/10.3% on the Switchboard/CallHome subsets of the Hub5 2000 evaluation, which - at least at the writing of this paper - is a new performance milestone (albeit not at what we measure to be human performance!). On the acoustic side, we use a score fusion of three models: one LSTM with multiple feature inputs, a second LSTM trained with speaker-adversarial multi-task learning and a third residual net (ResNet) with 25 convolutional layers and time-dilated convolutions. On the language modeling side, we use word and character LSTMs and convolutional WaveNet-style language models.


title: Upper Limit on the Milky Way Mass from the Orbit of the Sagittarius Dwarf Satellite

abstract: As one of the most massive Milky Way satellites, the Sagittarius dwarf galaxy has played an important role in shaping the Galactic disk and stellar halo morphologies. The disruption of Sagittarius over several close-in passages has populated the halo of our Galaxy with large-scale tidal streams and offers a unique diagnostic tool for measuring its gravitational potential. Here we test different progenitor mass models for the Milky Way and Sagittarius by modeling the full infall of the satellite. We constrain the mass of the Galaxy based on the observed orbital parameters and multiple tidal streams of Sagittarius. Our semi-analytic modeling of the orbital dynamics agrees with full $N$-body simulations, and favors low values for the Milky Way mass, $\lesssim 10^{12}M_\odot$. This conclusion eases the tension between $\Lambda$CDM and the observed parameters of the Milky Way satellites.