title: Features in the Spectrum of Cosmic-Ray Positrons from Pulsars

abstract: Pulsars have been invoked to explain the origin of recently observed high-energy Galactic cosmic-ray positrons. Since the positron propagation distance decreases with energy, the number of pulsars that can contribute to the observed positrons decreases from $O(10^3)$ for positron energies $E\gtrsim10$ GeV to only a few for $E \gtrsim 500$ GeV. Thus, if pulsars explain these positrons, the positron energy spectrum should become increasingly bumpy at higher energies. Here we present a power-spectrum analysis that can be applied to seek such spectral features in the energy spectrum for cosmic-ray positrons and for the energy spectrum of the combined electron/positron flux. We account for uncertainties in the pulsar distribution by generating hundreds of simulated spectra from pulsar distributions consistent with current observational constraints. Although the current AMS-02 data do not exhibit evidence for spectral features, we find that such features would be detectable in $\simeq 10\%$ of our simulations, with 20 years of AMS-02 data or three years of DAMPE measurements on the electron-plus-positron flux.


title: Generalized massive optimal data compression

abstract: Data compression has become one of the cornerstones of modern astronomical data analysis, with the vast majority of analyses compressing large raw datasets down to a manageable number of informative summaries. In this paper we provide a general procedure for optimally compressing $N$ data down to $n$ summary statistics, where $n$ is equal to the number of parameters of interest. We show that compression to the score function -- the gradient of the log-likelihood with respect to the parameters -- yields $n$ compressed statistics that are optimal in the sense that they preserve the Fisher information content of the data. Our method generalizes earlier work on linear Karhunen-Lo\'{e}ve compression for Gaussian data whilst recovering both lossless linear compression and quadratic estimation as special cases when they are optimal. We give a unified treatment that also includes the general non-Gaussian case as long as mild regularity conditions are satisfied, producing optimal non-linear summary statistics when appropriate. As a worked example, we derive explicitly the $n$ optimal compressed statistics for Gaussian data in the general case where both the mean and covariance depend on the parameters.


title: Cluster realization of positive representations of split real quantum Borel subalgebra

abstract: In our previous work, we studied the positive representations of split real quantum groups $\mathcal{U}_{q\tilde{q}}(\mathfrak{g}_\mathbb{R})$ restricted to its Borel part, and showed that they are closed under taking tensor products. However, the tensor product decomposition was only constructed abstractly using the GNS-representation of a $C^*$-algebraic version of the Drinfeld-Jimbo quantum groups. In this paper, using the recently discovered cluster realization of quantum groups, we write down the decomposition explicitly by realizing it as a sequence of cluster mutations in the corresponding quiver diagram representing the tensor product.


title: ALMACAL III: A combined ALMA and MUSE Survey for Neutral, Molecular, and Ionised Gas in an HI-Absorption-Selected System

abstract: Studying the flow of baryons into and out of galaxies is an important part of understanding the evolution of galaxies over time. We present a detailed case study of the environment around an intervening Ly $\alpha$ absorption line system at $z_{\rm abs} = 0.633$, seen towards the quasar J0423$-$0130 ($z_{\rm QSO} = 0.915$). We detect with ALMA the $^{12}$CO(2--1), $^{12}$CO(3--2) and $1.2$~mm continuum emission from a galaxy at the redshift of the Ly $\alpha$ absorber at a projected distance of $135$ kpc. From the ALMA detections, we infer ISM conditions similar to those in low redshift Luminous Infrared Galaxies. DDT MUSE integral field unit observations reveal the optical counterpart of the $^{12}$CO emission line source and three additional emission line galaxies at the absorber redshift, which together form a galaxy group. The $^{12}$CO emission line detections originate from the most massive galaxy in this group. While we cannot exclude that we miss a fainter host, we reach a dust-uncorrected star-formation rate (SFR) limit of > $0.3 \text{M}_{\odot} \text{ yr}^{-1}$ within $100$ kpc from the sightline to the background quasar. We measure the dust-corrected SFR (ranging from $3$ to $50$ M$_{\odot}$ yr$^{-1}$), the morpho-kinematics and the metallicities of the four group galaxies to understand the relation between the group and the neutral gas probed in absorption. We find that the Ly $\alpha$ absorber traces either an outflow from the most massive galaxy or intra-group gas. This case study illustrates the power of combining ALMA and MUSE to obtain a census of the cool baryons in a bounded structure at intermediate redshift.


title: Cavity quantum electrodynamics in the non-perturbative regime

abstract: We study a generic cavity-QED system where a set of (artificial) two-level dipoles is coupled to the electric field of a single-mode LC resonator. This setup is used to derive a minimal quantum mechanical model for cavity QED, which accounts for both dipole-field and direct dipole-dipole interactions. The model is applicable for arbitrary coupling strengths and allows us to extend the usual Dicke model into the non-perturbative regime of QED, where the dipole-field interaction can be associated with an effective finestructure constant of order unity. In this regime, we identify three distinct classes of normal, superradiant and subradiant vacuum states and discuss their characteristic properties and the transitions between them. Our findings reconcile many of the previous, often contradictory predictions in this field and establish a common theoretical framework to describe ultrastrong coupling phenomena in a diverse range of cavity-QED platforms.


title: The abundance, distribution, and physical nature of highly ionized oxygen OVI, OVII, and OVIII in IllustrisTNG

abstract: We explore the abundance, spatial distribution, and physical properties of the OVI, OVII, and OVIII ions of oxygen in circumgalactic and intergalactic media (the CGM, IGM, and WHIM). We use the TNG100 and TNG300 large volume cosmological magneto-hydrodynamical simulations. Modeling the ionization states of simulated oxygen, we find good agreement with observations of the low-redshift OVI column density distribution function (CDDF), and present its evolution for all three ions from z=0 to z=4. Producing mock quasar absorption line spectral surveys, we show that the IllustrisTNG simulations are fully consistent with constraints on the OVI content of the CGM from COS-Halos and other low redshift observations, producing columns as high as observed. We measure the total amount of mass and average column densities of each ion using hundreds of thousands of simulated galaxies spanning 10^11 < Mhalo/Msun < 10^15 corresponding to 10^9 < M*/Msun < 10^12 in stellar mass. Stacked radial profiles of OVI are computed in 3D number density and 2D projected column, decomposing into the 1-halo and 2-halo terms. Relating halo OVI to properties of the central galaxy, we find a correlation between the (g-r) color of a galaxy and the total amount of OVI in its CGM. In comparison to the COS-Halos finding, this leads to a dichotomy of columns around star-forming versus passive galaxies at fixed stellar (or halo) mass. We demonstrate that this correlation is a direct result of blackhole feedback associated with quenching and represents a causal consequence of galactic-scale baryonic feedback impacting the physical state of the circumgalactic medium.


title: Intensity-Corrected Herschel Observations of Nearby Isolated Low-Mass Clouds

abstract: We present intensity-corrected Herschel maps at 100 um, 160 um, 250 um, 350 um, and 500 um for 56 isolated low-mass clouds. We determine the zero-point corrections for Herschel PACS and SPIRE maps from the Herschel Science Archive (HSA) using Planck data. Since these HSA maps are small, we cannot correct them using typical methods. Here, we introduce a technique to measure the zero-point corrections for small Herschel maps. We use radial profiles to identify offsets between the observed HSA intensities and the expected intensities from Planck. Most clouds have reliable offset measurements with this technique. In addition, we find that roughly half of the clouds have underestimated HSA-SPIRE intensities in their outer envelopes relative to Planck, even though the HSA-SPIRE maps were previously zero-point corrected. Using our technique, we produce corrected Herschel intensity maps for all 56 clouds and determine their line-of-sight average dust temperatures and optical depths from modified black body fits. The clouds have typical temperatures of ~ 14-20 K and optical depths of ~ 1e-5 - 1e-3. Across the whole sample, we find an anti-correlation between temperature and optical depth. We also find lower temperatures than what was measured in previous Herschel studies, which subtracted out a background level from their intensity maps to circumvent the zero-point correction. Accurate Herschel observations of clouds are key to obtain accurate density and temperature profiles. To make such future analyses possible, intensity-corrected maps for all 56 clouds are publicly available in the electronic version.


title: First study of reionization in the Planck 2015 normalized closed $\Lambda$CDM inflation model

abstract: We study reionization in two non-flat $\Lambda$CDM inflation models that best fit the Planck 2015 cosmic microwave background anisotropy observations, ignoring or in conjunction with baryon acoustic oscillation distance measurements. We implement a principal component analysis (PCA) to estimate the uncertainties in the reionization history from a joint quasar-CMB dataset. A thorough Markov Chain Monte Carlo analysis is done over the parameter space of PCA modes for both non-flat $\Lambda$CDM inflation models as well as the original Planck 2016 tilted, spatially-flat $\Lambda$CDM inflation model. Although both flat and non-flat models can closely match the low-redshift ($z\lesssim6$) observations, we notice a possible tension between high-redshift ($z\sim8$) Lyman-$\alpha$ emitter data and the non-flat models. This is solely due to the fact that the closed models have a relatively higher reionization optical depth compared to the flat one, which in turn demands more high-redshift ionizing sources and favors an extended reionization starting as early as $z\approx14$. We conclude that as opposed to flat-cosmology, for the non-flat cosmology models (i) the escape fraction needs steep redshift evolution and even unrealistically high values at some redshifts and (ii) most of the physical parameters require to have non-monotonic redshift evolution, especially apparent when Lyman-$\alpha$ emitter data is included in the analysis.


title: The E-MOSAICS Project: simulating the formation and co-evolution of galaxies and their star cluster populations

