title: Exact solutions to the nonlinear Schr\"odinger equation with time-dependent coefficients

abstract: In this paper, the trial function method is employed to find the exact solutions for high-order nonlinear Schr\"odinger equations with time-dependent coefficients. This system describes the propagation of ultrashort light pulses in nonlinear fibers, with self-steepening and self-frequency shift effects. As a result, we derive a range of exact solutions which include Jacobi elliptic function solutions, solitary wave solutions, and rational function solutions.


title: Probabilistic simulation of partial differential equations

abstract: Computer simulations of differential equations require a time discretization, which inhibits to identify the exact solution with certainty. Probabilistic simulations take this into account via uncertainty quantification. The construction of a probabilistic simulation scheme can be regarded as Bayesian filtering by means of probabilistic numerics. Gaussian prior based filters, specifically Gauss-Markov priors, have successfully been applied to simulation of ordinary differential equations (ODEs) and give rise to filtering problems that can be solved efficiently. This work extends this approach to partial differential equations (PDEs) subject to periodic boundary conditions and utilizes continuous Gaussian processes in space and time to arrive at a Bayesian filtering problem structurally similar to the ODE setting. The usage of a process that is Markov in time and statistically homogeneous in space leads to a probabilistic spectral simulation method that allows for an efficient realization. Furthermore, the Bayesian perspective allows the incorporation of methods developed within the context of information field theory such as the estimation of the power spectrum associated with the prior distribution, to be jointly estimated along with the solution of the PDE.

categories: math.NA astro-ph.IM cs.NA physics.comp-ph stat.CO stat.ME


title: Jointly Optimizing Sensing Pipelines for Multimodal Mixed Reality Interaction

abstract: Natural human interactions for Mixed Reality Applications are overwhelmingly multimodal: humans communicate intent and instructions via a combination of visual, aural and gestural cues. However, supporting low-latency and accurate comprehension of such multimodal instructions (MMI), on resource-constrained wearable devices, remains an open challenge, especially as the state-of-the-art comprehension techniques for each individual modality increasingly utilize complex Deep Neural Network models. We demonstrate the possibility of overcoming the core limitation of latency--vs.--accuracy tradeoff by exploiting cross-modal dependencies -- i.e., by compensating for the inferior performance of one model with an increased accuracy of more complex model of a different modality. We present a sensor fusion architecture that performs MMI comprehension in a quasi-synchronous fashion, by fusing visual, speech and gestural input. The architecture is reconfigurable and supports dynamic modification of the complexity of the data processing pipeline for each individual modality in response to contextual changes. Using a representative "classroom" context and a set of four common interaction primitives, we then demonstrate how the choices between low and high complexity models for each individual modality are coupled. In particular, we show that (a) a judicious combination of low and high complexity models across modalities can offer a dramatic 3-fold decrease in comprehension latency together with an increase 10-15% in accuracy, and (b) the right collective choice of models is context dependent, with the performance of some model combinations being significantly more sensitive to changes in scene context or choice of interaction.


title: Non-commutative rational functions in the full Fock space

abstract: A rational function belongs to the Hardy space, $H^2$, of square-summable power series if and only if it is bounded in the complex unit disk. Any such rational function is necessarily analytic in a disk of radius greater than one. The inner-outer factorization of a rational function, $\mathfrak{r} \in H^2$ is particularly simple: The inner factor of $\mathfrak{r}$ is a (finite) Blaschke product and (hence) both the inner and outer factors are again rational. We extend these and other basic facts on rational functions in $H^2$ to the full Fock space over $\mathbb{C}^d$, identified as the \emph{non-commutative (NC) Hardy space} of square-summable power series in several NC variables. In particular, we characterize when an NC rational function belongs to the Fock space, we prove analogues of classical results for inner-outer factorizations of NC rational functions and NC polynomials, and we obtain spectral results for NC rational multipliers.


title: The explicit formula of Hankel determinant with Catalan elements

abstract: Applying Johann Cigler's Hankel determinant formula in terms of the binomial coefficient determinants, which is simplified from Christian Krattenthale's, we get an explicit formula of Hankel determinants for general. As far as I know, those are new results.


title: Measured Nondestructive Assay of $^{237}$Np Using Organic Scintillators and Active Neutron Multiplicity Counting

abstract: The purpose of nondestructive assay in the context of nuclear safeguards is to precisely verify the declared mass of a sample of nuclear material in a noninhibitive amount of time. 237Np is a proliferation concern, and the capacity to efficiently assay samples of it is a missing piece in the verification and safeguards toolbox. The material is subject to the same safeguards as 235U, is reportable in gram quantities, and is classified as "other nuclear material" according to the United States Department of Energy. Given that 3000 kg of 237Np is annually produced in the US and the bare sphere critical mass is 40-60 kg, it is desirable to augment the safeguards toolbox with a system capable of distinguishing 10 g of 237Np in a 20-minute measurement. One measurement modality is neutron multiplicity counting, which relates the detected multiplicity count rates to the amount of fissionable material. Prior simulation work shows that an organic scintillator-based multiplicity counter can achieve the design criteria, whereas the flagship 3He-based system, the Epithermal Neutron Multiplicity Counter, requires much longer measurement times to achieve the same precision. In this work, simultaneous measurements of a 6-kg sphere of 237Np by organic scintillator- and 3He-based systems are used to confirm the trends in the simulation study; the organic scintillator-based system achieves 1% uncertainty in the neutron double multiplicity rate on the order of minutes, while the 3He-based system requires days to reach the same precision. In conclusion, the International Atomic Energy Agency should consider the development and deployment of an organic scintillator-based multiplicity counter


title: Transportation Interventions Reshaping NYC Commute: the Probabilistic Simulation Framework Assessing the Impacts of Ridesharing and Manhattan Congestion Surcharge

abstract: Understanding holistic impact of planned transportation solutions and interventions on urban systems is challenged by their complexity but critical for decision making. The cornerstone for such impact assessments is estimating the transportation mode-shift resulting from the intervention. And while transportation planning has well-established models for the mode-choice assessment such as the nested multinomial logit model, an individual choice simulation could be better suited for addressing the mode-shift allowing to consistently account for individual preferences. In addition, no model perfectly represents the reality while the available ground truth data on the actual transportation choices needed to infer the model is often incomplete or inconsistent. The present paper addresses those challenges by offering an individual mode-choice and mode-shift simulation model and the Bayesian inference framework. It accounts for uncertainties in the data as well as the model estimate and translates them into uncertainties of the resulting mode-shift and the impacts. The framework is evaluated on the two intervention cases: introducing ride-sharing for-hire-vehicles in NYC as well the recent introduction of the Manhattan Congestion Surcharge. Being successfully evaluated on the cases above, the framework can be used for assessing mode-shift and resulting economic, social and environmental implications for any future urban transportation solutions and policies being considered by decision-makers or transportation companies.


title: Impact of Dust on Spectral Distortion Measurements of the Cosmic Microwave Background

abstract: Spectral distortions of the cosmic microwave background (CMB) are sensitive to energy injection by exotic physics in the early universe. The proposed Primordial Inflation Explorer (PIXIE) mission has the raw sensitivity to provide meaningful limits on new physics, but only if foreground emission can be adequately modeled. We quantify the impact of interstellar dust on Compton $y$ and $\mu$ measurements by considering a range of grain size distributions and compositions constrained by theoretical and observational priors (Zelko & Finkbeiner 2020). We find that PIXIE can marginalize over a modest number of dust parameters and still recover $y$ and $\mu$ estimates, though with increased uncertainty. As more foreground components are included (synchrotron, free-free), the estimates of $y$ degrade, and measurement of $\mu$ in the range sometimes considered for the standard $\Lambda$CDM of $2\times10^{-8}$ becomes infeasible without ancillary low-frequency foreground information. An additional concern is dust absorption of the CMB monopole, a subtle effect that must be included. We quantify one form of model discrepancy error, finding that the error introduced by fitting our interstellar medium dust model with a modified blackbody is too large for CMB spectral distortions to be detectable. The greatest challenge may be the cosmic infrared background (CIB). We find that $\mu$ and $y$ are extremely sensitive to modeling choices for the CIB, and quantify biases expected for a range of assumptions.


title: Globular clusters as tracers of the dark matter content of dwarfs in galaxy clusters

abstract: Globular clusters (GCs) are often used to estimate the dark matter content of galaxies, especially dwarfs, where other kinematic tracers are lacking. These estimates typically assume spherical symmetry and dynamical equilibrium, assumptions that may not hold for the sparse GC population of dwarfs in galaxy clusters. We use a catalog of GCs tagged onto the Illustris simulation to study the accuracy of GC-based mass estimates. We focus on galaxies in the stellar mass range 10$^{8} - 10^{11.8}$ M$_{\odot}$ identified in $9$ simulated Virgo-like clusters. Our results indicate that mass estimates are, on average, quite accurate in systems with GC numbers $N_{\rm GC} \geq 10$ and where the uncertainty of individual GC line-of-sight velocities is smaller than the inferred velocity dispersion, $\sigma_{\rm GC}$. In cases where $N_{\rm GC} \leq 10$, however, biases may result depending on how $\sigma_{\rm GC}$ is computed. We provide calibrations that may help alleviate these biases in methods widely used in the literature. As an application, we find a number of dwarfs with $M_{*} \sim 10^{8.5}\, M_{\odot}$ (comparable to the ultradiffuse galaxy DF2, notable for the low $\sigma_{GC}$ of its $10$ GCs) with $\sigma_{\rm GC} \sim 7$ - $15\; \rm km \rm s^{-1}$. These DF2 analogs correspond to relatively massive systems at their infall time ($M_{200} \sim 1$ - $3 \times 10^{11}$ $M_{\odot}$) which have retained only $3$-$17$ GCs and have been stripped of more than 95$\%$ of their dark matter. Our results suggest that extreme tidal mass loss in otherwise normal dwarf galaxies may be a possible formation channel for ultradiffuse objects like DF2.


title: Extended stellar systems in the solar neighborhood -- V. Discovery of coronae of nearby star clusters

