title: On locally conformally K\"ahler threefolds with algebraic dimension two

abstract: The paper is part of an attempt of understanding non-K\"ahler threefolds. We start by looking at compact complex non-K\"ahler threefolds with algebraic dimension two and admitting locally conformally K\"ahler metrics. Under certain assumptions, we prove that they are blown-up quasi-bundles over a projective surface.


title: Notes on Hamiltonian threshold and chain graphs

abstract: We revisit results obtained in [F. Harary, U. Peled, Hamiltonian threshold graphs, Discrete Appl.~Math., 16 (1987), 11--15], where several necessary and necessary and sufficient conditions for a connected threshold graph to be Hamiltonian were obtained. We present these results in new forms, now stated in terms of structural parameters that uniquely define the threshold graph and we extend them to chain graphs. We also identify the chain graph with minimum number of Hamilton cycles within the class of Hamiltonian chain graphs of a given order.


title: On the lightlike Lorentz violation

abstract: We consider the lightlike spontaneous Lorentz invariance violation (SLIV) appearing through the zero "length-fixing" constraint put on a gauge vector field, $A_{\mu }A^{\mu }=0$, and discuss its physical consequences in the framework of a conventional QED and beyond. Again, as in the timelike and spacelike SLIV cases, $A_{\mu }A^{\mu }=\pm M_{A}^{2}$ ($M_{A}$ is a scale of SLIV), while this constraint leads to an emergence of the Nambu-Goldstone modes collected in physical photon, the SLIV itself is still left unobservable unless gauge invariance in the theory is broken. At the same time, a crucial difference with the two former cases is that the asymmetrical vacuum corresponding to the lightlike Lorentz violation appears infinitely degenerated with all other vacua including the symmetrical one. We show that this degeneracy can be lifted out by introducing an extra gauge vector field being sterile with respect to an ordinary matter, though having some potential couplings with the basic $A_{\mu }$ field. A slight mixing of them makes the underlying gauge invariance to be partially broken due to which physical Lorentz invariance occurs broken as well. This may cause a variety of the Lorentz violating processes some of which are briefly discussed.


title: The singular Weinstein conjecture

abstract: In this article, we investigate Reeb dynamics on $b^m$-contact manifolds, previously introduced in [MiO], which are contact away from a hypersurface $Z$ but satisfy certain transversality conditions on $Z$. The study of these contact structures is motivated by that of contact manifolds with boundary. The search of periodic Reeb orbits on those manifolds thereby starts with a generalization of the well-known Weinstein conjecture. Contrary to the initial expectations, examples of compact $b^m$-contact manifolds without periodic Reeb orbits outside $Z$ are provided. Furthermore, we prove that in dimension $3$, there are always infinitely many periodic orbits on the critical set if it is compact. We prove that traps for the $b^m$-Reeb flow exist in any dimension. This investigation goes hand-in-hand with the Weinstein conjecture on non-compact manifolds having compact ends of convex type. In particular, we extend Hofer's arguments to open overtwisted contact manifolds that are $\mathbb R^+$-invariant in the open ends, obtaining as a corollary the existence of periodic $b^m$-Reeb orbits away from the critical set. The study of $b^m$-Reeb dynamics is motivated by well-known problems in fluid dynamics and celestial mechanics, where those geometric structures naturally appear. In particular, we prove that the dynamics on positive energy level-sets in the restricted planar circular three-body problem are described by the Reeb vector field of a $b^3$-contact form that admits an infinite number of periodic orbits at the critical set.


title: Towards Particle-Resolved Accuracy in Euler-Lagrange Simulations of Multiphase Flow Using Machine Learning and Pairwise Interaction Extended Point-particle (PIEP) Approximation

abstract: This study presents two different machine learning approaches for the modeling of hydrodynamic force on particles in a particle-laden multiphase flow. Results from particle-resolved direct numerical simulations (PR-DNS) of flow over a random array of stationary particles for eight combinations of particle Reynolds number ($Re$) and volume fraction ($\phi$) are used in the development of the models. The first approach follows a two step process. In the first flow prediction step, the perturbation flow due to a particle is obtained as an axisymmetric superposable wake using linear regression. In the second force prediction step, the force on a particle is evaluated in terms of the perturbation flow induced by all its neighbors using the generalized Fax\'en form of the force expression. In the second approach, the force data on all the particles from the PR-DNS simulations is used to develop an artificial neural network (ANN) model for direct prediction of force on a particle. Due to the unavoidable limitation on the number of fully resolved particles in the PR-DNS simulations, direct force prediction with the ANN model tends to over-fit the data and performs poorly in the prediction of test data. In contrast, due to the millions of grid points used in the PR-DNS simulations, accurate flow prediction is possible, which then allows accurate prediction of particle force. This hybridization of multiphase physics and machine learning is particularly important, since it blends the strength of each, and the resulting pairwise interaction extended point-particle (PIEP) model cannot be developed by either physics or machine learning alone.


title: On a Generalized Lam\'e-Navier system in $\mathbb{R}^3$

abstract: This paper is devoted to a fundamental system of equations in Linear Elasticity Theory: the famous Lam\'e-Navier system. The Clifford algebra language allows us to rewrite this system in terms of the euclidean Dirac operator, which at the same time suggests a very natural generalization involving the so-called structural sets. We are interested in finding some structures in the solutions of these generalized Lam\'e-Navier systems. Using MATLAB we also implement algorithms to compute with such partial differential operators as well as to verify some theoretical results obtained in the paper.


title: Toward Large-Scale Autonomous Monitoring and Sensing of Underwater Pollutants

abstract: Marine pollution is a growing worldwide concern, affecting health of marine ecosystems, human health, climate change, and weather patterns. To reduce underwater pollution, it is critical to have access to accurate information about the extent of marine pollutants as otherwise appropriate countermeasures and cleaning measures cannot be chosen. Currently such information is difficult to acquire as existing monitoring solutions are highly laborious or costly, limited to specific pollutants, and have limited spatial and temporal resolution. In this article, we present a research vision of large-scale autonomous marine pollution monitoring that uses coordinated groups of autonomous underwater vehicles (AUV)s to monitor extent and characteristics of marine pollutants. We highlight key requirements and reference technologies to establish a research roadmap for realizing this vision. We also address the feasibility of our vision, carrying out controlled experiments that address classification of pollutants and collaborative underwater processing, two key research challenges for our vision.


title: Ensemble Transfer Learning for the Prediction of Anti-Cancer Drug Response

abstract: Transfer learning has been shown to be effective in many applications in which training data for the target problem are limited but data for a related (source) problem are abundant. In this paper, we apply transfer learning to the prediction of anti-cancer drug response. Previous transfer learning studies for drug response prediction focused on building models that predict the response of tumor cells to a specific drug treatment. We target the more challenging task of building general prediction models that can make predictions for both new tumor cells and new drugs. We apply the classic transfer learning framework that trains a prediction model on the source dataset and refines it on the target dataset, and extends the framework through ensemble. The ensemble transfer learning pipeline is implemented using LightGBM and two deep neural network (DNN) models with different architectures. Uniquely, we investigate its power for three application settings including drug repurposing, precision oncology, and new drug development, through different data partition schemes in cross-validation. We test the proposed ensemble transfer learning on benchmark in vitro drug screening datasets, taking one dataset as the source domain and another dataset as the target domain. The analysis results demonstrate the benefit of applying ensemble transfer learning for predicting anti-cancer drug response in all three applications with both LightGBM and DNN models. Compared between the different prediction models, a DNN model with two subnetworks for the inputs of tumor features and drug features separately outperforms LightGBM and the other DNN model that concatenates tumor features and drug features for input in the drug repurposing and precision oncology applications. In the more challenging application of new drug development, LightGBM performs better than the other two DNN models.


title: Cytosolic flow induced symmetry breaking in a conceptual polarity model

abstract: Important cellular processes, such as cell motility and cell division, are coordinated by cell polarity, which is determined by the non-uniform distribution of certain proteins. Such protein patterns form via an interplay of protein reactions and protein transport. Since Turing's seminal work, the formation of protein patterns resulting from the interplay between reactions and diffusive transport has been widely studied. Over the last few years, increasing evidence shows that also advective transport, resulting from cytosolic and cortical flows, is present in many cells. However, it remains unclear how and whether these flows contribute to protein-pattern formation. To address this question, we use a minimal model that conserves the total protein mass to characterize the effects of cytosolic flow on pattern formation. Combining a linear stability analysis with numerical simulations, we find that membrane-bound protein patterns propagate against the direction of cytoplasmic flow with a speed that is maximal for intermediate flow speed. We show that the mechanism underlying this pattern propagation relies on a higher protein influx on the upstream side of the pattern compared to the downstream side. Furthermore, we find that cytosolic flow can change the membrane pattern qualitatively from a peak pattern to a mesa pattern. Finally, our study shows that a non-uniform flow profile can induce pattern formation by triggering a regional lateral instability.


title: Model independent measurement of the growth rate from the consistency relations of the LSS

abstract: The Consistency Relations for the Large Scale Structure provide a link between the amplitude of Baryonic Acoustic Oscillations in the squeezed bispectrum (BS) and in the power spectrum (PS). This relation depends on the large scale bias of the considered tracer, $b_\alpha$, and on the growth rate of structures, $f$. Remarkably, originating from basic symmetry principles, this relation is exact and independent on the underlying cosmological model. By analysing data from large volume simulations, both for dark matter and for haloes, we illustrate how BS and PS measurements can be used to extract $b_\alpha$ and $f$ without the need of any theoretical approximation scheme for the computation of the BS and the PS. We show that, combining measurements of the squeezed BS with the quadrupole to monopole ratios for the PS at large scales can successfully break the $b_\alpha -f$ degeneracy. We forecast that this method, applied to a Euclid-like survey, will be able to measure bias, and then the growth rate, at better than $10\%$ level, with no extra assumption.


title: Search for heavy neutral lepton production in $K^+$ decays to positrons