abstract: We introduce the MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE (E-MOSAICS) project. E-MOSAICS incorporates models describing the formation, evolution and disruption of star clusters into the EAGLE galaxy formation simulations, enabling the examination of the co-evolution of star clusters and their host galaxies in a fully cosmological context. A fraction of the star formation rate of dense gas is assumed to yield a cluster population; this fraction, and the population's initial properties, are governed by the physical properties of the natal gas. The subsequent evolution and disruption of the entire cluster population is followed accounting for two-body relaxation, stellar evolution, and gravitational shocks induced by the local tidal field. This introductory paper presents a detailed description of the model and initial results from a suite of 10 simulations of $\sim L^\star$ galaxies with disc-like morphologies at $z=0$. The simulations broadly reproduce key observed characteristics of young star clusters and globular clusters (GCs), without invoking separate formation mechanisms for each population. The simulated GCs are the surviving population of massive clusters formed at early epochs ($z\gtrsim1-2$), when the characteristic pressures and surface densities of star-forming gas were significantly higher than observed in local galaxies. We examine the influence of the star formation and assembly histories of galaxies on their cluster populations, finding that (at similar present-day mass) earlier-forming galaxies foster a more massive and disruption-resilient cluster population, while galaxies with late mergers are capable of forming massive clusters even at late cosmic epochs. (Abridged)


title: Three-dimensional dualities with bosons and fermions

abstract: We propose new infinite families of non-supersymmetric IR dualities in three space-time dimensions, between Chern-Simons gauge theories (with classical gauge groups) with both scalars and fermions in the fundamental representation. In all cases we study the phase diagram as we vary two relevant couplings, finding interesting lines of phase transitions. In various cases the dualities lead to predictions about multi-critical fixed points and the emergence of IR quantum symmetries. For unitary groups we also discuss the coupling to background gauge fields and the map of simple monopole operators.


title: Confronting Models of Massive Star Evolution and Explosions with Remnant Mass Measurements

abstract: The mass distribution of compact objects provides a fossil record that can be studied to uncover information on the late stages of massive star evolution, the supernova explosion mechanism, and the dense matter equation of state. Observations of neutron star masses indicate a bimodal Gaussian distribution, while the observed black hole mass distribution decays exponentially for stellar-mass black holes. We use these observed distributions to directly confront the predictions of stellar evolution models and the neutrino-driven supernova simulations of Sukhbold et al. (2016). We find excellent agreement between the black hole and low-mass neutron star distributions created by these simulations and the observations. We show that a large fraction of the stellar envelope must be ejected, either during the formation of stellar-mass black holes or prior to the implosion through tidal stripping due to a binary companion, in order to reproduce the observed black hole mass distribution. We also determine the origins of the bimodal peaks of the neutron star mass distribution, finding that the low-mass peak (centered at ~1.4 M_sun) originates from progenitors with M_zams ~ 9-18 M_sun. The simulations fail to reproduce the observed peak of high-mass neutron stars (centered at ~1.8 M_sun) and we explore several possible explanations. We argue that the close agreement between the observed and predicted black hole and low-mass neutron star mass distributions provides new promising evidence that these stellar evolution and explosion models are accurately capturing the relevant stellar, nuclear, and explosion physics involved in the formation of compact objects.


title: Finding binaries from phase modulation of pulsating stars with Kepler: V. Orbital parameters, with eccentricity and mass-ratio distributions of 341 new binaries

abstract: The orbital parameters of binaries at intermediate periods ($10^2$--$10^3$ d) are difficult to measure with conventional methods and are very incomplete. We have undertaken a new survey, applying our pulsation timing method to {\it Kepler} light curves of 2224 main-sequence A/F stars and found 341 non-eclipsing binaries. We calculate the orbital parameters for 317 PB1 systems (single-pulsator binaries) and 24 PB2s (double-pulsators). The method reaches down to small mass ratios $q=0.02$ and yields a highly homogeneous sample. We parametrize the mass-ratio distribution using both inversion and MCMC forward-modelling techniques, and find it to be skewed towards low-mass companions, peaking at $q=0.2$. While solar-type primaries show a brown dwarf desert across intermediate periods, we find a small but significant (2.6$\sigma$) population of extreme-mass-ratio companions ($q<0.1$) to our intermediate-mass primaries. Across periods of 100--1500 d and at $q>0.1$, we measure the binary fraction of current A/F primaries to be 15.4%$\pm$1.4%, though we find that a large fraction of the companions (21%$\pm$6%) are white dwarfs in post-mass-transfer systems with primaries that are now blue stragglers, some of which are the progenitors of Type Ia supernovae, barium stars, symbiotics, and related phenomena. Excluding these white dwarfs, we find the binary fraction of original A/F primaries to be 13.9%$\pm$2.1% over the same parameter space. Combining literature measurements with ours, we find the binary fraction across these periods is a constant 5% for primaries $M_1<0.8$ M$_{\odot}$, but then increases linearly with $\log M_1$, demonstrating that natal discs around more massive protostars $M_1>1$ M$_{\odot}$ become increasingly more prone to fragmentation. Finally, we find the eccentricity distribution of the main-sequence pairs to be much less eccentric than the thermal distribution.


title: Galactic nuclei evolution with spinning black holes: method and implementation

abstract: Supermassive black holes at the centre of galactic nuclei mostly grow in mass through gas accretion over cosmic time. This process also modifies the angular momentum (or spin) of black holes, both in magnitude and in orientation. Despite being often neglected in galaxy formation simulations, spin plays a crucial role in modulating accretion power, driving jet feedback, and determining recoil velocity of coalescing black hole binaries. We present a new accretion model for the moving-mesh code {\sc arepo} that incorporates (i) mass accretion through a thin $\alpha$-disc, and (ii) spin evolution through the Bardeen-Petterson effect. We use a diverse suite of idealised simulations to explore the physical connection between spin evolution and larger scale environment. We find that black holes with mass $\lesssim 10^{7}$ M$_{\odot}$ experience quick alignment with the accretion disc. This favours prolonged phases of spin-up, and the spin direction evolves according to the gas inflow on timescales as short as $\lesssim 100$ Myr, which might explain the observed jet direction distribution in Seyfert galaxies. Heavier black holes ($\gtrsim 10^{8}$ M$_{\odot}$) are instead more sensitive to the local gas kinematic. Here we find a wider distribution in spin magnitudes: spin-ups are favoured if gas inflow maintains a preferential direction, and spin-downs occur for nearly isotropic infall, while the spin direction does not change much over short timescales $\sim 100$ Myr. We therefore conclude that supermassive black holes with masses $\gtrsim 5 \times 10^{8}$ M$_{\odot}$ may be the ideal testbed to determine the main mode of black hole fuelling over cosmic time.


title: Using the CIFIST grid of CO5BOLD 3D model atmospheres to study the effects of stellar granulation on photometric colours. I. Grids of 3D corrections in the UBVRI, 2MASS, Hipparcos, Gaia, and SDSS systems

abstract: The atmospheres of cool stars are temporally and spatially inhomogeneous due to the effects of convection.The influence of this inhomogeneity, referred to as granulation, on colours has never been investigated over a large range of effective temperatures and gravities. We aim to study, in a quantitative way, the impact of granulation on colours. We use the CIFIST grid of CO5BOLD hydrodynamical models to compute emerging fluxes. These in turn are used to compute theoretical colours in the UBVRI,2MASS,Hipparcos,Gaia and SDSS systems. Every CO5BOLD model has a corresponding one dimensional (1D) plane-parallel LHD model computed for the same atmospheric parameters, which we used to define a "3D correction" that can be applied to colours computed from fluxes computed from any 1D model atmosphere code. The 3D corrections on colours are generally small, of the order of a few hundredths of a magnitude, yet they are far from negligible. We find that ignoring granulation effects can lead to underestimation of Teff by up to 200K and overestimation of gravity by up to 0.5dex, when using colours as diagnostics. We have identified a major shortcoming in how scattering is treated in the current version of the CIFIST grid, which could lead to offsets of the order 0.01mag, especially for colours involving blue and UV bands.We have investigated the Gaia and Hipparcos photometric systems and found that the (G-H_p),(BP-RP) diagram is immune to the effects of granulation. In addition, we point to the potential of the RVS photometry as a metallicity diagnostic. Our investigation shows that the effects of granulation should not be neglected if one wants to use colours as diagnostics of the stellar parameters of F,G,K stars. A limitation is that scattering is treated as true absorption in our current computations, thus our 3D corrections are likely an upper limit to the true effect. (Abridged)


title: A Spatially Resolved Radio Spectral Index Study of the Dwarf Irregular Galaxy NGC\,1569

abstract: We study the resolved radio-continuum spectral energy distribution of the dwarf irregular galaxy, NGC 1569, on a beam-by-beam basis to isolate and study its spatially resolved radio emission characteristics. Utilizing high quality NRAO Karl G. Jansky Very Large Array (VLA) observations that densely sample the 1--34\,GHz frequency range, we adopt a Bayesian fitting procedure, where we use H$\alpha$ emission that has not been corrected for extinction as a prior, to produce maps of how the separated thermal emission, non-thermal emission and non-thermal spectral index vary across NGC\,1569's main disk. We find a higher thermal fraction at 1\,GHz than is found in spiral galaxies ($26^{+2}_{-3}\%$) and find an average non-thermal spectral index $\alpha = -0.53\pm0.02$, suggesting that a young population of cosmic ray electrons is responsible for the observed non--thermal emission. By comparing our recovered map of the thermal radio emission with literature H$\alpha$ maps, we estimate the total reddening along the line of sight to NGC\,1569 to be $E(B-V) = 0.49 \pm 0.05$, which is in good agreement with other literature measurements. Spatial variations in the reddening indicate that a significant portion of the total reddening is due to internal extinction within NGC\,1569.


title: Eccentricity evolution during planet-disc interaction

abstract: During the process of planet formation, the planet-discs interactions might excite (or damp) the orbital eccentricity of the planet. In this paper, we present two long ($t\sim 3\times 10^5$ orbits) numerical simulations: (a) one (with a relatively light disc, $M_{\rm d}/M_{\rm p}=0.2$) where the eccentricity initially stalls before growing at later times and (b) one (with a more massive disc, $M_{\rm d}/M_{\rm p}=0.65$) with fast growth and a late decrease of the eccentricity. We recover the well-known result that a more massive disc promotes a faster initial growth of the planet eccentricity. However, at late times the planet eccentricity decreases in the massive disc case, but increases in the light disc case. Both simulations show periodic eccentricity oscillations superimposed on a growing/decreasing trend and a rapid transition between fast and slow pericentre precession. The peculiar and contrasting evolution of the eccentricity of both planet and disc in the two simulations can be understood by invoking a simple toy model where the disc is treated as a second point-like gravitating body, subject to secular planet-planet interaction and eccentricity pumping/damping provided by the disc. We show how the counterintuitive result that the more massive simulation produces a lower planet eccentricity at late times can be understood in terms of the different ratios of the disc-to-planet angular momentum in the two simulations. In our interpretation, at late times the planet eccentricity can increase more in low-mass discs rather than in high-mass discs, contrary to previous claims in the literature.