abstract: In this paper, we present a novel view on the morphology and the dynamical state of 10 prominent, nearby ($\leq$ 500 pc), and young ($\sim$30-300 Myr) open star clusters with Gaia DR2: $\alpha\,$Per, Blanco 1, IC 2602, IC 2391, Messier 39, NGC 2451A, NGC 2516, NGC 2547, Platais 9, and the Pleiades. We introduce a pioneering member identification method that is informed by cluster bulk velocities and deconvolves the spatial distribution with a mixture of Gaussians. Our approach enables inferring the clusters' true spatial distribution by effectively filtering field star contaminants while at the same time mitigating the impact of positional errors along the line of sight. This first application of the method reveals the existence of vast stellar coronae, extending for $\gtrsim\,$100 pc and surrounding the, by comparison tiny and compact, cluster cores. The coronae and cores form intertwined, co-eval, and co-moving extended cluster populations, each encompassing tens of thousands of cubic parsec and stretching across tens of degrees on the sky. Our analysis shows that the coronae are gravitationally unbound but largely comprise the bulk of the populations' stellar mass. Most systems are in a highly dynamic state, showing evidence of expansion and sometimes simultaneous contraction along different spatial axes. The velocity field of the extended populations for the cluster cores appears asymmetric but is aligned along a spatial axis unique to each cluster. The overall spatial distribution and the kinematic signature of the populations are largely consistent with the differential rotation pattern of the Milky Way. This finding underlines the important role of global Galactic dynamics to the fate of stellar systems. Our results highlight the complexity of the Milky Way's open cluster population and call for a new perspective on the characterization and dynamical state of open clusters.


title: Quantum operator growth bounds for kicked tops and semiclassical spin chains

abstract: We present a framework for understanding the dynamics of operator size, and bounding the growth of out-of-time-ordered correlators, in models of large-$S$ spins. Focusing on the dynamics of a single spin, we show the finiteness of the Lyapunov exponent in the large-$S$ limit; our bounds are tighter than the best known Lieb-Robinson-type bounds on these systems. We numerically find our upper bound on Lyapunov exponents is within an order of magnitude of numerically computed values in classical and quantum kicked top models. Generalizing our results to coupled large-$S$ spins on lattices, we show that the butterfly velocity, which characterizes the spatial speed of quantum information scrambling, is finite as $S\rightarrow\infty$. We emphasize qualitative differences between operator growth in semiclassical large-spin models, and quantum holographic systems including the Sachdev-Ye-Kitaev model.

categories: cond-mat.str-el math-ph math.MP quant-ph


title: Scalar Leptoquarks in Leptonic Processes

abstract: Leptoquarks are hypothetical new particles, which couple quarks directly to leptons. They experienced a renaissance in recent years as they are prime candidates to explain the so-called \textit{flavor anomalies}, i.e. the deviations between the Standard Model predictions and measurements in $b\to s\ell^{+}\ell^{-}$ and $b\to c\tau\nu$ processes and in the anomalous magnetic moment of the muon. At the one-loop level these particles unavoidably generate effects in the purely leptonic processes like $Z\to\ell^{+}\ell^{-}$, $Z\to\nu\bar\nu$, $W\to\ell\nu$ and $h\to\ell^{+}\ell^{-}$ and can even generate non-zero rates for lepton flavor violating processes such as $\ell\to \ell^{\prime}\gamma$, $Z\to\ell^+\ell^{\prime -}$, $h\to\ell^+\ell^{\prime -}$ and $\ell\to 3\ell^\prime$. In this article we calculate these processes for all five representations of scalar Leptoquarks. We include their most general interaction terms with the Standard Model Higgs boson, which leads to Leptoquark mixing after the former acquires a vacuum expectation value. In our phenomenological analysis we investigate the effects in modified lepton couplings to electroweak gauge bosons, we study the correlations of the anomalous magnetic moment of the muon with $h\to\mu^{+}\mu^{-}$ and $Z\to\mu^{+}\mu^{-}$ as well as the interplay between different lepton flavor violating decays.


title: Generalized symmetries and 2-groups via electromagnetic duality in AdS/CFT

abstract: We discuss how electromagnetically dualizing a 1-form to a 2-form in AdS$_5$ exchanges regular and alternate boundary conditions, and thus gauges the originally global $U(1)$ symmetry in the dual field theory. The generalized symmetry current dual to the 2-form in the bulk is identified as the dual field strength of the gauged $U(1)$, and the associated double-trace operator with a logarithmically running coupling is just the gauged $U(1)$ Maxwell action. Applying this dualization to an AdS Maxwell-Chern-Simons theory dual to a global $U(1) \times U(1)$ model with an 't Hooft anomaly results in a theory with a modified field strength that holographically realizes a 2-group symmetry. We explicitly carry out the holographic renormalization to verify this, and discuss the generalization to other rank fields in other dimensions.


title: With Little Power Comes Great Responsibility

abstract: Despite its importance to experimental design, statistical power (the probability that, given a real effect, an experiment will reject the null hypothesis) has largely been ignored by the NLP community. Underpowered experiments make it more difficult to discern the difference between statistical noise and meaningful model improvements, and increase the chances of exaggerated findings. By meta-analyzing a set of existing NLP papers and datasets, we characterize typical power for a variety of settings and conclude that underpowered experiments are common in the NLP literature. In particular, for several tasks in the popular GLUE benchmark, small test sets mean that most attempted comparisons to state of the art models will not be adequately powered. Similarly, based on reasonable assumptions, we find that the most typical experimental design for human rating studies will be underpowered to detect small model differences, of the sort that are frequently studied. For machine translation, we find that typical test sets of 2000 sentences have approximately 75% power to detect differences of 1 BLEU point. To improve the situation going forward, we give an overview of best practices for power analysis in NLP and release a series of notebooks to assist with future power analyses.


title: High-Cadence, Early-Time Observations of Core-Collapse Supernovae From the TESS Prime Mission

abstract: We present observations from the Transiting Exoplanet Survey Satellite (TESS) of twenty bright core-collapse supernovae with peak TESS-band magnitudes $\lesssim18$ mag. We reduce this data with an implementation of the image subtraction pipeline used by the All-Sky Automated Survey for Supernovae (ASAS-SN) optimized for use with the TESS images. In empirical fits to the rising light curves, we do not find strong correlations between the fit parameters and the peak luminosity. Existing semi-analytic models fit the light curves of the Type II supernovae well, but do not yield reasonable estimates of the progenitor radius or explosion energy, likely because they are derived for use with ultraviolet observations while TESS observes in the near-infrared. If we instead fit the data with numerically simulated light curves, the rising light curves of the Type~II SNe are consistent with the explosions of red supergiants. While we do not identify shock breakout emission for any individual event, when we combine the fit residuals of the Type II supernovae in our sample, we do find a $>5\sigma$ flux excess in the $\sim 0.5$~day before the start of the light curve rise. It is likely that this excess is due to shock breakout emission, and that during its extended mission TESS will observe a Type II supernova bright enough for this signal to be detected directly.


title: Flops, Gromov-Witten Invariants and Symmetries of Line Bundle Cohomology on Calabi-Yau Three-folds

abstract: The zeroth line bundle cohomology on Calabi-Yau three-folds encodes information about the existence of flop transitions and the genus zero Gromov-Witten invariants. We illustrate this claim by studying several Picard number 2 Calabi-Yau three-folds realised as complete intersections in products of projective spaces. Many of these manifolds exhibit certain symmetries on the Picard lattice which preserve the zeroth cohomology.


title: Elliptic Loci of SU(3) Vacua

abstract: The space of vacua of many four-dimensional, $\mathcal{N}=2$ supersymmetric gauge theories can famously be identified with a family of complex curves. For gauge group $SU(2)$, this gives a fully explicit description of the low-energy effective theory in terms of an elliptic curve and associated modular fundamental domain. The two-dimensional space of vacua for gauge group $SU(3)$ parametrizes an intricate family of genus two curves. We analyze this family using the so-called Rosenhain form for these curves. We demonstrate that two natural one-dimensional subloci of the space of $SU(3)$ vacua, $\mathcal{E}_u$ and $\mathcal{E}_v$, each parametrize a family of elliptic curves. For these elliptic loci, we describe the order parameters and fundamental domains explicitly. The locus $\mathcal{E}_u$ contains the points where mutually local dyons become massless, and is a fundamental domain for a classical congruence subgroup. Moreover, the locus $\mathcal{E}_v$ contains the superconformal Argyres-Douglas points, and is a fundamental domain for a Fricke group.


title: Quantum autoencoders with enhanced data encoding

abstract: We present the enhanced feature quantum autoencoder, or EF-QAE, a variational quantum algorithm capable of compressing quantum states of different models with higher fidelity. The key idea of the algorithm is to define a parameterized quantum circuit that depends upon adjustable parameters and a feature vector that characterizes such a model. We assess the validity of the method in simulations by compressing ground states of the Ising model and classical handwritten digits. The results show that EF-QAE improves the performance compared to the standard quantum autoencoder using the same amount of quantum resources, but at the expense of additional classical optimization. Therefore, EF-QAE makes the task of compressing quantum information better suited to be implemented in near-term quantum devices.


title: Stacking-engineered ferroelectricity in bilayer boron nitride

abstract: 2D ferroelectrics with robust polarization down to atomic thicknesses provide novel building blocks for functional heterostructures. Experimental reports, however, remain scarce because of the requirement of a layered polar crystal. Here, we demonstrate a rational design approach to engineering 2D ferroelectrics from a non-ferroelectric parent compound via employing van der Waals assembly. Parallel-stacked bilayer boron nitride is shown to exhibit out-of-plane electric polarization that reverses depending on the stacking order. The polarization switching is probed via the resistance of an adjacently-stacked graphene sheet. Furthermore, twisting the boron nitride sheets by a small-angle changes the dynamics of switching due to the formation of moir\'e ferroelectricity with staggered polarization. The ferroelectricity persists to room temperature while keeping the high mobility of graphene, paving the way for potential ultrathin nonvolatile memory applications.


title: Production of axion-like particles from photon conversions in large-scale solar magnetic fields

abstract: The Sun is a well-studied astrophysical source of axion-like particles (ALPs), produced mainly through the Primakoff process. Moreover, in the Sun there exist large-scale magnetic fields that catalyze an additional ALP production via a coherent conversion of thermal photons. We study this contribution to the solar ALP emissivity, typically neglected in previous investigations. Furthermore, we discuss additional bounds on the ALP-photon coupling from energy-loss arguments, and the detection perspectives of this new ALP flux at future helioscope and dark matter experiments.