abstract: A search for heavy neutral lepton ($N$) production in $K^+\to e^+N$ decays using the data sample collected by the NA62 experiment at CERN in 2017--2018 is reported. Upper limits of the extended neutrino mixing matrix element $|U_{e4}|^2$ are established at the level of $10^{-9}$ over most of the accessible heavy neutral lepton mass range 144--462 MeV/$c^2$, with the assumption that the lifetime exceeds 50 ns. These limits improve significantly upon those of previous production and decay searches. The $|U_{e4}|^2$ range favoured by Big Bang Nucleosynthesis is excluded up to a mass of about 340 MeV/$c^2$.


title: Nonsingular Black Holes From Charged Dust Collapse: A Concrete Mechanism to Evade Interior Singularities in General Relativity

abstract: In this essay we examine the gravitational collapse of a nonrelativistic charged perfect fluid interacting with a dark energy component. Given a simple factor for the energy transfer, we obtain a nonsingular interior solution which naturally matches the Reissner-Nordstr\"om-de Sitter exterior geometry. We also show that the interacting parameter is proportional to the overall charge of the final black hole thus formed. For the case of quasi-extremal configurations, we propose a statistical model for the entropy of the collapsed matter. This entropy extends Bekenstein's geometrical entropy by an additive constant proportional to the area of the extremal black hole.


title: Increasing Domain Infill Asymptotics for Stochastic Differential Equations Driven by Fractional Brownian Motion

abstract: Although statistical inference in stochastic differential equations (SDEs) driven by Wiener process has received significant attention in the literature, inference in those driven by fractional Brownian motion seem to have seen much less development in comparison, despite their importance in modeling long range dependence. In this article, we consider both classical and Bayesian inference in such fractional Brownian motion based SDEs. In particular, we consider asymptotic inference for two parameters in this regard; a multiplicative parameter associated with the drift function, and the so-called "Hurst parameter" of the fractional Brownian motion, when the time domain tends to infinity. For unknown Hurst parameter, the likelihood does not lend itself amenable to the popular Girsanov form, rendering usual asymptotic development difficult. As such, we develop increasing domain infill asymptotic theory, by discretizing the SDE. In this setup, we establish consistency and asymptotic normality of the maximum likelihood estimators, as well as consistency and asymptotic normality of the Bayesian posterior distributions. However, classical or Bayesian asymptotic normality with respect to the Hurst parameter could not be established. We supplement our theoretical investigations with simulation studies in a non-asymptotic setup, prescribing suitable methodologies for classical and Bayesian analyses of SDEs driven by fractional Brownian motion. Applications to a real, close price data, along with comparison with standard SDE driven by Wiener process, is also considered. As expected, it turned out that our Bayesian fractional SDE triumphed over the other model and methods, in both simulated and real data applications.


title: Stratification of polymer mixtures in drying droplets: hydrodynamics and diffusion

abstract: We study the evaporation-induced stratification of a mixture of short and long polymer chains in a drying droplet using molecular simulations. We systematically investigate the effects of hydrodynamic interactions (HI) on this process by comparing hybrid simulations accounting for HI between polymers through the multiparticle collision dynamics technique with free-draining Langevin dynamics simulations neglecting the same. We find that the dried supraparticle morphologies are homogeneous when HI are included but are stratified in core--shell structures (with the short polymers forming the shell) when HI are neglected. The simulation methodology unambiguously attributes this difference to the treatment of the solvent in the two models. We rationalize the presence (or absence) of stratification by measuring phenomenological multicomponent diffusion coefficients for the polymer mixtures. The diffusion coefficients show the importance of not only solvent backflow but also HI between polymers in controlling the dried supraparticle morphology.


title: Near-ideal spontaneous photon sources in silicon quantum photonics

abstract: While integrated photonics is a robust platform for quantum information processing, architectures for photonic quantum computing place stringent demands on high quality information carriers. Sources of single photons that are highly indistinguishable and pure, that are either near-deterministic or heralded with high efficiency, and that are suitable for mass-manufacture, have been elusive. Here, we demonstrate on-chip photon sources that simultaneously meet each of these requirements. Our photon sources are fabricated in silicon using mature processes, and exploit a novel dual-mode pump-delayed excitation scheme to engineer the emission of spectrally pure photon pairs through intermodal spontaneous four-wave mixing in low-loss spiralled multi-mode waveguides. We simultaneously measure a spectral purity of $0.9904 \pm 0.0006$, a mutual indistinguishably of $0.987 \pm 0.002$, and $>90\%$ intrinsic heralding efficiency. We measure on-chip quantum interference with a visibility of $0.96 \pm 0.02$ between heralded photons from different sources. These results represent a decisive step for scaling quantum information processing in integrated photonics.


title: A Ballistic Graphene Cooper Pair Splitter

abstract: We report an experimental study of Cooper pair splitting in an encapsulated graphene based multiterminal junction in the ballistic transport regime. Our device consists of two transverse junctions, namely the superconductor/graphene/superconductor and the normal metal/graphene/normal metal junctions. In this case, the electronic transport through one junction can be tuned by an applied bias along the other. We observe clear signatures of Cooper pair splitting in the local as well as nonlocal electronic transport measurements. Our experimental data can be very well described by using a modified Octavio-Tinkham-Blonder-Klapwijk model and a three-terminal beam splitter model.


title: Controlling the electrical and magnetic ground states by doping in the complete phase diagram of titanate Eu1-xLaxTiO3 thin films

abstract: EuTiO3, a band insulator, and LaTiO3, a Mott insulator, are both antiferromagnetic with transition temperatures ~ 5.5 K and ~ 160 K, respectively. Here, we report the synthesis of Eu1-xLaxTiO3 thin films with x = 0 to 1 by oxide molecular beam epitaxy. The films in the full range have high crystalline quality and show no phase segregation, allowing us carry out transport measurements to study their electrical and magnetic properties. From x = 0.03 to 0.95, Eu1-xLaxTiO3 films show conduction by electrons as charge carriers, with differences in carrier densities and mobilities, contrary to the insulating nature of pure EuTiO3 and LaTiO3. Following a rich phase diagram, the magnetic ground states of the films vary with increasing La-doping level, changing Eu1-xLaxTiO3 from an antiferromagnetic insulator to an antiferromagnetic metal, a ferromagnetic metal, a paramagnetic metal, and back to an antiferromagnetic insulator. These emergent properties reflect the evolutions of the band structure, mainly at the Ti t2g bands near the Fermi level, when Eu2+ are gradually replaced by La3+. This work sheds light on this method for designing the electrical and magnetic properties in strongly-correlated oxides and completes the phase diagram of the titanate Eu1-xLaxTiO3.


title: Preorders on Subharmonic Functions and Measures with Applications to the Distribution of Zeros of Holomorphic Functions

abstract: Let $X$ be a class of extended numerical functions on a domain $D$ of $d$-dimensional Euclidean space $\mathbb R^d$, $H\subset X$. Given $u,M\in X$, we write $u\prec_H M$ if there is a function $h\in H$ such that $u+h\leq M$ on $D$. We consider this special preorder $\prec_H$ for a pair of subharmonic unctions $u, M$ on $D$ in cases where $H$ is the space of all harmonic functions on $D$ or $H$ is the convex cone of all subharmonic functions $h \not\equiv -\infty$ on $D$. Main results are dual equivalent forms for this preorder $\prec_H$ in terms of balayage processes for Riesz measures of subharmonic functions $u$ and $M$, for Jensen and Arens-Singer (representing) measures, for potentials of these measures, and for special test functions generated by subharmonic functions on complements $D\setminus S$ of non-empty precompact subsets $S\Subset D$. Applications to holomorphic functions $f$ on a domain $D\subset \mathbb C^n$ relate to the distribution of zero sets of functions $f$ under upper restrictions $|f|\leq \exp M$ on $D$. If a domain $D\subset \mathbb C$ is a finitely connected domain with non-empty exterior or a simply connected domain with two different points on the boundary of $D$, then our conditions for the distribution of zeros of $f\neq 0$ with $|f|\leq \exp M$ on $D$ are both necessary and sufficient.


title: Instrumental Variables with Treatment-Induced Selection: Exact Bias Results

abstract: Instrumental variables (IV) estimation suffers selection bias when the analysis conditions on the treatment. Judea Pearl's early graphical definition of instrumental variables explicitly prohibited conditioning on the treatment. Nonetheless, the practice remains common. In this paper, we derive exact analytic expressions for IV selection bias across a range of data-generating models, and for various selection-inducing procedures. We present four sets of results for linear models. First, IV selection bias depends on the conditioning procedure (covariate adjustment vs. sample truncation). Second, IV selection bias due to covariate adjustment is the limiting case of IV selection bias due to sample truncation. Third, in certain models, the IV and OLS estimators under selection bound the true causal effect in large samples. Fourth, we characterize situations where IV remains preferred to OLS despite selection on the treatment. These results broaden the notion of IV selection bias beyond sample truncation, replace prior simulation findings with exact analytic formulas, and enable formal sensitivity analyses.


title: Angular dependence of the Wigner time delay upon tunnel ionization of $H_{2}$

abstract: More than 100 years after its discovery and its explanation in the energy domain, the duration of the photoelectric effect is still heavily studied. The emission time of a photoelectron can be quantified by the Wigner time delay. Experiments addressing this time delay for single-photon ionization became feasible during the last 10 years. A missing piece, which has not been studied, so far, is the Wigner time delay for strong-field ionization of molecules. Here we show experimental data on the Wigner time delay for tunnel ionization of $H_{2}$ molecules and demonstrate its dependence on the emission direction of the electron with respect to the molecular axis. We find, that the observed changes in the Wigner time delay can be quantitatively explained by elongated/shortened travel paths of the electrons that are due to spatial shifts of the electron's birth position after tunneling. This introduces an intuitive perspective towards the Wigner time delay in strong-field ionization.


title: A Frameless 2-Coloring of the Plane Lattice

abstract: A picture frame in two dimensions is a rectangular array of symbols, with at least two rows and columns, where the first and last rows are identical, and the first and last columns are identical. If a coloring of the plane lattice has no picture frames, we call it frameless. In this note we show how to create a simple 2-coloring of the plane lattice that is frameless.