title: SN 2017dio: a type-Ic supernova exploding in a hydrogen-rich circumstellar medium

abstract: SN 2017dio shows both spectral characteristics of a type-Ic supernova (SN) and signs of a hydrogen-rich circumstellar medium (CSM). Prominent, narrow emission lines of H and He are superposed on the continuum. Subsequent evolution revealed that the SN ejecta are interacting with the CSM. The initial SN Ic identification was confirmed by removing the CSM interaction component from the spectrum and comparing with known SNe Ic, and reversely, adding a CSM interaction component to the spectra of known SNe Ic and comparing them to SN 2017dio. Excellent agreement was obtained with both procedures, reinforcing the SN Ic classification. The light curve constrains the pre-interaction SN Ic peak absolute magnitude to be around $M_g = -17.6$ mag. No evidence of significant extinction is found, ruling out a brighter luminosity required by a SN Ia classification. These pieces of evidence support the view that SN 2017dio is a SN Ic, and therefore the first firm case of a SN Ic with signatures of hydrogen-rich CSM in the early spectrum. The CSM is unlikely to have been shaped by steady-state stellar winds. The mass loss of the progenitor star must have been intense, $\dot{M} \sim 0.02$ $(\epsilon_{H\alpha}/0.01)^{-1}$ $(v_\textrm{wind}/500$ km s$^{-1}$) $(v_\textrm{shock}/10 000$ km s$^{-1})^{-3}$ $M_\odot$~yr$^{-1}$, peaking at a few decades before the SN. Such a high mass loss rate might have been experienced by the progenitor through eruptions or binary stripping.


title: Feature discovery and visualization of robot mission data using convolutional autoencoders and Bayesian nonparametric topic models

abstract: The gap between our ability to collect interesting data and our ability to analyze these data is growing at an unprecedented rate. Recent algorithmic attempts to fill this gap have employed unsupervised tools to discover structure in data. Some of the most successful approaches have used probabilistic models to uncover latent thematic structure in discrete data. Despite the success of these models on textual data, they have not generalized as well to image data, in part because of the spatial and temporal structure that may exist in an image stream. We introduce a novel unsupervised machine learning framework that incorporates the ability of convolutional autoencoders to discover features from images that directly encode spatial information, within a Bayesian nonparametric topic model that discovers meaningful latent patterns within discrete data. By using this hybrid framework, we overcome the fundamental dependency of traditional topic models on rigidly hand-coded data representations, while simultaneously encoding spatial dependency in our topics without adding model complexity. We apply this model to the motivating application of high-level scene understanding and mission summarization for exploratory marine robots. Our experiments on a seafloor dataset collected by a marine robot show that the proposed hybrid framework outperforms current state-of-the-art approaches on the task of unsupervised seafloor terrain characterization.


title: The HectoMAP Cluster Survey - II. X-ray Clusters

abstract: We apply a friends-of-friends algorithm to the HectoMAP redshift survey and cross-identify associated X-ray emission in the ROSAT All-Sky Survey data (RASS). The resulting flux limited catalog of X-ray cluster survey is complete to a limiting flux of $\sim3 \times10^{-13}$ erg s$^{-1}$ cm$^{-2}$ and includes 15 clusters (7 newly discovered) with redshift $z \leq 0.4$. HectoMAP is a dense survey ($\sim1200$ galaxies deg$^{-2}$) that provides $\sim50$ members (median) in each X-ray cluster. We provide redshifts for the 1036 cluster members. Subaru/Hyper Suprime-Cam imaging covers three of the X-ray systems and confirms that they are impressive clusters. The HectoMAP X-ray clusters have an $L_{X} - {\sigma}_{cl}$ scaling relation similar to that of known massive X-ray clusters. The HectoMAP X-ray cluster sample predicts $\sim 12000 \pm3000$ detectable X-ray clusters in the RASS to the limiting flux, comparable with previous estimates.


title: Resource Sharing of a Computing Access Point for Multi-user Mobile Cloud Offloading with Delay Constraints

abstract: We consider a mobile cloud computing system with multiple users, a remote cloud server, and a computing access point (CAP). The CAP serves both as the network access gateway and a computation service provider to the mobile users. It can either process the received tasks from mobile users or offload them to the cloud. We jointly optimize the offloading decisions of all users, together with the allocation of computation and communication resources, to minimize the overall cost of energy consumption, computation, and maximum delay among users. The joint optimization problem is formulated as a mixed-integer program. We show that the problem can be reformulated and transformed into a non-convex quadratically constrained quadratic program, which is NP-hard in general. We then propose an efficient solution to this problem by semidefinite relaxation and a novel randomization mapping method. Furthermore, when there is a strict delay constraint for processing each user's task, we further propose a three-step algorithm to guarantee the feasibility and local optimality of the obtained solution. Our simulation results show that the proposed solutions give nearly optimal performance under a wide range of parameter settings, and the addition of a CAP can significantly reduce the cost of multi-user task offloading compared with conventional mobile cloud computing where only the remote cloud server is available.


title: Heavy X-ray obscuration in the most-luminous galaxies discovered by WISE

abstract: Hot Dust-Obscured Galaxies (Hot DOGs) are hyperluminous ($L_{\mathrm{8-1000\,\mu m}}>10^{13}\,\mathrm{L_\odot}$) infrared galaxies with extremely high (up to hundreds of K) dust temperatures. The sources powering both their extremely high luminosities and dust temperatures are thought to be deeply buried and rapidly accreting supermassive black holes (SMBHs). Hot DOGs could therefore represent a key evolutionary phase in which the SMBH growth peaks. X-ray observations can be used to study their obscuration levels and luminosities. In this work, we present the X-ray properties of the 20 most-luminous ($L_{\mathrm{bol}}\gtrsim10^{14}\, L_\odot$) known Hot DOGs at $z=2-4.6$. Five of them are covered by long-exposure ($10-70$ ks) Chandra and XMM-Newton observations, with three being X-ray detected, and we study their individual properties. One of these sources (W0116$-$0505) is a Compton-thick candidate, with column density $N_H=(1.0-1.5)\times10^{24}\,\mathrm{cm^{-2}}$ derived from X-ray spectral fitting. The remaining 15 Hot DOGs have been targeted by a Chandra snapshot (3.1 ks) survey. None of these 15 is individually detected; therefore we applied a stacking analysis to investigate their average emission. From hardness-ratio analysis, we constrained the average obscuring column density and intrinsic luminosity to be log$N_H\,\mathrm{[cm^{-2}]}>23.5$ and $L_X\gtrsim10^{44}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$, which are consistent with results for individually detected sources. We also investigated the $L_X-L_{6\mu\mathrm{m}}$ and $L_X-L_{bol}$ relations, finding hints that Hot DOGs are typically X-ray weaker than expected, although larger samples of luminous obscured QSOs are needed to derive solid conclusions.


title: Paris-Lille-3D: a large and high-quality ground truth urban point cloud dataset for automatic segmentation and classification

abstract: This paper introduces a new Urban Point Cloud Dataset for Automatic Segmentation and Classification acquired by Mobile Laser Scanning (MLS). We describe how the dataset is obtained from acquisition to post-processing and labeling. This dataset can be used to learn classification algorithm, however, given that a great attention has been paid to the split between the different objects, this dataset can also be used to learn the segmentation. The dataset consists of around 2km of MLS point cloud acquired in two cities. The number of points and range of classes make us consider that it can be used to train Deep-Learning methods. Besides we show some results of automatic segmentation and classification. The dataset is available at: http://caor-mines-paristech.fr/fr/paris-lille-3d-dataset/


title: Strongly Interacting Matter Phase Diagram in the presence of Magnetic Fields in an Extended Effective Lagrangian Approach with Explicit Chiral Symmetry Breaking Interactions

abstract: Extensions of the NJL model which go beyond the original 4-quark interaction, which drives the dynamical mass generation, have proven to be quite successful in describing low energy hadronic phenomenology. The inclusion of 8-quark interaction terms solved a metastability problem of the effective potential introduced by the inclusion of the 6-quark 't Hooft determinant term in the 3-flavor version of the model (needed to eliminate the unwanted U(1) axial symmetry) . This model, that has proven to be quite powerful and feature-rich, has been expanded to include all the spin-0 terms, without and with explicit chiral symmetry breaking, which are of the same order as the 't Hooft flavor determinant in a 1/Nc expansion resulting in an unprecedented success in reproducing the low lying scalar and pseudoscalar meson spectra. This success can be seen as a result of the inclusion of the full chiral symmetry breaking pattern. The two critical endpoints which are obtained in the temperature/chemical potential phase diagram are shifted to lower chemical potential and higher temperature when the effect of magnetic field is taken into account. For the studied magnetic field strengths (in the range $eH=0-0.4~GeV^2$) no significant extra transitions are seen to appear.


title: Thermodynamic and kinetic fragility of Freon113: the most fragile plastic crystal

abstract: We present a dynamic and thermodynamic study of the orientational glass former Freon113 (CCl2F-CClF2) in order to analyze its kinetic and thermodynamic fragilities. Freon113 displays internal molecular degrees of freedom which promote a complex energy landscape. Experimental specific heat and its microscopic origin, the vibrational density of states from inelastic neutron scattering, together with the orientational dynamics obtained by means of dielectric spectroscopy have revealed the highest fragility value, both thermodynamic and kinetic, found for this orientational glass former. The excess in both Debye-reduced specific heat and density of states (boson peak) evidences the existence of glassy low-energy excitations. We demonstrate that early proposed correlations between the boson peak and the Debye specific heat value are elusive as revealed by the clear counterexample of the studied case.


title: Glassy anomalies in the low-temperature thermal properties of a minimally disordered crystalline solid

