title: The relative contribution to heavy metals production from binary neutron star mergers and neutron star-black hole mergers

abstract: The origin of the heavy elements in the Universe is not fully determined. Neutron star-black hole (NSBH) and {binary neutron star} (BNS) mergers may both produce heavy elements via rapid neutron-capture (r-process). We use the recent detection of gravitational waves from NSBHs, improved measurements of the neutron star equation-of-state, and the most modern numerical simulations of ejected material from binary collisions to measure the relative contribution of NSBHs and BNSs to the production of heavy elements. As the amount of r-process ejecta depends on the mass and spin distribution of the compact objects, as well as on the equation-of-state of the neutron stars, we consider various models for these quantities, informed by gravitational-wave and pulsar data. We find that in most scenarios, BNSs have produced more r-process elements than NSBHs over the past 2.5 billion years. If black holes have preferentially small spins, BNSs can produce at least twice of the amount of r-process elements than NSBHs. If black hole spins are small and there is a dearth of low mass ($<5M_{\odot}$) black holes within NSBH binaries, BNSs can account for the near totality of the r-process elements from binaries. For NSBH to produce large fraction of r-process elements, black holes in NSBHs must have small masses and large aligned spins, which is disfavored by current data.


title: Geometrical Postural Optimisation of 7-DoF Limb-Like Manipulators

abstract: Robots are moving towards applications in less structured environments, but their model-based controllers are challenged by the tasks' complexity and intrinsic environmental unpredictability. Studying biological motor control can provide insights into overcoming these limitations due to the high dexterity and stability observable in humans and animals. This work presents a geometrical solution to the postural optimisation of 7-DoF limbs-like mechanisms, which are robust to singularities and computationally efficient. The theoretical formulation identified two separate decoupled optimisation strategies. The shoulder and elbow strategy align the plane of motion with the expected plane of motion and guarantee the reachability of the end-posture. The wrist strategy ensures the end-effector orientation, which is essential to retain manipulability when nearing a singular configuration. The numerical results confirmed the theoretical observations and allowed us to identify the effect of different grasp strategies on system manipulability. The geometrical method was numerically tested in thousands of configurations proving to be both robust and accurate. The tested scenarios include left and right arm postures, singular configurations, and walking scenarios. The proposed geometrical approach can find application in developing efficient and robust interaction controllers that could be applied in computational neuroscience and robotics.


title: Evaluating subgroup disparity using epistemic uncertainty in mammography

abstract: As machine learning (ML) continue to be integrated into healthcare systems that affect clinical decision making, new strategies will need to be incorporated in order to effectively detect and evaluate subgroup disparities to ensure accountability and generalizability in clinical workflows. In this paper, we explore how epistemic uncertainty can be used to evaluate disparity in patient demographics (race) and data acquisition (scanner) subgroups for breast density assessment on a dataset of 108,190 mammograms collected from 33 clinical sites. Our results show that even if aggregate performance is comparable, the choice of uncertainty quantification metric can significantly the subgroup level. We hope this analysis can promote further work on how uncertainty can be leveraged to increase transparency of machine learning applications for clinical deployment.


title: The rapid destruction of toroidal magnetic surfaces

abstract: The operation of ITER will require reliable simulations in order to avoid major damage to the device from disruptions. Disruptions are the sudden breakup of magnetic surfaces across the plasma volume -- a fast magnetic reconnection. This reconnection can be caused by the growth of perturbations outside of the plasma core causing an ideal perturbation to the core. This causes an increasing ratio of the maximum to the minimum separation, $\Delta_{max}/\Delta_{min}$, between neighboring magnetic surfaces. Magnetic reconnection becomes a dominant process when magnetic field lines can quickly interchange connections over a spatial scale $a_r$. This occurs when $\Delta_{max}/\Delta_{min}\gtrsim a_r/\Delta_d$, where $\Delta_d$ is the scale over which non-ideal effects make magnetic field lines indistinguishable. Traditional reconnection theory is fundamentally different. It is a study of the steady-state cancellation of oppositely directed magnetic field components across a thin layer. During more than sixty years, mathematical implications of Faraday's Law have been derived that clarify and constrain the physics of fast magnetic reconnection. These are reviewed because they are not commonly known but are needed to understand and to place in context how an ideal magnetic evolution can cause reconnection to quickly become a dominant process no matter how small $\Delta_d/a_r$ may be.


title: Foreground-Aware Stylization and Consensus Pseudo-Labeling for Domain Adaptation of First-Person Hand Segmentation

abstract: Hand segmentation is a crucial task in first-person vision. Since first-person images exhibit strong bias in appearance among different environments, adapting a pre-trained segmentation model to a new domain is required in hand segmentation. Here, we focus on appearance gaps for hand regions and backgrounds separately. We propose (i) foreground-aware image stylization and (ii) consensus pseudo-labeling for domain adaptation of hand segmentation. We stylize source images independently for the foreground and background using target images as style. To resolve the domain shift that the stylization has not addressed, we apply careful pseudo-labeling by taking a consensus between the models trained on the source and stylized source images. We validated our method on domain adaptation of hand segmentation from real and simulation images. Our method achieved state-of-the-art performance in both settings. We also demonstrated promising results in challenging multi-target domain adaptation and domain generalization settings. Code is available at https://github.com/ut-vision/FgSty-CPL.


title: Saturation-Aware Model Predictive Energy Management for Droop-Controlled Islanded Microgrids

abstract: In this paper, we propose a minimax model predictive control (MPC)-based energy management system that is robust with respect to uncertainties in renewable infeed and load. The MPC formulation includes a model of low-level droop control with saturation at the power and energy limits of the units. Robust MPC-based energy management systems tend to under-utilize the renewable energy sources to guarantee safe operation. In order to mitigate this effect, we further consider droop control of renewable energy sources. For a microgrid with droop-controlled units, we show that enhancing droop feedback with saturation enlarges the space of feasible control actions. However, the resulting controller requires to solve a mixed-integer problem with additional variables and equations representing saturation. We derive a computationally tractable formulation for this problem. Furthermore, we investigate the performance gained by using droop with saturation, renewable droop and combination of both in a case study.


title: Lexical Access Model for Italian -- Modeling human speech processing: identification of words in running speech toward lexical access based on the detection of landmarks and other acoustic cues to features

abstract: Modelling the process that a listener actuates in deriving the words intended by a speaker requires setting a hypothesis on how lexical items are stored in memory. This work aims at developing a system that imitates humans when identifying words in running speech and, in this way, provide a framework to better understand human speech processing. We build a speech recognizer for Italian based on the principles of Stevens' model of Lexical Access in which words are stored as hierarchical arrangements of distinctive features (Stevens, K. N. (2002). "Toward a model for lexical access based on acoustic landmarks and distinctive features," J. Acoust. Soc. Am., 111(4):1872-1891). Over the past few decades, the Speech Communication Group at the Massachusetts Institute of Technology (MIT) developed a speech recognition system for English based on this approach. Italian will be the first language beyond English to be explored; the extension to another language provides the opportunity to test the hypothesis that words are represented in memory as a set of hierarchically-arranged distinctive features, and reveal which of the underlying mechanisms may have a language-independent nature. This paper also introduces a new Lexical Access corpus, the LaMIT database, created and labeled specifically for this work, that will be provided freely to the speech research community. Future developments will test the hypothesis that specific acoustic discontinuities - called landmarks - that serve as cues to features, are language independent, while other cues may be language-dependent, with powerful implications for understanding how the human brain recognizes speech.


title: The topology of Gelfand-Zeitlin fibers

abstract: We prove several new results about the topology of fibers of Gelfand--Zeitlin systems on unitary and orthogonal coadjoint orbits, at the same time finding a unifying framework recovering and shedding light on essentially all known results. We find completely explicit descriptions of the diffeomorphism type of the fiber in many instances a direct factor decomposition of the fiber, and a torus factor corresponding to the action given by the Thimm trick. The new description also gives us a weak local normal form for a coadjoint orbit, which we use to define a topological toric degeneration, new in the orthogonal case. We also compute the first three homotopy groups (new in the orthogonal case) and cohomology rings of a fiber (new in both cases). All these descriptions can be read in a straightforward manner from the combinatorics of the associated Gelfand--Zeitlin pattern.


title: Cooling flows around cold clouds in the circumgalactic medium: steady-state models & comparison with TNG50