title: Development of a Confined Circular-cum-parallel Ribbon Flare and Associated Pre-flare Activity

abstract: We study a complex GOES M1.1 circular ribbon flare and related pre-flare activity on 26 January 2015 [SOL26-01-2015] in solar active region NOAA 12268. This flare activity was observed by the AIA on board SDO and the RHESSI. The examination of photospheric magnetograms during the extended period, prior to the event, suggests the successive development of a so-called 'anemone' type magnetic configuration. NLFFF extrapolation reveals a fan-spine magnetic configuration with the presence of a coronal null-point. We found that the pre-flare activity in the active region starts ~15 min prior to the main flare in the form of localized bright patches at two locations. A comparison of locations and spatial structures of the pre-flare activity with magnetic configuration of the corresponding region suggests onset of magnetic reconnection at the null-point along with the low-atmosphere magnetic reconnection caused by the emergence and the cancellation of the magnetic flux. The main flare of M1.1 class is characterized by the formation of a well-developed circular ribbon along with a region of remote brightening. Remarkably, a set of relatively compact parallel ribbons formed inside the periphery of the circular ribbon which developed lateral to the brightest part of the circular ribbon. During the peak phase of the flare, a coronal jet is observed at the north-east edge of the circular ribbon which suggests interchange reconnection between large-scale field lines and low-lying closed field lines. Our investigation suggests a combination of two distinct processes in which ongoing pre-flare null-point reconnection gets further intensified as the confined eruption along with jet activity proceeded from within the circular ribbon region which results to the formation of inner parallel ribbons and corresponding post-reconnection arcade.


title: GEV Beamforming Supported by DOA-based Masks Generated on Pairs of Microphones

abstract: Distant speech processing is a challenging task, especially when dealing with the cocktail party effect. Sound source separation is thus often required as a preprocessing step prior to speech recognition to improve the signal to distortion ratio (SDR). Recently, a combination of beamforming and speech separation networks have been proposed to improve the target source quality in the direction of arrival of interest. However, with this type of approach, the neural network needs to be trained in advance for a specific microphone array geometry, which limits versatility when adding/removing microphones, or changing the shape of the array. The solution presented in this paper is to train a neural network on pairs of microphones with different spacing and acoustic environmental conditions, and then use this network to estimate a time-frequency mask from all the pairs of microphones forming the array with an arbitrary shape. Using this mask, the target and noise covariance matrices can be estimated, and then used to perform generalized eigenvalue (GEV) beamforming. Results show that the proposed approach improves the SDR from 4.78 dB to 7.69 dB on average, for various microphone array geometries that correspond to commercially available hardware.


title: Faster Deterministic All Pairs Shortest Paths in Congest Model

abstract: We present a new deterministic algorithm for distributed weighted all pairs shortest paths (APSP) in both undirected and directed graphs. Our algorithm runs in $\tilde{O}(n^{4/3})$ rounds in the Congest models on graphs with arbitrary edge weights, and it improves on the previous $\tilde{O}(n^{3/2})$ bound of Agarwal et al. [ARKP18]. The main components of our new algorithm are a new faster technique for constructing blocker set deterministically and a new pipelined method for deterministically propagating distance values from source nodes to the blocker set nodes in the network. Both of these techniques have potential applications to other distributed algorithms. Our new deterministic algorithm for computing blocker set adapts the NC approximate hypergraph set cover algorithm in [BRS94] to the distributed construction of a blocker set. It follows the two-step process of first designing a randomized algorithm that uses only pairwise independence, and then derandomizes this algorithm using a sample space of linear size. This algorithm runs in almost the same number of rounds as the initial step in our APSP algorithm that computes $h$-hops shortest paths, and significantly improves on the deterministic blocker set algorithms in [ARKP18, AR19] by removing an additional $n\cdot |Q|$ term in the round bound, where Q is the blocker set. The other new component in our APSP algorithm is a deterministic pipelined approach to propagate distance values from source nodes to blocker nodes. We use a simple natural round-robin method for this step, and we show using a suitable progress measure that it achieve the $\tilde{O}(n^{4/3})$ bound on the number of rounds. It appears that the standard deterministic methods for efficiently broadcasting multiple values, and for sending or receiving messages using the routing schedule in [HPDG+19,LSP19] do not apply to this setting.


title: Towards portable muography with small-area, gas-tight glass Resistive Plate Chambers

abstract: Imaging techniques that use atmospheric muons, collectively named under the neologism "muography", have seen a tremendous growth in recent times, mainly due to their diverse range of applications. The most well-known ones include but are not limited to: volcanology, archaeology, civil engineering, nuclear reactor monitoring, nuclear waste characterization, underground mapping, etc. These methods are based on the attenuation or deviation of muons to image large and/or dense objects where conventional techniques cannot work or their use becomes challenging. In this context, we have constructed a muography telescope based on "mini glass-RPC planes" following a design similar to the glass-RPC detectors developed by the CALICE Collaboration and used by the TOMUVOL experiment in the context of volcano radiography, but with smaller active area (16 $\times$ 16 cm$^{2}$). The compact size makes it an attractive choice with respect to other detectors previously employed for imaging on similar scales. An important innovation in this design is that the detectors are sealed. This makes the detector more portable and solves the usual safety and logistic issues for gas detectors operated underground and/or inside small rooms. This paper provides an overview on our guiding principles, the detector development and our operational experiences. Drawing on the lessons learnt from the first prototype, we also discuss our future direction for an improved second prototype, focusing primarily on a recently adopted serigraphy technique for the resistive coating of the glass plates.


title: $\alpha$, $\beta$-expansions of the Riordan matrices of the associated subgroup

abstract: We consider the group of the matrices $\left( 1,g\left( x \right) \right)$ isomorphic to the group of formal power series $g\left( x \right)=x+{{g}_{2}}{{x}^{2}}+...$ under composition: $\left( 1,{{g}_{2}}\left( x \right) \right)\left( 1,{{g}_{1}}\left( x \right) \right)=\left( 1,{{g}_{1}}\left( {{g}_{2}}\left( x \right) \right) \right)$. Denote $P_{k}^{\alpha }=\left( 1,x{{\left( 1-k\alpha {{x}^{k}} \right)}^{{-1}/{k}\;}} \right)$. Matrix $\left( 1,g\left( x \right) \right)$is decomposed into an infinite product of the matrices $P_{k}^{\alpha }$ with suitable exponents in two ways: to left-handed and right-handed products with respect to the matrix $P_{1}^{{{\alpha }_{1}}={{\beta }_{1}}}$: $\left( 1,g\left( x \right) \right)=...P_{k}^{{{\alpha }_{k}}}...P_{2}^{{{\alpha }_{2}}}P_{1}^{{{\alpha }_{1}}}=P_{1}^{{{\beta }_{1}}}P_{2}^{{{\beta }_{2}}}...P_{k}^{{{\beta }_{k}}}...$. We obtain two formulas expressing the coefficients of the series ${{\left( {g\left( x \right)}/{x}\; \right)}^{z}}$ in terms of the expansion coefficients ${{\alpha }_{i}}$, ${{\beta }_{i}}$ and introduce two one-parameter families of series $g_{\alpha }^{\left( t \right)}\left( x \right)$ and $g_{\beta }^{\left( t \right)}\left( x \right)$ associated with these expansions.


title: Transferability of neural network potentials for varying stoichiometry: phonons and thermal conductivity of Mn$_x$Ge$_y$ compounds

abstract: Germanium manganese compounds exhibit a variety of stable and metastable phases with different stoichiometry. These materials entail interesting electronic, magnetic and thermal properties both in their bulk form and as heterostructures. Here we develop and validate a transferable machine learning potential, based on the high-dimensional neural network formalism, to enable the study of Mn$_x$Ge$_y$ materials over a wide range of compositions. We show that a neural network potential fitted on a minimal training set reproduces successfully the structural and vibrational properties and the thermal conductivity of systems with different local chemical environments, and it can be used to predict phononic effects in nanoscale heterostructures.


title: Minimal models of invasion and clonal selection in cancer

abstract: In this thesis we develop minimal models of the relationship between motility, growth, and evolution of cancer cells. We utilise simple simulations of a population of individual cells in space to examine how changes in mechanical properties of invasive cells and their surroundings can affect the speed of cell migration. We also find that the growth rate of large lesions depends weakly on the migration speed of escaping cells, and has stronger and more complex dependencies on the rates of other stochastic processes in the model, namely the rate at which cells transition to being motile and the reverse rate at which cells cease to be motile. To examine how the rates of growth and evolution of an ensemble of cancerous lesions depends on their geometry and underlying fitness landscape, we develop an analytical framework in which the spatial structure is coarse grained and the cancer treated as a continuously growing system with stochastic migration events. Both approaches conclude that the whole ensemble can undergo migration-driven exponential growth regardless of the dependence of size on time of individual lesions, and that the relationship between growth rate and rate of migration is determined by the geometrical constraints of individual lesions. We also find that linear fitness landscapes result in faster-than-exponential growth of the ensemble, and we can determine the expected number of driver mutations present in several important cases of the model. Finally, we study data from a clinical study of the effectiveness of a new low-dose combined chemotherapy. This enables us to test some important hypotheses about the growth rate of pancreatic cancers and the speed with which evolution occurs in reality. Despite this, we find that the frequency of resistant mutants is far too high to be explained without resorting to novel mechanisms of cross-resistance to multiple drugs.


title: A new family of infinitely braided Thompson's groups

abstract: We present a generalization of the Dehornoy-Brin braided Thompson group $BV_2$ that uses recursive braids. Our new groups are denoted by $BV_{n,r}(H)$, for all $n\geq 2,r\geq 1$ and $H \leq \mathcal{B}_n$, where $\mathcal{B}_n$ is the braid group on $n$ strands. We give a new approach to deal with braided Thompson groups by using strand diagrams. We show that $BV_{n,r}(H)$ is finitely generated if $H$ is finitely generated.


title: Novel flavour-changing neutral currents in the top quark sector