abstract: The low-temperature thermal and transport properties of an unusual kind of crystal exhibiting minimal molecular positional and tilting disorder have been measured. The material, namely, low-dimensional, highly anisotropic pentachloronitrobenzene (PCNB) has a layered structure of rhombohedral parallel planes in which the molecules execute large-amplitude in-plane as well as concurrent out-of-plane librational motions. Our study reveals that low-temperature glassy anomalies can be found in a system with minimal disorder due to the freezing of (mostly in-plane) reorientational jumps of molecules between equivalent crystallographic positions with partial site occupation. Our findings will pave the way to a deeper understanding of the origin of above-mentioned universal glassy properties at low temperature.


title: Robust Stereo Visual Inertial Odometry for Fast Autonomous Flight

abstract: In recent years, vision-aided inertial odometry for state estimation has matured significantly. However, we still encounter challenges in terms of improving the computational efficiency and robustness of the underlying algorithms for applications in autonomous flight with micro aerial vehicles in which it is difficult to use high quality sensors and pow- erful processors because of constraints on size and weight. In this paper, we present a filter-based stereo visual inertial odometry that uses the Multi-State Constraint Kalman Filter (MSCKF) [1]. Previous work on stereo visual inertial odometry has resulted in solutions that are computationally expensive. We demonstrate that our Stereo Multi-State Constraint Kalman Filter (S-MSCKF) is comparable to state-of-art monocular solutions in terms of computational cost, while providing signifi- cantly greater robustness. We evaluate our S-MSCKF algorithm and compare it with state-of-art methods including OKVIS, ROVIO, and VINS-MONO on both the EuRoC dataset, and our own experimental datasets demonstrating fast autonomous flight with maximum speed of 17.5m/s in indoor and outdoor environments. Our implementation of the S-MSCKF is available at https://github.com/KumarRobotics/msckf_vio.


title: Electron Flavored Dark Matter

abstract: In this paper we investigate the phenomenology of the electron flavored Dirac dark matter with two types of portal interactions. We analyze constraints from the electron magnetic moment anomaly, LHC searches of singly charged scalar, dark matter relic abundance as well as direct and indirect detections. Our study shows that the available parameter space is quite constrained, but there are parameter space that is compatible with the current data. We further show that the DAMPE cosmic ray electron excess, which indicates cosmic ray excess at around 1.5 TeV, can be interpreted as the annihilation of dark matter into electron positron pairs in this model.


title: Augmented Minimax Linear Estimation

abstract: Many statistical estimands can expressed as continuous linear functionals of a conditional expectation function. This includes the average treatment effect under unconfoundedness and generalizations for continuous-valued and personalized treatments. In this paper, we discuss a general approach to estimating such quantities: we begin with a simple plug-in estimator based on an estimate of the conditional expectation function, and then correct the plug-in estimator by subtracting a minimax linear estimate of its error. We show that our method is semiparametrically efficient under weak conditions and observe promising performance on both real and simulated data.


title: Effective interactions in a quantum Bose-Bose mixture

abstract: We generalize the Beliaev diagrammatic theory of an interacting spinless Bose-Einstein condensate to the case of a binary mixture. We derive a set of coupled Dyson equations and find analytically the Green's functions of the system. The elementary excitation spectrum consists of two branches, one of which takes the characteristic parabolic form in the limit of a spin-independent interaction. We observe renormalization of the magnon mass and the spin-wave velocity due to the Andreev-Bashkin entrainment effect. For a 3D weakly-interacting gas the spectrum can be obtained by applying the Bogoliubov transformation to the second-quantized Hamiltonian in which the microscopic two-body potentials in each channel are replaced by the corresponding off-shell scattering amplitudes. The superfluid drag density can be calculated by considering a mixture of phonons and magnons interacting via the effective potentials. We show that this problem is identical to the second-order perturbative treatment of a Bose polaron. In 2D the drag contributes to the magnon dispersion already in the first approximation. Our consideration provides a basis for systematic study of emergent phases in quantum degenerate Bose-Bose mixtures.


title: Manifestly Lorentz-covariant variables for the phase space of general relativity

abstract: We present a manifestly Lorentz-covariant description of the phase space of general relativity with the Immirzi parameter. This formulation emerges after solving the second-class constraints arising in the canonical analysis of the Holst action. We show that the new canonical variables give rise to other Lorentz-covariant parametrizations of the phase space via canonical transformations. The resulting form of the first-class constraints in terms of new variables is given. In the time gauge, these variables and the constraints become those found by Barbero.


title: Searching for Extragalactic Sources in the VISTA Variables in the V\'ia L\'actea Survey

abstract: We search for extragalactic sources in the VISTA Variables in the V\'ia L\'actea survey that are hidden by the Galaxy. Herein, we describe our photometric procedure to find and characterize extragalactic objects using a combination of SExtractor and PSFEx. It was applied in two tiles of the survey: d010 and d115, without previous extragalactic IR detections, in order to obtain photometric parameters of the detected sources. The adopted criteria to define extragalactic candidates include CLASS_STAR < 0.3; 1.0 < R1/2 < 5.0 arcsec; 2.1 < C < 5; and Phi > 0.002 and the colors: 0.5 < (J - K_s) < 2.0 mag; 0.0 < (J - H) < 1.0 mag; 0.0 < (H - K_s) < 2.0 mag and (J - H) + 0.9 (H - K_s) > 0.44 mag. We detected 345 and 185 extragalactic candidates in the d010 and d115 tiles, respectively. All of them were visually inspected and confirmed to be galaxies. In general, they are small and more circular objects, due to the near-IR sensitivity to select more compact objects with higher surface brightness. The procedure will be used to identify extragalactic objects in other tiles of the VVV disk, which will allow us to study the distribution of galaxies and filaments hidden by the Milky Way.


title: Regularization of non-normal matrices by Gaussian noise - the banded Toeplitz and twisted Toeplitz cases

abstract: We consider the spectrum of additive, polynomially vanishing random perturbations of deterministic matrices, as follows. Let $M_N$ be a deterministic $N\times N$ matrix, and let $G_N$ be a complex Ginibre matrix. We consider the matrix $\mathcal{M}_N=M_N+N^{-\gamma}G_N$, where $\gamma>1/2$. With $L_N$ the empirical measure of eigenvalues of $\mathcal{M}_N$, we provide a general deterministic equivalence theorem that ties $L_N$ to the singular values of $z-M_N$, with $z\in \mathbb{C}$. We then compute the limit of $L_N$ when $M_N$ is an upper triangular Toeplitz matrix of finite symbol: if $M_N=\sum_{i=0}^{\mathfrak{d}} a_i J^i$ where $\mathfrak{d}$ is fixed, $a_i\in\mathcal{C}$ are deterministic scalars and $J$ is the nilpotent matrix $J(i,j)={\bf 1}_{j=i+1}$, then $L_N$ converges, as $N\to\infty$, to the law of $\sum_{i=0}^{\mathfrak{d}} a_i U^i$ where $U$ is a uniform random variable on the unit circle in the complex plane. We also consider the case of slowly varying diagonals (twisted Toeplitz matrices), and, when $\mathfrak{d}=1$, also of i.i.d.~entries on the diagonals in $M_N$.


title: A Color Intensity Invariant Low Level Feature Optimization Framework for Image Quality Assessment

abstract: Image Quality Assessment (IQA) algorithms evaluate the perceptual quality of an image using evaluation scores that assess the similarity or difference between two images. We propose a new low-level feature based IQA technique, which applies filter-bank decomposition and center-surround methodology. Differing from existing methods, our model incorporates color intensity adaptation and frequency scaling optimization at each filter-bank level and spatial orientation to extract and enhance perceptually significant features. Our computational model exploits the concept of object detection and encapsulates characteristics proposed in other IQA algorithms in a unified architecture. We also propose a systematic approach to review the evolution of IQA algorithms using unbiased test datasets, instead of looking at individual scores in isolation. Experimental results demonstrate the feasibility of our approach.


title: Graph Centrality Measures for Boosting Popularity-Based Entity Linking

abstract: Many Entity Linking systems use collective graph-based methods to disambiguate the entity mentions within a document. Most of them have focused on graph construction and initial weighting of the candidate entities, less attention has been devoted to compare the graph ranking algorithms. In this work, we focus on the graph-based ranking algorithms, therefore we propose to apply five centrality measures: Degree, HITS, PageRank, Betweenness and Closeness. A disambiguation graph of candidate entities is constructed for each document using the popularity method, then centrality measures are applied to choose the most relevant candidate to boost the results of entity popularity method. We investigate the effectiveness of each centrality measure on the performance across different domains and datasets. Our experiments show that a simple and fast centrality measure such as Degree centrality can outperform other more time-consuming measures.


title: Categorical Logarithmic Hodge Theory, I

abstract: We write down a new "logarithmic" quasicoherent category $\operatorname{Qcoh}_{log}(U, X, D)$ attached to a smooth open algebraic variety $U$ with toroidal compactification $X$ and boundary divisor $D$. This is a (large) symmetric monoidal Abelian category, which we argue can be thought of as the categorical substrate for logarithmic Hodge theory of $U$. We show that its Hochschild homology theory coincides with the theory of log-forms on $X$ with logarithmic structure induced by $D$, and in particular, that the noncommutative Hodge-to de Rham sequence on $\operatorname{Qcoh}_{log}(U, X, D)$ recovers known log Hodge structure on the de Rham cohomology of the open variety $U$. As an application, we compute the Hochschild homology of the category of coherent sheaves on the infinite root stack of Talpo and Vistoli in the toroidal setting. We prove a derived invariance result for this theory: namely, that strictly toroidal changes of compactification do not change the derived category of $\operatorname{Qcoh}_{log}(U, X, D)$. The definition is motivated by the coherent object appearing in the author's microlocal mirror symmetry result [20]. In this paper, the first in a series, we work over an algebraically closed field of characteristic zero. The next installment will develop the characteristic p and mixed-characteristic theories.


title: Scattering using real-time path integrals

abstract: Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. Purpose: To test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals. Methods: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on a space of continuous paths with respect to a complex probability distribution. The method has the advantage that it can be applied to arbitrary short-range potentials. Results: The new method is tested by applying it to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half shell transition operator matrix elements are in agreement with a numerical solution of the Lippmann-Schwinger equation. The computational method has a straightforward generalization to more complicated systems.