title: Observations of Hawking radiation: the Page curve and baby universes

abstract: We reformulate recent insights into black hole information in a manner emphasizing operationally-defined notions of entropy, Lorentz-signature descriptions, and asymptotically flat spacetimes. With the help of replica wormholes, we find that experiments of asymptotic observers are consistent with black holes as unitary quantum systems, with density of states given by the Bekenstein-Hawking formula. However, this comes at the cost of superselection sectors associated with the state of baby universes. Spacetimes studied by Polchinski and Strominger in 1994 provide a simple illustration of the associated concepts and techniques, and we argue them to be a natural late-time extrapolation of replica wormholes. The work aims to be self-contained and, in particular, to be accessible to readers who have not yet mastered earlier formulations of the ideas above.


title: The RR Lyrae projected density distribution from the Galactic centre to the halo

abstract: The projected density distribution of type ab RR Lyrae (RRab) stars was characterised from the innermost regions of the Milky Way to the halo, with the aim of placing constraints on the Galaxys evolution. The compiled sample (N_RRab = 64,850) stems from fundamental mode RR Lyrae variables identified by the VVV, OGLE, and Gaia surveys. The distribution is well fitted by three power laws over three radial intervals. In the innermost region (R < 2.2 deg) the distribution follows Sigma_RRab[1] propto R ^(-0.94 +- 0.051), while in the external region the distribution adheres to Sigma_RRab[2] propto R^(-1.50 +- 0.019) for 2.2 deg< R <8.0 deg and Sigma_RRab[3] propto R ^(-2.43 +- 0.043) for 8.0 deg < R <30.0 deg. Conversely, the cumulative distribution of red clump (RC) giants exhibits a more concentrated distribution in the mean, but in the central R < 2.2 deg the RRab population is more peaked, whereas globular clusters (GCs) follow a density power law (Sigma_GCs propto R ^(-1.59 +- 0.060) for R<30.0 deg) similar to that of RRab stars, especially when considering a more metal-poor subsample ([Fe/H]<-1.1 dex). The main conclusion emerging from the analysis is that the RRab distribution favours the star cluster infall and merger scenario for creating an important fraction (>18 %) of the central Galactic region. The radii containing half of the populations (half populations radii) are R_H=6.8 deg (0.99 kpc), R_H =4.2 deg (0.61 kpc), and R_H =11.9 deg (1.75 kpc) for the RRab stars, RC giants, and GCs, respectively. Finally, merely 1% of the stars have been actually discovered in the innermost region (R < 35 pc) out of the expected (based on our considerations) total number of RRab therein: N sim 1,562. That deficit will be substantially ameliorated with future space missions like the Nancy Grace Roman Space Telescope (formerly WFIRST).


title: Markovian Dynamics in de Sitter

abstract: The equilibrium state of fields in the causal wedge of a dS observer is thermal, though realistic observers have only partial access to the state. To them, out-of-equilibrium states of a light scalar field appear to thermalize in a Markovian fashion. We show this by formulating a systematic expansion for tracing out the environment. As an example, we calculate the $O(\lambda)$ correction to the result of Starobinsky and Yokoyama for the relaxation exponents of $\lambda \phi^4$ theory.


title: Tracing the evolution of dust-obscured activity using sub-millimetre galaxy populations from STUDIES and AS2UDS

abstract: We analyse the physical properties of 121 SNR $\geq$ 5 sub-millimetre galaxies (SMGs) from the STUDIES 450-$\mu$m survey. We model their UV-to-radio spectral energy distributions using MAGPHYS+photo-$z$ and compare the results to similar modelling of 850-$\mu$m-selected SMG sample from AS2UDS, to understand the fundamental physical differences between the two populations at the observed depths. The redshift distribution of the 450-$\mu$m sample has a median of $z$ = 1.85 $\pm$ 0.12 and can be described by strong evolution of the far-infrared luminosity function. The fainter 450-$\mu$m sample has $\sim$14 times higher space density than the brighter 850-$\mu$m sample at $z$ $\lesssim$2, and a comparable space density at $z$ = 2-3, before rapidly declining, suggesting LIRGs are the main obscured population at $z$ $\sim$ 1-2, while ULIRGs dominate at higher redshifts. We construct rest-frame $\sim$ 180-$\mu$m-selected and dust-mass-matched samples at $z$ = 1-2 and $z$ = 3-4 from the 450-$\mu$m and 850-$\mu$m samples, respectively, to probe the evolution of a uniform sample of galaxies spanning the cosmic noon era. Using far-infrared luminosity, dust masses and an optically-thick dust model, we suggest that higher-redshift sources have higher dust densities due to inferred dust continuum sizes which are roughly half of those for the lower-redshift population at a given dust mass, leading to higher dust attenuation. We track the evolution in the cosmic dust mass density and suggest that the dust content of galaxies is governed by a combination of both the variation of gas content and dust destruction timescale.


title: A Pareto Dominance Principle for Data-Driven Optimization

abstract: We propose a statistically optimal approach to construct data-driven decisions for stochastic optimization problems. Fundamentally, a data-driven decision is simply a function that maps the available training data to a feasible action. It can always be expressed as the minimizer of a surrogate optimization model constructed from the data. The quality of a data-driven decision is measured by its out-of-sample risk. An additional quality measure is its out-of-sample disappointment, which we define as the probability that the out-of-sample risk exceeds the optimal value of the surrogate optimization model. An ideal data-driven decision should minimize the out-of-sample risk simultaneously with respect to every conceivable probability measure as the true measure is unkown. Unfortunately, such ideal data-driven decisions are generally unavailable. This prompts us to seek data-driven decisions that minimize the in-sample risk subject to an upper bound on the out-of-sample disappointment. We prove that such Pareto-dominant data-driven decisions exist under conditions that allow for interesting applications: the unknown data-generating probability measure must belong to a parametric ambiguity set, and the corresponding parameters must admit a sufficient statistic that satisfies a large deviation principle. We can further prove that the surrogate optimization model must be a distributionally robust optimization problem constructed from the sufficient statistic and the rate function of its large deviation principle. Hence the optimal method for mapping data to decisions is to solve a distributionally robust optimization model. Maybe surprisingly, this result holds even when the training data is non-i.i.d. Our analysis reveals how the structural properties of the data-generating stochastic process impact the shape of the ambiguity set underlying the optimal distributionally robust model.


title: Towards Runtime Verification via Event Stream Processing in Cloud Computing Infrastructures

abstract: Software bugs in cloud management systems often cause erratic behavior, hindering detection, and recovery of failures. As a consequence, the failures are not timely detected and notified, and can silently propagate through the system. To face these issues, we propose a lightweight approach to runtime verification, for monitoring and failure detection of cloud computing systems. We performed a preliminary evaluation of the proposed approach in the OpenStack cloud management platform, an "off-the-shelf" distributed system, showing that the approach can be applied with high failure detection coverage.


title: AI-assisted super-resolution cosmological simulations

abstract: Cosmological simulations of galaxy formation are limited by finite computational resources. We draw from the ongoing rapid advances in Artificial Intelligence (specifically Deep Learning) to address this problem. Neural networks have been developed to learn from high-resolution (HR) image data, and then make accurate super-resolution (SR) versions of different low-resolution (LR) images. We apply such techniques to LR cosmological N-body simulations, generating SR versions. Specifically, we are able to enhance the simulation resolution by generating 512 times more particles and predicting their displacements from the initial positions. Therefore our results can be viewed as new simulation realizations themselves rather than projections, e.g., to their density fields. Furthermore, the generation process is stochastic, enabling us to sample the small-scale modes conditioning on the large-scale environment. Our model learns from only 16 pairs of small-volume LR-HR simulations, and is then able to generate SR simulations that successfully reproduce the HR matter power spectrum to percent level up to $16\,h^{-1}\mathrm{Mpc}$, and the HR halo mass function to within $10 \%$ down to $10^{11} \, M_\odot$. We successfully deploy the model in a box 1000 times larger than the training simulation box, showing that high-resolution mock surveys can be generated rapidly. We conclude that AI assistance has the potential to revolutionize modeling of small-scale galaxy formation physics in large cosmological volumes.


title: Hubble Space Telescope Astrometry of the Metal-Poor Visual Binary $\mu$ Cassiopeiae: Dynamical Masses, Helium Content, and Age

abstract: $\mu$ Cassiopeiae is a nearby, high-velocity, metal-poor ($\rm[Fe/H]=-0.81$) visual binary. We have used high-resolution imaging with the Hubble Space Telescope (HST), obtained over nearly two decades, to determine the period (21.568 yr) and precise orbital elements. Combining these with published ground- and space-based astrometry, we determined dynamical masses for both components of $\mu$ Cas: $0.7440\pm0.0122\,M_\odot$ for the G5 V primary, and $0.1728\pm0.0035\,M_\odot$ for its faint dM companion. We detect no significant perturbations in the HST astrometry due to a third body in the system. The primary aim of our program was to determine, with the aid of stellar models, the helium content and age of the metal-deficient primary star, $\mu$ Cas A. Although we now have a precise mass, there remain uncertainties about other parameters, including its effective temperature. Moreover, a re-examination of archival interferometric observations leads to a suspicion that the angular diameter was overestimated by a few percent. In the absolute magnitude versus color plane, $\mu$ Cas A lies slightly cooler and more luminous than the main sequence of the globular cluster 47 Tucanae; this may imply that the star has a lower helium content, and/or is older, and/or has a higher metallicity, than the cluster. Our best estimates for the helium content and age of $\mu$ Cas A are $Y=0.255\pm0.014$ and $12.7\pm2.7$ Gyr--making $\mu$ Cas possibly the oldest star in the sky visible to the naked eye. Improved measurements of the absolute parallax of the system, the effective temperature of $\mu$ Cas A, and its angular diameter would provide tighter constraints.


title: Training independent subnetworks for robust prediction