abstract: We demonstrate that flavour-changing neutral currents in the top sector, mediated by leptophilic scalars at the electroweak scale, can easily arise in scenarios of new physics, and in particular in composite Higgs models. We moreover show that such interactions are poorly constrained by current experiments, while they can be searched for at the LHC in rare top decays and, more generally, in the channels $pp\to tS(S)+j$, with $S\to\ell^+\ell^-$. We provide dedicated analyses in this respect, obtaining that cut-off scales as large as $\Lambda\sim$ 90 TeV can be probed with an integrated luminosity of $\mathcal{L} = 150$ fb$^{-1}$.


title: Continuous LWE

abstract: We introduce a continuous analogue of the Learning with Errors (LWE) problem, which we name CLWE. We give a polynomial-time quantum reduction from worst-case lattice problems to CLWE, showing that CLWE enjoys similar hardness guarantees to those of LWE. Alternatively, our result can also be seen as opening new avenues of (quantum) attacks on lattice problems. Our work resolves an open problem regarding the computational complexity of learning mixtures of Gaussians without separability assumptions (Diakonikolas 2016, Moitra 2018). As an additional motivation, (a slight variant of) CLWE was considered in the context of robust machine learning (Diakonikolas et al.~FOCS 2017), where hardness in the statistical query (SQ) model was shown; our work addresses the open question regarding its computational hardness (Bubeck et al.~ICML 2019).


title: Spatial analysis and prediction of COVID-19 spread in South Africa after lockdown

abstract: What is the impact of COVID-19 on South Africa? This paper envisages assisting researchers and decision-makers in battling the COVID-19 pandemic focusing on South Africa. This paper focuses on the spread of the disease by applying heatmap retrieval of hotspot areas and spatial analysis is carried out using the Moran index. For capturing spatial autocorrelation between the provinces of South Africa, the adjacent, as well as the geographical distance measures, are used as a weight matrix for both absolute and relative counts. Furthermore, generalized logistic growth curve modeling is used for the prediction of the COVID-19 spread. We expect this data-driven modeling to provide some insights into hotspot identification and timeous action controlling the spread of the virus.


title: Dissipative polarization domain walls in a passive driven Kerr resonator

abstract: Using a passive driven nonlinear optical fiber ring resonator, we report the experimental realization of dissipative polarization domain walls. The domain walls arise through a symmetry breaking bifurcation and consist of temporally localized structures where the amplitudes of the two polarization modes of the resonator interchange, segregating domains of orthogonal polarization states. We show that dissipative polarization domain walls can persist in the resonator without changing shape. We also demonstrate on-demand excitation, as well as pinning of domain walls at specific positions for arbitrary long times. Our results could prove useful for the analog simulation of ubiquitous domain-wall related phenomena, and pave the way to an all-optical buffer adapted to the transmission of topological bits.


title: Animal social networks: an introduction for complex systems scientists

abstract: Many animals live in societies where individuals frequently interact socially with each other. The social structures of these systems can be studied in depth by means of network analysis. A large number of studies on animal social networks in many species have in recent years been carried out in the biological research field of animal behaviour and have provided new insights into behaviour, ecology, and social evolution. This line of research is currently not so well connected to the field of complex systems as could be expected. The purpose of this paper is to provide an introduction to animal social networks for complex systems scientists and highlight areas of synergy. We believe that an increased integration of animal social networks with the interdisciplinary field of complex systems and networks would be beneficial for various reasons. Increased collaboration between researchers in this field and biologists studying animal social systems could be valuable in solving challenges that are of importance to animal social network research. Furthermore, animal social networks provide the opportunity to investigate hypotheses about complex systems across a range of natural real-world social systems. In this paper, we describe what animal social networks are and main research themes where they are studied; we give an overview of the methods commonly used to study animal social networks; we highlight challenges in the study of animal social networks where complex systems expertise may be particularly valuable; and we consider aspects of animal social networks that may be of particular interest to complex systems researchers. We hope that this will help to facilitate further interdisciplinary collaborations involving animal social networks, and further integration of these networks into the field of complex systems.


title: Asymmetric scale functions for $t$-digests

abstract: The $t$-digest is a data structure that can be queried for approximate quantiles, with greater accuracy near the minimum and maximum of the distribution. We develop a $t$-digest variant with accuracy asymmetric about the median, thereby making possible alternative tradeoffs between computational resources and accuracy which may be of particular interest for distributions with significant skew. After establishing some theoretical properties of scale functions for $t$-digests, we show that a tangent line construction on the familiar scale functions preserves the crucial properties that allow $t$-digests to operate online and be mergeable. We conclude with an empirical study demonstrating the asymmetric variant preserves accuracy on one side of the distribution with a much smaller memory footprint.


title: Generalised regression estimation given imperfectly matched auxiliary data

abstract: Generalised regression estimation allows one to make use of available auxiliary information in survey sampling. We develop three types of generalised regression estimator when the auxiliary data cannot be matched perfectly to the sample units, so that the standard estimator is inapplicable. The inference remains design-based. Consistency of the proposed estimators is either given by construction or else can be tested given the observed sample and links. Mean square errors can be estimated. A simulation study is used to explore the potentials of the proposed estimators.


title: Restricted classes of Veronese type ideals and algebras

abstract: We study ideals which are generated by monomials of degree $d$ in the polynomial ring in $n$ variables and which satisfy certain numerical side conditions regarding their exponents. Typical examples of such ideals are the ideals of Veronese type, squarefree Veronese ideals or $t$-spread Veronese ideals. In this paper we focus on $c$-bounded $t$-spread Veronese ideals and on Veronese ideals of bounded support. Their powers as well as their fiber cone is also be considered.


title: Prospects for New Physics searches with $\Lambda_b \to \Lambda(1520)\ell^+\ell^-$ decays

abstract: We present the prospects of an angular analysis of the $\Lambda_b \to \Lambda(1520)\ell^+\ell^-$ decay. Using the expected yield in the current dataset collected at the LHCb experiment, as well as the foreseen ones after the LHCb upgrades, sensitivity studies are presented to determine the experimental precision on angular observables related to the lepton distribution and their potential to identify New Physics. The forward-backward lepton asymmetry at low dilepton invariant mass is particularly promising. NP scenarios favoured by the current anomalies in $b\to s\ell^+\ell^-$ decays can be distinguished from the SM case with the data collected between the Run 3 and the Upgrade 2 of the LHCb experiment.


title: On hyperspherical associated Legendre functions: the extension of spherical harmonics to $N$ dimensions

abstract: The solution in hyperspherical coordinates for $N$ dimensions is given for a general class of partial differential equations of mathematical physics including the Laplace, wave, heat and Helmholtz, Schr\"{o}dinger, Klein-Gordon and telegraph equations and their combinations. The starting point is the Laplacian operator specified by the scale factors of hyperspherical coordinates. The general equation of mathematical physics is solved by separation of variables leading to the dependencies: (i) on time by the usual exponential function; (ii) on longitude by the usual sinusoidal function; (iii) on radius by Bessel functions of order generally distinct from cylindrical or spherical Bessel functions; (iv) on one latitude by associated Legendre functions; (v) on the remaining latitudes by an extension, namely the hyperspherical associated Legendre functions. The original associated Legendre functions are a particular case of the Gaussian hypergeometric functions, and the hyperspherical associated Legendre functions are also a more general particular case of the Gaussian hypergeometric functions so that it is not necessary to consider extended Gaussian hypergeometric functions.


title: Characterizing quantum correlations in spin chains

abstract: The growth in the demand for precisely crafted many-body systems of spin-$1/2$ particles/qubits is due to their top-notch versatility in application-oriented quantum-enhanced protocols and the fundamental tests of quantum theory. Here we address the question: how quantum is a chain of spins? We demonstrate that a single element of the density matrix carries the answer. Properly analyzed it brings information about the extent of the many-body entanglement and the non-locality. This method can be used to tailor and witness highly non-classical effects in many-body systems with possible applications to quantum computing, ultra-precise metrology or large-scale tests of quantum mechanics. As a proof of principle, we investigate the extend of non-locality and entanglement in ground states and thermal states of experimentally accessible spin chains.


title: Evaluating Policies Early in a Pandemic: Bounding Policy Effects with Nonrandomly Missing Data

abstract: During the early part of the Covid-19 pandemic, national and local governments introduced a number of policies to combat the spread of Covid-19. In this paper, we propose a new approach to bound the effects of such early-pandemic policies on Covid-19 cases and other outcomes while dealing with complications arising from (i) limited availability of Covid-19 tests, (ii) differential availability of Covid-19 tests across locations, and (iii) eligibility requirements for individuals to be tested. We use our approach study the effects of Tennessee's expansion of Covid-19 testing early in the pandemic and find that the policy decreased Covid-19 cases.


title: A Recipe for Creating Multimodal Aligned Datasets for Sequential Tasks

abstract: Many high-level procedural tasks can be decomposed into sequences of instructions that vary in their order and choice of tools. In the cooking domain, the web offers many partially-overlapping text and video recipes (i.e. procedures) that describe how to make the same dish (i.e. high-level task). Aligning instructions for the same dish across different sources can yield descriptive visual explanations that are far richer semantically than conventional textual instructions, providing commonsense insight into how real-world procedures are structured. Learning to align these different instruction sets is challenging because: a) different recipes vary in their order of instructions and use of ingredients; and b) video instructions can be noisy and tend to contain far more information than text instructions. To address these challenges, we first use an unsupervised alignment algorithm that learns pairwise alignments between instructions of different recipes for the same dish. We then use a graph algorithm to derive a joint alignment between multiple text and multiple video recipes for the same dish. We release the Microsoft Research Multimodal Aligned Recipe Corpus containing 150K pairwise alignments between recipes across 4,262 dishes with rich commonsense information.


title: Multiparameter tests of general relativity using multiband gravitational-wave observations