title: On Optimal Stochastic Ballistic Transports

abstract: For a given Lagrangian $L:[0,T]\times M\times M^\ast\rightarrow \mathbb{R}_+$ and probability measures $\mu\in\mathcal{P}(M^\ast)$, $\nu\in \mathcal{P}(M)$, we introduce the stochastic ballistic transportation problems \begin{align}\tag{$\star$} \underline{B}(\mu,\nu):=\inf\left\{\mathbb{E}\left[\langle V,X_0\rangle +\int_0^T L(t,X,\beta(t,X))\,dt\right]\middle\rvert V\sim\mu,X_T\sim \nu\right\}\\\tag{$\star\star$} \overline{B}(\nu,\mu):=\sup\left\{\mathbb{E}\left[\langle V,X_T\rangle -\int_0^T L(t,X,\beta(t,X))\,dt\right]\middle\rvert V\sim\mu,X_0\sim \nu\right\} \end{align} where $X$ is a diffusion process with drift $\beta$. This cost is based on the stochastic optimal transport problem presented by Mikami and the deterministic ballistic transport introduced by Ghoussoub. We obtain a Kantorovich-style duality result that reformulates this problem in terms of solutions to the Hamilton-Jacobi-Bellman equation \begin{equation*} \frac{\partial\phi}{\partial t}+\frac{1}{2}\Delta \phi+H(t,x,\nabla\phi)=0, \end{equation*} and show how optimal processes may be thereby attained.


title: Investigation of Gaze Patterns in Multi View Laparoscopic Surgery

abstract: Laparoscopic Surgery (LS) is a modern surgical technique whereby the surgery is performed through an incision with tools and camera as opposed to conventional open surgery. This promises minimal recovery times and less hemorrhaging. Multi view LS is the latest development in the field, where the system uses multiple cameras to give the surgeon more information about the surgical site, potentially making the surgery easier. In this publication, we study the gaze patterns of a high performing subject in a multi-view LS environment and compare it with that of a novice to detect the differences between the gaze behavior. This was done by conducting a user study with 20 university students with varying levels of expertise in Multi-view LS. The subjects performed an laparoscopic task in simulation with three cameras (front/top/side). The subjects were then separated as high and low performers depending on the performance times and their data was analyzed. Our results show statistically significant differences between the two behaviors. This opens up new areas from of training novices to Multi-view LS to making smart displays that guide your shows the optimum view depending on the situation.


title: Integrated Nanophotonics Architecture for Residue Number System Arithmetic