abstract: Recent approaches to efficiently ensemble neural networks have shown that strong robustness and uncertainty performance can be achieved with a negligible gain in parameters over the original network. However, these methods still require multiple forward passes for prediction, leading to a significant computational cost. In this work, we show a surprising result: the benefits of using multiple predictions can be achieved `for free' under a single model's forward pass. In particular, we show that, using a multi-input multi-output (MIMO) configuration, one can utilize a single model's capacity to train multiple subnetworks that independently learn the task at hand. By ensembling the predictions made by the subnetworks, we improve model robustness without increasing compute. We observe a significant improvement in negative log-likelihood, accuracy, and calibration error on CIFAR10, CIFAR100, ImageNet, and their out-of-distribution variants compared to previous methods.


title: Listening to the quantum vacuum: a perspective on the dynamical Casimir effect

abstract: Modern quantum field theory has offered us a very intriguing picture of empty space. The vacuum state is no longer an inert, motionless state. We are instead dealing with an entity teeming with fluctuations that continuously produce virtual particles popping in and out of existence. The dynamical Casimir effect is a paradigmatic phenomenon, whereby these particles are converted into real particles (photons) by changing the boundary conditions of the field. It was predicted 50 years ago by Gerald T. Moore and it took more than 40 years until the first experimental verification.


title: Late-Time Radio and Millimeter Observations of Superluminous Supernovae and Long Gamma Ray Bursts: Implications for Obscured Star Formation, Central Engines, and Fast Radio Bursts

abstract: We present the largest and deepest late-time radio and millimeter survey to date of superluminous supernovae (SLSNe) and long duration gamma-ray bursts (LGRBs) to search for associated non-thermal synchrotron emission. Using the Karl G. Jansky Very Large Array (VLA) and the Atacama Large Millimeter/submillimeter Array (ALMA), we observed 43 sources at 6 and 100 GHz on a timescale of $\sim 1 - 19$ yr post-explosion. We do not detect radio/mm emission from any of the sources, with the exception of a 6 GHz detection of PTF10hgi (Eftekhari et al. 2019), as well as the detection of 6 GHz emission near the location of the SLSN PTF12dam, which we associate with its host galaxy. We use our data to place constraints on central engine emission due to magnetar wind nebulae and off-axis relativistic jets. We also explore non-relativistic emission from the SN ejecta, and place constraints on obscured star formation in the host galaxies. In addition, we conduct a search for fast radio bursts (FRBs) from some of the sources using VLA Phased-Array observations; no FRBs are detected to a limit of $16$ mJy ($7\sigma$; 10 ms duration) in about 40 min on source per event. A comparison to theoretical models suggests that continued radio monitoring may lead to detections of persistent radio emission on timescales of $\gtrsim {\rm decade}$.


title: Intelligent Reflecting Surface-Assisted Wireless Key Generation for Low-Entropy Environments

abstract: Physical layer key generation is a promising candidate for cryptographic key establishment between two wireless communication parties. It offers information-theoretic security and is an attractive alternative to public-key techniques. Here, the inherent randomness of wireless radio channels is used as a shared entropy source to generate cryptographic key material. However, practical implementations often suffer from static channel conditions which exhibit a limited amount of randomness. In the past, considerable research efforts have been made to address this fundamental limitation. However, current solutions are not generic or require dedicated hardware extensions such as reconfigurable antennas. In this paper, we propose a novel wireless key generation architecture based on randomized channel responses from an intelligent reflecting surface (IRS). Due to its passive nature, a cooperative IRS is well-suited to provide randomness for conventional resource-constrained radios. We conduct the first practical studies to successfully demonstrate IRS-based physical-layer key generation with an OFDM system. In a static environment, using a single subcarrier only, our IRS-assisted prototype system achieves a key generation rate (KGR) of 97.39 bps with 6.5% key disagreement rate (KDR) after quantization, while passing standard randomness tests.


title: Mass ratio, the hills mechanism, and the galactic centre S-stars

abstract: The Galactic centre contains several young populations within its central parsec: a disk between $\sim$0.05 and 0.5 pc from the centre, and the isotropic S-star cluster extending an order of magnitude further inwards in radius. Recent observations (i.e. spectroscopy and hypervelocity stars) suggest that some S-stars originate in the disk. In particular, the S-stars may be remnants of tidally disrupted disk binaries. However, there is an apparent inconsistency in this scenario: the disk contains massive O and Wolf--Rayet stars while the S-stars are lower mass, B stars. We explore two different explanations for this apparent discrepancy: (i) a built-in bias in binary disruptions, where the primary star remains closer in energy to the centre-of-mass orbit than the secondary and (ii) selective tidal disruption of massive stars within the S-star cluster. The first explanation is plausible. On the other hand, tidal disruptions have not strongly affected the mass distribution of the S-stars over the last several Myr.


title: Effects of spin-orbit torque on the ferromagnetic and exchange spin wave modes in ferrimagnetic CoGd alloy

abstract: We use micro-focus Brillouin light scattering spectroscopy to study the effects of spin-orbit torque on thermal spin waves in almost angular-momentum compensated ferrimagnetic CoGd alloy films. The spin-orbit torque is produced by the electric current flowing in the Pt layer adjacent to CoGd. Both the ferromagnetic and the exchange modes are detected in our measurements. The intensity and the linewidth of the ferromagnetic mode are modified by the spin-orbit torque. In contrast, the properties of the exchange mode are unaffected by the spin-orbit torque. We also find that the frequencies and the linewidths of both modes are significantly modified by Joule heating, due to the strong temperature dependence of the magnetic properties of CoGd in the vicinity of angular momentum compensation point. Our results provide insight into the mechanisms that can enable the implementation of sub-THz magnetic nano-oscillators based on ferrimagnetic materials, as well as related effects in antiferromagnets.


title: Finite-Time Model Inference From A Single Noisy Trajectory

abstract: This paper proposes a novel model inference procedure to identify system matrix from a single noisy trajectory over a finite-time interval. The proposed inference procedure comprises an observation data processor, a redundant data processor and an ordinary least-square estimator, wherein the data processors mitigate the influence of observation noise on inference error. We first systematically investigate the comparisons with naive least-square-regression based model inference and uncover that 1) the same observation data has identical influence on the feasibility of the proposed and the naive model inferences, 2) the naive model inference uses all of the redundant data, while the proposed model inference optimally uses the basis and the redundant data. We then study the sample complexity of the proposed model inference in the presence of observation noise, which leads to the dependence of the processed bias in the observed system trajectory on time and coordinates. Particularly, we derive the sample-complexity upper bound (on the number of observations sufficient to infer a model with prescribed levels of accuracy and confidence) and the sample-complexity lower bound (high-probability lower bound on model error). Finally, the proposed model inference is numerically validated and analyzed.


title: HI 21-centimetre emission from an ensemble of galaxies at an average redshift of one

abstract: The baryonic processes in galaxy evolution include gas infall onto galaxies to form neutral atomic hydrogen (HI), the conversion of HI to the molecular state (H$_2$), and, finally, the conversion of H$_2$ to stars. Understanding galaxy evolution thus requires understanding the evolution of both the stars, and the neutral atomic and molecular gas, the primary fuel for star-formation, in galaxies. For the stars, the cosmic star-formation rate density is known to peak in the redshift range $z \approx 1-3$, and to decline by an order of magnitude over the next $\approx 10$ billion years; the causes of this decline are not known. For the gas, the weakness of the hyperfine HI 21cm transition, the main tracer of the HI content of galaxies, has meant that it has not hitherto been possible to measure the atomic gas mass of galaxies at redshifts higher than $\approx 0.4$; this is a critical lacuna in our understanding of galaxy evolution. Here, we report a measurement of the average HI mass of star-forming galaxies at a redshift $z \approx 1$, by stacking their individual HI 21 cm emission signals. We obtain an average HI mass similar to the average stellar mass of the sample. We also estimate the average star-formation rate of the same galaxies from the 1.4 GHz radio continuum, and find that the HI mass can fuel the observed star-formation rates for only $\approx 1-2$ billion years in the absence of fresh gas infall. This suggests that gas accretion onto galaxies at $z < 1$ may have been insufficient to sustain high star-formation rates in star-forming galaxies. This is likely to be the cause of the decline in the cosmic star-formation rate density at redshifts below 1.


title: Scaling Laws for Three-Dimensional Combined Heaving and Pitching Propulsors

abstract: We present new scaling laws for the thrust production and power consumption of three-dimensional combined heaving and pitching hydrofoils by extending the three-dimensional pitching scaling laws introduced by Ayancik et al. (2019). New self-propelled inviscid simulations and previously published experimental data are used to validate the scaling laws over a wide range of motion amplitudes, Strouhal numbers, heave ratios, aspect ratios, and pitching axis locations. The developed scaling laws are shown to predict inviscid numerical data and experimental data well, within $\pm$ 25% and $\pm$ 16% of the thrust and power data, respectively. The scaling laws reveal that both the circulatory and added mass forces are important when considering a wide range of motion amplitudes and that nonlinear corrections to the classic linear theory are essential to modeling the performance across this wide amplitude range. By using the scaling laws as a guide, it is determined that peak efficiencies occur when $A^*$ > 1 and for these large-amplitude motions, there is an optimal $h^*$ that maximizes the efficiency in the narrow range of 0.75 < $h^*$ < 0.94. Finally, the scaling laws show that to further improve efficiency in this high-efficiency regime, the aspect ratio, and dimensionless amplitude should be increased, while the Lighthill number should be decreased (lower drag and/or a larger propulsor planform area to wetted surface area ratio), and the pitch axis should be located behind the leading edge. This scaling model can be used to guide the design of the next generation of high-efficiency bio-inspired machines.


title: Trigonal Warping, Satellite Dirac Points and Multiple Field Tuned Topological Transitions in Twisted Double Bilayer Graphene

abstract: We show that the valley Chern number of the low energy band in twisted double bilayer graphene can be tuned through two successive topological transitions, where the direct bandgap closes, by changing the electric field perpendicular to the plane of the graphene layers. The two transitions with Chern number changes of -3 and +1 can be explained by the formation of three satellite Dirac points around the central Dirac cone in the moir\'e Brillouin zone due to the presence of trigonal warping. The satellite cones have opposite chirality to the central Dirac cone. Considering the overlap of the bands in energy, which lead to metallic states, we construct the experimentally observable phase diagram of the system in terms of the indirect bandgap and the anomalous valley Hall conductivity. We show that while most of the intermediate phase becomes metallic, there is a narrow parameter regime where the transition through three insulating phases with different quantized valley Hall conductivity can be seen. We systematically study the effects of variations in the model parameters on the phase diagram of the system to reveal the importance of particle-hole asymmetry and trigonal warping in constructing the phase diagram. We also study the effect of changes in interlayer tunneling on this phase diagram.


title: Improved limits for violations of local position invariance from atomic clock comparisons

abstract: We compare two optical clocks based on the $^2$S$_{1/2}(F=0)\to {}^2$D$_{3/2}(F=2)$ electric quadrupole (E2) and the $^2$S$_{1/2}(F=0)\to {}^2$F$_{7/2}(F=3)$ electric octupole (E3) transition of $^{171}$Yb$^{+}$ and measure the frequency ratio $\nu_{\mathrm{E3}}/\nu_{\mathrm{E2}}=0.932\,829\,404\,530\,965\,376(32)$. We determine the transition frequency $\nu_{E3}=642\,121\,496\,772\,645.10(8)$ Hz using two caesium fountain clocks. Repeated measurements of both quantities over several years are analyzed for potential violations of local position invariance. We improve by factors of about 20 and 2 the limits for fractional temporal variations of the fine structure constant $\alpha$ to $1.0(1.1)\times10^{-18}/\mathrm{yr}$ and of the proton-to-electron mass ratio $\mu$ to $-8(36)\times10^{-18}/\mathrm{yr}$. Using the annual variation of the Sun's gravitational potential at Earth $\Phi$, we improve limits for a potential coupling of both constants to gravity, $(c^2/\alpha) (d\alpha/d\Phi)=14(11)\times 10^{-9}$ and $(c^2/\mu) (d\mu/d\Phi)=7(45)\times 10^{-8}$.


title: Monte Carlo simulations for the ANTARES underwater neutrino telescope

abstract: Monte Carlo simulations are a unique tool to check the response of a detector and to monitor its performance. For a deep-sea neutrino telescope, the variability of the environmental conditions that can affect the behaviour of the data acquisition system must be considered, in addition to a reliable description of the active parts of the detector and of the features of physics events, in order to produce a realistic set of simulated events. In this paper, the software tools used to produce neutrino and cosmic ray signatures in the telescope and the strategy developed to represent the time evolution of the natural environment and of the detector efficiency are described.


title: CISE3: Verifying Weakly Consistent Applications with Why3

abstract: In this paper we present a tool for the formal analysis of applications built on top of replicated databases, where data integrity can be at stake. To address this issue, one can introduce synchronization in the system. Introducing synchronization in too many places can hurt the system's availability but if introduced in too few places, then data integrity can be compromised. The goal of our tool is to aid the programmer reason about the correct balance of synchronization in the system. Our tool analyses a sequential specification and deduces which operations require synchronization in order for the program to safely execute in a distributed environment. Our prototype is built on top of the deductive verification platform Why3, which provides a friendly and integrated user experience. Several case studies have been successfully verified using our tool.


title: Renormalon Effects in Quasi Parton Distributions

abstract: We investigate the renormalon ambiguity from bubble-chain diagrams in the isovector unpolarized quasi-parton distribution function (PDF) of a hadron. We confirm the assertion by Braun, Vladimirov and Zhang (arXiv:1810.00048) that the leading IR renormalon ambiguity is formally an $\mathcal{O}({\Lambda}^2_{QCD}/x^2P_z^2)$ effect, with $x$ the parton momentum fraction and $P_z$ the hadron momentum, together with a new $\mathcal{O}(\delta(x)\Lambda^2_{QCD}/P_z^2)$ contribution such that the quark number is conserved. This implies the convergence of the perturbative matching kernel between a quasi-PDF and a PDF would eventually fail for small $x$. However, in both the R-scheme designed to cancel the leading IR renormalon and the typically used RI/MOM scheme in lattice QCD for the same quasi-PDF, we find good convergence in the kernel based on three-loop bubble-chain diagram analyses. These results are encouraging for the quasi-PDF program. However, firm conclusions can only be drawn after the complete higher loop QCD calculations are carried out.


title: A study of the stellar photosphere-hydrogen ionisation front interaction in pulsating variables using period-color relations

abstract: Period-color (PC) relations may be used to study the interaction of the stellar photosphere and the hydrogen ionization front (HIF). RR Lyraes (RRLs) and long period classical Cepheids (P > 10d) have been found to exhibit different PC behavior at minimum and maximum light which can be explained by the HIF-photosphere interaction based on their location on the HR diagram. In this work, we extend the study to include type II Cepheids (T2Cs) with an aim to test the HIF-photosphere interaction theory across a broad spectrum of variable star types. We find W Vir stars and BL Her stars to have similar PC relations as those from long period and short period classical Cepheids, respectively. We also use MESA to compute RRL, BL Her and classical Cepheid models to study the theoretical HIF-photosphere distance and find the results to be fairly consistent with the HIF-photosphere interaction theory.


title: Full-field stress computation from measured deformation fields: a hyperbolic formulation

abstract: Recent developments in imaging techniques and correlation algorithms enable measurement of strain fields on a deforming material at high spatial and temporal resolution. In such cases, the computation of the stress field from the known deformation field becomes an interesting possibility. This is known as an inverse problem. Current approaches to this problem, such as the finite element update method, are generally over-determined and must rely on statistical approaches to minimize error. This provides approximate solutions in some cases, however, implementation difficulties, computational requirements, and accuracy are still significant challenges. Here, we show how the inverse problem can be formulated deterministically and solved exactly in two or three dimensions for large classes of materials including isotropic elastic solids, Newtonian fluids, non-Newtonian fluids, granular materials and isotropic plastic solids. This solution is based on a single assumption of the alignment of the principal directions of stress and strain or strain rate. No further assumptions regarding incompressibility, pressure independence, yield surface shape or the hardening law are necessary. This assumption leads to a closed, first order, linear system of hyperbolic partial differential equations with variable coefficients and can be solved to give the solution for any geometry and loading condition. We provide a numerical proof-of-principle study of the plastic deformation of a two-dimensional bar with spatially varying yield stress and strain hardening coefficient. The results are validated indicating the solution is exact up to numerical error. No model calibration or material parameters are required. Interestingly, this solution procedure lends itself to a simple physical interpretation of stress propagation through the material.


title: On Deep Learning Techniques to Boost Monocular Depth Estimation for Autonomous Navigation

abstract: Inferring the depth of images is a fundamental inverse problem within the field of Computer Vision since depth information is obtained through 2D images, which can be generated from infinite possibilities of observed real scenes. Benefiting from the progress of Convolutional Neural Networks (CNNs) to explore structural features and spatial image information, Single Image Depth Estimation (SIDE) is often highlighted in scopes of scientific and technological innovation, as this concept provides advantages related to its low implementation cost and robustness to environmental conditions. In the context of autonomous vehicles, state-of-the-art CNNs optimize the SIDE task by producing high-quality depth maps, which are essential during the autonomous navigation process in different locations. However, such networks are usually supervised by sparse and noisy depth data, from Light Detection and Ranging (LiDAR) laser scans, and are carried out at high computational cost, requiring high-performance Graphic Processing Units (GPUs). Therefore, we propose a new lightweight and fast supervised CNN architecture combined with novel feature extraction models which are designed for real-world autonomous navigation. We also introduce an efficient surface normals module, jointly with a simple geometric 2.5D loss function, to solve SIDE problems. We also innovate by incorporating multiple Deep Learning techniques, such as the use of densification algorithms and additional semantic, surface normals and depth information to train our framework. The method introduced in this work focuses on robotic applications in indoor and outdoor environments and its results are evaluated on the competitive and publicly available NYU Depth V2 and KITTI Depth datasets.


title: Video Game Level Repair via Mixed Integer Linear Programming

abstract: Recent advancements in procedural content generation via machine learning enable the generation of video-game levels that are aesthetically similar to human-authored examples. However, the generated levels are often unplayable without additional editing. We propose a generate-then-repair framework for automatic generation of playable levels adhering to specific styles. The framework constructs levels using a generative adversarial network (GAN) trained with human-authored examples and repairs them using a mixed-integer linear program (MIP) with playability constraints. A key component of the framework is computing minimum cost edits between the GAN generated level and the solution of the MIP solver, which we cast as a minimum cost network flow problem. Results show that the proposed framework generates a diverse range of playable levels, that capture the spatial relationships between objects exhibited in the human-authored levels.


title: The surface code on the rhombic dodecahedron

abstract: I present a jaunty little $[[14, 3, 3]]$ non-CSS surface code that can be described using a rhombic dodecahedron. Do with it what you will.


title: Interferometric geometry from symmetry-broken Uhlmann gauge group with applications to topological phase transitions

abstract: We provide a natural generalization of a Riemannian structure, i.e., a metric, recently introduced by Sj\"{o}qvist for the space of non degenerate density matrices, to the degenerate case, i.e., the case in which the eigenspaces have dimension greater than or equal to 1. We present a physical interpretation of the metric in terms of an interferometric measurement. We apply this metric, physically interpreted as an interferometric susceptibility, to the study of topological phase transitions at finite temperatures for band insulators. We compare the behaviors of this susceptibility and the one coming from the well-known Bures metric, showing them to be dramatically different. While both infer zero temperature phase transitions, only the former predicts finite temperature phase transitions as well. The difference in behaviors can be traced back to a symmetry breaking mechanism, akin to Landau-Ginzburg theory, by which the Uhlmann gauge group is broken down to a subgroup determined by the type of the system's density matrix (i.e., the ranks of its spectral projectors).

categories: quant-ph cond-mat.str-el math-ph math.MP


title: Motivations and Preliminary Design for Mid-Air Deployment of a Science Rotorcraft on Mars

abstract: Mid-Air Deployment (MAD) of a rotorcraft during Entry, Descent and Landing (EDL) on Mars eliminates the need to carry a propulsion or airbag landing system. This reduces the total mass inside the aeroshell by more than 100 kg and simplifies the aeroshell architecture. MAD's lighter and simpler design is likely to bring the risk and cost associated with the mission down. Moreover, the lighter entry mass enables landing in the Martian highlands, at elevations inaccessible to current EDL technologies. This paper proposes a novel MAD concept for a Mars helicopter. We suggest a minimum science payload package to perform relevant science in the highlands. A variant of the Ingenuity helicopter is proposed to provide increased deceleration during MAD, and enough lift to fly the science payload in the highlands. We show in simulation that the lighter aeroshell results in a lower terminal velocity (30 m/s) at the end of the parachute phase of the EDL, and at higher altitudes than other approaches. After discussing the aerodynamics, controls, guidance, and mechanical challenges associated with deploying at such speed, we propose a backshell architecture that addresses them to release the helicopter in the safest conditions. Finally, we implemented the helicopter model and aerodynamic descent perturbations in the JPL Dynamics and Real-Time Simulation (DARTS)framework. Preliminary performance evaluation indicates landing and helicopter operation scan be achieved up to 5 km MOLA (Mars Orbiter Laser Altimeter reference).


title: Succinct Explanations With Cascading Decision Trees

abstract: The decision tree is one of the most popular and classical machine learning models from the 1980s. However, in many practical applications, decision trees tend to generate decision paths with excessive depth. Long decision paths often cause overfitting problems, and make models difficult to interpret. With longer decision paths, inference is also more likely to fail when the data contain missing values. In this work, we propose a new tree model called Cascading Decision Trees to alleviate this problem. The key insight of Cascading Decision Trees is to separate the decision path and the explanation path. Our experiments show that on average, Cascading Decision Trees generate 63.38% shorter explanation paths, avoiding overfitting and thus achieve higher test accuracy. We also empirically demonstrate that Cascading Decision Trees have advantages in the robustness against missing values.


title: A Fully General, Non-Perturbative Treatment of Impulsive Heating

abstract: Impulsive encounters between astrophysical objects are usually treated using the distant tide approximation (DTA) for which the impact parameter, $b$, is assumed to be significantly larger than the characteristic radii of the subject, $r_{\mathrm{S}}$, and the perturber, $r_{\mathrm{P}}$. The perturber potential is then expanded as a multipole series and truncated at the quadrupole term. When the perturber is more extended than the subject, this standard approach can be extended to the case where $r_{\mathrm{S}} \ll b < r_{\mathrm{P}}$. However, for encounters with $b$ of order $r_{\mathrm{S}}$ or smaller, the DTA typically overpredicts the impulse, $\Delta \mathbf{v}$, and hence the internal energy change of the subject, $\Delta E_{\mathrm{int}}$. This is unfortunate, as these close encounters are the most interesting, potentially leading to tidal capture, mass stripping, or tidal disruption. Another drawback of the DTA is that $\Delta E_{\mathrm{int}}$ is proportional to the moment of inertia, which diverges unless the subject is truncated or has a density profile that falls off faster than $r^{-5}$. To overcome these shortcomings, this paper presents a fully general, non-perturbative treatment of impulsive encounters which is valid for any impact parameter, and not hampered by divergence issues, thereby negating the necessity to truncate the subject. We present analytical expressions for $\Delta \mathbf{v}$ for a variety of perturber profiles, apply our formalism to both straight-path encounters and eccentric orbits, and discuss the mass loss due to tidal shocks in gravitational encounters between equal mass galaxies.


title: Large cross-polarized Raman signal in CrI$_3$: A first-principles study

abstract: We find unusually large cross-polarized (and anti-symmetric) Raman signature of A$_{\rm g}$ phonon mode in CrI$_3$, in agreement with experiments. The signal is present only when the following three effects are considered in concert: ferromagnetism on Cr atoms, spin-orbit interaction, and resonant effects. Somewhat surprisingly, we find that the relevant spin-orbit interaction potential originates from iodine atoms, despite magnetism being mostly on chromium atoms. We analyze the Raman signature as a function of magnetic order, the direction of the magnetic moment, energy and polarization of light used for Raman scattering, as well as carrier lifetime. In addition to a strong cross-polarized Raman signal, we also find unusually strong phonon modulated magneto-optical Kerr effect (MOKE) in CrI$_3$.


title: Functional Theory for Bose-Einstein Condensates

abstract: One-particle reduced density matrix functional theory would potentially be the ideal approach for describing Bose-Einstein condensates. It namely replaces the macroscopically complex wavefunction by the simple one-particle reduced density matrix, therefore provides direct access to the degree of condensation and still recovers quantum correlations in an exact manner. We eventually initiate and establish this novel theory by deriving the respective universal functional $\mathcal{F}$ for general homogeneous Bose-Einstein condensates with arbitrary pair interaction. Most importantly, the successful derivation necessitates a particle-number conserving modification of Bogoliubov theory and a solution of the common phase dilemma of functional theories. We then illustrate this novel approach in several bosonic systems such as homogeneous Bose gases and the Bose-Hubbard model. Remarkably, the general form of $\mathcal{F}$ reveals the existence of a universal Bose-Einstein condensation force which provides an alternative and more fundamental explanation for quantum depletion.


title: Measuring the cosmic homogeneity scale with SDSS-IV DR16 Quasars

abstract: We report measurements of the scale of cosmic homogeneity ($r_{h}$) using the recently released quasar sample of the sixteenth data release of the Sloan Digital Sky Survey (SDSS-IV DR16). We perform our analysis in 2 redshift bins lying in the redshift interval $2.2 < z < 3.2$ by means of the fractal dimension $D_2$. By adopting the usual assumption that $r_{h}$ is obtained when $D_2 \sim 2.97$, that is, within 1% of $D_2=3$, we find the cosmic homogeneity scale with a decreasing trend with redshift, and in good agreement with the $\Lambda$CDM prediction. Our results confirm the presence of a homogeneity scale in the spatial distribution of quasars as predicted by the fundamental assumptions of the standard cosmological model.


title: Electromagnetic emission due to nonlinear interaction of laser wakefields colliding in plasma at an oblique angle

abstract: Head-on collision of laser-induced plasma wakefields with differing profiles of electrostatic potential has been recently found to be an efficient mechanism for generating high-power electromagnetic emission at the second harmonic of the plasma frequency [I.V.Timofeev et al. Phys. Plasmas 24, 103106 (2017)]. This mechanism is attractive for creating a source of tunable narrow-band coherent radiation in the terahertz frequency range. In this paper, we generalize the theory of electromagnetic emission produced by nonlinear interaction of two plasma wakes to the case of an arbitrary collision angle. Such a theory is used to evaluate the angular distribution of the second harmonic radiation as well as its total generation efficiency for parameters of the proof-of-principle experiment in which laser axes will be aligned with a small finite angle. Theoretical predictions are qualitatively confirmed by particle-in-cell simulations.


title: Type 0 $\mathbb{Z}_2\times \mathbb{Z}_2$ Heterotic String Orbifolds and Misaligned Supersymmetry

abstract: The $\mathbb{Z}_2\times \mathbb{Z}_2$ heterotic string orbifold yielded a large space of phenomenological three generation models and serves as a testing ground to explore how the Standard Model of particle physics may be incorporated in a theory of quantum gravity. In this paper we explore the existence of type 0 models in this class of string compactifications. We demonstrate the existence of type 0 $\mathbb{Z}_2\times \mathbb{Z}_2$ heterotic string orbifolds, and show that there exist a large degree of redundancy in the space of GGSO projection coefficients when the type 0 restrictions are implemented. We explore the existence of such configurations in several constructions. The first correspond to essentially a unique configuration out of a priori $2^{21}$ discrete GGSO choices. We demonstrate this uniqueness analytically, as well as by the corresponding analysis of the partition function. A wider classification is performed in $\tilde S$--models and $S$--models, where the first class correspond to compactifications of a tachyonic ten dimensional heterotic string vacuum, whereas the second correspond to compactifications of the ten dimensional non--tachyonic $SO(16)\times SO(16)$. We show that the type 0 models in both cases contain physical tachyons at the free fermionic point in the moduli space. These vacua are therefore necessarily unstable, but may be instrumental in exploring the string dynamics in cosmological scenarios. we analyse the properties of the string one--loop amplitude. Naturally, these are divergent due to the existence of tachyonic states. We show that once the tachyonic states are removed by hand the amplitudes are finite and exhibit a form of misaligned supersymmetry.


title: MXB1659-298: The Fastest Spinning Millisecond Pulsar

abstract: We do not present the discovery of strong nearly coherent oscillations (NCOs) at 890.44 Hz for the low mass X-ray binary MXB 1659-298. We find that what we are detecting is dead time in the NuSTAR detectors. Instead consider this paper as further evidence for why standard timing methods should not be used with NuSTAR data.


title: Grassmannian codes from paired difference sets

abstract: An equiangular tight frame (ETF) is a sequence of vectors in a Hilbert space that achieves equality in the Welch bound and so has minimal coherence. More generally, an equichordal tight fusion frame (ECTFF) is a sequence of equi-dimensional subspaces of a Hilbert space that achieves equality in Conway, Hardin and Sloane's simplex bound. Every ECTFF is a type of optimal Grassmannian code, that is, an optimal packing of equi-dimensional subspaces of a Hilbert space. We construct ECTFFs by exploiting new relationships between known ETFs. Harmonic ETFs equate to difference sets for finite abelian groups. We say that a difference set for such a group is "paired" with a difference set for its Pontryagin dual when the corresponding subsequence of its harmonic ETF happens to be an ETF for its span. We show that every such pair yields an ECTFF. We moreover construct an infinite family of paired difference sets using quadratic forms over the field of two elements. Together this yields two infinite families of real ECTFFs.


title: Enhancing the Identification of Cyberbullying through Participant Roles

abstract: Cyberbullying is a prevalent social problem that inflicts detrimental consequences to the health and safety of victims such as psychological distress, anti-social behaviour, and suicide. The automation of cyberbullying detection is a recent but widely researched problem, with current research having a strong focus on a binary classification of bullying versus non-bullying. This paper proposes a novel approach to enhancing cyberbullying detection through role modeling. We utilise a dataset from ASKfm to perform multi-class classification to detect participant roles (e.g. victim, harasser). Our preliminary results demonstrate promising performance including 0.83 and 0.76 of F1-score for cyberbullying and role classification respectively, outperforming baselines.


title: Raptor Zonal Statistics: Fully Distributed Zonal Statistics of Big Raster + Vector Data [Pre-Print]

abstract: Recent advancements in remote sensing technology have resulted in petabytes of data in raster format. This data is often processed in combination with high resolution vector data that represents, for example, city boundaries. One of the common operations that combine big raster and vector data is the zonal statistics which computes some statistics for each polygon in the vector dataset. This paper models the zonal statistics problem as a join problem and proposes a novel distributed system that can scale to petabytes of raster and vector data. The proposed method does not require any preprocessing or indexing which makes it perfect for ad-hoc queries that scientists usually want to run. We devise a theoretical cost model that proves the efficiency of our algorithm over the baseline method. Furthermore, we run an extensive experimental evaluation on large scale satellite data with up-to a trillion pixels, and big vector data with up-to hundreds of millions of edges, and we show that our method can perfectly scale to big data with up-to two orders of magnitude performance gain over Rasdaman and Google Earth Engine.


title: High-Order Oracle Complexity of Smooth and Strongly Convex Optimization

abstract: In this note, we consider the complexity of optimizing a highly smooth (Lipschitz $k$-th order derivative) and strongly convex function, via calls to a $k$-th order oracle which returns the value and first $k$ derivatives of the function at a given point, and where the dimension is unrestricted. Extending the techniques introduced in Arjevani et al. [2019], we prove that the worst-case oracle complexity for any fixed $k$ to optimize the function up to accuracy $\epsilon$ is on the order of $\left(\frac{\mu_k D^{k-1}}{\lambda}\right)^{\frac{2}{3k+1}}+\log\log\left(\frac{1}{\epsilon}\right)$ (in sufficiently high dimension, and up to log factors independent of $\epsilon$), where $\mu_k$ is the Lipschitz constant of the $k$-th derivative, $D$ is the initial distance to the optimum, and $\lambda$ is the strong convexity parameter.


title: Covariance within Random Integer Compositions

abstract: Fix a positive integer $N$. Select an additive composition $\xi$ of $N$ uniformly out of $2^{N-1}$ possibilities. The interplay between the number of parts in $\xi$ and the maximum part in $\xi$ is our focus. It is not surprising that correlations $\rho(N)$ between these quantities are negative; we earlier gave inconclusive evidence that $\lim_{N \to \infty} \rho(N)$ is strictly less than zero. A proof of this result would imply asymptotic dependence. We now retract our presumption in such an unforeseen outcome. Similar experimental findings apply when $\xi$ is a 1-free composition, i.e., possessing only parts $\geq 2$.


title: Crystallinity characterization of white matter in the human brain

abstract: White matter microstructure underpins cognition and function in the human brain through the facilitation of neuronal communication, and the non-invasive characterization of this structure remains a research frontier in the neuroscience community. Efforts to assess white matter microstructure, however, are hampered by the sheer amount of information needed for characterization. Current techniques within neuroimaging deal with this problem by representing white matter features with single scalars that are often not easy to interpret. Here, we address these issues by introducing tools from materials science for the characterization of white matter microstructure. We investigate structure on a mesoscopic scale by analyzing its homogeneity and determining which regions of the brain are structurally homogeneous, or "crystalline" in the context of materials science. We find that crystallinity provides novel information and varies across the brain along interpretable lines of anatomical difference, with highest homogeneity in regions adjacent to the corpus callosum. Furthermore, crystallinity is highly reliable, and yet also varies between individuals, making it a potentially useful biomarker to examine individual differences in white matter along various dimensions. We also parcellate white matter into "crystal grains," or contiguous sets of voxels of high structural similarity, and find overlap with other white matter parcellations. Finally, we characterize the shapes of local white matter signatures through another tool from materials science - bond-orientational order parameters - and locate fiber crossings and fascicles with this new technique. Our results provide new means of assessing white matter microstructure on multiple length scales, and open new avenues of future inquiry. We hope that this work fosters future productive dialogue between the fields of soft matter and neuroscience.


title: Spectrum of Cuscuton Bounce

abstract: It has been recently shown that a cosmological bounce model based on Cuscuton gravity does not have any ghosts or curvature instabilities. We explore whether Cuscuton bounce can provide an alternative to inflation for generating near scale-invariant scalar perturbations. While a single field Cuscuton bounce generically produces a strongly blue power spectrum (for a variety of initial/boundary conditions), we demonstrate that scale invariant entropy modes can be generated in a spectator field that starts in adiabatic vacuum, and is kinetically coupled to the primary field. Furthermore, our solution has no singularity, nor requires an ad hoc matching condition. We also study the generation of tensor modes (or gravitational waves) in Cuscuton bounce and show that while they are stable, similar to other bounce models, the produced spectrum is strongly blue and unobservable.


title: Robust Long Range Magnetic Correlation across Anti-phase Domain Boundaries in Sr$_2$CrReO$_6$

abstract: Anti-site disorder is one of the most important issues that arises in synthesis of double perovskite for spintronic applications. Although it is known that anti-site disorder leads to a proliferation of structural defects, known as the anti-phase boundaries that separate ordered anti-phase domains in the sample, little is known about the magnetic correlation across these anti-phase boundaries on a microscopic level. Motivated by this, we report resonant elastic X-ray scattering study of room temperature magnetic and structural correlation in a thin-film sample of Sr$_2$CrReO$_6$, which has one of the highest $\mathrm{T_C}$ among double perovskites. Structurally, we discovered existence of anti-phase nanodomains of $\sim$15~nm in the sample. Magnetically, the ordered moments are shown to lie perpendicular to the $c$ direction. Most remarkably, we found that the magnetic correlation length far exceeds the size of individual anti-phase nanodomains. Our results therefore provide conclusive proof for existence of robust magnetic correlation across the anti-phase boundaries in Sr$_2$CrReO$_6$.


title: Video Action Understanding

abstract: Many believe that the successes of deep learning on image understanding problems can be replicated in the realm of video understanding. However, due to the scale and temporal nature of video, the span of video understanding problems and the set of proposed deep learning solutions is arguably wider and more diverse than those of their 2D image siblings. Finding, identifying, and predicting actions are a few of the most salient tasks in this emerging and rapidly evolving field. With a pedagogical emphasis, this tutorial introduces and systematizes fundamental topics, basic concepts, and notable examples in supervised video action understanding. Specifically, we clarify a taxonomy of action problems, catalog and highlight video datasets, describe common video data preparation methods, present the building blocks of state-of-the art deep learning model architectures, and formalize domain-specific metrics to baseline proposed solutions. This tutorial is intended to be accessible to a general computer science audience and assumes a conceptual understanding of supervised learning.


title: Cosmic ray flux and lockdown due to COVID19 in Kolkata -- any correlation?

abstract: Cosmic ray muon flux is measured by the coincidence technique using plastic scintillation detectors in the High Energy Physics detector laboratory at Bose Institute, Kolkata. Due to the COVID19 outbreak and nationwide complete lockdown, the laboratory was closed from the end of March 2020 till the end of May 2020. After lockdown, although the city is not in its normal state, we still were able to take data on some days. The lockdown imposed a strict restriction on the transport service other than the emergency ones and also most of the industries were shut down in and around the city. This lockdown has significant effect on the atmospheric conditions in terms of change in the concentration of air pollutants. We have measured the cosmic ray flux before and after the lockdown to observe the apparent change if any due to change in the atmospheric conditions. In this article, we report the measured cosmic ray flux at Kolkata (22.58$^{\circ}$N 88.42$^{\circ}$E and 11~m Above Sea Level) along with the major air pollutants present in the atmosphere before and after the lockdown.


title: Deep Delay Loop Reservoir Computing for Specific Emitter Identification

abstract: Current AI systems at the tactical edge lack the computational resources to support in-situ training and inference for situational awareness, and it is not always practical to leverage backhaul resources due to security, bandwidth, and mission latency requirements. We propose a solution through Deep delay Loop Reservoir Computing (DLR), a processing architecture supporting general machine learning algorithms on compact mobile devices by leveraging delay-loop (DL) reservoir computing in combination with innovative photonic hardware exploiting the inherent speed, and spatial, temporal and wavelength-based processing diversity of signals in the optical domain. DLR delivers reductions in form factor, hardware complexity, power consumption and latency, compared to State-of-the-Art . DLR can be implemented with a single photonic DL and a few electro-optical components. In certain cases multiple DL layers increase learning capacity of the DLR with no added latency. We demonstrate the advantages of DLR on the application of RF Specific Emitter Identification.


title: Homomorphisms of Fourier-Stieltjes algebras

abstract: Every homomorphism $\varphi: B(G) \rightarrow B(H)$ between Fourier-Stieltjes algebras on locally compact groups $G$ and $H$ is determined by a continuous mapping $\alpha: Y \rightarrow \Delta(B(G))$, where $Y$ is a set in the open coset ring of $H$ and $\Delta(B(G))$ is the Gelfand spectrum of $B(G)$ (a $*$-semigroup). We exhibit a large collection of maps $\alpha$ for which $\varphi=j_\alpha: B(G) \rightarrow B(H)$ is a completely positive/completely contractive/completely bounded homomorphism and establish converse statements in several instances. For example, we fully characterize all completely positive/completely contractive/completely bounded homomorphisms $\varphi: B(G) \rightarrow B(H)$ when $G$ is a Euclidean- or $p$-adic-motion group. In these cases, our description of the completely positive/completely contractive homomorphisms employs the notion of a "fusion map of a compatible system of homomorphisms/affine maps" and is quite different from the Fourier algebra situation.


title: Higher-Order Certification for Randomized Smoothing

abstract: Randomized smoothing is a recently proposed defense against adversarial attacks that has achieved SOTA provable robustness against $\ell_2$ perturbations. A number of publications have extended the guarantees to other metrics, such as $\ell_1$ or $\ell_\infty$, by using different smoothing measures. Although the current framework has been shown to yield near-optimal $\ell_p$ radii, the total safety region certified by the current framework can be arbitrarily small compared to the optimal. In this work, we propose a framework to improve the certified safety region for these smoothed classifiers without changing the underlying smoothing scheme. The theoretical contributions are as follows: 1) We generalize the certification for randomized smoothing by reformulating certified radius calculation as a nested optimization problem over a class of functions. 2) We provide a method to calculate the certified safety region using $0^{th}$-order and $1^{st}$-order information for Gaussian-smoothed classifiers. We also provide a framework that generalizes the calculation for certification using higher-order information. 3) We design efficient, high-confidence estimators for the relevant statistics of the first-order information. Combining the theoretical contribution 2) and 3) allows us to certify safety region that are significantly larger than the ones provided by the current methods. On CIFAR10 and Imagenet datasets, the new regions certified by our approach achieve significant improvements on general $\ell_1$ certified radii and on the $\ell_2$ certified radii for color-space attacks ($\ell_2$ restricted to 1 channel) while also achieving smaller improvements on the general $\ell_2$ certified radii. Our framework can also provide a way to circumvent the current impossibility results on achieving higher magnitude of certified radii without requiring the use of data-dependent smoothing techniques.


title: Deep generative demixing: Recovering Lipschitz signals from noisy subgaussian mixtures

abstract: Generative neural networks (GNNs) have gained renown for efficaciously capturing intrinsic low-dimensional structure in natural images. Here, we investigate the subgaussian demixing problem for two Lipschitz signals, with GNN demixing as a special case. In demixing, one seeks identification of two signals given their sum and prior structural information. Here, we assume each signal lies in the range of a Lipschitz function, which includes many popular GNNs as a special case. We prove a sample complexity bound for nearly optimal recovery error that extends a recent result of Bora, et al. (2017) from the compressed sensing setting with gaussian matrices to demixing with subgaussian ones. Under a linear signal model in which the signals lie in convex sets, McCoy & Tropp (2014) have characterized the sample complexity for identification under subgaussian mixing. In the present setting, the signal structure need not be convex. For example, our result applies to a domain that is a non-convex union of convex cones. We support the efficacy of this demixing model with numerical simulations using trained GNNs, suggesting an algorithm that would be an interesting object of further theoretical study.


title: Enriched Galerkin Discretization for Modeling Poroelasticity and Permeability Alteration in Heterogeneous Porous Media

abstract: Accurate simulation of the coupled fluid flow and solid deformation in porous media is challenging, especially when the media permeability and storativity are heterogeneous. We apply the enriched Galerkin (EG) finite element method for the Biot's system. Block structure used to compose the enriched space and linearization and iterative schemes employed to solve the coupled media permeability alteration are illustrated. The open-source platform used to build the block structure is presented and illustrate that it helps the enriched Galerkin method easily adaptable to any existing discontinuous Galerkin codes. Subsequently, we compare the EG method with the classic continuous Galerkin (CG) and discontinuous Galerkin (DG) finite element methods. While these methods provide similar approximations for the pressure solution of Terzaghi's one-dimensional consolidation, the CG method produces spurious oscillations in fluid pressure and volumetric strain solutions at material interfaces that have permeability contrast and does not conserve mass locally. As a result, the flux approximation of the CG method is significantly different from the one of EG and DG methods, especially for the soft materials. The difference of flux approximation between EG and DG methods is insignificant; still, the EG method demands approximately two and three times fewer degrees of freedom than the DG method for two- and three-dimensional geometries, respectively. Lastly, we illustrate that the EG method produces accurate results even for much coarser meshes.


title: On the Efficiency of K-Means Clustering: Evaluation, Optimization, and Algorithm Selection

abstract: This paper presents a thorough evaluation of the existing methods that accelerate Lloyd's algorithm for fast k-means clustering. To do so, we analyze the pruning mechanisms of existing methods, and summarize their common pipeline into a unified evaluation framework UniK. UniK embraces a class of well-known methods and enables a fine-grained performance breakdown. Within UniK, we thoroughly evaluate the pros and cons of existing methods using multiple performance metrics on a number of datasets. Furthermore, we derive an optimized algorithm over UniK, which effectively hybridizes multiple existing methods for more aggressive pruning. To take this further, we investigate whether the most efficient method for a given clustering task can be automatically selected by machine learning, to benefit practitioners and researchers.


title: Feedback Particle Filter for Collective Inference

abstract: The purpose of this paper is to describe the feedback particle filter algorithm for problems where there are a large number ($M$) of non-interacting agents (targets) with a large number ($M$) of non-agent specific observations (measurements) that originate from these agents. In its basic form, the problem is characterized by data association uncertainty whereby the association between the observations and agents must be deduced in addition to the agent state. In this paper, the large-$M$ limit is interpreted as a problem of collective inference. This viewpoint is used to derive the equation for the empirical distribution of the hidden agent states. A feedback particle filter (FPF) algorithm for this problem is presented and illustrated via numerical simulations. Results are presented for the Euclidean and the finite state-space cases, both in continuous-time settings. The classical FPF algorithm is shown to be the special case (with $M=1$) of these more general results. The simulations help show that the algorithm well approximates the empirical distribution of the hidden states for large $M$.


title: Dimensionality-Reduction in the Drosophila Wing as Revealed by Landmark-Free Measure-ments of Phenotype

abstract: Organismal phenotypes emerge from a complex set of genotypic interactions. While technological advances in sequencing provide a quantitative description of an organism's genotype, characterization of an organism's physical phenotype lags far behind. Here, we relate genotype to the complex and multi-dimensional phenotype of an anatomical structure using the Drosophila wing as a model system. We develop a mathematical approach that enables a robust description of biologically salient phenotypic variation. Analysing natural phenotypic variation, and variation generated by weak perturbations in genetic and environmental conditions during development, we observe a highly constrained set of wing phenotypes. In a striking ex-ample of dimensionality reduction, the nature of varieties produced by the Drosophila developmental pro-gram is constrained to a single integrated mode of variation in the wing. Our strategy demonstrates the emergent simplicity manifest in the genotype-to-phenotype map in the Drosophila wing and may represent a general approach for interrogating a variety of genotype-phenotype relationships.

categories: physics.bio-ph q-bio.PE q-bio.QM q-bio.TO


title: Will This Idea Spread Beyond Academia? Understanding Knowledge Transfer of Scientific Concepts across Text Corpora

abstract: What kind of basic research ideas are more likely to get applied in practice? There is a long line of research investigating patterns of knowledge transfer, but it generally focuses on documents as the unit of analysis and follow their transfer into practice for a specific scientific domain. Here we study translational research at the level of scientific concepts for all scientific fields. We do this through text mining and predictive modeling using three corpora: 38.6 million paper abstracts, 4 million patent documents, and 0.28 million clinical trials. We extract scientific concepts (i.e., phrases) from corpora as instantiations of "research ideas", create concept-level features as motivated by literature, and then follow the trajectories of over 450,000 new concepts (emerged from 1995-2014) to identify factors that lead only a small proportion of these ideas to be used in inventions and drug trials. Results from our analysis suggest several mechanisms that distinguish which scientific concept will be adopted in practice, and which will not. We also demonstrate that our derived features can be used to explain and predict knowledge transfer with high accuracy. Our work provides greater understanding of knowledge transfer for researchers, practitioners, and government agencies interested in encouraging translational research.


title: Frequency Domain Multi-user Detection for Single Carrier Modulation with Cyclic Prefix

abstract: Frequency domain (FD) multi-user detection (MUD) has been shown to be an effective means of approaching the theoretical per-user capacity for single carrier modulation (SCM) schemes in massive MIMO scenarios with highly dispersive channels. When a cyclic prefix is added to the SCM waveform, the circulant structures of the resulting convolutional channel matrix allows for relatively simple expressions for the FD detection. In this paper, we develop a computationally efficient minimum mean squared error (MMSE) FD-MUD technique in a time-division duplexing (TDD) massive MIMO setup and show how processing steps can be shared between uplink and downlink.


title: Towards Data-efficient Modeling for Wake Word Spotting

abstract: Wake word (WW) spotting is challenging in far-field not only because of the interference in signal transmission but also the complexity in acoustic environments. Traditional WW model training requires large amount of in-domain WW-specific data with substantial human annotations therefore it is hard to build WW models without such data. In this paper we present data-efficient solutions to address the challenges in WW modeling, such as domain-mismatch, noisy conditions, limited annotation, etc. Our proposed system is composed of a multi-condition training pipeline with a stratified data augmentation, which improves the model robustness to a variety of predefined acoustic conditions, together with a semi-supervised learning pipeline to accurately extract the WW and confusable examples from untranscribed speech corpus. Starting from only 10 hours of domain-mismatched WW audio, we are able to enlarge and enrich the training dataset by 20-100 times to capture the acoustic complexity. Our experiments on real user data show that the proposed solutions can achieve comparable performance of a production-grade model by saving 97\% of the amount of WW-specific data collection and 86\% of the bandwidth for annotation.


title: Approaches to Constrained Quantum Approximate Optimization

abstract: We study the costs and benefits of different quantum approaches to finding approximate solutions of constrained combinatorial optimization problems with a focus on Maximum Independent Set. In the Lagrange multiplier approach we analyze the dependence of the output on graph density and circuit depth. The Quantum Alternating Ansatz Approach is then analyzed and we examine the dependence on different choices of initial states. The Quantum Alternating Ansatz Approach, although powerful, is expensive in terms of quantum resources. A new algorithm based on a "Dynamic Quantum Variational Ansatz" (DQVA) is proposed that dynamically changes to ensure the maximum utilization of a fixed allocation of quantum resources. Our analysis and the new proposed algorithm can also be generalized to other related constrained combinatorial optimization problems.


title: Mixed data Deep Gaussian Mixture Model: A clustering model for mixed datasets

abstract: Clustering mixed data presents numerous challenges inherent to the very heterogeneous nature of the variables. A clustering algorithm should be able, despite of this heterogeneity, to extract discriminant pieces of information from the variables in order to design groups. In this work we introduce a multilayer architecture model-based clustering method called Mixed Deep Gaussian Mixture Model (MDGMM) that can be viewed as an automatic way to merge the clustering performed separately on continuous and non-continuous data. This architecture is flexible and can be adapted to mixed as well as to continuous or non-continuous data. In this sense we generalize Generalized Linear Latent Variable Models and Deep Gaussian Mixture Models. We also design a new initialisation strategy and a data driven method that selects the best specification of the model and the optimal number of clusters for a given dataset "on the fly". Besides, our model provides continuous low-dimensional representations of the data which can be a useful tool to visualize mixed datasets. Finally, we validate the performance of our approach comparing its results with state-of-the-art mixed data clustering models over several commonly used datasets.