abstract: Cold, non-self-gravitating clumps occur in various astrophysical systems, ranging from the interstellar and circumgalactic medium (CGM), to AGN outflows and solar coronal loops. Cold gas has diverse origins such as turbulent mixing or precipitation from hotter phases. We obtain the analytic solution for a steady pressure-driven 1-D cooling flow around cold, local over-densities, irrespective of their origin. Our solutions describe the slow and steady radiative cooling-driven gas inflow in the saturated regime of nonlinear thermal instability in clouds, sheets and filaments. Such a cooling flow develops when the gas around small clumps undergoes radiative cooling. These small-scale, cold `seeds' are embedded in a large volume-filling hot CGM maintained by feedback. We use a simple two-fluid treatment to include magnetic fields as an additional polytropic fluid. To test the limits of applicability of these analytic solutions, we compare with the gas structure found in and around small-scale cold clouds in the CGM of massive halos in the TNG50 cosmological MHD simulation from the IllustrisTNG suite. Despite qualitative resemblance of the gas structure, we find deviations from steady state profiles generated by our model. Complex geometries and turbulence all add complexity beyond our analytic solutions. We derive an exact relation between the mass cooling rate ($\dot{\rm M}_{\rm cool}$) and the radiative cooling rate ($\dot{\rm E}_{\rm cool}$) for a steady cooling flow. A comparison with the TNG50 clouds shows that this cooling flow relation only applies in a narrow temperature range around $\rm \sim 10^{4.5}$ K where the isobaric cooling time is the shortest. In general, turbulence and mixing, instead of radiative cooling, may dominate the transition of gas between different temperature phases.


title: Non-monotonic dynamics in the onset of frictional slip

abstract: The transition from static to dynamic friction is often described as a fracture-like instantaneous slip. However, studies on slow sliding processes aimed at understanding frictional instabilities and earthquakes report slow friction transients that are usually explained by empirical rate-and-state formulations. We perform very slow ($\sim nm/s$) macroscopic-scale sliding experiments and show that the onset of frictional slip is governed by continuous non-monotonic dynamics originating from a competition between contact aging and shear-induced rejuvenation. This allows to describe both our non-monotonic dynamics and the simpler rate-and-state transients with a single evolution equation.


title: The proportion of derangements characterizes the symmetric and alternating groups

abstract: Let $G$ be a subgroup of the symmetric group $S_n$. If the proportion of fixed-point-free elements in $G$ (or a coset) equals the proportion of fixed-point-free elements in $S_n$, then $G=S_n$. The analogue for $A_n$ holds if $n \ge 7$. We give an application to monodromy groups.


title: Families of extensions of the Kantorovich-Rubinstein and Lipschitz norms

abstract: We propose a family of extensions of the Kantorovich-Rubinstein norm from the space of zero-charge countably additive measures on a compact metric space to the space of all countably additive measures, and a family of extensions of the Lipschitz norm from the quotient space of Lipschitz functions on a compact metric space to the space of all Lipschitz functions. These families are parameterized by $p,q \in [1,\infty]$, and if $p,q$ are H\"older conjugates, then the dual of the resulting $p$-Kantorovich space is isometrically isomorphic to the resulting $q$-Lipschitz space.


title: Distributed Adaptive Huber Regression

abstract: Distributed data naturally arise in scenarios involving multiple sources of observations, each stored at a different location. Directly pooling all the data together is often prohibited due to limited bandwidth and storage, or due to privacy protocols. This paper introduces a new robust distributed algorithm for fitting linear regressions when data are subject to heavy-tailed and/or asymmetric errors with finite second moments. The algorithm only communicates gradient information at each iteration and therefore is communication-efficient. Statistically, the resulting estimator achieves the centralized nonasymptotic error bound as if all the data were pooled together and came from a distribution with sub-Gaussian tails. Under a finite $(2+\delta)$-th moment condition, we derive a Berry-Esseen bound for the distributed estimator, based on which we construct robust confidence intervals. Numerical studies further confirm that compared with extant distributed methods, the proposed methods achieve near-optimal accuracy with low variability and better coverage with tighter confidence width.


title: Dependent finitely homogneneous rosy structures

abstract: We study finitely homogeneous dependent rosy structures, adapting results of Cherlin, Harrington, and Lachlan proved for $\omega$-stable $\omega$-categorical structures. In particular, we prove that such structures have finite {\th}-rank and are coordinatized by a {\th}-rank 1 set. We show that they admit a distal, finitely axiomatizable, expansion. These results show that there are, up to inter-definability, at most countably many dependent rosy structures M which are homogeneous in a finite relational language.


title: A quantum prediction as a collection of epistemically restricted classical predictions

abstract: Spekkens has introduced an epistemically restricted classical theory of discrete systems, based on discrete phase space. The theory manifests a number of quantum-like properties but cannot fully imitate quantum theory because it is noncontextual. In this paper we show how, for a certain class of quantum systems, the quantum description of an experiment can be decomposed into classical descriptions that are epistemically restricted, though in a different sense than in Spekkens' work. For each aspect of the experiment -- the preparation, the transformations, and the measurement -- the epistemic restriction limits the form of the probability distribution an imagined classical observer may use. There are also global constraints that the whole collection of classical descriptions must satisfy. Each classical description generates its own prediction regarding the outcome of the experiment. One recovers the quantum prediction via a simple but highly nonclassical rule: the "nonrandom part" of the predicted quantum probabilities is obtained by summing the nonrandom parts of the classically predicted probabilities. By "nonrandom part" we mean the deviation from complete randomness, that is, from what one would expect upon measuring the fully mixed state.


title: AdaRL: What, Where, and How to Adapt in Transfer Reinforcement Learning

abstract: One practical challenge in reinforcement learning (RL) is how to make quick adaptations when faced with new environments. In this paper, we propose a principled framework for adaptive RL, called \textit{AdaRL}, that adapts reliably and efficiently to changes across domains with a few samples from the target domain, even in partially observable environments. Specifically, we leverage a parsimonious graphical representation that characterizes structural relationships over variables in the RL system. Such graphical representations provide a compact way to encode what and where the changes across domains are, and furthermore inform us with a minimal set of changes that one has to consider for the purpose of policy adaptation. We show that by explicitly leveraging this compact representation to encode changes, we can efficiently adapt the policy to the target domain, in which only a few samples are needed and further policy optimization is avoided. We illustrate the efficacy of AdaRL through a series of experiments that vary factors in the observation, transition, and reward functions for Cartpole and Atari games.


title: A provable two-stage algorithm for penalized hazards regression

abstract: From an optimizer's perspective, achieving the global optimum for a general nonconvex problem is often provably NP-hard using the classical worst-case analysis. In the case of Cox's proportional hazards model, by taking its statistical model structures into account, we identify local strong convexity near the global optimum, motivated by which we propose to use two convex programs to optimize the folded-concave penalized Cox's proportional hazards regression. Theoretically, we investigate the statistical and computational tradeoffs of the proposed algorithm and establish the strong oracle property of the resulting estimators. Numerical studies and real data analysis lend further support to our algorithm and theory.


title: A two-mechanism and multiscale compatible approach for solid state electrolytes of (Li-ion) batteries

abstract: All solid state batteries are claimed to be the next-generation battery system, in view of their safety accompanied by high energy densities. A new advanced, multiscale compatible, and fully three dimensional model for solid electrolytes is presented in this note. The response of the electrolyte is profoundly studied theoretically and numerically, analyzing the equilibrium and steady state behaviors, the limiting factors, as well as the most relevant constitutive parameters according to the sensitivity analysis of the model.


title: Provable Lipschitz Certification for Generative Models

abstract: We present a scalable technique for upper bounding the Lipschitz constant of generative models. We relate this quantity to the maximal norm over the set of attainable vector-Jacobian products of a given generative model. We approximate this set by layerwise convex approximations using zonotopes. Our approach generalizes and improves upon prior work using zonotope transformers and we extend to Lipschitz estimation of neural networks with large output dimension. This provides efficient and tight bounds on small networks and can scale to generative models on VAE and DCGAN architectures.


title: The observed cosmic star formation rate density has an evolution which resembles a {\Gamma}(a, bt) distribution and can be described successfully by only 2 parameters

abstract: A debate is emerging regarding the recent inconsistent results of different studies for the Cosmic Star Formation Rate Density (CSFRD) at high-z. We employ UV and IR datasets to investigate the star formation rate function (SFRF) at ${\rm z \sim 0-9}$. We find that the SFRFs derived from the dust corrected ${\rm UV}$ (${\rm UV_{corr}}$) data contradict those from IR on some key issues since they are described by different distributions (Schechter vs double-power law), imply different physics for galaxy formation (${\rm UV_{corr}}$ data suggest a SFR limit/strong mechanism that diminish the number density of high star forming systems with respect IR) and compare differently with the stellar mass density evolution obtained from SED fitting (${\rm UV_{corr}}$ is in agreement, while IR in tension up to 0.5 dex). However, both tracers agree on a constant CSFRD evolution at ${\rm z \sim 1-4}$ and point to a plateau instead of a peak. In addition, using both indicators we demonstrate that the evolution of the {\it observed} CSFRD can be described by only {\bf 2} parameters and a function that has the form of a Gamma distribution (${\bf \Gamma(a,bt)}$). In contrast to previous parameterizations used in the literature our framework connects the parameters to physical properties like the star formation rate depletion time and cosmic baryonic gas density. The build up of stellar mass occurs in $\Gamma(a,bt)$ distributed steps and is the result of gas consumption up to the limit that there is no eligible gas for SF at t = ${\rm \infty}$, resulting to a final cosmic stellar mass density of $\sim 0.5 \times 10^9 \, {\rm \frac{M_{\odot}}{Mpc^3}}$.


title: Phase-Aberration Correction in Shear-wave Elastography Imaging Using Local Speed-of-Sound Adaptive Beamforming

abstract: Shear-wave Elastography Imaging (SWEI) is a noninvasive imaging modality that provides tissue elasticity information by measuring the travelling speed of an induced shear-wave. It is commercially available on clinical ultrasound scanners and popularly used in the diagnosis and staging of liver disease and breast cancer. In conventional SWEI methods, a sequence of acoustic radiation force (ARF) pushes are used for inducing a shear-wave, which is tracked using high frame-rate multi-angle plane wave imaging (MA-PWI) to estimate the shear-wave speed (SWS). Conventionally, these plane waves are beamformed using a constant speed-of-sound (SoS), assuming an a-priori known and homogeneous tissue medium. However, soft tissues are inhomogeneous, with intrinsic SoS variations. In this work, study the SoS effects and inhomogeneities on SWS estimation, using simulation and phantoms experiments with porcine muscle as an abbarator, and show how these aberrations can be corrected using local speed-of-sound adaptive beamforming. For shear-wave tracking, we compare standard beamform with spatially constant SoS values to software beamforming with locally varying SoS maps. We show that, given SoS aberrations, traditional beamforming using a constant SoS, regardless of the utilized SoS value, introduces a substantial bias in resulting SWS estimations. Average SWS estimation disparity for the same material was observed over 4.3 times worse when a constant SoS value is used compared to that when a known SoS map is used for beamforming. Such biases are shown to be corrected by using a local SoS map in beamforming, indicating the importance of and a need for local SoS reconstruction techniques.


title: Strong coupling constants of charmed and bottom mesons with light vector mesons in QCD sum rules

abstract: We estimate the strong coupling constants of charmed and bottom mesons $ D_{(s)}^* $, $ D_{(s)1} $, $ B_{(s)}^* $, and $ B_{(s)1} $ with light vector mesons $ \rho $, $ \omega $, $ K^* $, and $ \phi $ within the framework of light-cone QCD sum rules. We compare our estimations to the ones predicted by other approaches.


title: DEANN: Speeding up Kernel-Density Estimation using Approximate Nearest Neighbor Search

abstract: Kernel Density Estimation (KDE) is a nonparametric method for estimating the shape of a density function, given a set of samples from the distribution. Recently, locality-sensitive hashing, originally proposed as a tool for nearest neighbor search, has been shown to enable fast KDE data structures. However, these approaches do not take advantage of the many other advances that have been made in algorithms for nearest neighbor algorithms. We present an algorithm called Density Estimation from Approximate Nearest Neighbors (DEANN) where we apply Approximate Nearest Neighbor (ANN) algorithms as a black box subroutine to compute an unbiased KDE. The idea is to find points that have a large contribution to the KDE using ANN, compute their contribution exactly, and approximate the remainder with Random Sampling (RS). We present a theoretical argument that supports the idea that an ANN subroutine can speed up the evaluation. Furthermore, we provide a C++ implementation with a Python interface that can make use of an arbitrary ANN implementation as a subroutine for kernel density estimation. We show empirically that our implementation outperforms state of the art implementations in all high dimensional datasets we considered, and matches the performance of RS in cases where the ANN yield no gains in performance.


title: Quantum-based Molecular Dynamics Simulations Using Tensor Cores

abstract: Tensor cores, along with tensor processing units, represent a new form of hardware acceleration specifically designed for deep neural network calculations in artificial intelligence applications. Tensor cores provide extraordinary computational speed and energy efficiency, but with the caveat that they were designed for tensor contractions (matrix-matrix multiplications) using only low-precision floating point operations. In spite of this, we demonstrate how tensor cores can be applied with high efficiency to the challenging and numerically sensitive problem of quantum-based Born-Oppenheimer molecular dynamics, which requires highly accurate electronic structure optimizations and conservative force evaluations. The interatomic forces are calculated on-the-fly from an electronic structure that is obtained from a generalized deep neural network, where the computational structure naturally takes advantage of the exceptional processing power of the tensor cores and allows for high performance in excess of 100 Tflops on the tensor cores of a single Nvidia A100 GPU. Stable molecular dynamics trajectories are generated using the framework of extended Lagrangian Born-Oppenheimer molecular dynamics, which combines computational efficiency with long-term stability, even when using approximate charge relaxations and force evaluations that are limited in accuracy by the numerically noisy conditions caused by the low precision tensor core floating-point operations. A canonical ensemble simulation scheme is also presented, where the additional numerical noise in the calculated forces is absorbed into a Langevin-like dynamics.


title: Dueling Bandits with Team Comparisons

abstract: We introduce the dueling teams problem, a new online-learning setting in which the learner observes noisy comparisons of disjoint pairs of $k$-sized teams from a universe of $n$ players. The goal of the learner is to minimize the number of duels required to identify, with high probability, a Condorcet winning team, i.e., a team which wins against any other disjoint team (with probability at least $1/2$). Noisy comparisons are linked to a total order on the teams. We formalize our model by building upon the dueling bandits setting (Yue et al.2012) and provide several algorithms, both for stochastic and deterministic settings. For the stochastic setting, we provide a reduction to the classical dueling bandits setting, yielding an algorithm that identifies a Condorcet winning team within $\mathcal{O}((n + k \log (k)) \frac{\max(\log\log n, \log k)}{\Delta^2})$ duels, where $\Delta$ is a gap parameter. For deterministic feedback, we additionally present a gap-independent algorithm that identifies a Condorcet winning team within $\mathcal{O}(nk\log(k)+k^5)$ duels.


title: Shapes as Product Differentiation: Neural Network Embedding in the Analysis of Markets for Fonts

abstract: Many differentiated products have key attributes that are unstructured and thus high-dimensional (e.g., design, text). Instead of treating unstructured attributes as unobservables in economic models, quantifying them can be important to answer interesting economic questions. To propose an analytical framework for these types of products, this paper considers one of the simplest design products-fonts-and investigates merger and product differentiation using an original dataset from the world's largest online marketplace for fonts. We quantify font shapes by constructing embeddings from a deep convolutional neural network. Each embedding maps a font's shape onto a low-dimensional vector. In the resulting product space, designers are assumed to engage in Hotelling-type spatial competition. From the image embeddings, we construct two alternative measures that capture the degree of design differentiation. We then study the causal effects of a merger on the merging firm's creative decisions using the constructed measures in a synthetic control method. We find that the merger causes the merging firm to increase the visual variety of font design. Notably, such effects are not captured when using traditional measures for product offerings (e.g., specifications and the number of products) constructed from structured data.


title: High-speed imaging of spatiotemporal correlations in Hong-Ou-Mandel interference

abstract: The Hong-Ou-Mandel interference effect lies at the heart of many emerging quantum technologies whose performance can be significantly enhanced with increasing numbers of entangled modes one could measure and thus utilize. Photon pairs generated through the process of spontaneous parametric down conversion are known to be entangled in a vast number of modes in the various degrees of freedom (DOF) the photons possess such as time, energy, and momentum, etc. Due to limitations in detection technology and techniques, often only one such DOFs can be effectively measured at a time, resulting in much lost potential. Here, we experimentally demonstrate, with the aid of a time tagging camera, high speed measurement and characterization of two-photon interference. {With a data acquisition time of only a few seconds, we observe a bi-photon interference and coalescence visibility of $\sim64\%$ with potentially up to $\sim2\times10^3$ spatial modes}. These results open up a route for practical applications of using the high dimensionality of spatiotemporal DOF in two-photon interference, and in particular, for quantum sensing and communication.


title: Invariant Variational Schemes for Ordinary Differential Equations

abstract: We propose a novel algorithmic method for constructing invariant variational schemes of systems of ordinary differential equations that are the Euler-Lagrange equations of a variational principle. The method is based on the invariantization of standard, non-invariant discrete Lagrangian functionals using equivariant moving frames. The invariant variational schemes are given by the Euler-Lagrange equations of the corresponding invariantized discrete Lagrangian functionals. We showcase this general method by constructing invariant variational schemes of ordinary differential equations that preserve variational and divergence symmetries of the associated continuous Lagrangians. Noether's theorem automatically implies that the resulting schemes are exactly conservative. Numerical simulations are carried out and show that these invariant variational schemes outperform standard numerical discretizations.


title: How Big Should Your Data Really Be? Data-Driven Newsvendor: Learning One Sample at a Time

abstract: We study the classical newsvendor problem in which the decision-maker must trade-off underage and overage costs. In contrast to the typical setting, we assume that the decision-maker does not know the underlying distribution driving uncertainty but has only access to historical data. In turn, the key questions are how to map existing data to a decision and what type of performance to expect as a function of the data size. We analyze the classical setting with access to past samples drawn from the distribution (e.g., past demand), focusing not only on asymptotic performance but also on what we call the transient regime of learning, i.e., performance for arbitrary data sizes. We evaluate the performance of any algorithm through its worst-case relative expected regret, compared to an oracle with knowledge of the distribution. We provide the first finite sample exact analysis of the classical Sample Average Approximation (SAA) algorithm for this class of problems across all data sizes. This allows to uncover novel fundamental insights on the value of data: it reveals that tens of samples are sufficient to perform very efficiently but also that more data can lead to worse out-of-sample performance for SAA. We then focus on the general class of mappings from data to decisions without any restriction on the set of policies and derive an optimal algorithm (in the minimax sense) as well as characterize its associated performance. This leads to significant improvements for limited data sizes, and allows to exactly quantify the value of historical information.


title: Submodular Order Functions and Assortment Optimization

abstract: We define a new class of set functions that in addition to being monotone and subadditive, also admit a very limited form of submodularity defined over a permutation of the ground set. We refer to this permutation as a submodular order. This class of functions includes monotone submodular functions as a sub-family. We give fast algorithms with strong approximation guarantees for maximizing submodular order functions under a variety of constraints and show a nearly tight upper bound on the highest approximation guarantee achievable by algorithms with polynomial query complexity. Applying this new notion to the problem of constrained assortment optimization in fundamental choice models, we obtain new algorithms that are both faster and have stronger approximation guarantees (in some cases, first algorithm with constant factor guarantee). We also show an intriguing connection to the maximization of monotone submodular functions in the streaming model, where we recover best known approximation guarantees as a corollary of our results.


title: Neural Computing

abstract: This chapter aims to provide next-level understanding of the problems of the world and the solutions available to those problems, which lie very well within the domain of neural computing, and at the same time are intelligent in their approach, to invoke a sense of innovation among the educationalists, researchers, academic professionals, students and people concerned, by highlighting the work done by major researchers and innovators in this field and thus, encouraging the readers to develop newer and more advanced techniques for the same. By means of this chapter, the societal problems are discussed and various solutions are also given by means of the theories presented and researches done so far. Different types of neural networks discovered so far and applications of some of those neural networks are focused on, apart from their theoretical understanding, the working and core concepts involved in the applications.


title: Thermal monopole condensation in QCD with physical quark masses

abstract: Thermal monopoles, identified after Abelian projection as magnetic currents wrapping non-trivially around the thermal circle, are studied in $N_f = 2+1$ QCD at the physical point. The distribution in the number of wrappings, which in pure gauge theories points to a condensation temperature coinciding with deconfinement, points in this case to around 275 MeV, almost twice the QCD crossover temperature $T_c$; similar indications emerge looking for the formation of a percolating current cluster. The possible relation with other non-perturbative phenomena observed above $T_c$ is discussed.


title: Kepler K2 Campaign 9: I. Candidate short-duration events from the first space-based survey for planetary microlensing

abstract: We present the first short-duration candidate microlensing events from the Kepler K2 mission. From late April to early July 2016, Campaign 9 of K2 obtained high temporal cadence observations over a 3.7 square degree region of the Galactic bulge. Its primary objectives were to look for evidence of a free-floating planet (FFP) population using microlensing, and demonstrate the feasibility of space-based planetary microlensing surveys. Though Kepler K2 is far from optimal for microlensing, the recently developed MCPM photometric pipeline enables us to identify and model microlensing events. We describe our blind event-selection pipeline in detail and use it to recover 22 short-duration events with effective timescales of less than 10 days previously announced by the OGLE and KMTNet ground-based surveys. We also announce five new candidate events. One of these is a caustic-crossing binary event, consistent with a bound planet and modelled as such in a companion study. The other four have very short durations (effective timescales less than 0.1 days) typical of an Earth-mass FFP population. Whilst Kepler was not designed for crowded-field photometry, the K2C9 dataset clearly demonstrates the feasibility of conducting blind space-based microlensing surveys towards the Galactic bulge.


title: Nonlinear mirror image method for nonlinear Schr\"odinger equation: Absorption/emission of one soliton by a boundary

abstract: We perform the analysis of the focusing nonlinear Schr\"odinger equation on the half-line with time-dependent boundary conditions along the lines of the nonlinear method of images with the help of B\"acklund transformations. The difficulty arising from having such time-dependent boundary conditions at $x=0$ is overcome by changing the viewpoint of the method and fixing the B\"acklund transformation at infinity as well as relating its value at $x=0$ to a time-dependent reflection matrix. The interplay between the various aspects of integrable boundary conditions is reviewed in detail to brush a picture of the area. We find two possible classes of solutions. One is very similar to the case of Robin boundary conditions whereby solitons are reflected at the boundary, as a result of an effective interaction with their images on the other half-line. The new regime of solutions supports the existence of one soliton that is not reflected at the boundary but can be either absorbed or emitted by it. We demonstrate that this is a unique feature of time-dependent integrable boundary conditions.


title: MAJORITY-3SAT (and Related Problems) in Polynomial Time

abstract: Majority-SAT is the problem of determining whether an input $n$-variable formula in conjunctive normal form (CNF) has at least $2^{n-1}$ satisfying assignments. Majority-SAT and related problems have been studied extensively in various AI communities interested in the complexity of probabilistic planning and inference. Although Majority-SAT has been known to be PP-complete for over 40 years, the complexity of a natural variant has remained open: Majority-$k$SAT, where the input CNF formula is restricted to have clause width at most $k$. We prove that for every $k$, Majority-$k$SAT is in P. In fact, for any positive integer $k$ and rational $\rho \in (0,1)$ with bounded denominator, we give an algorithm that can determine whether a given $k$-CNF has at least $\rho \cdot 2^n$ satisfying assignments, in deterministic linear time (whereas the previous best-known algorithm ran in exponential time). Our algorithms have interesting positive implications for counting complexity and the complexity of inference, significantly reducing the known complexities of related problems such as E-MAJ-$k$SAT and MAJ-MAJ-$k$SAT. At the heart of our approach is an efficient method for solving threshold counting problems by extracting sunflowers found in the corresponding set system of a $k$-CNF. We also show that the tractability of Majority-$k$SAT is somewhat fragile. For the closely related GtMajority-SAT problem (where we ask whether a given formula has greater than $2^{n-1}$ satisfying assignments) which is known to be PP-complete, we show that GtMajority-$k$SAT is in P for $k\le 3$, but becomes NP-complete for $k\geq 4$. These results are counterintuitive, because the ``natural'' classifications of these problems would have been PP-completeness, and because there is a stark difference in the complexity of GtMajority-$k$SAT and Majority-$k$SAT for all $k\ge 4$.


title: Spatial Modulation Search Applied to the Search and Confirmation of Highly Scintillated Pulsars at FAST with A Pulsar Discovered in M3

abstract: We present a pulsar candidate identification and confirmation procedure based on a position-switch mode during the pulsar search observations. This method enables the simultaneous search and confirmation of a pulsar in a single observation, by utilizing the different spatial features of a pulsar signal and a radio frequency interference (RFI). Based on this method, we performed test pulsar search observations in globular clusters M3, M15, and M92. We discovered and confirmed a new pulsar, M3F, and detected the known pulsars M3B, M15 A to G (except C), and M92A.


title: Benchmarking (multi)wavelet-based dynamic and static non-uniform grid solvers for flood inundation modelling

abstract: This paper explores static non-uniform grid solvers that adapt three raster-based flood models on an optimised non-uniform grid: the second-order discontinuous Galerkin (DG2) model representing the modelled data as piecewise-planar fields, the first-order finite volume (FV1) model using piecewise-constant fields, and the local inertial (ACC) model only evolving piecewise-constant water depth fields. The optimised grid is generated by applying the multiresolution analysis (MRA) of multiwavelets (MWs) to piecewise-planar representation of raster-formatted topography data, for more sensible grid coarsening based on one user-specified parameter. Two adaptive solvers are also explored that apply the MRA of MWs and of Haar wavelets (HWs) to, respectively, scale and adapt the DG2 (MWDG2) and FV1 (HWFV1) modelled data dynamically in time. The performance of the non-uniform grid and adaptive solvers is assessed in terms of flood depth and extent, velocities, and CPU runtimes, with reference to the raster-based DG2 model predictions on their finest resolution grid. The assessments considered three large-scale flooding scenarios, involving rapid and slow-to-gradual flows. MWDG2 is found to be the most favourable choice when modelling rapid flows, where it excels in capturing small velocity variations. For slow-to-gradual flows, the adaptive solvers deliver less accurate outcomes, and their efficiency can be hampered by overhead costs of the dynamic MRA. Instead, non-uniform DG2 is recommended to capture urban flow interactions more accurately. Non-uniform ACC is 5 times faster to run than non-uniform DG2 but delivers close flooding depth and extent predictions, thus is more attractive for fluvial/pluvial flood simulation over large areas.


title: Quantum Annealing Formulation for Binary Neural Networks

abstract: Quantum annealing is a promising paradigm for building practical quantum computers. Compared to other approaches, quantum annealing technology has been scaled up to a larger number of qubits. On the other hand, deep learning has been profoundly successful in pushing the boundaries of AI. It is thus natural to investigate potentially game changing technologies such as quantum annealers to augment the capabilities of deep learning. In this work, we explore binary neural networks, which are lightweight yet powerful models typically intended for resource constrained devices. Departing from current training regimes for binary networks that smooth/approximate the activation functions to make the network differentiable, we devise a quadratic unconstrained binary optimization formulation for the training problem. While the problem is intractable, i.e., the cost to estimate the binary weights scales exponentially with network size, we show how the problem can be optimized directly on a quantum annealer, thereby opening up to the potential gains of quantum computing. We experimentally validated our formulation via simulation and testing on an actual quantum annealer (D-Wave Advantage), the latter to the extent allowable by the capacity of current technology.


title: Theory of Probability and Statistics as Exemplified in Short Dictums

abstract: I am presenting a first-ever scientific collection of short sayings on probability and statistics expressed by most various men of science, many classics included, from antiquity to Kepler to our time. Quite understandably, the reader will find here no mathematical formulas and in some instances he will miss a worthy subject. Markov chains provide a good example: their inventor had not said anything about them suitable for my goal. Nevertheless, the scope of the collected dictums is amazingly broad which reflects both the discussions concerning the lack of a solid foundation of probability theory until the 1930s and the great extent of applications of probability and statistics. And I have also included two related and most important topics, the treatment of observations and randomness.


title: Machine Learning for Network-based Intrusion Detection Systems: an Analysis of the CIDDS-001 Dataset

abstract: With the increasing amount of reliance on digital data and computer networks by corporations and the public in general, the occurrence of cyber attacks has become a great threat to the normal functioning of our society. Intrusion detection systems seek to address this threat by preemptively detecting attacks in real time while attempting to block them or minimizing their damage. These systems can function in many ways being some of them based on artificial intelligence methods. Datasets containing both normal network traffic and cyber attacks are used for training these algorithms so that they can learn the underlying patterns of network-based data. The CIDDS-001 is one of the most used datasets for network-based intrusion detection research. Regarding this dataset, in the majority of works published so far, the Class label was used for training machine learning algorithms. However, there is another label in the CIDDS-001, AttackType, that seems very promising for this purpose and remains considerably unexplored. This work seeks to make a comparison between two machine learning models, K-Nearest Neighbours and Random Forest, which were trained with both these labels in order to ascertain whether AttackType can produce reliable results in comparison with the Class label.


title: Frequency Shift Algorithm: Design of a Baseband Phase Locked Loop for Frequency-Domain Multiplexing Readout of X-ray Transition-Edge Sensor Microcalorimeters

abstract: The Transition-Edge Sensor (TES) is an extremely sensitive device which is used to measure the energy of individual X-ray photons. For astronomical spectrometry applications, SRON develops a Frequency Domain Multiplexing (FDM) read-out system for kilopixel arrays of such TESs. Each TES is voltage biased at a specific frequency in the range 1 to 5 MHz. Isolation between the individual pixels is obtained through very narrow-band (high-Q) lithographic LC resonators. To prevent energy resolution degradation due to intermodulation line noise, the bias frequencies are distributed on a regular grid. The requirements on the accuracy of the LC resonance frequency are very high. The deviation of the resonance frequencies due to production tolerances is significant with respect to the bandwidth, and a controller is necessary to compensate for the LC series impedance. We present two such controllers: a simple orthogonal proportional-integrating (PI) controller and a more complex impedance estimator. Both controllers operate in baseband and try to make the TES current in-phase with the bias voltage, effectively operating as phase-locked loops (PLL). They allow off-LC-resonance operation of the TES pixels, while preserving TES thermal response and energy resolution. Extensive experimental results -- published in a companion paper recently -- with the proposed methods, show that these controllers allow the preservation of single pixel energy resolution in multiplexed operation.

categories: physics.ins-det astro-ph.IM cs.SY eess.SY


title: FedFog: Network-Aware Optimization of Federated Learning over Wireless Fog-Cloud Systems

abstract: Federated learning (FL) is capable of performing large distributed machine learning tasks across multiple edge users by periodically aggregating trained local parameters. To address key challenges of enabling FL over a wireless fog-cloud system (e.g., non-i.i.d. data, users' heterogeneity), we first propose an efficient FL algorithm based on Federated Averaging (called FedFog) to perform the local aggregation of gradient parameters at fog servers and global training update at the cloud. Next, we employ FedFog in wireless fog-cloud systems by investigating a novel network-aware FL optimization problem that strikes the balance between the global loss and completion time. An iterative algorithm is then developed to obtain a precise measurement of the system performance, which helps design an efficient stopping criteria to output an appropriate number of global rounds. To mitigate the straggler effect, we propose a flexible user aggregation strategy that trains fast users first to obtain a certain level of accuracy before allowing slow users to join the global training updates. Extensive numerical results using several real-world FL tasks are provided to verify the theoretical convergence of FedFog. We also show that the proposed co-design of FL and communication is essential to substantially improve resource utilization while achieving comparable accuracy of the learning model.


title: Classification in chains of three-dimensional real evolution algebras

abstract: A chain of evolution algebras (CEA) is an uncountable family (depending on time) of evolution algebras on the field of real numbers. The matrix of structural constants of a CEA satisfies Kolmogorov-Chapman equation. In this paper, we consider three CEAs of three-dimensional real evolution algebras. These CEAs depend on several (non-zero) functions defined on the set of time. For each chain we give full classification (up to isomorphism) of the algebras depending on the time-parameter. We find concrete functions ensuring that the corresponding CEA contains all possible three-dimensional evolution algebras.


title: Sawtooth Factorial Topic Embeddings Guided Gamma Belief Network

abstract: Hierarchical topic models such as the gamma belief network (GBN) have delivered promising results in mining multi-layer document representations and discovering interpretable topic taxonomies. However, they often assume in the prior that the topics at each layer are independently drawn from the Dirichlet distribution, ignoring the dependencies between the topics both at the same layer and across different layers. To relax this assumption, we propose sawtooth factorial topic embedding guided GBN, a deep generative model of documents that captures the dependencies and semantic similarities between the topics in the embedding space. Specifically, both the words and topics are represented as embedding vectors of the same dimension. The topic matrix at a layer is factorized into the product of a factor loading matrix and a topic embedding matrix, the transpose of which is set as the factor loading matrix of the layer above. Repeating this particular type of factorization, which shares components between adjacent layers, leads to a structure referred to as sawtooth factorization. An auto-encoding variational inference network is constructed to optimize the model parameter via stochastic gradient descent. Experiments on big corpora show that our models outperform other neural topic models on extracting deeper interpretable topics and deriving better document representations.


title: $a_1(1260)$-meson longitudinal twist-2 distribution amplitude and the $D\to a_1(1260)\ell^+\nu_\ell$ decay processes

abstract: In the paper, we investigate the moments $\langle\xi_{2;a_1}^{\|;n}\rangle$ of the axial-vector $a_1(1260)$-meson distribution amplitude by using the QCD sum rules approach under the background field theory. By considering the vacuum condensates up to dimension-six and the perturbative part up to next-to-leading order QCD corrections, its first five moments at an initial scale $\mu_0=1~{\rm GeV}$ are $\langle\xi_{2;a_1}^{\|;2}\rangle|_{\mu_0} = 0.223 \pm 0.029$, $\langle\xi_{2;a_1}^{\|;4}\rangle|_{\mu_0} = 0.098 \pm 0.008$, $\langle\xi_{2;a_1}^{\|;6}\rangle|_{\mu_0} = 0.056 \pm 0.006$, $\langle\xi_{2;a_1}^{\|;8}\rangle|_{\mu_0} = 0.039 \pm 0.004$ and $\langle\xi_{2;a_1}^{\|;10}\rangle|_{\mu_0} = 0.028 \pm 0.003$, respectively. We then construct a light-cone harmonic oscillator model for $a_1(1260)$-meson longitudinal twist-2 distribution amplitude $\phi_{2;a_1}^{\|}(x,\mu)$, whose model parameters are fitted by using the least squares method. As an application of $\phi_{2;a_1}^{\|}(x,\mu)$, we calculate the transition form factors (TFFs) of $D\to a_1(1260)$ in large and intermediate momentum transfers by using the QCD light-cone sum rules approach. At the largest recoil point ($q^2=0$), we obtain $ A(0) = 0.130_{ - 0.013}^{ + 0.015}$, $V_1(0) = 1.898_{-0.121}^{+0.128}$, $V_2(0) = 0.228_{-0.021}^{ + 0.020}$, and $V_0(0) = 0.217_{ - 0.025}^{ + 0.023}$. By applying the extrapolated TFFs to the semi-leptonic decay $D^{0(+)} \to a_1^{-(0)}(1260)\ell^+\nu_\ell$, we obtain ${\cal B}(D^0\to a_1^-(1260) e^+\nu_e) = (5.261_{-0.639}^{+0.745}) \times 10^{-5}$, ${\cal B}(D^+\to a_1^0(1260) e^+\nu_e) = (6.673_{-0.811}^{+0.947}) \times 10^{-5}$, ${\cal B}(D^0\to a_1^-(1260) \mu^+ \nu_\mu)=(4.732_{-0.590}^{+0.685}) \times 10^{-5}$, ${\cal B}(D^+ \to a_1^0(1260) \mu^+ \nu_\mu)=(6.002_{-0.748}^{+0.796}) \times 10^{-5}$.


title: Gender Recognition in Informal and Formal Language Scenarios via Transfer Learning

abstract: The interest in demographic information retrieval based on text data has increased in the research community because applications have shown success in different sectors such as security, marketing, heath-care, and others. Recognition and identification of demographic traits such as gender, age, location, or personality based on text data can help to improve different marketing strategies. For instance it makes it possible to segment and to personalize offers, thus products and services are exposed to the group of greatest interest. This type of technology has been discussed widely in documents from social media. However, the methods have been poorly studied in data with a more formal structure, where there is no access to emoticons, mentions, and other linguistic phenomena that are only present in social media. This paper proposes the use of recurrent and convolutional neural networks, and a transfer learning strategy for gender recognition in documents that are written in informal and formal languages. Models are tested in two different databases consisting of Tweets and call-center conversations. Accuracies of up to 75\% are achieved for both databases. The results also indicate that it is possible to transfer the knowledge from a system trained on a specific type of expressions or idioms such as those typically used in social media into a more formal type of text data, where the amount of data is more scarce and its structure is completely different.


title: Additivity of $n$-Multiplicative Mappings of Gamma Rings

abstract: In this paper, we address the additivity of $n$-multiplicative isomorphisms and $n$-multiplicative derivations on Gamma rings. We proved that, if $\M$ is a $\Gamma$-ring satisfying the some conditions, then any $n$-multiplicative isomorphism $\left(\varphi, \phi\right)$ of $\M$ onto an arbitrary gamma ring is additive.


title: Spin optics for gravitational waves

abstract: We present the geometric optics expansion for circularly polarized gravitational waves on a curved spacetime background, to subleading order. We call spin optics to the subleading order geometric optics expansion, which involves modifying the standard eikonal function by including a specially chosen helicity-dependent correction. We show that the techniques developed for the propagation of electromagnetic waves can also be applied to gravitational waves in the limit of spin optics. However, one needs to account for the difference in the photon and graviton helicity, which we do here.


title: Area-Delay-Efficeint FPGA Design of 32-bit Euclid's GCD based on Sum of Absolute Difference

abstract: Euclids algorithm is widely used in calculating of GCD (Greatest Common Divisor) of two positive numbers. There are various fields where this division is used such as channel coding, cryptography, and error correction codes. This makes the GCD a fundamental algorithm in number theory, so a number of methods have been discovered to efficiently compute it. The main contribution of this paper is to investigate a method that computes the GCD of two 32-bit numbers based on Euclidean algorithm which targets six different Xilinx chips. The complexity of this method that we call Optimized_GCDSAD is achieved by utilizing Sum of Absolute Difference (SAD) block which is based on a fast carry-out generation function. The efficiency of the proposed architecture is evaluated based on criteria such as time (latency), area delay product (ADP) and space (slice number) complexity. The VHDL codes of these architectures have been implemented and synthesized through ISE 14.7. A detailed comparative analysis indicates that the proposed Optimized_GCDSAD method based on SAD block outperforms previously known results.


title: Predicting Surface Heat Flux on Complex Systems via Conv-LSTM

abstract: Existing algorithms with iterations as the principle for 3D inverse heat conduction problems (IHCPs) are usually time-consuming. With the recent advancements in deep learning techniques, it is possible to apply the neural network to compute IHCPs. In this paper, a new framework based on Convolutional-LSTM is introduced to predict the transient heat flux via measured temperature. The inverse heat conduction models concerned in this work have 3D complex structures with non-linear boundary conditions and thermophysical parameters. In order to reach high precision, a forward solver based on the finite element method is utilized to generate sufficient data for training. The fully trained framework can provide accurate predictions efficiently once the measured temperature and models are acquired. It is believed that the proposed framework offers a new pattern for real-time heat flux inversion.


title: Collaborative Insurance Sustainability and Network Structure

abstract: The peer-to-peer (P2P) economy has been growing with the advent of the Internet, with well known brands such as Uber or Airbnb being examples thereof. In the insurance sector the approach is still in its infancy, but some companies have started to explore P2P-based collaborative insurance products (eg. Lemonade in the U.S. or Inspeer in France). The actuarial literature only recently started to consider those risk sharing mechanisms, as in Denuit and Robert (2021) or Feng et al. (2021). In this paper, describe and analyse such a P2P product, with some reciprocal risk sharing contracts. Here, we consider the case where policyholders still have an insurance contract, but the first self-insurance layer, below the deductible, can be shared with friends. We study the impact of the shape of the network (through the distribution of degrees) on the risk reduction. We consider also some optimal setting of the reciprocal commitments, and discuss the introduction of contracts with friends of friends to mitigate some possible drawbacks of having people without enough connections to exchange risks.

categories: q-fin.RM cs.SI econ.GN q-fin.CP q-fin.EC


title: Exploring the Scope of Using News Articles to Understand Development Patterns of Districts in India

abstract: Understanding what factors bring about socio-economic development may often suffer from the streetlight effect, of analyzing the effect of only those variables that have been measured and are therefore available for analysis. How do we check whether all worthwhile variables have been instrumented and considered when building an econometric development model? We attempt to address this question by building unsupervised learning methods to identify and rank news articles about diverse events occurring in different districts of India, that can provide insights about what may have transpired in the districts. This can help determine whether variables related to these events are indeed available or not to model the development of these districts. We also describe several other applications that emerge from this approach, such as to use news articles to understand why pairs of districts that may have had similar socio-economic indicators approximately ten years back ended up at different levels of development currently, and another application that generates a newsfeed of unusual news articles that do not conform to news articles about typical districts with a similar socio-economic profile. These applications outline the need for qualitative data to augment models based on quantitative data, and are meant to open up research on new ways to mine information from unstructured qualitative data to understand development.


title: Magnetic fields in the accretion disks for various inner boundary conditions

abstract: The magnetic fields of accretion disks play an important role in studying their evolution. We may assume that its generation is connected to the dynamo mechanism, which is similar with that in the galactic disks. Here, we propose a model of the magnetic field of the accretion disk that uses the same approaches that have been used for galaxies. It is necessary to obtain the field, which is expected to be less than the equipartition value, and without destroying the disk. To do so, it is necessary to formulate the basic properties of the ionized medium and to estimate the parameters governing the dynamo. We used the no-z approximation that has been developed for thin disks. We also take different boundary conditions that can change the value of the field significantly. We show that the magnetic field strictly depends on the boundary conditions. Taking zero conditions and the fixed magnetic field condition on the inner boundary, which are connected to the physical properties of the accretion disk, we can avoid solutions that are greater than the equipartition field.


title: Inferring the instability of a dynamical system from the skill of data assimilation exercises

abstract: Data assimilation (DA) aims at optimally merging observational data and model outputs to create a coherent statistical and dynamical picture of the system under investigation. Indeed, DA aims at minimizing the effect of observational and model error, and at distilling the correct ingredients of its dynamics. DA is of critical importance for the analysis of systems featuring sensitive dependence on the initial conditions, as chaos wins over any finitely accurate knowledge of the state of the system, even in absence of model error. Clearly, the skill of DA is guided by the properties of dynamical system under investigation, as merging optimally observational data and model outputs is harder when strong instabilities are present. In this paper we reverse the usual angle on the problem and show that it is indeed possible to use the skill of DA to infer some basic properties of the tangent space of the system, which may be hard to compute in very high-dimensional systems. Here, we focus our attention on the first Lyapunov exponent and the Kolmogorov-Sinai entropy, and perform numerical experiments on the Vissio-Lucarini 2020 model, a recently proposed generalisation of the Lorenz 1996 model that is able to describe in a simple yet meaningful way the interplay between dynamical and thermodynamical variables.

categories: physics.data-an nlin.CD physics.ao-ph physics.comp-ph


title: Equi-Lipschitz minimizing trajectories for non coercive, discontinuous, non convex Bolza controlled-linear optimal control problems

abstract: This article deals with the Lipschitz regularity of the ''approximate`` minimizers for the Bolza type control functional of the form \[J_t(y,u):=\int_t^T\Lambda(s,y(s), u(s))\,ds+g(y(T))\] among the pairs $(y,u)$ satisfying a prescribed initial condition $y(t)=x$, where the state $y$ is absolutely continuous, the control $u$ is summable and the dynamic is controlled-linear of the form $y'=b(y)u$. For $b\equiv 1$ the above becomes a problem of the calculus of variations. The Lagrangian $\Lambda(s,y,u)$ is assumed to be either convex in the variable $u$ on every half-line from the origin (radial convexity in $u$), or partial differentiable in the control variable and satisfies a local Lipschitz regularity on the time variable, named Condition (S). It is allowed to be extended valued, discontinuous in $y$ or in $u$, and non convex in $u$.\\ We assume a very mild growth condition, that is fulfilled if the Lagrangian is coercive as well as in some almost linear cases. The main result states that, given any admissible pair $(y,u)$, there exists a more convenient admissible pair $(\overline y, \overline u)$ for $J_t$ where $\overline u$ is bounded, $\overline y$ is Lipschitz, with bounds and ranks that are uniform with respect to $t,x$ in the compact subsets of $[0,T[\times\mathbb R^n$. The result is new even in the superlinear case. As a consequence, there are minimizing sequences that are formed by pairs of equi-Lipschitz trajectories and equi--bounded controls.\\ A new existence and regularity result follows without assuming any kind of Lipschitzianity in the state variable.\\ We deduce, without any need of growth conditions, the nonoccurrence of the Lavrentiev phenomenon for a wide class of Lagrangians containing those that satisfy Condition (S), are bounded on bounded sets ``well'' inside the effective domain and are radially convex in the control variable.


title: Exploring a Dynamic Ring without Landmark

abstract: Consider a group of autonomous mobile computational entities, called agents, arbitrarily placed at some nodes of a dynamic but always connected ring. The agents neither have any knowledge about the size of the ring nor have a common notion of orientation. We consider the \textsc{Exploration} problem where the agents have to collaboratively to explore the graph and terminate, with the requirement that each node has to be visited by at least one agent. It has been shown by Di Luna et al. [Distrib. Comput. 2020] that the problem is solvable by two anonymous agents if there is a single observably different node in the ring called landmark node. The problem is unsolvable by any number of anonymous agents in absence of a landmark node. We consider the problem with non-anonymous agents (agents with distinct identifiers) in a ring with no landmark node. The assumption of agents with distinct identifiers is strictly weaker than having a landmark node as the problem is unsolvable by two agents with distinct identifiers in absence of a landmark node. This setting has been recently studied by Mandal et al. [ALGOSENSORS 2020]. There it is shown that the problem is solvable in this setting by three agents assuming that they have edge crossing detection capability. Edge crossing detection capability is a strong assumption which enables two agents moving in opposite directions through an edge in the same round to detect each other and also exchange information. In this paper we give an algorithm that solves the problem with three agents without the edge crossing detection capability.


title: Neutron Imaging of Paramagnetic Ions: Electrosorption by Carbon Aerogels and Macroscopic Magnetic Forces

abstract: The electrosorption of Gd$^{3+}$ ions from aqueous 70 mM Gd(NO$_3$)$_3$ solution in monolithic carbon aerogel electrodes was recorded by dynamic neutron imaging. The aerogels have a bimodal pore size distribution consisting of macropores and mesopores centered at 115 nm and 15 nm, respectively. After the uptake of Gd$^{3+}$ ions by the negatively charged surface of the porous structure, an inhomogeneous magnetic field was applied to the system of discharging electrodes. This led to a convective flow and confinement of Gd(NO$_3$)$_3$ solution in the magnetic field gradient. Thus, a way to desalt and capture paramagnetic ions from an initially homogeneous solution is established.


title: Astrometric and photometric observations of six brightest trans-Neptunian objects at the Kyiv comet station

abstract: In this work, we focused on observations of six trans-Neptunian objects (TNOs) whose apparent magnitudes are brighter than 20 mag. We present the results of astrometric and photometric observations of (134340) Pluto, (136108) Haumea, (136472) Makemake, (136199) Eris, (90482) Orcus and (20000) Varuna obtained at the Kyiv comet station (Code MPC 585) in 2017-2019. For observations, we used the 0.7-m reflector AZT-8 with the FLI PL4710 CCD camera and filters of the Johnson-Cousins photometric system. From our images, we measured the objects astrometric positions, calculated apparent magnitudes in the BVRI (mostly R) bands using the aperture photometry method, and found the absolute magnitudes together with the colour indices in several bands. Analysing our results, we investigate the limitation on the astrometry and photometry of faint objects with the 0.7-m telescope.


title: A Causal Bandit Approach to Learning Good Atomic Interventions in Presence of Unobserved Confounders

abstract: We study the problem of determining the best intervention in a Causal Bayesian Network (CBN) specified only by its causal graph. We model this as a stochastic multi-armed bandit (MAB) problem with side-information, where the interventions correspond to the arms of the bandit instance. First, we propose a simple regret minimization algorithm that takes as input a semi-Markovian causal graph with atomic interventions and possibly unobservable variables, and achieves $\tilde{O}(\sqrt{M/T})$ expected simple regret, where $M$ is dependent on the input CBN and could be very small compared to the number of arms. We also show that this is almost optimal for CBNs described by causal graphs having an $n$-ary tree structure. Our simple regret minimization results, both upper and lower bound, subsume previous results in the literature, which assumed additional structural restrictions on the input causal graph. In particular, our results indicate that the simple regret guarantee of our proposed algorithm can only be improved by considering more nuanced structural restrictions on the causal graph. Next, we propose a cumulative regret minimization algorithm that takes as input a general causal graph with all observable nodes and atomic interventions and performs better than the optimal MAB algorithm that does not take causal side-information into account. We also experimentally compare both our algorithms with the best known algorithms in the literature. To the best of our knowledge, this work gives the first simple and cumulative regret minimization algorithms for CBNs with general causal graphs under atomic interventions and having unobserved confounders.


title: Deformations of JT Gravity via Topological Gravity and Applications

abstract: We study the duality between JT gravity and the double-scaled matrix model including their respective deformations. For these deformed theories we relate the thermal partition function to the generating function of topological gravity correlators that are determined as solutions to the KdV hierarchy. We specialise to those deformations of JT gravity coupled to a gas of defects, which conforms with known results in the literature. We express the (asymptotic) thermal partition functions in a low temperature limit, in which non-perturbative corrections are suppressed and the thermal partition function becomes exact. In this limit we demonstrate that there is a Hawking-Page phase transition between connected and disconnected surfaces for this instance of JT gravity with a transition temperature affected by the presence of defects. Furthermore, the calculated spectral form factors show the qualitative behaviour expected for a Hawking-Page phase transition. The considered deformations cause the ramp to be shifted along the real time axis. Finally, we comment on recent results related to conical Weil-Petersson volumes and the analytic continuation to two-dimensional de Sitter space.


title: Quantum illumination with noisy probes: Conditional advantages of non-Gaussianity

abstract: Entangled states, like the two-mode squeezed vacuum state, are known to give quantum advantage in the illumination protocol, a method to detect a weakly reflecting target submerged in a thermal background. We use non-Gaussian photon-added and -subtracted states, affected by local Gaussian noise on top of the omnipresent thermal noise, as probes in the illumination protocol. Based on the difference between the Chernoff bounds obtained with the coherent state and the non-Gaussian state having equal signal strengths, whose positive values denote quantum advantage in illumination, we highlight the hierarchy among non-Gaussian states, which is compatible with correlations per unit signal strength, although the Gaussian states offer the best performance. Interestingly, such hierarchy is different when comparisons are made using the Chernoff bounds. The entire analysis is performed in the presence of different imperfect apparatus like faulty twin-beam generator, imperfect photon addition (subtraction) as well as with noisy non-Gaussian probe states.


title: Countering Misinformation on Social Media Through Educational Interventions: Evidence from a Randomized Experiment in Pakistan

abstract: Fake news is a growing problem in developing countries with potentially far-reaching consequences. We conduct a randomized experiment in urban Pakistan to evaluate the effectiveness of two educational interventions to counter misinformation among low-digital literacy populations. We do not find a significant effect of video-based general educational messages about misinformation. However, when such messages are augmented with personalized feedback based on individuals' past engagement with fake news, we find an improvement of 0.14 standard deviations in identifying fake news. We also find negative but insignificant effects on identifying true news, driven by female respondents. Our results suggest that educational interventions can enable information discernment but their effectiveness critically depends on how well their features and delivery are customized for the population of interest.


title: Counterfactual Explanations in Sequential Decision Making Under Uncertainty

abstract: Methods to find counterfactual explanations have predominantly focused on one step decision making processes. In this work, we initiate the development of methods to find counterfactual explanations for decision making processes in which multiple, dependent actions are taken sequentially over time. We start by formally characterizing a sequence of actions and states using finite horizon Markov decision processes and the Gumbel-Max structural causal model. Building upon this characterization, we formally state the problem of finding counterfactual explanations for sequential decision making processes. In our problem formulation, the counterfactual explanation specifies an alternative sequence of actions differing in at most k actions from the observed sequence that could have led the observed process realization to a better outcome. Then, we introduce a polynomial time algorithm based on dynamic programming to build a counterfactual policy that is guaranteed to always provide the optimal counterfactual explanation on every possible realization of the counterfactual environment dynamics. We validate our algorithm using both synthetic and real data from cognitive behavioral therapy and show that the counterfactual explanations our algorithm finds can provide valuable insights to enhance sequential decision making under uncertainty.


title: Developing and delivering a remote experiment based on the experiential learning framework during COVID-19 pandemic

abstract: The students following Engineering disciplines should not only acquire the conceptual understanding of the concepts but also the processors and attitudes. There are two recognizable learning environments for students, namely, classroom environment and laboratory environment. With the COVID-19 pandemic, both environments merged to online environments, impacting students' development of processes and characteristic attitudes. This paper introduces a theoretical framework based on experiential learning to plan and deliver processes through an online environment. A case study based on the power factor correction experiment was presented. The traditional experiment that runs for 3 hours was broken into smaller tasks such as a pre-lab activity, a simulation exercise, a PowerPoint presentation, a remote laboratory activity, and a final report based on the experiential learning approach. A questionnaire that carries close and open-ended questions were administered to obtain students' reflections about developing the processes through an online-friendly experiential learning approach. The majority of the students like the approach followed and praise for providing them with an opportunity to perform the experiment in a novel way during the COVID-19 situation.


title: Anomaly Detection using Edge Computing in Video Surveillance System: Review

abstract: The current concept of Smart Cities influences urban planners and researchers to provide modern, secured and sustainable infrastructure and give a decent quality of life to its residents. To fulfill this need video surveillance cameras have been deployed to enhance the safety and well-being of the citizens. Despite technical developments in modern science, abnormal event detection in surveillance video systems is challenging and requires exhaustive human efforts. In this paper, we surveyed various methodologies developed to detect anomalies in intelligent video surveillance. Firstly, we revisit the surveys on anomaly detection in the last decade. We then present a systematic categorization of methodologies developed for ease of understanding. Considering the notion of anomaly depends on context, we identify different objects-of-interest and publicly available datasets in anomaly detection. Since anomaly detection is considered a time-critical application of computer vision, our emphasis is on anomaly detection using edge devices and approaches explicitly designed for them. Further, we discuss the challenges and opportunities involved in anomaly detection at the edge.


title: Pile-Up Mitigation using Attention

abstract: Particle production from secondary proton-proton collisions, commonly referred to as pile-up, impair the sensitivity of both new physics searches and precision measurements at LHC experiments. We propose a novel algorithm, PUMA, for identifying pile-up objects with the help of deep neural networks based on sparse transformers. These attention mechanisms were developed for natural language processing but have become popular in other applications. In a realistic detector simulation, our method outperforms classical benchmark algorithms for pile-up mitigation in key observables. It provides a perspective for mitigating the effects of pile-up in the high luminosity era of the LHC, where up to 200 proton-proton collisions are expected to occur simultaneously.


title: Causal Inference with Corrupted Data: Measurement Error, Missing Values, Discretization, and Differential Privacy

abstract: The US Census Bureau will deliberately corrupt data sets derived from the 2020 US Census, enhancing the privacy of respondents while potentially reducing the precision of economic analysis. To investigate whether this trade-off is inevitable, we formulate a semiparametric model of causal inference with high dimensional corrupted data. We propose a procedure for data cleaning, estimation, and inference with data cleaning-adjusted confidence intervals. We prove consistency and Gaussian approximation by finite sample arguments, with a rate of $n^{ 1/2}$ for semiparametric estimands that degrades gracefully for nonparametric estimands. Our key assumption is that the true covariates are approximately low rank, which we interpret as approximate repeated measurements and empirically validate. Our analysis provides nonasymptotic theoretical contributions to matrix completion, statistical learning, and semiparametric statistics. Calibrated simulations verify the coverage of our data cleaning adjusted confidence intervals and demonstrate the relevance of our results for Census-derived data.


title: Inclusive Astronomy in Peru: Contribution of Astronomy Teaching for Visually Impaired People

abstract: Everything we know about the environment around us is thanks to light, a kind of electromagnetic radiation. Astronomy takes advantage of it and all the electromagnetic spectrum with the help of many devices to record them and determine from which places in our Universe they come. These signals must be processed to obtain the images that are will be then exposed to the public. Astronomers know that these images will inspire and generate curiosity in each person who sees them. This Science is inclusive and wants to transmit the knowledge and the beautiful events that happen in the universe to all people. We try to do that by developing and showing other forms of teaching, taking advantage of new technologies that are available today to bring this knowledge closer to the minorities in our country like visually impaired people. Because of this problem, the AstroBVI project arrived in Peru for the first time, thanks to the distribution of tactile images of 3D galaxies that were delivered to us from the Centro de Astronom\'ia de la Universidad de Antofagasta of the and financed by the International Astronomical Union Office of Astronomy for Development (IAU-OAD). This allowed the holding of seven workshops during 2019, visiting various institutions and benefiting more than 160 participants with blindness and low vision, identifying in them a lot of interest, which shows the enormous potential presented by these 3D tactile materials.


title: Sangrahaka: A Tool for Annotating and Querying Knowledge Graphs

abstract: In this work, we present a web-based annotation and querying tool Sangrahaka. It annotates entities and relationships from text corpora and constructs a knowledge graph (KG). The KG is queried using templatized natural language queries. The application is language and corpus agnostic, but can be tuned for special needs of a specific language or a corpus. A customized version of the framework has been used in two annotation tasks. The application is available for download and installation. Besides having a user-friendly interface, it is fast, supports customization, and is fault tolerant on both client and server side. The code is available at https://github.com/hrishikeshrt/sangrahaka and the presentation with a demo is available at https://youtu.be/nw9GFLVZMMo.


title: SAGE: Intrusion Alert-driven Attack Graph Extractor

abstract: Attack graphs (AG) are used to assess pathways availed by cyber adversaries to penetrate a network. State-of-the-art approaches for AG generation focus mostly on deriving dependencies between system vulnerabilities based on network scans and expert knowledge. In real-world operations however, it is costly and ineffective to rely on constant vulnerability scanning and expert-crafted AGs. We propose to automatically learn AGs based on actions observed through intrusion alerts, without prior expert knowledge. Specifically, we develop an unsupervised sequence learning system, SAGE, that leverages the temporal and probabilistic dependence between alerts in a suffix-based probabilistic deterministic finite automaton (S-PDFA) -- a model that accentuates infrequent severe alerts and summarizes paths leading to them. AGs are then derived from the S-PDFA on a per-objective, per-victim basis. Tested with intrusion alerts collected through Collegiate Penetration Testing Competition, SAGE compresses over 330k alerts into 93 AGs. These AGs reflect the strategies used by the participating teams. The AGs are succinct, interpretable, and capture behavioral dynamics, e.g., that attackers will often follow shorter paths to re-exploit objectives.


title: Data-driven reduced order modeling of environmental hydrodynamics using deep autoencoders and neural ODEs

abstract: Model reduction for fluid flow simulation continues to be of great interest across a number of scientific and engineering fields. In a previous work [arXiv:2104.13962], we explored the use of Neural Ordinary Differential Equations (NODE) as a non-intrusive method for propagating the latent-space dynamics in reduced order models. Here, we investigate employing deep autoencoders for discovering the reduced basis representation, the dynamics of which are then approximated by NODE. The ability of deep autoencoders to represent the latent-space is compared to the traditional proper orthogonal decomposition (POD) approach, again in conjunction with NODE for capturing the dynamics. Additionally, we compare their behavior with two classical non-intrusive methods based on POD and radial basis function interpolation as well as dynamic mode decomposition. The test problems we consider include incompressible flow around a cylinder as well as a real-world application of shallow water hydrodynamics in an estuarine system. Our findings indicate that deep autoencoders can leverage nonlinear manifold learning to achieve a highly efficient compression of spatial information and define a latent-space that appears to be more suitable for capturing the temporal dynamics through the NODE framework.


title: Quasiconformal Flows on non-Conformally Flat Spheres

abstract: We study integral curvature conditions for a Riemannian metric $g$ on $S^4$ that quantify the best bilipschitz constant between $(S^4,g)$ and the standard metric on $S^4$. Our results show that the best bilipschitz constant is controlled by the $L^2$-norm of the Weyl tensor and the $L^1$-norm of the $Q$-curvature, under the conditions that those quantities are sufficiently small, $g$ has a positive Yamabe constant and the $Q$-curvature is mean-positive. The proof of the result is achieved in two steps. Firstly, we construct a quasiconformal map between two conformally related metrics in a positive Yamabe class. Secondly, we apply the Ricci flow to establish the bilipschitz equivalence from such a conformal class to the standard conformal class on $S^4$.


title: On Quanta of Information and Electromagnetic Fields

abstract: A preliminary model of meaningful quanta of information arising from the quantum vacuum, and giving rise to the elementary quantum fields, is presented. We attempt to establish an equivalence of information and physical action, with the quanta (or packets of words) of information as a source of quantum fields (with a possible role in governing the interactions among those), similar to the quantum fields being the source of classical information. The information and quantum fields emerge in an entangled form out of the quantum vacuum, with their common underlying basis bearing intrinsically both a well-defined deterministic (geometric) and probabilistic (stochastic) structure. The description of equivalent energy of a quantum of the electromagnetic field, in terms of its information content, as well as an expression of the signal power, is provided. Various aspects of this problem and its possible implications and repercussions on our understanding of elementary quantum fields and their interactions are discussed. It is suggested that information, in a deterministic form, as well as chaotic processes within the quantum vacuum, give rise to a possible condensation of information into quantum fields.


title: Space Efficient Two-Dimensional Orthogonal Colored Range Counting

abstract: In the two-dimensional orthogonal colored range counting problem, we preprocess a set, $P$, of $n$ colored points on the plane, such that given an orthogonal query rectangle, the number of distinct colors of the points contained in this rectangle can be computed efficiently. For this problem, we design three new solutions, and the bounds of each can be expressed in some form of time-space tradeoff. By setting appropriate parameter values for these solutions, we can achieve new specific results with (the space are in words and $\epsilon$ is an arbitrary constant in $(0,1)$): ** $O(n\lg^3 n)$ space and $O(\sqrt{n}\lg^{5/2} n \lg \lg n)$ query time; ** $O(n\lg^2 n)$ space and $O(\sqrt{n}\lg^{4+\epsilon} n)$ query time; ** $O(n\frac{\lg^2 n}{\lg \lg n})$ space and $O(\sqrt{n}\lg^{5+\epsilon} n)$ query time; ** $O(n\lg n)$ space and $O(n^{1/2+\epsilon})$ query time. A known conditional lower bound to this problem based on Boolean matrix multiplication gives some evidence on the difficulty of achieving near-linear space solutions with query time better than $\sqrt{n}$ by more than a polylogarithmic factor using purely combinatorial approaches. Thus the time and space bounds in all these results are efficient. Previously, among solutions with similar query times, the most space-efficient solution uses $O(n\lg^4 n)$ space to answer queries in $O(\sqrt{n}\lg^8 n)$ time (SIAM. J. Comp.~2008). Thus the new results listed above all achieve improvements in space efficiency, while all but the last result achieve speed-up in query time as well.


title: Gravitational Wave Timing Array

abstract: We describe the design of a gravitational wave timing array, a novel scheme that can be used to search for low-frequency gravitational waves by monitoring continuous gravitational waves at higher frequencies. We show that observations of gravitational waves produced by Galactic binaries using a space-based detector like LISA provide sensitivity in the nanohertz to microhertz band. While the expected sensitivity of this proposal is not competitive with other methods, it fills a gap in frequency space around the microhertz regime, which is above the range probed by current pulsar timing arrays and below the expected direct frequency coverage of LISA. The low-frequency extension of sensitivity does not require any experimental design change to space-based gravitational wave detectors, and can be achieved with the data products that would already be collected by them.


title: Digitized-counterdiabatic quantum approximate optimization algorithm

abstract: The quantum approximate optimization algorithm (QAOA) has proved to be an effective classical-quantum algorithm serving multiple purposes, from solving combinatorial optimization problems to finding the ground state of many-body quantum systems. Since QAOA is an ansatz-dependent algorithm, there is always a need to design ansatz for better optimization. To this end, we propose a digitized version of QAOA enhanced via the use of shortcuts to adiabaticity. Specifically, we use a counterdiabatic (CD) driving term to design a better ansatz, along with the Hamiltonian and mixing terms, enhancing the global performance. We apply our digitized-counterdiabatic QAOA to Ising models, classical optimization problems, and the P-spin model, demonstrating that it outperforms standard QAOA in all cases we study.


title: iPOKE: Poking a Still Image for Controlled Stochastic Video Synthesis

abstract: How would a static scene react to a local poke? What are the effects on other parts of an object if you could locally push it? There will be distinctive movement, despite evident variations caused by the stochastic nature of our world. These outcomes are governed by the characteristic kinematics of objects that dictate their overall motion caused by a local interaction. Conversely, the movement of an object provides crucial information about its underlying distinctive kinematics and the interdependencies between its parts. This two-way relation motivates learning a bijective mapping between object kinematics and plausible future image sequences. Therefore, we propose iPOKE -- invertible Prediction of Object Kinematics -- that, conditioned on an initial frame and a local poke, allows to sample object kinematics and establishes a one-to-one correspondence to the corresponding plausible videos, thereby providing a controlled stochastic video synthesis. In contrast to previous works, we do not generate arbitrary realistic videos, but provide efficient control of movements, while still capturing the stochastic nature of our environment and the diversity of plausible outcomes it entails. Moreover, our approach can transfer kinematics onto novel object instances and is not confined to particular object classes. Our project page is available at https://bit.ly/3dJN4Lf.


title: Depth-supervised NeRF: Fewer Views and Faster Training for Free

abstract: A commonly observed failure mode of Neural Radiance Field (NeRF) is fitting incorrect geometries when given an insufficient number of input views. One potential reason is that standard volumetric rendering does not enforce the constraint that most of a scene's geometry consist of empty space and opaque surfaces. We formalize the above assumption through DS-NeRF (Depth-supervised Neural Radiance Fields), a loss for learning radiance fields that takes advantage of readily-available depth supervision. We leverage the fact that current NeRF pipelines require images with known camera poses that are typically estimated by running structure-from-motion (SFM). Crucially, SFM also produces sparse 3D points that can be used as "free" depth supervision during training: we add a loss to encourage the distribution of a ray's terminating depth matches a given 3D keypoint, incorporating depth uncertainty. DS-NeRF can render better images given fewer training views while training 2-3x faster. Further, we show that our loss is compatible with other recently proposed NeRF methods, demonstrating that depth is a cheap and easily digestible supervisory signal. And finally, we find that DS-NeRF can support other types of depth supervision such as scanned depth sensors and RGB-D reconstruction outputs.


title: Learned Visual Navigation for Under-Canopy Agricultural Robots

abstract: We describe a system for visually guided autonomous navigation of under-canopy farm robots. Low-cost under-canopy robots can drive between crop rows under the plant canopy and accomplish tasks that are infeasible for over-the-canopy drones or larger agricultural equipment. However, autonomously navigating them under the canopy presents a number of challenges: unreliable GPS and LiDAR, high cost of sensing, challenging farm terrain, clutter due to leaves and weeds, and large variability in appearance over the season and across crop types. We address these challenges by building a modular system that leverages machine learning for robust and generalizable perception from monocular RGB images from low-cost cameras, and model predictive control for accurate control in challenging terrain. Our system, CropFollow, is able to autonomously drive 485 meters per intervention on average, outperforming a state-of-the-art LiDAR based system (286 meters per intervention) in extensive field testing spanning over 25 km.


title: The COMBS Survey -- III. The Chemodynamical Origins of Metal-Poor Bulge Stars

abstract: The characteristics of the stellar populations in the Galactic Bulge inform and constrain the Milky Way's formation and evolution. The metal-poor population is particularly important in light of cosmological simulations, which predict that some of the oldest stars in the Galaxy now reside in its center. The metal-poor bulge appears to consist of multiple stellar populations that require dynamical analyses to disentangle. In this work, we undertake a detailed chemodynamical study of the metal-poor stars in the inner Galaxy. Using R$\sim$ 20,000 VLT/GIRAFFE spectra of 319 metal-poor (-2.55 dex$\leq$[Fe/H]$\leq$0.83 dex, with $\overline{\rm{[Fe/H]}}$=-0.84 dex) stars, we perform stellar parameter analysis and report 12 elemental abundances (C, Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Mn, Zn, Ba, and Ce) with precisions of $\approx$0.10 dex. Based on kinematic and spatial properties, we categorise the stars into four groups, associated with the following Galactic structures: the inner bulge, the outer bulge, the halo, and the disk. We find evidence that the inner and outer bulge population is more chemically complex (i.e., higher chemical dimensionality and less correlated abundances) than the halo population. This result suggests that the older bulge population was enriched by a larger diversity of nucleosynthetic events. We also find one inner bulge star with a [Ca/Mg] ratio consistent with theoretical pair-instability supernova yields and two stars that have chemistry consistent with globular cluster stars.