abstract: Residue number system (RNS) enables dimensionality reduction of an arithmetic problem by representing a large number as a set of smaller integers, where the number is decomposed by prime number factorization using the moduli as basic functions. These reduced problem sets can then be processed independently and in parallel, thus improving computational efficiency and speed. Here we show an optical RNS hardware representation based on integrated nanophotonics. The digit-wise shifting in RNS arithmetic is expressed as spatial routing of an optical signal in 2x2 hybrid photonic-plasmonic switches. Here the residue is represented by spatially shifting the input waveguides relative to the routers outputs, where the moduli are represented by the number of waveguides. By cascading the photonic 2x2 switches, we design a photonic RNS adder and a multiplier forming an all-to-all sparse directional network. The advantage of this photonic arithmetic processor is the short (10's ps) computational execution time given by the optical propagation delay through the integrated nanophotonic router. Furthermore, we show how photonic processing in-the-network leverages the natural parallelism of optics such as wavelength-division-multiplexing or optical angular momentum in this RNS processor. A key application for photonic RNS is the functional analysis convolution with widespread usage in numerical linear algebra, computer vision, language- image- and signal processing, and neural networks.


title: Fluctuation theory for level-dependent L\'evy risk processes

abstract: A level-dependent L\'evy process solves the stochastic differential equation $dU(t) = dX(t)-{\phi}(U(t)) dt$, where $X$ is a spectrally negative L\'evy process. A special case is a multi-refracted L\'evy process with $\phi_k(x)=\sum_{j=1}^k\delta_j1_{\{x\geq b_j\}}$. A general rate function $\phi$ that is non-decreasing and continuously differentiable is also considered. We discuss solutions of the above stochastic differential equation and investigate the so-called scale functions, which are counterparts of the scale functions from the theory of L\'evy processes. We show how fluctuation identities for $U$ can be expressed via these scale functions. We demonstrate that the derivatives of the scale functions are solutions of Volterra integral equations.


title: The pion holographic wavefunction with dynamical spin effects

abstract: We report that the inclusion of dynamical spin effects in the pion holographic light-front wavefunction leads to a remarkable improvement in describing pion observables (pion mean charge radius, decay constant, spacelike electromagnetic form factor) without the need to invoke higher Fock state contributions and/or a special AdS/QCD mass scale for the pion.


title: Redshift-space distortions from vector perturbations

abstract: We compute a general expression for the contribution of vector perturbations to the redshift-space distortion of galaxy surveys. We show that they contribute to the same multipoles of the correlation function as scalar perturbations and should thus in principle be taken into account in data analysis. We derive constraints for next-generation surveys on the amplitude of two sources of vector perturbations, namely non-linear clustering and topological defects. While topological defects leave a very small imprint on redshift-space distortions, we show that the multipoles of the correlation function are sensitive to vorticity induced by non-linear clustering. Therefore future redshift surveys such as DESI or the SKA should be capable of measuring such vector modes, especially with the hexadecapole which appears to be the most sensitive to the presence of vorticity.


title: Irreversibility in two-body system. Toy model

abstract: It is shown that the irreversible behavior of classical systems of interacting particles is a common property for both few-body and many-body systems. It is due to the delay in the interactions between the particles. PACS: 45.50.Jf Few- and many-body systems; 82.40.Bj Oscillations, chaos, and bifurcations; 05.70.-a Thermodynamics. Keywords: irreversibility, interactions retardation.


title: Multi-scale methods in quantum field theory

abstract: Daubechies wavelets are used to make an exact multi-scale decomposition of quantum fields. For reactions that involve a finite energy that take place in a finite volume, the number of relevant quantum mechanical degrees of freedom is finite. The wavelet decomposition has natural resolution and volume truncations that can be used to isolate the relevant degrees of freedom. The application of flow equation methods to construct effective theories that decouple coarse and fine scale degrees of freedom is examined.


title: Topological Hopf and chain link semimetal states and their application to Co2MnGa (Theory and Materials Prediction)

abstract: Topological semimetals can be classified by the connectivity and dimensionality of the band cross- ing in momentum space. The band crossings of a Dirac, Weyl, or an unconventional fermion semimet- al are 0D points, whereas the band crossings of a nodal-line semimetal are 1D closed loops. Here we propose that the presence of perpendicular crystalline mirror planes can protect 3D band crossings characterized by nontrivial links such as a Hopf link or a coupled-chain, giving rise to a variety of new types of topological semimetals. We show that the nontrivial winding number protects topolog- ical surface states distinct from those in previously known topological semimetals with a vanishing spin-orbit interaction. We also show that these nontrivial links can be engineered by tuning the mirror eigenvalues associated with the perpendicular mirror planes. Using first-principles band structure calculations, we predict the ferromagnetic full Heusler compound Co2MnGa as a candidate. Both Hopf link and chain-like bulk band crossings and unconventional topological surface states are identified.


title: Models of the Mass-Ejection Histories of pre Planetary Nebulae. II. The Formation of the Butterfly and its Proboscis in M2-9

abstract: M2-9, or the "Butterfly Nebula" is one of the most iconic outflow sources from an evolved star. In this paper we present a hydrodynamic model of M2-9 in which the nebula is formed and shaped by a steady, low-density ("light"), mildly collimated "spray" of gas injected at 200 km s^-1 that interacts with a far denser, intrinsically simple pre-existing AGB wind has slowly formed all of the complex features within M2-9's lobes (including the knot pairs N3/S3 and N4/S4 at their respective leading edges, and the radial gradient of Doppler shifts within 20" of the nucleus). We emphasize that the knot pairs are not ejected from the star but formed in situ. In addition, the observed radial speed of the knots is only indirectly related to the speed of the gas injected by the star. The model allows us to probe the early history of the wind geometry and lobe formation. We also formulate a new estimate of the nebular distance D = 1.3 kpc. The physical mechanism that accounts for the linear radial speed gradient in M2-9 applies generally to many other pre planetary nebulae whose hollow lobes exhibit similar gradients along their edges.


title: Madness in vector spaces

abstract: We consider maximal almost disjoint families of block subspaces of countable vector spaces, focusing on questions of their size and definability. We prove that the minimum infinite cardinality of such a family cannot be decided in ZFC and that the "spectrum" of cardinalities of mad families of subspaces can be made arbitrarily large, in analogy to results for mad families on $\omega$ . We apply the author's local Ramsey theory for vector spaces to give partial results concerning their definability.


title: The dynamics of stellar disks in live dark-matter halos

abstract: Recent developments in computer hardware and software enable researchers to simulate the self-gravitating evolution of galaxies at a resolution comparable to the actual number of stars. Here we present the results of a series of such simulations. We performed $N$-body simulations of disk galaxies with between 100 and 500 million particles over a wide range of initial conditions. Our calculations include a live bulge, disk, and dark matter halo, each of which is represented by self-gravitating particles in the $N$-body code. The simulations are performed using the gravitational $N$-body tree-code Bonsai running on the Piz Daint supercomputer. We find that the time scale over which the bar forms increases exponentially with decreasing disk-mass fraction and that the bar formation epoch exceeds a Hubble time when the disk-mass fraction is $\sim0.35$. These results can be explained with the swing-amplification theory. The condition for the formation of $m=2$ spirals is consistent with that for the formation of the bar, which is also an $m=2$ phenomenon. We further argue that the non-barred grand-design spiral galaxies are transitional, and that they evolve to barred galaxies on a dynamical timescale. We also confirm that the disk-mass fraction and shear rate are important parameters for the morphology of disk galaxies. The former affects the number of spiral arms and the bar formation epoch, and the latter determines the pitch angle of the spiral arms.


title: Efficient Third Harmonic Generation in Composite Aluminum Nitride / Silicon Nitride Microrings

abstract: Aluminum nitride and silicon nitride have recently emerged as important nonlinear optical materials in integrated photonics respectively for their quadratic and cubic optical nonlinearity. A composite aluminum nitride and silicon nitride waveguide structure, if realized, will simultaneously allow highly efficient second and third harmonic generation on the same chip platform and therefore assists 2f-3f self-referenced frequency combs. On-chip third harmonic generation, being a higher order nonlinear optics effect, is more demanding than second harmonic generation due to the large frequency difference between the fundamental and third harmonic frequencies which implies large change of refractive indices and more strigent requirements on phase matching. In this work we demonstrate high-efficiency third harmonic generation in a high Q composite aluminum nitride / silicon nitide ring cavity. By carefully engineering the microring resonator geometry of the bilayer structure to optimize the quality factor, mode volume and modal overlap of the optical fields, we report a maximum conversion efficiency of $180\%\,\mathrm{W^{-2}}$. This composite photonic chip design provides a solution for efficient frequency conversion over a large wavelength span, broadband comb generation and self-referenced frequency combs.


title: Precession and glitches in the framework of three-component model of neutron star

abstract: We consider the pulsar rotation assuming that the neutron star consists of crust component (which rotation is observed) and two core components. One of the core components contains pinned superfluid which can, for some reasons, suddenly inject small fraction of stored angular momentum in it. In the framework of this simple model the star can demonstrate glitch-like events together with long period precession (with period $1-10^{4}$ years).


title: Agile Software Engineering and Systems Engineering at SKA Scale

abstract: Systems Engineering (SE) is the set of processes and documentation required for successfully realising large-scale engineering projects, but the classical approach is not a good fit for software-intensive projects, especially when the needs of the different stakeholders are not fully known from the beginning, and requirement priorities might change. The SKA is the ultimate software-enabled telescope, with enormous amounts of computing hardware and software required to perform its data reduction. We give an overview of the system and software engineering processes in the SKA1 development, and the tension between classical and agile SE.


title: Adaptive fast gradient method in stochastic optimization tasks

abstract: In this paper, we describe a stochastic adaptive fast gradient descent method based on the mirror variant of similar triangles method. To our knowledge, this is the first attempt to use adaptivity in stochastic method. Additionally, a main result was proved in terms of probabilities of large deviations.


title: Probabilities of causation of climate changes

abstract: Multiple changes in Earth's climate system have been observed over the past decades. Determining how likely each of these changes are to have been caused by human influence, is important for decision making on mitigation and adaptation policy. Here we describe an approach for deriving the probability that anthropogenic forcings have caused a given observed change. The proposed approach is anchored into causal counterfactual theory (Pearl 2009) which has been introduced recently, and was in fact partly used already, in the context of extreme weather event attribution (EA). We argue that these concepts are also relevant, and can be straightforwardly extended to, the context of detection and attribution of long term trends associated to climate change (D&A). For this purpose, and in agreement with the principle of "fingerprinting" applied in the conventional D&A framework, a trajectory of change is converted into an event occurrence defined by maximizing the causal evidence associated to the forcing under scrutiny. Other key assumptions used in the conventional D&A framework, in particular those related to numerical models error, can also be adapted conveniently to this approach. Our proposal thus allows to bridge the conventional framework with the standard causal theory, in an attempt to improve the quantification of causal probabilities. An illustration suggests that our approach is prone to yield a significantly higher estimate of the probability that anthropogenic forcings have caused the observed temperature change, thus supporting more assertive causal claims.


title: A Short-term Intervention for Long-term Fairness in the Labor Market

abstract: The persistence of racial inequality in the U.S. labor market against a general backdrop of formal equality of opportunity is a troubling phenomenon that has significant ramifications on the design of hiring policies. In this paper, we show that current group disparate outcomes may be immovable even when hiring decisions are bound by an input-output notion of "individual fairness." Instead, we construct a dynamic reputational model of the labor market that illustrates the reinforcing nature of asymmetric outcomes resulting from groups' divergent accesses to resources and as a result, investment choices. To address these disparities, we adopt a dual labor market composed of a Temporary Labor Market (TLM), in which firms' hiring strategies are constrained to ensure statistical parity of workers granted entry into the pipeline, and a Permanent Labor Market (PLM), in which firms hire top performers as desired. Individual worker reputations produce externalities for their group; the corresponding feedback loop raises the collective reputation of the initially disadvantaged group via a TLM fairness intervention that need not be permanent. We show that such a restriction on hiring practices induces an equilibrium that, under particular market conditions, Pareto-dominates those arising from strategies that statistically discriminate or employ a "group-blind" criterion. The enduring nature of equilibria that are both inequitable and Pareto suboptimal suggests that fairness interventions beyond procedural checks of hiring decisions will be of critical importance in a world where machines play a greater role in the employment process.


title: Multiscale modeling of a rectifying bipolar nanopore: explicit-water versus implicit-water simulations

abstract: In a multiscale modeling approach, we present computer simulation results for a rectifying bipolar nanopore on two modeling levels. In an all-atom model, we use explicit water to simulate ion transport directly with the molecular dynamics technique. In a reduced model, we use implicit water and apply the Local Equilibrium Monte Carlo method together with the Nernst-Planck transport equation. This hybrid method makes the fast calculation of ion transport possible at the price of lost details. We show that the implicit-water model is an appropriate representation of the explicit-water model when we look at the system at the device (i.e., input vs.\ output) level. The two models produce qualitatively similar behavior of the electrical current for different voltages and model parameters. Looking at details of concentration and potential profiles, we find profound differences between the two models. These differences, however, do not influence the basic behavior of the model as a device because they do not influence the $z$-dependence of the concentration profiles which are the main determinants of current. These results then address an old paradox: how do reduced models, whose assumptions should break down in a nanoscale device, predict experimental data? Our simulations show that reduced models can still capture the overall device physics correctly, even though they get some important aspects of the molecular-scale physics quite wrong; reduced models work because they include the physics that is necessary from the point of view of device function. Therefore, reduced models can suffice for general device understanding and device design, but more detailed models might be needed for molecular level understanding.


title: Spaces of orders of some one-relator groups

abstract: We show that certain orderable groups admit no isolated left orders. The groups we consider are cyclic amalgamations of a free group with a general orderable group, the HNN extensions of free groups over cyclic subgroups, and a particular class of one-relator groups. In order to prove the results about orders, we develop perturbation techniques for actions of these groups on the line.


title: Inference of Dynamic Regimes in the Microbiome

abstract: Many studies have been performed to characterize the dynamics and stability of the microbiome across a range of environmental contexts [Costello et al., 2012, Faust et al., 2015]. For example, it is often of interest to identify time intervals within which certain subsets of taxa have an interesting pattern of behavior. Viewed abstractly, these problems often have a flavor not just of time series modeling but also of regime detection, a problem with a rich history across a variety of applications, including speech recognition [Fox et al., 2011], finance [Lee, 2009], EEG analysis [Camilleri et al., 2014], and geophysics [Weatherley and Mora, 2002]. In spite of the parallels, regime detection methods are rarely used in microbiome analysis, most likely due to the fact that references for these methods are scattered across several literatures, descriptions are inaccessible outside limited research communities, and implementations are difficult to come across. We distill the core ideas of different regime detection methods, provide example applications, and share reproducible code, making these techniques more accessible to microbiome researchers. We re-analyze data of Dethlefsen and Relman [2011], a study of the effects of antibiotics on the microbiome, using Classification and Regression Trees (CART) [Breiman et al., 1984], Hidden Markov Models (HMMs) [Rabiner and Juang, 1986], Bayesian nonparametric HMMs [Teh and Jordan, 2010, Fox et al., 2008], mixtures of Gaussian Processes (GPs) [Rasmussen and Ghahramani, 2002], switching dynamical systems [Linderman et al., 2016], and multiple changepoint detection [Fan and Mackey, 2015]. Along the way, we summarize each method, their relevance to the microbiome, and tradeoffs associated with using them. Ultimately, our goal is to describe types of temporal or regime switching structure that can be incorporated into studies of microbiome dynamics.


title: Measurements and calculations of air activation in the NuMI neutrino production facility at Fermilab with the 120-GeV proton beam on target

abstract: Measurements and calculations of the air activation at a high-energy proton accelerator are described. The quantity of radionuclides released outdoors depends on operation scenarios including details of the air exchange inside the facility. To improve the prediction of the air activation levels, the MARS15 Monte Carlo code radionuclide production model was modified to be used for these studies. Measurements were done to benchmark the new model and verify its use in optimization studies for the new DUNE experiment at the Long Baseline Neutrino Facility (LBNF) at Fermilab. The measured production rates for the most important radionuclides - $^{11}$C, $^{13}$N, $^{15}$O and $^{41}$Ar - are in a good agreement with those calculated with the improved MARS15 code.


title: On the importance of normative data in speech-based assessment

abstract: Data sets for identifying Alzheimer's disease (AD) are often relatively sparse, which limits their ability to train generalizable models. Here, we augment such a data set, DementiaBank, with each of two normative data sets, the Wisconsin Longitudinal Study and Talk2Me, each of which employs a speech-based picture-description assessment. Through minority class oversampling with ADASYN, we outperform state-of-the-art results in binary classification of people with and without AD in DementiaBank. This work highlights the effectiveness of combining sparse and difficult-to-acquire patient data with relatively large and easily accessible normative datasets.


title: When and how much the altruism impacts your privileged information? Proposing a new paradigm in game theory: The boxers game

abstract: In this work, we proposed a new $N$-person game in which the players can bet on two options, for example represented by two boxers. Some of the players have privileged information about the boxers and part of them can provide this information to uninformed players. However, this information may be true if the informed player is altruist or false if he is selfish. So, in this game, the players are divided in three categories: informed and altruist players, informed and selfish players, and uninformed players. By considering the matchings ($N/2$ distinct pairs of randomly chosen players) and that the payoff of the winning group follows aspects captured from two important games, the public goods game and minority game, we showed quantitatively and qualitatively how the altruism can impact on the privileged information. We localized analytically the regions of positive payoffs which were corroborated by numerical simulations performed for all values of information and altruism densities given that we know the information level of the informed players. Finally, in an evolutionary version of the game ,we showed that the gain of the informed players can get worse if we adopted the following procedure: the players increase their investment for situations of positive payoffs, and decrease their investment when negative payoffs occur.


title: Factorization statistics and the twisted Grothendieck-Lefschetz formula

abstract: We announce recent results on a connection between factorization statistics of polynomials over a finite field and the structure of the cohomology of configurations in $\mathbb{R}^3$ as a representation of the symmetric group. This connection parallels a result of Church, Ellenberg, and Farb relating factorization statistics of squarefree polynomials and the cohomology of configurations in $\mathbb{R}^2$.


title: A mathematical approach to aspects of LET, micro- and nano-dosimetry in radiation therapy with photons and charged particles

abstract: The radiation effects induced by Co60 serve as a reference system for the consideration of LET and RBE in normal and tumor tissue dose-effect relations are usually handled by the linear-quadratic model (LQ). This approximation excellently works up to the shoulder domain. In particle therapy we have strictly differ between RBE in the initial plateau and environment of the Bragg peak. Thus for protons LET and RBE of the initial plateau agree with Co60, whereas in the Bragg peak domain both properties are increased,, but RBE of SOBP only varies between 1.1 and 1.17. The RBE of carbon ions is increased once again Their dose-effect curves are much steeper with a rather small shoulder domain due to dense ionizing radiation effect. Thus protons are also dense ionizing in the Bragg peak region, but with rather smaller magnitude compared to carbon ions. A generalization of the LQ-model based on the nonlinear reaction-diffusion model is proposed to describe LET and RBE of dense ionizing particles, which accounts for properties of micro- and nano-dosimetry. A linear term of a reaction diffusion formula describes the destroy of cells, the nonlinear term is related to repair and the diffusion term accounts for the density of the radiation damages. Based on dose-effect properties of Co60 the parameters of dense ionizing particles can be determined and compared with measurement data. The local dense of radiation effects and their consequences in RBE and dose effect curves provide a key of understanding modern therapy planning with different modalities and properties of nano-dosimetry are interpreted by mathematical descriptions. The irradiation of spheroids is a feature of micro-dosimetry, whereas intracellular exposure refers to nano-dosimetry.


title: Johnson-Levine Homomorphisms and the tree reduction of the LMO functor

abstract: Let M denote the mapping class group of S, a compact connected oriented surface with one boundary component. The action of M on the nilpotent quotients of the fundamental group of S allows to define the so-called Johnson filtration and the Johnson homomorphisms. J. Levine introduced a new filtration of M, called the Lagrangian filtration. He also introduced a version of the Johnson homomorphisms for this new filtration. The first term of the Lagrangian filtration is the Lagrangian mapping class group, whose definition involves a handlebody bounded by S, and which contains the Torelli group. These constructions extend in a natural way to the monoid of homology cobordisms. Besides, D. Cheptea, K. Habiro and G. Massuyeau constructed a functorial extension of the LMO invariant, called the LMO functor, which takes values in a category of diagrams. In this paper we give a topological interpretation of the upper part of the tree reduction of the LMO functor in terms of the homomorphisms defined by J. Levine for the Lagrangian mapping class group. We also compare the Johnson filtration with the filtration introduced by J. Levine.


title: Recent advances on chiral partners and patterns

abstract: We review recent results on chiral $SU(2)_L\times SU(2)_R\approx O(4)$ and $U(1)_A$ symmetry restoration in QCD. In particular, we discuss how Ward Identities allow one to derive general results on partner degeneration, which shed light on the distinction between the $O(4)$ and $O(4)\times U(1)_A$ patterns of the chiral transition. For that purpose, susceptibilities associated with the $O(4)$ and $U(1)_A$ symmetries are studied. From this analysis we conclude that in the ideal regime of exact $O(4)$ restoration (formally achieved in the limit of two massless flavours), $U(1)_A$ partners degenerate as well. We also discuss the role of the thermal $f_0(500)$ state to describe thermodynamic observables sensitive to chiral restoration, such as the scalar susceptibility. We pay special attention to the consistency of our results with recent lattice analysis.


title: Multi-Channel CNN-based Object Detection for Enhanced Situation Awareness

abstract: Object Detection is critical for automatic military operations. However, the performance of current object detection algorithms is deficient in terms of the requirements in military scenarios. This is mainly because the object presence is hard to detect due to the indistinguishable appearance and dramatic changes of object's size which is determined by the distance to the detection sensors. Recent advances in deep learning have achieved promising results in many challenging tasks. The state-of-the-art in object detection is represented by convolutional neural networks (CNNs), such as the fast R-CNN algorithm. These CNN-based methods improve the detection performance significantly on several public generic object detection datasets. However, their performance on detecting small objects or undistinguishable objects in visible spectrum images is still insufficient. In this study, we propose a novel detection algorithm for military objects by fusing multi-channel CNNs. We combine spatial, temporal and thermal information by generating a three-channel image, and they will be fused as CNN feature maps in an unsupervised manner. The backbone of our object detection framework is from the fast R-CNN algorithm, and we utilize cross-domain transfer learning technique to fine-tune the CNN model on generated multi-channel images. In the experiments, we validated the proposed method with the images from SENSIAC (Military Sensing Information Analysis Centre) database and compared it with the state-of-the-art. The experimental results demonstrated the effectiveness of the proposed method on both accuracy and computational efficiency.


title: Machine Learning and Manycore Systems Design: A Serendipitous Symbiosis

abstract: Tight collaboration between experts of machine learning and manycore system design is necessary to create a data-driven manycore design framework that integrates both learning and expert knowledge. Such a framework will be necessary to address the rising complexity of designing large-scale manycore systems and machine learning techniques.


title: Fluctuation identities with continuous monitoring and their application to price barrier options

abstract: We present a numerical scheme to calculate fluctuation identities for exponential L\'evy processes in the continuous monitoring case. This includes the Spitzer identities for touching a single upper or lower barrier, and the more difficult case of the two-barriers exit problem. These identities are given in the Fourier-Laplace domain and require numerical inverse transforms. Thus we cover a gap in the literature that has mainly studied the discrete monitoring case; indeed, there are no existing numerical methods that deal with the continuous case. As a motivating application we price continuously monitored barrier options with the underlying asset modelled by an exponential L\'evy process. We perform a detailed error analysis of the method and develop error bounds to show how the performance is limited by the truncation error of the sinc-based fast Hilbert transform used for the Wiener-Hopf factorisation. By comparing the results for our new technique with those for the discretely monitored case (which is in the Fourier-$z$ domain) as the monitoring time step approaches zero, we show that the error convergence with continuous monitoring represents a limit for the discretely monitored scheme.


title: Mining Precision Interfaces From Query Logs

abstract: Interactive tools make data analysis both more efficient and more accessible to a broad population. Simple interfaces such as Google Finance as well as complex visual exploration interfaces such as Tableau are effective because they are tailored to the desired user tasks. Yet, designing interactive interfaces requires technical expertise and domain knowledge. Experts are scarce and expensive, and therefore it is currently infeasible to provide tailored (or precise) interfaces for every user and every task. We envision a data-driven approach to generate tailored interactive interfaces. We observe that interactive interfaces are designed to express sets of programs; thus, samples of programs-increasingly collected by data systems-may help us build interactive interfaces. Based on this idea, Precision Interfaces is a language-agnostic system that examines an input query log, identifies how the queries structurally change, and generates interactive web interfaces to express these changes. The focus of this paper is on applying this idea towards logs of structured queries. Our experiments show that Precision Interfaces can support multiple query languages (SQL and SPARQL), derive Tableau's salient interaction components from OLAP queries, analyze <75k queries in <12 minutes, and generate interaction designs that improve upon existing interfaces and are comparable to human-crafted interfaces.


title: Dynamical density functional theory for the evaporation of droplets of nanoparticle suspension

abstract: We develop a lattice gas model for the drying of droplets of a nanoparticle suspension on a planar surface, using dynamical density functional theory (DDFT) to describe the time evolution of the solvent and nanoparticle density profiles. The DDFT assumes a diffusive dynamics but does not include the advective hydrodynamics of the solvent, so the model is relevant to highly viscous or near to equilibrium systems. Nonetheless, we see an equivalent of the coffee-ring stain effect, but in the present model it occurs for thermodynamic rather the fluid-mechanical reasons. The model incorporates the effect of phase separation and vertical density variations within the droplet and the consequence of these on the nanoparticle deposition pattern on the surface. We show how to include the effect of slip or no-slip at the surface and how this is related to the receding contact angle. We also determine how the equilibrium contact angle depends on the microscopic interaction parameters.


title: Super SloMo: High Quality Estimation of Multiple Intermediate Frames for Video Interpolation

abstract: Given two consecutive frames, video interpolation aims at generating intermediate frame(s) to form both spatially and temporally coherent video sequences. While most existing methods focus on single-frame interpolation, we propose an end-to-end convolutional neural network for variable-length multi-frame video interpolation, where the motion interpretation and occlusion reasoning are jointly modeled. We start by computing bi-directional optical flow between the input images using a U-Net architecture. These flows are then linearly combined at each time step to approximate the intermediate bi-directional optical flows. These approximate flows, however, only work well in locally smooth regions and produce artifacts around motion boundaries. To address this shortcoming, we employ another U-Net to refine the approximated flow and also predict soft visibility maps. Finally, the two input images are warped and linearly fused to form each intermediate frame. By applying the visibility maps to the warped images before fusion, we exclude the contribution of occluded pixels to the interpolated intermediate frame to avoid artifacts. Since none of our learned network parameters are time-dependent, our approach is able to produce as many intermediate frames as needed. We use 1,132 video clips with 240-fps, containing 300K individual video frames, to train our network. Experimental results on several datasets, predicting different numbers of interpolated frames, demonstrate that our approach performs consistently better than existing methods.


title: Quantizing time: Interacting clocks and systems

abstract: This article generalizes the conditional probability interpretation of time in which time evolution is realized through entanglement between a clock and a system of interest. This formalism is based upon conditioning a solution to the Wheeler-DeWitt equation on a subsystem of the Universe, serving as a clock, being in a state corresponding to a time $t$. Doing so assigns a conditional state to the rest of the Universe $|\psi_S(t)\rangle$, referred to as the system. We demonstrate that when the total Hamiltonian appearing in the Wheeler-DeWitt equation contains an interaction term coupling the clock and system, the conditional state $|\psi_S(t)\rangle$ satisfies a time-nonlocal Schr\"{o}dinger equation in which the system Hamiltonian is replaced with a self-adjoint integral operator. This time-nonlocal Schr\"{o}dinger equation is solved perturbatively and three examples of clock-system interactions are examined. One example considered supposes that the clock and system interact via Newtonian gravity, which leads to the system's Hamiltonian developing corrections on the order of $G/c^4$ and inversely proportional to the distance between the clock and system.


title: A Quantized Inter-level Character in Quantum Systems

abstract: For a quantum system subject to external parameters, the Berry phase is an intra-level property, which is gauge invariant module $2\pi$ for a closed loop in the parameter space and generally is non-quantized. In contrast, we define a inter-band character $\Theta$ for a closed loop, which is gauge invariant and quantized as integer values. It is a quantum mechanical analogy of the Euler character based on the Gauss-Bonnet theorem for a manifold with a boundary. The role of the Gaussian curvature is mimicked by the difference between the Berry curvatures of the two levels, and the counterpart of the geodesic curvature is the quantum geometric potential which was proposed to improve the quantum adiabatic condition. This quantized inter-band character is also generalized to quantum degenerate systems.


title: Free Energy of a Large-$N$ Pion Gas and Chiral Symmetry Restoration

abstract: We study thermal properties of a large-N massless pion gas using a low-energy QCD approach given by an $O(N +1)/O(N)$ Nonlinear Sigma Model. We build diagrammatically the associated finite free energy to $O(TM^{3})$ in the pion mass expansion through an effective vertex that considers all the contributions coming from closed diagrams. Subsequently, we calculate finite order parameters such as the quark condensate and its respective derivative, i.e., the scalar susceptibility, in the chiral limit, along with their associated critical exponents. These results are compared with our previous unitarized scattering analyses for the chiral transition universality class, thus showing a reasonable agreement both with lattice simulations and these resonant studies.


title: Initial Demonstration Of 9-MHz Framing Camera Rates On The FAST Drive Laser Pulse Trains

abstract: We report the configuration of a Hamamatsu C5680 streak camera as a framing camera to record transverse spatial information of green-component laser micropulses at 3- and 9-MHz rates for the first time. The latter is near the time scale of the ~7.5-MHz revolution frequency of the Integrable Optics Test Accelerator (IOTA) ring and its expected synchroton radiation source temporal structure. The 2-D images are recorded with a Gig-E readout CCD camera. We also report a first proof of principle with an OTR source using the linac streak camera in a semi-framing mode.


title: Decay angular distributions of $K^*$ and $D^*$ mesons as a tool for the dynamics of open strange and charm production

abstract: We analyze decay spin-density matrix elements of $K^*$ and $D^*$ vector mesons and find that the corresponding vector meson decay distributions are sensitive to the production mechanisms. This suggests that the measurement of these quantities can be used to unravel the production mechanisms.


title: AdS/QCD Modified Soft Wall Model and Light Meson Spectra

abstract: We analyze here the mass spectrum of light vector and scalar mesons applying a novel approach where a modified soft wall model that includes a UV-cutoff at a finite $z$-position in the AdS space is used, thus introducing an extra energy scale. For this model, we found that the masses for the scalar and vector spectra are well fitted within $\delta_{RMS}=7.64\%$ for these states, with non-linear trajectories given by two common parameters, the UV locus $z_{0}$ and the quadratic dilaton profile slope $\kappa$. We also conclude that in this model, the $f_{0}(500)$ scalar resonance cannot be fitted holographycally as a $q\overline{q}$ state since we cannot find a trajectory that include this pole. This result is in agreement with the most recent phenomenological and theoretical methods.


title: Gravitational waves from phase transition in split NMSSM

abstract: We discuss gravitational wave signal from the strongly first order electroweak phase transition in the split NMSSM. We find that for sets of parameters predicting successful electroweak baryogenesis the gravitational wave signal can be within the reach of future experiments LISA, BBO and Ultimate DECIGO.


title: The Effect of Focal Distance, Age, and Brightness on Near-Field Augmented Reality Depth Matching

abstract: Many augmented reality (AR) applications operate within near-field reaching distances, and require matching the depth of a virtual object with a real object. The accuracy of this matching was measured in three experiments, which examined the effect of focal distance, age, and brightness, within distances of 33.3 to 50 cm, using a custom-built AR haploscope. Experiment I examined the effect of focal demand, at the levels of collimated (infinite focal distance), consistent with other depth cues, and at the midpoint of reaching distance. Observers were too young to exhibit age-related reductions in accommodative ability. The depth matches of collimated targets were increasingly overestimated with increasing distance, consistent targets were slightly underestimated, and midpoint targets were accurately estimated. Experiment II replicated Experiment I, with older observers. Results were similar to Experiment I. Experiment III replicated Experiment I with dimmer targets, using young observers. Results were again consistent with Experiment I, except that both consistent and midpoint targets were accurately estimated. In all cases, collimated results were explained by a model, where the collimation biases the eyes' vergence angle outwards by a constant amount. Focal demand and brightness affect near-field AR depth matching, while age-related reductions in accommodative ability have no effect.


title: One spectrum to rule them all?

abstract: We show that in absence of entropy or effective anisotropic stress the freedom in the choice of the initial energy scale of inflation implies the existence of an infinite family of dual slow-roll parameters histories which can produce the same spectrum of comoving curvature perturbations. This implies that in general there is no one-to-one correspondence between the spectrum and higher order correlation functions. We give some numerical examples of expansion histories corresponding to different initial energy scales, with the same spectrum of curvature perturbations, the same squeezed limit bispectrum, in agreement with the squeezed limit consistency condition, but with different bispectra in other configurations and different spectra of primordial gravitational waves. The combined analysis of data from future CMB and gravitational wave experiments could allow to distinguish between dual models.


title: A remark on the two dimensional water wave problem with surface tension

abstract: We consider the motion of a two-dimensional interface between air (above) and an irrotational, incompressible, inviscid, infinitely deep water (below), with surface tension present. We propose a new way to reduce the original problem into an equivalent quasilinear system which are related to the interface's tangent angle and a quantity related to the difference of tangential velocities of the interface in the Lagrangian and the arc-length coordinates. The new way is relatively simple because it involves only taking differentiation and the real and the imaginary parts. Then if assuming that waves are periodic, we establish a priori energy inequality.


title: Radiative effects in deep virtual Compton scattering

abstract: Radiative corrections to the cross sections of photon electroproduction and the single spin asymmetries induced by the interference between the Bethe Heitler and deep virtual Compton scattering amplitudes are calculated within the leading log approximation. The deep virtual Compton scattering amplitude is presented in the Belitsky, M\"uller, and Kirchner approximation for the polarized initial particles. The Fortran code for estimation of the radiative effects in a given kinematic point and Monte Carlo generator for simulation of one or two photons are developed. Numerical results are performed for beam-spin asymmetries in kinematical conditions of current experiments in the Jefferson Laboratory.


title: A local asymptotic expansion for a local solution of the Stokes system

abstract: We consider solutions of the Stokes system in a neighborhood of a point in which the velocity $u$ vanishes of order $d$. We prove that there exists a divergence-free polynomial $P$ in $x$ with $t$-dependent coefficients of degree $d$ which approximates the solution $u$ of order $d+\alpha$ for certain $\alpha>0$. The polynomial $P$ satisfies a Stokes equation with a forcing term which is a sum of two polynomials in $x$ of degrees $d-1$ and $d$. The results extend to Oseen systems and to the Navier-Stokes equation.


title: Residential Energy Storage Management with Bidirectional Energy Control

abstract: We consider the residential energy storage management system with integrated renewable generation, with the availability of bidirectional energy flow from and to the grid thorough buying and selling. We propose a real-time bidirectional energy control algorithm, aiming to minimize the net system cost, due to energy buying and selling and battery deterioration and inefficiency from storage activities, within a given time period, subject to the battery operational constraints and energy buying and selling constraints. We formulate the problem as a stochastic control optimization problem. We then modify and transform this difficult problem into one that enables us to develop the real-time energy control algorithm through Lyapunov optimization. Our developed algorithm is applicable to arbitrary and unknown statistics of renewable generation, load, and electricity prices. It provides a simple closed-form control solution only based on current system states with minimum complexity for real-time implementation. Furthermore, the solution structure reveals how the battery energy level and energy prices affect the decision on energy flow and storage. The proposed algorithm possesses a bounded performance guarantee to that of the optimal non-causal T-slot look-ahead control policy. Simulation shows the effectiveness of our proposed algorithm as compared with alternative real-time and non-causal algorithms, as well as the effect of selling-to-buying price ratio and battery inefficiency on the storage behavior and system cost.


title: Strong orientation dependence of surface mass density profiles of dark haloes at large scales

abstract: We study the dependence of surface mass density profiles, which can be directly measured by weak gravitational lensing, on the orientation of haloes with respect to the line-of-sight direction, using a suite of $N$-body simulations. We find that, when major axes of haloes are aligned with the line-of-sight direction, surface mass density profiles have higher amplitudes than those averaged over all halo orientations, over all scales from $0.1$ to $100\,\mathrm{Mpc}/h$ we studied. While the orientation dependence at small scales is ascribed to the halo triaxiality, our results indicate even stronger orientation dependence in the so-called two-halo regime, up to $100\,\mathrm{Mpc}/h$. The orientation dependence for the two-halo term is well approximated by a multiplicative shift of the amplitude and therefore a shift in the halo bias parameter value. The halo bias from the two-halo term can be overestimated or underestimated by up to $\sim 30 \%$ depending on the viewing angle, which translates into the bias in estimated halo masses by up to a factor of two from halo bias measurements. The orientation dependence at large scales originates from the anisotropic halo-matter correlation function, which has an elliptical shape with the axis ratio of $\sim 0.55$ up to $100\, \mathrm{Mpc}/h$. We discuss potential impacts of halo orientation bias on other observables such as optically selected cluster samples and a clustering analysis of large-scale structure tracers such as quasars.