abstract: In this Letter we show that multiband observations of stellar-mass binary black holes by the next generation of ground-based observatories (3G) and the space-based Laser Interferometer Space Antenna (LISA) would facilitate a comprehensive test of general relativity by simultaneously measuring all the post-Newtonian (PN) coefficients. Multiband observations would measure most of the known PN phasing coefficients to an accuracy below a few percent---two orders-of-magnitude better than the best bounds achievable from even `golden' binaries in the 3G or LISA bands. Such multiparameter bounds would play a pivotal role in constraining the parameter space of modified theories of gravity beyond general relativity.


title: A moment inequality and positivity for signed graph Laplacians

abstract: A number of recent papers have considered signed graph Laplacians, a generalization of the classical graph Laplacian, where the edge weights are allowed to take either sign. In the classical case, where the edge weights are all positive, the Laplacian is positive semi-definite with the dimension of the kernel representing the number of connected components of the graph. In many applications one is interested in establishing conditions which guarantee the positive semi-definiteness of the matrix. In this paper we present an inequality on the eigenvalues of a weighted graph Laplacian (where the weights need not have any particular sign) in terms of the first two moments of the edge weights. This bound involves the eigenvalues of the equally weighted Laplacian on the graph as well as the eigenvalues of the adjacency matrix of the line graph (the edge-to-vertex dual graph). For a regular graph the bound can be expressed entirely in terms of the second eigenvalue of the equally weighted Laplacian, an object that has been extensively studied in connection with expander graphs and spectral measures of graph connectivity. We present several examples including Erd\H{o}s-R\'enyi random graphs in the critical and subcritical regimes, random large $d$-regular graphs, and the complete graph, for which the inequalities here are tight.


title: Deep learning classification of chest x-ray images

abstract: We propose a deep learning based method for classification of commonly occurring pathologies in chest X-ray images. The vast number of publicly available chest X-ray images provides the data necessary for successfully employing deep learning methodologies to reduce the misdiagnosis of thoracic diseases. We applied our method to the classification of two example pathologies, pulmonary nodules and cardiomegaly, and we compared the performance of our method to three existing methods. The results show an improvement in AUC for detection of nodules and cardiomegaly compared to the existing methods.


title: Free2Shard: Adaptive-adversary-resistant sharding via Dynamic Self Allocation

abstract: Propelled by the growth of large-scale blockchain deployments, much recent progress has been made in designing sharding protocols that achieve throughput scaling linearly in the number of nodes. However, existing protocols are not robust to an adversary adaptively corrupting a fixed fraction of nodes. In this paper, we propose Free2Shard -- a new architecture that achieves near-linear scaling while being secure against a fully adaptive adversary. The focal point of this architecture is a dynamic self-allocation algorithm that lets users allocate themselves to shards in response to adversarial action, without requiring a central or cryptographic proof. This architecture has several attractive features unusual for sharding protocols, including: (a) the ability to handle the regime of large number of shards (relative to the number of nodes); (b) heterogeneous shard demands; (c) requiring only a small minority to follow the self-allocation; (d) asynchronous shard rotation; (e) operation in a purely identity-free proof-of-work setting. The key technical contribution is a deep mathematical connection to the classical work of Blackwell in dynamic game theory.


title: Information-Theoretic Abstractions for Planning in Agents with Computational Constraints

abstract: In this paper, we develop a framework for path-planning on abstractions that are not provided to the agent a priori but instead emerge as a function of the available computational resources. We show how a path-planning problem in an environment can be systematically approximated by solving a sequence of easier to solve problems on abstractions of the original space. The properties of the problem are analyzed, and a number of theoretical results are presented and discussed. A numerical example is presented to show the utility of the approach and to corroborate the theoretical findings. We conclude by providing a discussion detailing the connections of the proposed approach to anytime algorithms and bounded rationality.


title: Quantum resource covariance

abstract: The developments of special relativity and quantum mechanics marked the beginning of the modern physics age. The former has taught us that while space and time are frame dependent notions, there is a quantity -- the space-time interval -- whose value all inertial observers agree upon. This reveals, so to speak, a genuine "fact" of the universe, a relativistic invariant. On the other hand, since the dawn of quantum mechanics, there is no consensus on what the theory is all about. The situation is admittedly subtler: quantum theory is grounded on a complex vector space and the very notions of observer and reference frame are controversial. Here we construct a theoretical framework within which a given combination of quantum resources is shown to be a Galilean invariant. To this end, we postulate a principle of relational symmetry between "the observer" and "the observed" and employ the notion of quantum reference frame. Unitary transformations then follow that allow us to perceive the physical resources seen from the viewpoint of any quantum system. Interestingly, we find that one needs more than quantum coherence and quantum correlations to prove quantum resources covariance. Finally, we show that the notion of physical reality implied by quantum mechanics is not absolute.


title: Imaging and spectral study on the null point of a fan-spine structure during a solar flare

abstract: Using the multi-instrument observations, we make the first simultaneous imaging and spectral study on the null point of a fan-spine magnetic topology during a solar flare. When magnetic reconnection occurs at the null point, the fan-spine configuration brightens in the (extreme-)ultraviolet channels. In the H$\alpha$ images, the fan-spine structure is partly filled and outlined by the bi-directional material flows ejected from the reconnection site. The extrapolated coronal magnetic field confirms the existence of the fan-spine topology. Before and after the flare peak, the total velocity of the outflows is estimated to be about 60 km s$^{-1}$. During the flare, the Si IV line profile at the reconnection region is enhanced both in the blue-wing and red-wing. At the flare peak time, the total velocity of the outflows is found to be 144 km s$^{-1}$. Superposed on the Si IV profile, there are several deep absorption lines with the blueshift of several tens of km s$^{-1}$. The reason is inferred to be that the bright reconnection region observed in Si IV channel is located under the cooler material appearing as dark features in the H$\alpha$ line. The blueshifted absorption lines indicate the movement of the cooler material toward the observer. The depth of the absorption lines also depends on the amount of cooler material. These results imply that this kind of spectral profiles can be used as a tool to diagnose the properties of cooler material above reconnection site.


title: Every complete Pick space satisfies the column-row property

abstract: In the theory of complete Pick spaces, the column-row property has appeared in a variety of contexts. We show that it is satisfied by every complete Pick space in the following strong form: each sequence of multipliers that induces a contractive column multiplication operator also induces a contractive row multiplication operator. In combination with known results, this yields a number of consequences. Firstly, we obtain multiple applications to the theory of weak product spaces, including factorization, multipliers and invariant subspaces. Secondly, there is a short proof of the characterization of interpolating sequences in terms of separation and Carleson measure conditions, independent of the solution of the Kadison-Singer problem. Thirdly, we find that in the theory of de Branges-Rovnyak spaces on the ball, the column-extreme multipliers of Jury and Martin are precisely the extreme points of the unit ball of the multiplier algebra.


title: ARES II: Characterising the Hot Jupiters WASP-127 b, WASP-79 b and WASP-62 b with HST

abstract: This paper presents the atmospheric characterisation of three large, gaseous planets: WASP-127b, WASP-79b and WASP-62b. We analysed spectroscopic data obtained with the G141 grism (1.088 - 1.68 $\mu$m) of the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) using the Iraclis pipeline and the TauREx3 retrieval code, both of which are publicly available. For WASP-127 b, which is the least dense planet discovered so far and is located in the short-period Neptune desert, our retrieval results found strong water absorption corresponding to an abundance of log(H$_2$O) = -2.71$^{+0.78}_{-1.05}$, and absorption compatible with an iron hydride abundance of log(FeH)=$-5.25^{+0.88}_{-1.10}$, with an extended cloudy atmosphere. We also detected water vapour in the atmospheres of WASP-79 b and WASP-62 b, with best-fit models indicating the presence of iron hydride, too. We used the Atmospheric Detectability Index (ADI) as well as Bayesian log evidence to quantify the strength of the detection and compared our results to the hot Jupiter population study by Tsiaras et al. 2018. While all the planets studied here are suitable targets for characterisation with upcoming facilities such as the James Webb Space Telescope (JWST) and Ariel, WASP-127 b is of particular interest due to its low density, and a thorough atmospheric study would develop our understanding of planet formation and migration.


title: Group segmentation and heterogeneity in the choice of cooking fuels in post-earthquake Nepal

abstract: Segmenting population into subgroups with higher intergroup, but lower intragroup, heterogeneity can be useful in enhancing the effectiveness of many socio-economic policy interventions; yet it has received little attention in promoting clean cooking. Here, we use PERMANOVA, a distance-based multivariate analysis, to identify the factor that captures the highest intergroup heterogeneity in the choice of cooking fuels. Applying this approach to the post-earthquake data on 747,137 households from Nepal, we find that ethnicity explains 39.12% of variation in fuel choice, followed by income (26.30%), education (12.62%), and location (4.05%). This finding indicates that ethnicity, rather than income or other factors, as a basis of policy interventions may be more effective in promoting clean cooking. We also find that, among the ethnic groups in Nepal, the most marginalized Chepang/Thami community exhibits the lowest intragroup diversity (Shannon index = 0.101) while Newars the highest (0.667). This information on intra-ethnic diversity in fuel choice can have important policy implications for reducing ethnic gap in clean cooking.


title: Unlocking New York City Crime Insights using Relational Database Embeddings

abstract: This version withdrawn by arXiv administrators because the author did not have the right to agree to our license at the time of submission.


title: Transport in neutron star mergers

abstract: After an introduction to the QCD phase diagram, the nuclear equations of state, and neutron star mergers, I discuss three projects related to transport and nuclear matter in neutron star mergers. The first is the nature of beta equilibrium in the portion of a merger that is transparent to neutrinos. We calculate the weak interaction (Urca) rates and find that the beta equilibrium condition needs to be modified by adding an additional chemical potential, which changes slightly the particle content in neutrino-transparent beta equilibrium. Secondly, we calculate the bulk viscosity in neutrino-transparent nuclear matter in conditions encountered in neutron star mergers. Bulk viscosity arises from a phase lag between the pressure and density in the nuclear matter, which is due to the finite rate of beta equilibration. When bulk viscosity is sufficiently strong, which happens when the equilibration rate nearly matches the frequency of the density oscillation, it can noticeably dampen the oscillation. We find that in certain thermodynamic conditions likely encountered in mergers, oscillations in nuclear matter can be damped on timescales on the order of 10 milliseconds, so we conclude that bulk viscosity should be included in merger simulations. Finally, we study thermal transport due to axions in neutron star mergers. We conclude that axions are never trapped in mergers, but instead escape, carrying energy away from the merger. We calculate the cooling time due to the energy carried away by axions and find that within current constraints on the axion-nucleon coupling, axions could cool fluid elements in mergers on timescales which could affect the dynamics of the merger.


title: Identifying Statistical Bias in Dataset Replication

abstract: Dataset replication is a useful tool for assessing whether improvements in test accuracy on a specific benchmark correspond to improvements in models' ability to generalize reliably. In this work, we present unintuitive yet significant ways in which standard approaches to dataset replication introduce statistical bias, skewing the resulting observations. We study ImageNet-v2, a replication of the ImageNet dataset on which models exhibit a significant (11-14%) drop in accuracy, even after controlling for a standard human-in-the-loop measure of data quality. We show that after correcting for the identified statistical bias, only an estimated $3.6\% \pm 1.5\%$ of the original $11.7\% \pm 1.0\%$ accuracy drop remains unaccounted for. We conclude with concrete recommendations for recognizing and avoiding bias in dataset replication. Code for our study is publicly available at http://github.com/MadryLab/dataset-replication-analysis .


title: Parametric Hamilton's equations for Lagrangian model for passive scalar gradients

abstract: In the context of instanton method for stochastic system this paper purposes a modification of the arclength parametrization of the Hamilton's equations allowing for an arbitrary instanton speed. The main results of the paper are: (i) it generalizes the parametrized Hamilton's equations to any speed required. (ii) corrects the parametric action on the occasion that the Hamiltonian is small but finite and how it adjusts to the probability density function (pdf). (iii) Improves instanton approximation to pdf by noise and propagator renormalization. As an application of the above set up we predict the statistics of passive scalar gradients in a Lagrangian model for turbulence, namely the scalar gradient Recent Fluid Deformation Closure (sgRFD).


title: Reconfigurable electro-optic frequency shifter

abstract: Here we demonstrate an on-chip electro-optic frequency shifter that is precisely controlled using only a single-tone microwave signal. This is accomplished by engineering the density of states of, and coupling between, optical modes in ultra-low loss electro-optic waveguides and resonators realized in lithium niobate nanophotonics. Our device provides frequency shifts as high as 28 GHz with measured shift efficiencies of ~99% and insertion loss of <0.5 dB. Importantly, the device can be reconfigured as a tunable frequency-domain beam splitter, in which the splitting ratio and splitting frequency are controlled by microwave power and frequency, respectively. Using the device, we also demonstrate (non-blocking) frequency routing through an efficient exchange of information between two distinct frequency channels, i.e. swap operation. Finally, we show that our scheme can be scaled to achieve cascaded frequency shifts beyond 100 GHz. Our device could become an essential building-block for future high-speed and large-scale classical information processors as well as emerging frequency-domain photonic quantum computers.


title: Learning Equations from Biological Data with Limited Time Samples

abstract: Equation learning methods present a promising tool to aid scientists in the modeling process for biological data. Previous equation learning studies have demonstrated that these methods can infer models from rich datasets, however, the performance of these methods in the presence of common challenges from biological data has not been thoroughly explored. We present an equation learning methodology comprised of data denoising, equation learning, model selection and post-processing steps that infers a dynamical systems model from noisy spatiotemporal data. The performance of this methodology is thoroughly investigated in the face of several common challenges presented by biological data, namely, sparse data sampling, large noise levels, and heterogeneity between datasets. We find that this methodology can accurately infer the correct underlying equation and predict unobserved system dynamics from a small number of time samples when the data is sampled over a time interval exhibiting both linear and nonlinear dynamics. Our findings suggest that equation learning methods can be used for model discovery and selection in many areas of biology when an informative dataset is used. We focus on glioblastoma multiforme modeling as a case study in this work to highlight how these results are informative for data-driven modeling-based tumor invasion predictions.


title: Focus on defocus: bridging the synthetic to real domain gap for depth estimation

abstract: Data-driven depth estimation methods struggle with the generalization outside their training scenes due to the immense variability of the real-world scenes. This problem can be partially addressed by utilising synthetically generated images, but closing the synthetic-real domain gap is far from trivial. In this paper, we tackle this issue by using domain invariant defocus blur as direct supervision. We leverage defocus cues by using a permutation invariant convolutional neural network that encourages the network to learn from the differences between images with a different point of focus. Our proposed network uses the defocus map as an intermediate supervisory signal. We are able to train our model completely on synthetic data and directly apply it to a wide range of real-world images. We evaluate our model on synthetic and real datasets, showing compelling generalization results and state-of-the-art depth prediction.


title: Batch-Augmented Multi-Agent Reinforcement Learning for Efficient Traffic Signal Optimization

abstract: The goal of this work is to provide a viable solution based on reinforcement learning for traffic signal control problems. Although the state-of-the-art reinforcement learning approaches have yielded great success in a variety of domains, directly applying it to alleviate traffic congestion can be challenging, considering the requirement of high sample efficiency and how training data is gathered. In this work, we address several challenges that we encountered when we attempted to mitigate serious traffic congestion occurring in a metropolitan area. Specifically, we are required to provide a solution that is able to (1) handle the traffic signal control when certain surveillance cameras that retrieve information for reinforcement learning are down, (2) learn from batch data without a traffic simulator, and (3) make control decisions without shared information across intersections. We present a two-stage framework to deal with the above-mentioned situations. The framework can be decomposed into an Evolution Strategies approach that gives a fixed-time traffic signal control schedule and a multi-agent off-policy reinforcement learning that is capable of learning from batch data with the aid of three proposed components, bounded action, batch augmentation, and surrogate reward clipping. Our experiments show that the proposed framework reduces traffic congestion by 36% in terms of waiting time compared with the currently used fixed-time traffic signal plan. Furthermore, the framework requires only 600 queries to a simulator to achieve the result.


title: Inference, prediction and optimization of non-pharmaceutical interventions using compartment models: the PyRoss library

abstract: PyRoss is an open-source Python library that offers an integrated platform for inference, prediction and optimisation of NPIs in age- and contact-structured epidemiological compartment models. This report outlines the rationale and functionality of the PyRoss library, with various illustrations and examples focusing on well-mixed, age-structured populations. The PyRoss library supports arbitrary structured models formulated stochastically (as master equations) or deterministically (as ODEs) and allows mid-run transitioning from one to the other. By supporting additional compartmental subdivision ad libitum, PyRoss can emulate time-since-infection models and allows medical stages such as hospitalization or quarantine to be modelled and forecast. The PyRoss library enables fitting to epidemiological data, as available, using Bayesian parameter inference, so that competing models can be weighed by their evidence. PyRoss allows fully Bayesian forecasts of the impact of idealized NPIs by convolving uncertainties arising from epidemiological data, model choice, parameters, and intrinsic stochasticity. Algorithms to optimize time-dependent NPI scenarios against user-defined cost functions are included. PyRoss's current age-structured compartment framework for well-mixed populations will in future reports be extended to include compartments structured by location, occupation, use of travel networks and other attributes relevant to assessing disease spread and the impact of NPIs. We argue that such compartment models, by allowing social data of arbitrary granularity to be combined with Bayesian parameter estimation for poorly-known disease variables, could enable more powerful and robust prediction than other approaches to detailed epidemic modelling. We invite others to use the PyRoss library for research to address today's COVID-19 crisis, and to plan for future pandemics.


title: On boundedness of divisors computing minimal log discrepancies for surfaces

abstract: Let $\Gamma$ be a finite set, and $X\ni x$ a fixed klt germ. For any lc germ $(X\ni x,B:=\sum_{i} b_iB_i)$ such that $b_i\in \Gamma$, Nakamura's conjecture, which is equivalent to the ACC conjecture for minimal log discrepancies for fixed germs, predicts that there always exists a prime divisor $E$ over $X\ni x$, such that $a(E,X,B)={\rm{mld}}(X\ni x,B)$, and $a(E,X,0)$ is bounded from above. We extend Nakamura's conjecture to the setting that $X\ni x$ is not necessarily fixed and $\Gamma$ satisfies the DCC, and show it holds for surfaces. We also find some sufficient conditions for the boundedness of $a(E,X,0)$ for any such $E$.


title: One Size Fits All: Can We Train One Denoiser for All Noise Levels?

abstract: When training an estimator such as a neural network for tasks like image denoising, it is often preferred to train one estimator and apply it to all noise levels. The de facto training protocol to achieve this goal is to train the estimator with noisy samples whose noise levels are uniformly distributed across the range of interest. However, why should we allocate the samples uniformly? Can we have more training samples that are less noisy, and fewer samples that are more noisy? What is the optimal distribution? How do we obtain such a distribution? The goal of this paper is to address this training sample distribution problem from a minimax risk optimization perspective. We derive a dual ascent algorithm to determine the optimal sampling distribution of which the convergence is guaranteed as long as the set of admissible estimators is closed and convex. For estimators with non-convex admissible sets such as deep neural networks, our dual formulation converges to a solution of the convex relaxation. We discuss how the algorithm can be implemented in practice. We evaluate the algorithm on linear estimators and deep networks.


title: Lattice polytopes from Schur and symmetric Grothendieck polynomials

abstract: Given a family of lattice polytopes, two common questions in Ehrhart Theory are determining when a polytope has the integer decomposition property and determining when a polytope is reflexive. While these properties are of independent interest, the confluence of these properties is a source of active investigation due to conjectures regarding the unimodality of the $h^\ast$-polynomial. In this paper, we consider the Newton polytopes arising from two families of polynomials in algebraic combinatorics: Schur polynomials and inflated symmetric Grothendieck polynomials. In both cases, we prove that these polytopes have the integer decomposition property by using the fact that both families of polynomials have saturated Newton polytope. Furthermore, in both cases, we provide a complete characterization of when these polytopes are reflexive. We conclude with some explicit formulas and unimodality implications of the $h^\ast$-vector in the case of Schur polynomials.


title: Improved Noisy Student Training for Automatic Speech Recognition

abstract: Recently, a semi-supervised learning method known as "noisy student training" has been shown to improve image classification performance of deep networks significantly. Noisy student training is an iterative self-training method that leverages augmentation to improve network performance. In this work, we adapt and improve noisy student training for automatic speech recognition, employing (adaptive) SpecAugment as the augmentation method. We find effective methods to filter, balance and augment the data generated in between self-training iterations. By doing so, we are able to obtain word error rates (WERs) 4.2%/8.6% on the clean/noisy LibriSpeech test sets by only using the clean 100h subset of LibriSpeech as the supervised set and the rest (860h) as the unlabeled set. Furthermore, we are able to achieve WERs 1.7%/3.4% on the clean/noisy LibriSpeech test sets by using the unlab-60k subset of LibriLight as the unlabeled set for LibriSpeech 960h. We are thus able to improve upon the previous state-of-the-art clean/noisy test WERs achieved on LibriSpeech 100h (4.74%/12.20%) and LibriSpeech (1.9%/4.1%).


title: A systems biology approach to COVID-19 progression in a population

abstract: A number of models in mathematical epidemiology have been developed to account for control measures such as vaccination or quarantine. However, COVID-19 has brought unprecedented social distancing measures, with a challenge on how to include these in a manner that can explain the data but avoid overfitting in parameter inference. We here develop a simple time-dependent model, where social distancing effects are introduced analogous to coarse-grained models of gene expression control in systems biology. We apply our approach to understand drastic differences in COVID-19 infection and fatality counts, observed between Hubei (Wuhan) and other Mainland China provinces. We find that these unintuitive data may be explained through an interplay of differences in transmissibility, effective protection, and detection efficiencies between Hubei and other provinces. More generally, our results demonstrate that regional differences may drastically shape infection outbursts. The obtained results demonstrate the applicability of our developed method to extract key infection parameters directly from publically available data so that it can be globally applied to outbreaks of COVID-19 in a number of countries. Overall, we show that applications of uncommon strategies, such as methods and approaches from molecular systems biology research to mathematical epidemiology, may significantly advance our understanding of COVID-19 and other infectious diseases.


title: Confirmation of water emission in the dayside spectrum of the ultrahot Jupiter WASP-121b

abstract: We present four new secondary eclipse observations for the ultrahot Jupiter WASP-121b acquired using the Hubble Space Telescope Wide Field Camera 3. The eclipse depth is measured to a median precision of 60ppm across 28 spectroscopic channels spanning the 1.12-1.64 micron wavelength range. This is a considerable improvement to the 90ppm precision we achieved previously for a single eclipse observation using the same observing setup. Combining these data with those reported at other wavelengths, a blackbody spectrum for WASP-121b is ruled out at >6-sigma confidence and we confirm the interpretation of previous retrieval analyses that found the data is best explained by a dayside thermal inversion. The updated spectrum clearly resolves the water emission band at 1.3-1.6 micron, with higher signal-to-noise than before. It also fails to reproduce a bump in the spectrum at 1.25 micron derived from the first eclipse observation, which had tentatively been attributed to VO emission. We conclude the latter was either a statistical fluctuation or a systematic artefact specific to the first eclipse dataset.


title: Markov Chains for Horizons (MARCH). I. Identifying Biases in Fitting Theoretical Models to Event Horizon Telescope Observations

abstract: We introduce a new Markov Chain Monte Carlo (MCMC) algorithm with parallel tempering for fitting theoretical models of horizon-scale images of black holes to the interferometric data from the Event Horizon Telescope (EHT). The algorithm implements forms of the noise distribution in the data that are accurate for all signal-to-noise ratios. In addition to being trivially parallelizable, the algorithm is optimized for high performance, achieving 1 million MCMC chain steps in under 20 seconds on a single processor. We use synthetic data for the 2017 EHT coverage of M87 that are generated based on analytic as well as General Relativistic Magnetohydrodynamic (GRMHD) model images to explore several potential sources of biases in fitting models to sparse interferometric data. We demonstrate that a very small number of data points that lie near salient features of the interferometric data exert disproportionate influence on the inferred model parameters. We also show that the preferred orientations of the EHT baselines introduce significant biases in the inference of the orientation of the model images. Finally, we discuss strategies that help identify the presence and severity of such biases in realistic applications.


title: A New Metal-poor Globular Cluster and Resolved Stars in the Outer Disk of the Black Eye Galaxy M64: Implication for the Origin of the Type III Disk Break

abstract: M64 is a nearby spiral galaxy with a Type III anti-truncation component. To trace the origin of the Type III component, we present Hubble Space Telescope/Advanced Camera for Surveys $F606W/F814W$ photometry of resolved stars in a field located in the outer disk ($2.5' \lesssim r \lesssim 6.5'$) of M64. At $r\approx 5.5'$ (7 kpc) to the east, we discover a new metal-poor globular cluster ($R_{\rm eff}=5.73\pm0.02$ pc and $M_V=-9.54\pm0.09$ mag), M64-GC1. This is the first globular cluster found in M64. The color-magnitude diagram (CMD) of the resolved stars in M64-GC1 is well matched by 12 Gyr isochrones with [Fe/H] $=-1.5\pm0.2$, showing that this cluster belongs to a halo. The CMD of the resolved stars in the entire ACS field shows two distinguishable red giant branches (RGBs): a curved metal-rich RGB and a vertical metal-poor RGB. The metal-rich RGB represents an old metal-rich ([Fe/H] $\approx -0.4$) disk population. In contrast, the CMD of the metal-poor RGB stars is very similar to the CMD of M64-GC1, showing that the metal-poor RGB represents a halo population. The radial number density profile of the metal-rich RGB stars is described by an exponential disk law, while the profile of the metal-poor RGB stars is described by a de Vaucouleurs's law. From these, we conclude that the origin of the Type III component in M64 is a halo which has a much lower metallicity than a disk or bulge population.


title: Automated Copper Alloy Grain Size Evaluation Using a Deep-learning CNN

abstract: Moog Inc. has automated the evaluation of copper (Cu) alloy grain size using a deep-learning convolutional neural network (CNN). The proof-of-concept automated image acquisition and batch-wise image processing offers the potential for significantly reduced labor, improved accuracy of grain evaluation, and decreased overall turnaround times for approving Cu alloy bar stock for use in flight critical aircraft hardware. A classification accuracy of 91.1% on individual sub-images of the Cu alloy coupons was achieved. Process development included minimizing the variation in acquired image color, brightness, and resolution to create a dataset with 12300 sub-images, and then optimizing the CNN hyperparameters on this dataset using statistical design of experiments (DoE). Over the development of the automated Cu alloy grain size evaluation, a degree of "explainability" in the artificial intelligence (XAI) output was realized, based on the decomposition of the large raw images into many smaller dataset sub-images, through the ability to explain the CNN ensemble image output via inspection of the classification results from the individual smaller sub-images.


title: InterFaceGAN: Interpreting the Disentangled Face Representation Learned by GANs

abstract: Although Generative Adversarial Networks (GANs) have made significant progress in face synthesis, there lacks enough understanding of what GANs have learned in the latent representation to map a random code to a photo-realistic image. In this work, we propose a framework called InterFaceGAN to interpret the disentangled face representation learned by the state-of-the-art GAN models and study the properties of the facial semantics encoded in the latent space. We first find that GANs learn various semantics in some linear subspaces of the latent space. After identifying these subspaces, we can realistically manipulate the corresponding facial attributes without retraining the model. We then conduct a detailed study on the correlation between different semantics and manage to better disentangle them via subspace projection, resulting in more precise control of the attribute manipulation. Besides manipulating the gender, age, expression, and presence of eyeglasses, we can even alter the face pose and fix the artifacts accidentally made by GANs. Furthermore, we perform an in-depth face identity analysis and a layer-wise analysis to evaluate the editing results quantitatively. Finally, we apply our approach to real face editing by employing GAN inversion approaches and explicitly training feed-forward models based on the synthetic data established by InterFaceGAN. Extensive experimental results suggest that learning to synthesize faces spontaneously brings a disentangled and controllable face representation.


title: Molecular-field-theory fits to magnetic susceptibilities of antiferromagnetic GdCu2Si2, CuO, LiCrO2, and alpha-CaCr2O4 single crystals below their Neel temperatures

abstract: A recently-developed molecular field theory (MFT) has been used to fit single-crystal magnetic susceptibility chi versus temperature T data below the respective antiferromagnetic ordering temperatures TN for a variety of collinear and coplanar noncollinear Heisenberg antiferromagnets. The spins in the system are assumed to interact by Heisenberg exchange and to be identical and crystallographically equivalent. The fitting parameters for chi(T) of collinear antiferromagnets are measurable quantities: the Weiss temperature theta_p in the Curie-Weiss law, TN, chi(TN), and the spin S. For coplanar noncollinear helix and cycloid structures, an additional fitting parameter is the turn angle between layers of ferromagnetically-aligned spins. Here MFT fits to anisotropic chi(T) data from the literature for single crystals of the collinear antiferromagnets GdCu2Si2 and CuO and the noncollinear antiferromagnets LiCr2As2 with a 120 degree cycloidal structure and alpha-CaCr2O4 with a 120 degree helical structure below their respective Neel temperatures are presented. The MFT fit to the anisotropic chi(T < TN) data for CuO is poor, whereas the fits to the data for GdCu2Si2, LiCrO2, and alpha-CaCr2O4 are quite good. The poor fit for CuO is attributed to the influence of strong quantum fluctuations associated with the small Cu spin and the quasi-one-dimensional magnetism that are not taken into account by the MFT. The magnetic contribution to the zero-field heat capacity of the collinear antiferromagnet GdNiGe3 at T < TN is also fitted by the MFT.


title: On the inverse source identification problem in $L^\infty$ for fully nonlinear elliptic PDE

abstract: In this paper we generalise the results proved in [N. Katzourakis, An $L^\infty$ regularisation strategy to the inverse source identification problem for elliptic equations, SIAM J. Math. Anal. 51:2, 1349-1370 (2019)] by studying the ill-posed problem of identifying the source of a fully nonlinear elliptic equation. We assume Dirichlet data and some partial noisy information for the solution on a compact set through a fully nonlinear observation operator. We deal with the highly nonlinear nonconvex nature of the problem and the lack of weak continuity by introducing a two-parameter Tykhonov regularisation with a higher order $L^2$ "viscosity term" for the $L^\infty$ minimisation problem which allows to approximate by weakly lower semicontinuous cost functionals.


title: Physics-informed Neural Networks for Solving Inverse Problems of Nonlinear Biot's Equations: Batch Training

abstract: In biomedical engineering, earthquake prediction, and underground energy harvesting, it is crucial to indirectly estimate the physical properties of porous media since the direct measurement of those are usually impractical/prohibitive. Here we apply the physics-informed neural networks to solve the inverse problem with regard to the nonlinear Biot's equations. Specifically, we consider batch training and explore the effect of different batch sizes. The results show that training with small batch sizes, i.e., a few examples per batch, provides better approximations (lower percentage error) of the physical parameters than using large batches or the full batch. The increased accuracy of the physical parameters, comes at the cost of longer training time. Specifically, we find the size should not be too small since a very small batch size requires a very long training time without a corresponding improvement in estimation accuracy. We find that a batch size of 8 or 32 is a good compromise, which is also robust to additive noise in the data. The learning rate also plays an important role and should be used as a hyperparameter.


title: Webpage Segmentation for Extracting Images and Their Surrounding Contextual Information

abstract: Web images come in hand with valuable contextual information. Although this information has long been mined for various uses such as image annotation, clustering of images, inference of image semantic content, etc., insufficient attention has been given to address issues in mining this contextual information. In this paper, we propose a webpage segmentation algorithm targeting the extraction of web images and their contextual information based on their characteristics as they appear on webpages. We conducted a user study to obtain a human-labeled dataset to validate the effectiveness of our method and experiments demonstrated that our method can achieve better results compared to an existing segmentation algorithm.


title: Torus-breakdown near a Bykov attractor: a case study

abstract: There are few explicit examples in the literature of vector fields exhibiting complex dynamics that may be proved analytically. This paper reports numerical experiments performed for an explicit two-parameter family of vector fields unfolding an attracting heteroclinic network, linking two saddle-foci with $ (\mathbb{SO}(2) \oplus \mathbb{Z}_2)$-symmetry. The vector field is the restriction to $\mathbb{S}^3$ of a polynomial vector field in $\mathbb{R}^4$. We investigate global bifurcations due to symmetry-breaking and we detect strange attractors via a phenomenon called Torus-Breakdown theory. We explain how an attracting torus gets destroyed by following the changes in the invariant manifolds of the saddle-foci. Although a complete understanding of the corresponding bifurcation diagram and the mechanisms underlying the dynamical changes is still out of reach, using a combination of theoretical tools and computer simulations, we have uncovered some complex patterns for the symmetric family under analysis. This also suggests a route to obtain rotational horseshoes; additionally, we give an attempt to elucidate some of the bifurcations involved in an Arnold wedge.


title: Bessel Type Orthogonality For Hermite Polynomials

abstract: It is shown that Hermite polynomials satisfy a Bessel type orthogonality relation, based on the zeros of a single index Hermite polynomial and with a finite integration interval. Because of the role of non-symmetric zeros in the final relation, its applicability covers Hermite polynomials $P_n(x)$ with $n\ge 3$.


title: COMPASS: VLBI Beacons In Support of Lunar Science and Exploration

abstract: The large constellations of spacecraft planned for use in cislunar space (on the Lunar surface, in Lunar orbit, and in the vicinity of the Lunar Gateway) require new solutions for positioning, navigation and timing (PNT). Here, I describe COMPASS (Combined Observational Methods for Positional Awareness in the Solar System), a spacecraft navigation system to provide cost-effective techniques for the positioning of large numbers of spacecraft in cislunar space. COMPASS will use beacons that emit coherent ultra-wideband signals designed to be interoperable with existing and future Very Long Baseline Interferometry (VLBI) networks. Using differential VLBI, COMPASS will provide rapid determination of the interferometric phase delay with picosecond level accuracy during routine VLBI observing sessions. Multi-baseline phase-referenced COMPASS-VLBI observations with simultaneous calibrator observations should thus enable sub-meter accuracy transverse positioning and meter level lunar orbit determination using with small femtospacecraft beacons and a few seconds of observation per position determination.


title: An Innovative Approach to Determine Rebar Depth and Size by Comparing GPR Data with a Theoretical Database

abstract: Ground penetrating radar (GPR) is an efficient technique used for rapidly recognizing embedded rebar in concrete structures. However, due to the difficulty in extracting signals from GPR data and the intrinsic coupling between the rebar depth and size showing in the data, simultaneously determining rebar depth and size is challenging. This paper proposes an innovative algorithm to address this issue. First, the hyperbola signal from the GPR data is identified by direct wave removal, signal reconstruction and separation. Subsequently, a database is developed from a series of theoretical hyperbolas and then compared with the extracted hyperbola outlines. Finally, the rebar depth and size are determined by searching for the closest counterpart in the database. The obtained results are very promising and indicate that: (1) implementing the method presented in this paper can completely remove the direct wave noise from the GPR data, and can successfully extract the outlines from the interlaced hyperbolas; and (2) the proposed method can simultaneously determine the rebar depth and size with the accuracy of 100% and 95.11%, respectively.


title: A Statistical Model for Imaging Systems

abstract: The behavior of photons is controlled by quantum mechanics, not as deterministic as classical optics shows. To this end, we defined a new statistic $Z$, which is equal to the variance minus the expectation or mean. Then, we established a statistical model for imaging systems and obtained three fundamental imaging formulas. Among them, the first formula is entirely consistent with the classic convolution equation. The second and third ones link the $Z$ quantities of the object and noise elegantly with the $Z$ image and covariance image, revealing new laws. Consequently, besides the flux density, the $Z$ quantity of an object is also imageable, which opens a new realm for imaging systems to explore the physical world.

categories: eess.IV physics.data-an physics.optics stat.ME


title: The Second Type of Uncertainty in Monte Carlo Tree Search

abstract: Monte Carlo Tree Search (MCTS) efficiently balances exploration and exploitation in tree search based on count-derived uncertainty. However, these local visit counts ignore a second type of uncertainty induced by the size of the subtree below an action. We first show how, due to the lack of this second uncertainty type, MCTS may completely fail in well-known sparse exploration problems, known from the reinforcement learning community. We then introduce a new algorithm, which estimates the size of the subtree below an action, and leverages this information in the UCB formula to better direct exploration. Subsequently, we generalize these ideas by showing that loops, i.e., the repeated occurrence of (approximately) the same state in the same trace, are actually a special case of subtree depth variation. Testing on a variety of tasks shows that our algorithms increase sample efficiency, especially when the planning budget per timestep is small.


title: Study of the hidden charm pentaquarks $\Lambda_{c\bar{c}}$ above $4$ GeV

abstract: A dynamical calculation of the strange hidden charm pentaquark systems $\Lambda_{c\bar{c}}$ is performed in the framework of the quark delocalization color screening model. The effective potentials between two clusters are calculated to investigate the interaction between the baryons and mesons. The dynamic calculations indicate that the strange hidden charm pentaquark system with $IJ^{P}=0\frac{1}{2}^{-}$ and $IJ^{P}=0\frac{3}{2}^{-}$ can form bound states with the help of the channel-coupling. The energies of these two system are 4094.3 MeV and 4207.4 Mev, respectively. These two $\Lambda_{c\bar{c}}$ states are possible to be intriguing pentaquark candidates which are worth searching in experiments. Whereas, for the systems with isospin $I=1$, the effective potentials of all channels are repulsive, and there is no any bound state for these systems.


title: Von Neumann Entropy from Mean Spin Vector

abstract: We show for a general pure entangled state of two two-level atoms, the von Neumann entropy of the partial traces can be directly measured from the magnitude of the mean spin vector of a single atom of the pair. We emphasize the fact that the von Neumann entropy of the partial traces for such a system can be obtained without knowing the exact form of the quantum state of the two atoms, if we have the value of the magnitude of the mean spin vector of a single atom of the pair. Mean spin vector, used in the context of spin squeezing and spectroscopic squeezing in population spectroscopy, is experimentally measurable and provides an exact measure of von Neumann entropy of the partial traces for such a system. The idea developed in this paper can be used in the context of other quantum mechanical two level systems as the algebra of two level systems can be described by that of spin- $\frac{1}{2}$ particles.


title: A comparison of limited-stretch models of rubber elasticity

abstract: In this paper we describe various limited-stretch models of nonlinear rubber elasticity, each dependent on only the first invariant of the left Cauchy-Green strain tensor and having only two independent material constants. The models are described as limited-stretch, or restricted elastic, because the strain energy and stress response become infinite at a finite value of the first invariant. These models describe well the limited stretch of the polymer chains of which rubber is composed. We discuss Gent's model which is the simplest limited-stretch model and agrees well with experiment. Various statistical models are then described: the one-chain, three-chain, four-chain and Arruda-Boyce eight-chain models, all of which involve the inverse Langevin function. A numerical comparison between the three-chain and eight-chain models is provided. Next, we compare various models which involve approximations to the inverse Langevin function with the exact inverse Langevin function of the eight-chain model. A new approximate model is proposed that is as simple as Cohen's original model but significantly more accurate. We show that effectively the eight-chain model may be regarded as a linear combination of the neo-Hookean and Gent models. Treloar's model is shown to have about half the percentage error of our new model but it is much more complicated. For completeness a modified Treloar model is introduced but this is only slightly more accurate than Treloar's original model. For the deformations of uniaxial tension, biaxial tension, pure shear and simple shear we compare the accuracy of these models, and that of Puso, with the eight-chain model by means of graphs and a table. Our approximations compare extremely well with models frequently used and described in the literature, having the smallest mean percentage error over most of the range of the argument.