title: The X-shooter Spectral Library (XSL): Data Release 3

abstract: We present the third data release (DR3) of the X-shooter Spectral Library (XSL). This moderate-to-high resolution, near-ultraviolet-to-near-infrared ($350-2480$ nm, R $\sim$ 10 000) spectral library is composed of 830 stellar spectra of 683 stars. DR3 improves upon the previous data release by providing the combined de-reddened spectra of the three X-shooter segments over the full $350-2480$ nm wavelength range. It also includes additional 20 M-dwarf spectra from the ESO archive. We provide detailed comparisons between this library and Gaia EDR3, MILES, NGSL, CaT library, and (E-)IRTF. The normalised rms deviation is better than $D=0.05$ or 5$\%$ for the majority of spectra in common between MILES (144 spectra of 180), NGSL (112$/$116), and (E-)IRTF (55$/$77) libraries. Comparing synthetic colours of those spectra reveals only negligible offsets and small rms scatter, such as the median offset(rms) 0.001$\pm$0.040 mag in the (box1-box2) colour of the UVB arm,-0.004$\pm$0.028 mag in (box3-box4) of the VIS arm, and -0.001$\pm$0.045 mag in (box2-box3) colour between the UVB and VIS arms, when comparing stars in common with MILES. We also find an excellent agreement between the Gaia published (BP-RP) colours and those measured from the XSL DR3 spectra, with a zero median offset and an rms scatter of 0.037 mag for 449 non-variable stars. The unmatched characteristics of this library, which combine a relatively high resolution, a large number of stars, and an extended wavelength coverage, will help us to bridge the gap between the optical and the near-IR studies of intermediate and old stellar populations, and to probe low-mass stellar systems.


title: StructFormer: Learning Spatial Structure for Language-Guided Semantic Rearrangement of Novel Objects

abstract: Geometric organization of objects into semantically meaningful arrangements pervades the built world. As such, assistive robots operating in warehouses, offices, and homes would greatly benefit from the ability to recognize and rearrange objects into these semantically meaningful structures. To be useful, these robots must contend with previously unseen objects and receive instructions without significant programming. While previous works have examined recognizing pairwise semantic relations and sequential manipulation to change these simple relations none have shown the ability to arrange objects into complex structures such as circles or table settings. To address this problem we propose a novel transformer-based neural network, StructFormer, which takes as input a partial-view point cloud of the current object arrangement and a structured language command encoding the desired object configuration. We show through rigorous experiments that StructFormer enables a physical robot to rearrange novel objects into semantically meaningful structures with multi-object relational constraints inferred from the language command.


title: Modelling simple stellar populations in the near-ultraviolet to near-infrared with the X-shooter Spectral Library (XSL)

abstract: We present simple stellar population models based on the empirical X-shooter Spectral Library (XSL) from NUV to NIR wavelengths. The unmatched characteristics of relatively high resolution and extended wavelength coverage ($350-2480$ nm, $R\sim10\,000$) of the XSL population models bring us closer to bridging optical and NIR studies of intermediate and old stellar populations. It is now common to find good agreement between observed and predicted NUV and optical properties of stellar clusters due to our good understanding of the main-sequence and early giant phases of stars. However, NIR spectra of intermediate-age and old stellar populations are sensitive to cool K and M giants. The asymptotic giant branch, especially the thermally pulsing asymptotic giant branch, shapes the NIR spectra of $0.5-2$ Gyr old stellar populations; the tip of the red giant branch defines the NIR spectra of populations with ages larger than that. We construct sequences of the average spectra of static giants, variable-rich giants, and C-rich giants to include in the models separately. The models span the metallicity range $-2.2<[Fe/H]<+0.2$ and ages above 50 Myr, a broader range in the NIR than in other models based on empirical spectral libraries. Our models can reproduce the integrated optical colours of the Coma cluster galaxies at the same level as other semi-empirical models found in the literature. In the NIR, there are notable differences between the colours of the models and Coma cluster galaxies. The XSL models expand the range of predicted values of NIR indices compared to other models based on empirical libraries. Our models make it possible to perform in-depth studies of colours and spectral features consistently throughout the optical and the NIR range to clarify the role of evolved cool stars in stellar populations.


title: Iharm3D: Vectorized General Relativistic Magnetohydrodynamics

abstract: Iharm3D is an open-source C code for simulating black hole accretion systems in arbitrary stationary spacetimes using ideal general-relativistic magnetohydrodynamics (GRMHD). It is an implementation of the HARM ("High Accuracy Relativistic Magnetohydrodynamics") algorithm outlined in Gammie et al. (2003) with updates as outlined in McKinney & Gammie (2004) and Noble et al. (2006). The code is most directly derived from Ryan et al. (2015) but with radiative transfer portions removed. HARM is a conservative finite-volume scheme for solving the equations of ideal GRMHD, a hyperbolic system of partial differential equations, on a logically Cartesian mesh in arbitrary coordinates.


title: Difference-in-Differences with Geocoded Microdata

abstract: This paper formalizes a common approach for estimating effects of treatment at a specific location using geocoded microdata. This estimator compares units immediately next to treatment (an inner-ring) to units just slightly further away (an outer-ring). I introduce intuitive assumptions needed to identify the average treatment effect among the affected units and illustrates pitfalls that occur when these assumptions fail. Since one of these assumptions requires knowledge of exactly how far treatment effects are experienced, I propose a new method that relaxes this assumption and allows for nonparametric estimation using partitioning-based least squares developed in Cattaneo et. al. (2019). Since treatment effects typically decay/change over distance, this estimator improves analysis by estimating a treatment effect curve as a function of distance from treatment. This is contrast to the traditional method which, at best, identifies the average effect of treatment. To illustrate the advantages of this method, I show that Linden and Rockoff (2008) under estimate the effects of increased crime risk on home values closest to the treatment and overestimate how far the effects extend by selecting a treatment ring that is too wide.


title: On the local limit theorems for lower psi-mixing Markov chains

abstract: In this paper we investigate the local limit theorem for additive functionals of nonstationary Markov chains that converge in distribution. We consider both the lattice and the non-lattice cases. The results are also new in the stationary setting and lead to local limit theorems linked to convergence to stable distributions. The conditions are imposed to individual summands and are expressed in terms of lower psi-mixing coefficients.


title: CoFi: Coarse-to-Fine ICP for LiDAR Localization in an Efficient Long-lasting Point Cloud Map

abstract: LiDAR odometry and localization has attracted increasing research interest in recent years. In the existing works, iterative closest point (ICP) is widely used since it is precise and efficient. Due to its non-convexity and its local iterative strategy, however, ICP-based method easily falls into local optima, which in turn calls for a precise initialization. In this paper, we propose CoFi, a Coarse-to-Fine ICP algorithm for LiDAR localization. Specifically, the proposed algorithm down-samples the input point sets under multiple voxel resolution, and gradually refines the transformation from the coarse point sets to the fine-grained point sets. In addition, we propose a map based LiDAR localization algorithm that extracts semantic feature points from the LiDAR frames and apply CoFi to estimate the pose on an efficient point cloud map. With the help of the Cylinder3D algorithm for LiDAR scan semantic segmentation, the proposed CoFi localization algorithm demonstrates the state-of-the-art performance on the KITTI odometry benchmark, with significant improvement over the literature.


title: Operator-induced structural variable selection for identifying materials genes

abstract: In the emerging field of materials informatics, a fundamental task is to identify physicochemically meaningful descriptors, or materials genes, which are engineered from primary features and a set of elementary algebraic operators through compositions. Standard practice directly analyzes the high-dimensional candidate predictor space in a linear model; statistical analyses are then substantially hampered by the daunting challenge posed by the astronomically large number of correlated predictors with limited sample size. We formulate this problem as variable selection with operator-induced structure (OIS) and propose a new method to achieve unconventional dimension reduction by utilizing the geometry embedded in OIS. Although the model remains linear, we iterate nonparametric variable selection for effective dimension reduction. This enables variable selection based on ab initio primary features, leading to a method that is orders of magnitude faster than existing methods, with improved accuracy. To select the nonparametric module, we discuss a desired performance criterion that is uniquely induced by variable selection with OIS; in particular, we propose to employ a Bayesian Additive Regression Trees (BART)-based variable selection method. Numerical studies show superiority of the proposed method, which continues to exhibit robust performance when the input dimension is out of reach of existing methods. Our analysis of single-atom catalysis identifies physical descriptors that explain the binding energy of metal-support pairs with high explanatory power, leading to interpretable insights to guide the prevention of a notorious problem called sintering and aid catalysis design.


title: Diversity metric for evaluation of quantum annealing

abstract: Solving discrete NP-hard problems is an important part of scientific discoveries and operations research as well as many commercial applications. A commonly used metric to compare meta-heuristic solvers is the time required to obtain an optimal solution, known as time to solution. However, for some applications it is desirable to have a set of high-quality and diverse solutions, instead of a single optimal one. For these applications, time to solution may not be informative of the performance of a solver, and another metric would be necessary. In particular, it is not known how well quantum solvers sample the configuration space in comparison to their classical counterparts. Here, we apply a recently introduced collective distance measure in solution space to quantify diversity by Mohseni et. al. and, based on that, we employ time-to-diversity as a metric for evaluation of meta-heuristics solvers. We use this measure to compare the performance of the D-Wave quantum annealing processor with a few classical heuristic solvers on a set of synthetic problems and show that D-Wave quantum annealing processor is indeed a competitive heuristic, and on many instances outperforms state-of-the-art classical solvers, while it remains on par on other instances. This suggests that a portfolio solver that combines quantum and classical solutions may win over all solvers.


title: Ghost-imaging-enhanced non-invasive spectral characterization of stochastic x-ray free-electron-laser pulses

abstract: High-intensity ultrashort X-ray free-electron laser (XFEL) pulses are revolutionizing the study of fundamental nonlinear x-ray matter interactions and coupled electronic and nuclear dynamics. To fully exploit the potential of this powerful tool for advanced x-ray spectroscopies, a noninvasive spectral characterization of incident stochastic XFEL pulses with high resolution is a key requirement. Here we present a methodology that combines high-acceptance angle-resolved photoelectron time-of-flight spectroscopy and ghost imaging to enhance the quality of spectral characterization of x-ray free-electron laser pulses. Implementation of this non-invasive high-resolution x-ray diagnostic can greatly benefit the ultrafast x-ray spectroscopy community by functioning as a transparent beamsplitter for applications such as transient absorption spectroscopy in averaging mode as well as covariance-based x-ray nonlinear spectroscopies in single-shot mode where the shot-to-shot fluctuations inherent to a self-amplified spontaneous emission (SASE) XFEL pulse are a powerful asset.


title: The Rayleigh Benard problem for compressible fluid flows

abstract: We consider the physically relevant fully compressible setting of the Rayleigh Benard problem of a fluid confined between two parallel plates, heated from the bottom, and subjected to the gravitational force. Under suitable restrictions imposed on the constitutive relations we show that this open system is dissipative in the sense of Levinson, meaning there exists a bounded absorbing set for any global-in-time weak solution. In addition, global-in-time trajectories are asymptotically compact in suitable topologies and the system possesses a global compact trajectory attractor A. The standard technique of Krylov and Bogolyubov then yields the existence of an invariant measure - a stationary statistical solution sitting on A. In addition, the Birkhoff - Khinchin ergodic theorem provides convergence of ergodic averages of solutions belonging to A a.s. with respect to the invariant measure.


title: Theoretical Advances in Current Estimation and Navigation from a Glider-Based Acoustic Doppler Current Profiler (ADCP)

abstract: We examine acoustic Doppler current profiler (ADCP) measurements from underwater gliders to determine glider position, glider velocity, and subsurface current. ADCPs, however, do not directly observe the quantities of interest; instead, they measure the relative motion of the vehicle and the water column. We examine the lineage of mathematical innovations that have previously been applied to this problem, discovering an unstated but incorrect assumption of independence. We reframe a recent method to form a joint probability model of current and vehicle navigation, which allows us to correct this assumption and extend the classic Kalman smoothing method. Detailed simulations affirm the efficacy of our approach for computing estimates and their uncertainty. The joint model developed here sets the stage for future work to incorporate constraints, range measurements, and robust statistical modeling.


title: fairadapt: Causal Reasoning for Fair Data Pre-processing

abstract: Machine learning algorithms are useful for various predictions tasks, but they can also learn how to discriminate, based on gender, race or other sensitive attributes. This realization gave rise to the field of fair machine learning, which aims to measure and mitigate such algorithmic bias. This manuscript describes the R-package fairadapt, which implements a causal inference pre-processing method. By making use of a causal graphical model and the observed data, the method can be used to address hypothetical questions of the form "What would my salary have been, had I been of a different gender/race?". Such individual level counterfactual reasoning can help eliminate discrimination and help justify fair decisions. We also discuss appropriate relaxations which assume certain causal pathways from the sensitive attribute to the outcome are not discriminatory.


title: Ab initio Calculations in Atoms, Molecules, and Solids, Treating Spin-Orbit Coupling and Electron Interaction on Equal Footing

abstract: We incorporate explicit, non-perturbative treatment of spin-orbit coupling into ab initio auxiliary-field quantum Monte Carlo (AFQMC) calculations. The approach allows a general computational framework for molecular and bulk systems in which materials specificity, electron correlation, and spin-orbit coupling effects can be captured accurately and on equal footing, with favorable computational scaling versus system size. We adopt relativistic effective-core potentials which have been obtained by fitting to fully relativistic data and which have demonstrated a high degree of reliability and transferability in molecular systems. This results in a 2-component spin-coupled Hamiltonian, which is then treated by generalizing the ab initio AFQMC approach. We demonstrate the method by computing the electron affinity in Pb, the bond dissociation energy in Br$_2$ and I$_2$, and solid Bi.


title: $m_b$ at $m_H$: the running bottom quark mass and the Higgs boson

abstract: We present a new measurement of the bottom quark mass in the $\bar{MS}$ scheme at the renormalization scale of the Higgs boson mass from measurements of Higgs boson decay rates at the LHC: $m_b(m_H)= 2.60^{+0.36}_{-0.31}$ GeV. The measurement has a negligible theory uncertainty and excellent prospects to improve at the HL-LHC and a future Higgs factory. Confronting this result, and $m_b(m_b)$ from low-energy measurements and $m_b(m_Z)$ from $Z-$pole data, with the prediction of the scale evolution of the renormalization group equations, we find strong evidence for the "running" of the bottom quark mass.


title: Countable-configurations and paradoxical decompositions

abstract: In this paper we define countable-configuration of groups and prove that two Hopfian groups with the same set of countable-configurations are isomorphic and vice versa. We also study the countable paradoxical decomposition of groups. It is proved that a group G admits a countable paradoxical decomposition if and only if it is infinite.


title: Rational Ehrhart Theory

abstract: The Ehrhart quasipolynomial of a rational polytope $\mathsf{P}$ encodes fundamental arithmetic data of $\mathsf{P}$, namely, the number of integer lattice points in positive integral dilates of $\mathsf{P}$. Ehrhart quasipolynomials were introduced in the 1960s, satisfy several fundamental structural results and have applications in many areas of mathematics and beyond. The enumerative theory of lattice points in rational (equivalently, real) dilates of rational polytopes is much younger, starting with work by Linke (2011), Baldoni-Berline-Koeppe-Vergne (2013), and Stapledon (2017). We introduce a generating-function ansatz for rational Ehrhart quasipolynomials, which unifies several known results in classical and rational Ehrhart theory. In particular, we define $\gamma$-rational Gorenstein polytopes, which extend the classical notion to the rational setting and encompass the generalized reflexive polytopes studied by Fiset-Kasprzyk (2008) and Kasprzyk-Nill (2012).


title: MultiHead MultiModal Deep Interest Recommendation Network

abstract: With the development of information technology, human beings are constantly producing a large amount of information at all times. How to obtain the information that users are interested in from the large amount of information has become an issue of great concern to users and even business managers. In order to solve this problem, from traditional machine learning to deep learning recommendation systems, researchers continue to improve optimization models and explore solutions. Because researchers have optimized more on the recommendation model network structure, they have less research on enriching recommendation model features, and there is still room for in-depth recommendation model optimization. Based on the DIN\cite{Authors01} model, this paper adds multi-head and multi-modal modules, which enriches the feature sets that the model can use, and at the same time strengthens the cross-combination and fitting capabilities of the model. Experiments show that the multi-head multi-modal DIN improves the recommendation prediction effect, and outperforms current state-of-the-art methods on various comprehensive indicators.


title: Evaluating and Improving Interactions with Hazy Oracles

abstract: Many AI systems integrate sensor inputs, world knowledge, and human-provided information to perform inference. While such systems often treat the human input as flawless, humans are better thought of as hazy oracles whose input may be ambiguous or outside of the AI system's understanding. In such situations it makes sense for the AI system to defer its inference while it disambiguates the human-provided information by, for example, asking the human to rephrase the query. Though this approach has been considered in the past, current work is typically limited to application-specific methods and non-standardized human experiments. We instead introduce and formalize a general notion of deferred inference. Using this formulation, we then propose a novel evaluation centered around the Deferred Error Volume (DEV) metric, which explicitly considers the tradeoff between error reduction and the additional human effort required to achieve it. We demonstrate this new formalization and an innovative deferred inference method on the disparate tasks of Single-Target Video Object Tracking and Referring Expression Comprehension, ultimately reducing error by up to 48% without any change to the underlying model or its parameters.


title: Radioactive isotopes in the interstellar medium

abstract: Radioactive components of the interstellar medium provide an entirely-different and new aspect to the studies of the interstellar medium. Injected from sources of nucleosynthesis, unstable nuclei decay along their trajectories. Measurements can occur through characteristic gamma rays that are emitted with the decay, or in cosmic material samples through abundances of parent and daughter isotopes as they change with decay. The dynamics and material flows within interstellar medium are thus accessible to measurement, making use of the intrinsic clock that radioactive decay provides. We describe how measurements of radioactive decay have obtained a break-through in studies of the interstellar medium, after first summarizing the characteristics of radioactivity and the sources of unstable nuclei.


title: Waveguiding valley excitons in monolayer transition metal dicalcogenides by dielectric interfaces in the substrate

abstract: In monolayers of the semiconducting transition metal dichalcogenides, the electron-hole exchange interaction splits the exciton dispersion into a massive transverse branch, and a longitudinal branch that has very light or even zero mass depending on the form of screened Coulomb interaction. The group velocity of the longitudinal branch is sensitive to the strength of electron-hole exchange, which can be engineered through the dielectric environment. Here we show that dielectric patterning of the substrate can be exploited to realize waveguide of the exciton in the longitudinal branch in a homogeneous monolayer, leaving the massive transverse branch unaffected. At a lateral interface of different dielectric constant in the substrate, the transmission and reflection of exciton in the longitudinal branch obey the Snell-Descartes law of optical system, and total reflection can be exploited to realize excitonic waveguide using two parallel interfaces. The same dielectric pattern of the substrate appears to be completely transparent for the massive transverse branch exciton, which has no interface scattering. When the monolayer is placed on a one-dimensional dielectric superlattice, the dispersion of the longitudinal branch is strongly renormalized, and the wavefunctions exhibit one-dimensional features, confined to either the low-dielectric or high-dielectric regions. In contrast, the massive transverse branch excitons are not affected by the substrate dielectric pattern, exhibiting pristine properties as in a freestanding monolayer.


title: Batalin-Vilkovisky formality for Chern-Simons theory

abstract: We prove that the differential graded Lie algebra of functionals associated to the Chern-Simons theory of a semisimple Lie algebra is homotopy abelian. For a general field theory, we show that the variational complex in the Batalin-Vilkovisky formalism is a differential graded Lie algebra.


title: Long Random Matrices and Tensor Unfolding

abstract: In this paper, we consider the singular values and singular vectors of low rank perturbations of large rectangular random matrices, in the regime the matrix is "long": we allow the number of rows (columns) to grow polynomially in the number of columns (rows). We prove there exists a critical signal-to-noise ratio (depending on the dimensions of the matrix), and the extreme singular values and singular vectors exhibit a BBP type phase transition. As a main application, we investigate the tensor unfolding algorithm for the asymmetric rank-one spiked tensor model, and obtain an exact threshold, which is independent of the procedure of tensor unfolding. If the signal-to-noise ratio is above the threshold, tensor unfolding detects the signals; otherwise, it fails to capture the signals.


title: Learning Partial Equivariances from Data

abstract: Group Convolutional Neural Networks (G-CNNs) constrain learned features to respect the symmetries in the selected group, and lead to better generalization when these symmetries appear in the data. If this is not the case, however, equivariance leads to overly constrained models and worse performance. Frequently, transformations occurring in data can be better represented by a subset of a group than by a group as a whole, e.g., rotations in $[-90^{\circ}, 90^{\circ}]$. In such cases, a model that respects equivariance $\textit{partially}$ is better suited to represent the data. In addition, relevant transformations may differ for low and high-level features. For instance, full rotation equivariance is useful to describe edge orientations in a face, but partial rotation equivariance is better suited to describe face poses relative to the camera. In other words, the optimal level of equivariance may differ per layer. In this work, we introduce $\textit{Partial G-CNNs}$: G-CNNs able to learn layer-wise levels of partial and full equivariance to discrete, continuous groups and combinations thereof as part of training. Partial G-CNNs retain full equivariance when beneficial, e.g., for rotated MNIST, but adjust it whenever it becomes harmful, e.g., for classification of 6 / 9 digits or natural images. We empirically show that partial G-CNNs pair G-CNNs when full equivariance is advantageous, and outperform them otherwise.


title: The six operations in topology

abstract: In this paper we show that the six functor formalism for sheaves on locally compact Hausdorff topological spaces, as developed for example in Kashiwara and Schapira's book Sheaves on Manifolds, can be extended to sheaves with values in any closed symmetric monoidal $\infty$-category which is stable and bicomplete. Notice that, since we do not assume that our coefficients are presentable or restrict to hypercomplete sheaves, our arguments are not obvious and are substantially different from the ones explained by Kashiwara and Schapira. Along the way we also study locally contractible geometric morphisms and prove that, if $f:X\rightarrow Y$ is a continuous map which induces a locally contractible geometric morphism, then the exceptional pullback functor $f^!$ preserves colimits and can be related to the pullback $f^*$. At the end of our paper we also show how one can express Atiyah duality by means of the six functor formalism.


title: Neural Medication Extraction: A Comparison of Recent Models in Supervised and Semi-supervised Learning Settings

abstract: Drug prescriptions are essential information that must be encoded in electronic medical records. However, much of this information is hidden within free-text reports. This is why the medication extraction task has emerged. To date, most of the research effort has focused on small amount of data and has only recently considered deep learning methods. In this paper, we present an independent and comprehensive evaluation of state-of-the-art neural architectures on the I2B2 medical prescription extraction task both in the supervised and semi-supervised settings. The study shows the very competitive performance of simple DNN models on the task as well as the high interest of pre-trained models. Adapting the latter models on the I2B2 dataset enables to push medication extraction performances above the state-of-the-art. Finally, the study also confirms that semi-supervised techniques are promising to leverage large amounts of unlabeled data in particular in low resource setting when labeled data is too costly to acquire.


title: Persistent homology of quantum entanglement

abstract: Structure in quantum entanglement entropy is often leveraged to focus on a small corner of the exponentially large Hilbert space and efficiently parameterize the problem of finding ground states. A typical example is the use of matrix product states for local and gapped Hamiltonians. We study the structure of entanglement entropy using persistent homology, a relatively new method from the field of topological data analysis. The inverse quantum mutual information between pairs of sites is used as a distance metric to form a filtered simplicial complex. Both ground states and excited states of common spin models are analyzed as an example. Furthermore, the effect of homology with different coefficients and boundary conditions is also explored. Beyond these basic examples, we also discuss the promising future applications of this modern computational approach, including its connection to the question of how spacetime could emerge from entanglement.


title: TOI-2285b: A 1.7 Earth-radius Planet Near the Habitable Zone around a Nearby M Dwarf

abstract: We report the discovery of TOI-2285b, a sub-Neptune-sized planet transiting a nearby (42 pc) M dwarf with a period of 27.3 days. We identified the transit signal from the Transiting Exoplanet Survey Satellite photometric data, which we confirmed with ground-based photometric observations using the multiband imagers MuSCAT2 and MuSCAT3. Combining these data with other follow-up observations including high resolution spectroscopy with the Tillinghast Reflector Echelle Spectrograph, high resolution imaging with the SPeckle Polarimeter, and radial velocity (RV) measurements with the InfraRed Doppler instrument, we find that the planet has a radius of 1.74 $\pm$ 0.08 $R_\oplus$, a mass of $<$ 19.5 $M_\oplus$ (95\% c.l.), and an insolation flux of 1.54 $\pm$ 0.14 times that of the Earth. Although the planet resides just outside the habitable zone for a rocky planet, if the planet harbors an H$_2$O layer under a hydrogen-rich atmosphere, then liquid water could exist on the surface of the H$_2$O layer depending on the planetary mass and water mass fraction. The bright host star in the near infrared ($K_s=9.0$) makes this planet an excellent target for further RV and atmospheric observations to improve our understanding on the composition, formation, and habitability of sub-Neptune-sized planets.


title: A Bayesian Analysis of Two-Stage Randomized Experiments in the Presence of Interference, Treatment Nonadherence, and Missing Outcomes

abstract: Three critical issues for causal inference that often occur in modern, complicated experiments are interference, treatment nonadherence, and missing outcomes. A great deal of research efforts has been dedicated to developing causal inferential methodologies that address these issues separately. However, methodologies that can address these issues simultaneously are lacking. We propose a Bayesian causal inference methodology to address this gap. Our methodology extends existing causal frameworks and methods, specifically, two-staged randomized experiments and the principal stratification framework. In contrast to existing methods that invoke strong structural assumptions to identify principal causal effects, our Bayesian approach uses flexible distributional models that can accommodate the complexities of interference and missing outcomes, and that ensure that principal causal effects are weakly identifiable. We illustrate our methodology via simulation studies and a re-analysis of real-life data from an evaluation of India's National Health Insurance Program. Our methodology enables us to identify new active causal effects that were not identified in past analyses. Ultimately, our simulation studies and case study demonstrate how our methodology can yield more informative analyses in modern experiments with interference, treatment nonadherence, missing outcomes, and complicated outcome generation mechanisms.


title: An Adaptive Sampling and Edge Detection Approach for Encoding Static Images for Spiking Neural Networks

abstract: Current state-of-the-art methods of image classification using convolutional neural networks are often constrained by both latency and power consumption. This places a limit on the devices, particularly low-power edge devices, that can employ these methods. Spiking neural networks (SNNs) are considered to be the third generation of artificial neural networks which aim to address these latency and power constraints by taking inspiration from biological neuronal communication processes. Before data such as images can be input into an SNN, however, they must be first encoded into spike trains. Herein, we propose a method for encoding static images into temporal spike trains using edge detection and an adaptive signal sampling method for use in SNNs. The edge detection process consists of first performing Canny edge detection on the 2D static images and then converting the edge detected images into two X and Y signals using an image-to-signal conversion method. The adaptive signaling approach consists of sampling the signals such that the signals maintain enough detail and are sensitive to abrupt changes in the signal. Temporal encoding mechanisms such as threshold-based representation (TBR) and step-forward (SF) are then able to be used to convert the sampled signals into spike trains. We use various error and indicator metrics to optimize and evaluate the efficiency and precision of the proposed image encoding approach. Comparison results between the original and reconstructed signals from spike trains generated using edge-detection and adaptive temporal encoding mechanism exhibit 18x and 7x reduction in average root mean square error (RMSE) compared to the conventional SF and TBR encoding, respectively, while used for encoding MNIST dataset.


title: Gravitational wave background discovered by NANOGrav as evidence of a cyclic universe

abstract: The consortium NANOGrav discovered the isotropic gravitational wave background (GWB) with an amplitude of $h\sim 10^{-15}$ and a frequency of $f\sim 10^{-8}$ Hz using observations of millisecond pulsars. We hypothesize that the GWB is relic radiation left over from the merging stellar mass black holes (SBHs) during Big Crunch at $z\sim 10^{10}$. The relic gravitational waves are similar to the gravitational waves with $f\sim 10^2$ Hz discovered by LIGO in 2015, taking into account a decrease in frequency by a factor $\sim 10^{10}$ due to the expansion of the Universe. We take as a basis the observed spectrum of 139 SBHs, discovered by the LIGO observatory. Our model explains well all the observed features of the GWB. Unlike all other GWB models, our model predicts a sharp decrease in the GWB amplitude at frequencies $f>3.5*10^{-8}$ Hz, reflecting the deficit of SBHs with masses $<4M_\odot$. The SBMH mergers at Big Crunch should generate yet undiscovered GWB with a frequency of $f\sim 10^{-(14-17)}$ Hz.


title: Power Line Communication and Sensing Using Time Series Forecasting

abstract: Smart electrical grids rely on data communication to support their operation and on sensing for diagnostics and maintenance. Usually, the roles of communication and sensing equipment are different, i.e., communication equipment does not participate in sensing tasks and vice versa. Power line communication (PLC) offers a cost-effective solution for joint communication and sensing for smart grids. This is because the high-frequency PLC signals used for data communication also reveal detailed information regarding the health of the power lines that they travel through. Traditional PLC-based power line or cable diagnostic solutions are dependent on prior knowledge of the cable type, network topology, and/or characteristics of the anomalies. In this paper, we develop a power line sensing technique that can detect various types of cable anomalies without any prior domain knowledge. To this end, we design a solution that first uses time-series forecasting to predict the PLC channel state information at any given point in time based on its historical data. Under the approximation that the prediction error follows a Gaussian distribution, we then perform chi-squared statistical test to build an anomaly detector which identifies the occurrence of a cable fault. We demonstrate the effectiveness and universality of our sensing solution via evaluations conducted using both synthetic and real-world data extracted from low- and medium-voltage distribution networks.


title: Patch Based Transformation for Minimum Variance Beamformer Image Approximation Using Delay and Sum Pipeline

abstract: In the recent past, there have been several efforts in accelerating computationally heavy beamforming algorithms such as minimum variance distortionless response (MVDR) beamforming to achieve real-time performance comparable to the popular delay and sum (DAS) beamforming. This has been achieved using a variety of neural network architectures ranging from fully connected neural networks (FCNNs), convolutional neural networks (CNNs) and general adversarial networks (GANs). However most of these approaches are working with optimizations considering image level losses and hence require a significant amount of dataset to ensure that the process of beamforming is learned. In this work, a patch level U-Net based neural network is proposed, where the delay compensated radio frequency (RF) patch for a fixed region in space (e.g. 32x32) is transformed through a U-Net architecture and multiplied with DAS apodization weights and optimized for similarity with MVDR image of the patch. Instead of framing the beamforming problem as a regression problem to estimate the apodization weights, the proposed approach treats the non-linear transformation of the RF data space that can account for the data driven weight adaptation done by the MVDR approach in the parameters of the network. In this way, it is also observed that by restricting the input to a patch the model will learn the beamforming pipeline as an image non-linear transformation problem.


title: The CoRa Tensor Compiler: Compilation for Ragged Tensors with Minimal Padding

abstract: There is often variation in the shape and size of input data used for deep learning. In many cases, such data can be represented using tensors with non-uniform shapes, or ragged tensors. Due to limited and non-portable support for efficient execution on ragged tensors, current deep learning frameworks generally use techniques such as padding and masking to make the data shapes uniform and then offload the computations to optimized kernels for dense tensor algebra. Such techniques can, however, lead to a lot of wasted computation and therefore, a loss in performance. This paper presents CoRa, a tensor compiler that allows users to easily generate efficient code for ragged tensor operators targeting a wide range of CPUs and GPUs. Evaluating CoRa on a variety of operators on ragged tensors as well as on an encoder layer of the transformer model, we find that CoRa (i)performs competitively with hand-optimized implementations of the operators and the transformer encoder and (ii) achieves, over PyTorch, a 1.6X geomean speedup for the encoder on an Nvidia GPU and a 1.86X geomean speedup for the multi-head attention module used in transformers on an ARM CPU.


title: Periods of the Long-Term Variability of the Blazar 0716+714 and Their Inter-Correlations in a Helical Jet Model

abstract: Various quasi-periods for the long-term variability of the radio emission, optical emission, and structural position angle of the inner part of the parsec-scale jet in the blazar 0716+714 have been detected. The relationships between these quasi-periods are interpreted assuming that the variability arises due to a helical structure of the jet, which is preserved from regions near the jet base to at least 1 milliarcsecond from the core observed in radio interferometric observations. The radiating jet components should display radial motions with Lorentz factors of $\approx3$, and decelerate with distance from the jet base. The best agreement with the data is given in the case of non-radial motions of these components with a constant physical speed. It is also shown that the helical shape of the jet strongly influences correlations both between fluxes observed in different spectral ranges and between the flux and position angle of the inner part of the parsec-scale jet.


title: A Federated Learning Aggregation Algorithm for Pervasive Computing: Evaluation and Comparison

abstract: Pervasive computing promotes the installation of connected devices in our living spaces in order to provide services. Two major developments have gained significant momentum recently: an advanced use of edge resources and the integration of machine learning techniques for engineering applications. This evolution raises major challenges, in particular related to the appropriate distribution of computing elements along an edge-to-cloud continuum. About this, Federated Learning has been recently proposed for distributed model training in the edge. The principle of this approach is to aggregate models learned on distributed clients in order to obtain a new, more general model. The resulting model is then redistributed to clients for further training. To date, the most popular federated learning algorithm uses coordinate-wise averaging of the model parameters for aggregation. However, it has been shown that this method is not adapted in heterogeneous environments where data is not identically and independently distributed (non-iid). This corresponds directly to some pervasive computing scenarios where heterogeneity of devices and users challenges machine learning with the double objective of generalization and personalization. In this paper, we propose a novel aggregation algorithm, termed FedDist, which is able to modify its model architecture (here, deep neural network) by identifying dissimilarities between specific neurons amongst the clients. This permits to account for clients' specificity without impairing generalization. Furthermore, we define a complete method to evaluate federated learning in a realistic way taking generalization and personalization into account. Using this method, FedDist is extensively tested and compared with three state-of-the-art federated learning algorithms on the pervasive domain of Human Activity Recognition with smartphones.


abstract: A graph $G$ is fully reconstructible in $\mathbb{C}^d$ if the graph is determined from its $d$-dimensional measurement variety. The full reconstructibility problem has been solved for $d=1$ and $d=2$. For $d=3$, some necessary and some sufficient conditions are known and $K_{5,5}$ falls squarely within the gap in the theory. In this paper, we show that $K_{5,5}$ is fully reconstructible in $\mathbb{C}^3$.


title: What Averages Do Not Tell -- Predicting Real Life Processes with Sequential Deep Learning

abstract: Deep Learning is proven to be an effective tool for modeling sequential data as shown by the success in Natural Language, Computer Vision and Signal Processing. Process Mining concerns discovering insights on business processes from their execution data that are logged by supporting information systems. The logged data (event log) is formed of event sequences (traces) that correspond to executions of a process. Many Deep Learning techniques have been successfully adapted for predictive Process Mining that aims to predict process outcomes, remaining time, the next event, or even the suffix of running traces. Traces in Process Mining are multimodal sequences and very differently structured than natural language sentences or images. This may require a different approach to processing. So far, there has been little focus on these differences and the challenges introduced. Looking at suffix prediction as the most challenging of these tasks, the performance of Deep Learning models was evaluated only on average measures and for a small number of real-life event logs. Comparing the results between papers is difficult due to different pre-processing and evaluation strategies. Challenges that may be relevant are the skewness of trace-length distribution and the skewness of the activity distribution in real-life event logs. We provide an end-to-end framework which enables to compare the performance of seven state-of-the-art sequential architectures in common settings. Results show that sequence modeling still has a lot of room for improvement for majority of the more complex datasets. Further research and insights are required to get consistent performance not just in average measures but additionally over all the prefixes.


title: A magnetic excitation linking quasi-1D Chevrel-type selenide and arsenide superconductors

abstract: The quasi-one-dimensional Chevrel phases, A$_2$Mo$_6$Se$_6$ (A = Tl, In, K, Rb, Cs), are of interest due to their atypical electronic properties. The Tl and In analogues undergo a superconducting transition whereas the alkali metal analogues show charge gapping of another, not well understood type. We report the results of inelastic neutron scattering on polycrystalline In$_2$Mo$_6$Se$_6$ ($T_c=2.85\,$K) and Rb$_2$Mo$_6$Se$_6$ (non-superconducting) samples, which reveal a column of intensity with linear dispersion from [0 0 1/2] to [0 0 1] in both compounds. The observed temperature and |$Q$| independence together suggest the presence of unconventional carriers with a spin contribution to the excitation. This is contrary to the prevailing model for these materials, which is that they are non-magnetic. The excitation has similar dispersion and $S(Q,E,T)$ behavior as one observed in the structurally related superconducting compounds A$_2$Cr$_3$As$_3$ and A$_2$Mo$_3$As$_3$ (A = K, Rb, Cs), which has been interpreted as magnetic in origin and related to Fermi surface nesting. The connection is unexpected because the calculated Fermi surface of the arsenides differs substantially from the A$_2$Mo$_6$Se$_6$ compounds, and many consider them distinct classes of materials. The new observation suggests a hidden link in the physics between both classes of superconductors, perhaps originating from their quasi-low-dimensional character.


title: On the uniform Besov regularity of local times of general processes

abstract: Our main purpose is to use a new condition, $\alpha$-local nondeterminism, which is an alternative to the classical local nondeterminism usually utilized in the Gaussian framework, in order to investigate Besov regularity, in the time variable $t$ uniformly in the space variable $x$, for local times $L(x, t)$ of a class of continuous processes. We also extend the classical Adler's theorem [1, Theorem 8.7.1] to the Besov spaces case. These results are then exploited to study the Besov irregularity of the sample paths of the underlying processes. Based on similar known results in the case of the bifractional Brownian motion, we believe that our results are sharp. As applications, we get sharp Besov regularity results for some classical Gaussian processes and the solutions of systems of non-linear stochastic heat equations. The Besov regularity of their corresponding local times is also obtained.


title: A Measurement of Proton, Deuteron, Triton and Alpha Particle Emission after Nuclear Muon Capture on Al, Si and Ti with the AlCap Experiment

abstract: Heavy charged particles after nuclear muon capture are an important nuclear physics background to the muon-to-electron conversion experiments Mu2e and COMET, which will search for charged lepton flavor violation at an unprecedented level of sensitivity. The AlCap experiment measured the yield and energy spectra of protons, deuterons, tritons, and alpha particles emitted after the nuclear capture of muons stopped in Al, Si, and Ti in the low energy range relevant for the muon-to-electron conversion experiments. Individual charged particle types were identified in layered silicon detector packages and their initial energy distributions were unfolded from the observed energy spectra. Detailed information on yields and energy spectra for all observed nuclei are presented in the paper.


title: Mathematical model of subthalamic nucleus neuron: characteristic activity patterns and bifurcation analysis

abstract: The subthalamic nucleus (STN) has an important role in the pathophysiology of the basal ganglia in Parkinson's disease. The ability of STN cells to generate bursting rhythms under either transient or sustained hyperpolarization may underlie the excessively synchronous beta rhythms observed in Parkinson's disease. In this study, we developed a conductance-based single compartment model of an STN neuron, which is able to generate characteristic activity patterns observed in experiments including hyperpolarization-induced bursts and post-inhibitory rebound bursts. This study focused on the role of three currents in rhythm generation: T-type calcium (CaT) current, L-type calcium (CaL) current, and hyperpolarization-activated cyclic nucleotide-gated (HCN) current. To investigate the effects of these currents in rhythm generation, we performed a bifurcation analysis using slow variables in these currents. Bifurcation analysis showed that the HCN current promotes single-spike activity patterns rather than bursting in agreement with experimental results. It also showed that the CaT current is necessary for characteristic bursting activity patterns. In particular, the CaT current enables STN neurons to generate these activity patterns under hyperpolarizing stimuli. The CaL current enriches and reinforces these characteristic activity patterns. In hyperpolarization-induced bursts or post-inhibitory rebound bursts, the CaL current allows STN neurons to generate long bursting patterns. Thus, bifurcation analysis explained the synergistic interaction of the CaT and CaL currents, which enables STN neurons to respond to hyperpolarizing stimuli in a salient way. The results of this study implicate the importance of CaT and CaL currents in the pathophysiology of the basal ganglia in Parkinson's disease.


title: Optimally Targeting Interventions in Networks during a Pandemic: Theory and Evidence from the Networks of Nursing Homes in the United States

abstract: This study develops an economic model for a social planner who prioritizes health over short-term wealth accumulation during a pandemic. Agents are connected through a weighted undirected network of contacts, and the planner's objective is to determine the policy that contains the spread of infection below a tolerable incidence level, and that maximizes the present discounted value of real income, in that order of priority. The optimal unique policy depends both on the configuration of the contact network and the tolerable infection incidence. Comparative statics analyses are conducted: (i) they reveal the tradeoff between the economic cost of the pandemic and the infection incidence allowed; and (ii) they suggest a correlation between different measures of network centrality and individual lockdown probability with the correlation increasing with the tolerable infection incidence level. Using unique data on the networks of nursing and long-term homes in the U.S., we calibrate our model at the state level and estimate the tolerable COVID-19 infection incidence level. We find that laissez-faire (more tolerance to the virus spread) pandemic policy is associated with an increased number of deaths in nursing homes and higher state GDP growth. In terms of the death count, laissez-faire is more harmful to nursing homes than more peripheral in the networks, those located in deprived counties, and those who work for a profit. We also find that U.S. states with a Republican governor have a higher level of tolerable incidence, but policies tend to converge with high death count.


title: Essential Surfaces in the Exterior of K13n586

abstract: We count the number of isotopy classes of closed, connected, orientable, essential surfaces embedded in the exterior B of the knot K13n586.The main result is that the count of surfaces by genus is equal to the Euler totent function. This is the first manifold for which we know the number of surfaces for any genus. The main argument is to show when normal surfaces in B are connected by counting their number of components. We implement tools from Agol, Hass and Thurston to convert the problem of counting components of surfaces into counting the number of orbits in a set of integers under a collection of bijections defined on its subsets.


title: Test time Adaptation through Perturbation Robustness

abstract: Data samples generated by several real world processes are dynamic in nature \textit{i.e.}, their characteristics vary with time. Thus it is not possible to train and tackle all possible distributional shifts between training and inference, using the host of transfer learning methods in literature. In this paper, we tackle this problem of adapting to domain shift at inference time \textit{i.e.}, we do not change the training process, but quickly adapt the model at test-time to handle any domain shift. For this, we propose to enforce consistency of predictions of data sampled in the vicinity of test sample on the image manifold. On a host of test scenarios like dealing with corruptions (CIFAR-10-C and CIFAR-100-C), and domain adaptation (VisDA-C), our method is at par or significantly outperforms previous methods.


title: Forecasting Market Prices using DL with Data Augmentation and Meta-learning: ARIMA still wins!

abstract: Deep-learning techniques have been successfully used for time-series forecasting and have often shown superior performance on many standard benchmark datasets as compared to traditional techniques. Here we present a comprehensive and comparative study of performance of deep-learning techniques for forecasting prices in financial markets. We benchmark state-of-the-art deep-learning baselines, such as NBeats, etc., on data from currency as well as stock markets. We also generate synthetic data using a fuzzy-logic based model of demand driven by technical rules such as moving averages, which are often used by traders. We benchmark the baseline techniques on this synthetic data as well as use it for data augmentation. We also apply gradient-based meta-learning to account for non-stationarity of financial time-series. Our extensive experiments notwithstanding, the surprising result is that the standard ARIMA models outperforms deep-learning even using data augmentation or meta-learning. We conclude by speculating as to why this might be the case.


title: Collaboration Challenges in Building ML-Enabled Systems: Communication, Documentation, Engineering, and Process

abstract: The introduction of machine learning (ML) components in software projects has created the need for software engineers to collaborate with data scientists and other specialists. While collaboration can always be challenging, ML introduces additional challenges with its exploratory model development process, additional skills and knowledge needed, difficulties testing ML systems, need for continuous evolution and monitoring, and non-traditional quality requirements such as fairness and explainability. Through interviews with 45 practitioners from 28 organizations, we identified key collaboration challenges that teams face when building and deploying ML systems into production. We report on common collaboration points in the development of production ML systems for requirements, data, and integration, as well as corresponding team patterns and challenges. We find that most of these challenges center around communication, documentation, engineering, and process and collect recommendations to address these challenges.


title: Entanglement Thresholds of Doubly-Parametric Quantum Transducers

abstract: Doubly-parametric quantum transducers, such as electro-opto-mechanical devices, are quickly approaching quantum operation as decoherence mechanisms such as thermal noise, loss, and limited cooperativities are improved. These devices show potential as the critical link between quantum information contained at frequencies as disparate as those in the optical and microwave domains, thus enabling applications such as long distance networking of superconducting quantum computers. However, the requirements on the operating parameters of the transducers necessary to achieve quantum operation have yet to be characterized. In this work we find simple, explicit expressions for the necessary and sufficient conditions under which doubly-parametric transducers in the resolved-sideband, steady-state limit are capable of entangling optical and microwave modes. Our analysis treats the transducer as a two-mode bosonic Gaussian channel capable of both beamsplitter-type and two-mode squeezing-type interactions between optical and microwave modes. For the beamsplitter-type interaction, we find parameter thresholds which distinguish regions of the channel's separability, capacity for bound entanglement, and capacity for distillable entanglement. By contrast, the two-mode squeezing-type interaction always produces distillable entanglement with no restrictions on temperature, cooperativities, or losses. Finally, we find the entanglement breaking conditions on the reduced one-mode upconversion and downconversion channels obtained by initializing the unused input with vacuum while ignoring the unused output. These differences between the entanglement thresholds of the beamsplitter-type and two-mode squeezing-type interactions are then important considerations in the construction of larger quantum networks that integrate multiple transducers.


title: Optimal Sequential Stochastic Deployment of Multiple Passenger Robots

abstract: We present a new algorithm for deploying passenger robots in marsupial robot systems. A marsupial robot system consists of a carrier robot (e.g., a ground vehicle), which is highly capable and has a long mission duration, and at least one passenger robot (e.g., a short-duration aerial vehicle) transported by the carrier. We optimize the performance of passenger robot deployment by proposing an algorithm that reasons over uncertainty by exploiting information about the prior probability distribution of features of interest in the environment. Our algorithm is formulated as a solution to a sequential stochastic assignment problem (SSAP). The key feature of the algorithm is a recurrence relationship that defines a set of observation thresholds that are used to decide when to deploy passenger robots. Our algorithm computes the optimal policy in $O(NR)$ time, where $N$ is the number of deployment decision points and $R$ is the number of passenger robots to be deployed. We conducted drone deployment exploration experiments on real-world data from the DARPA Subterranean challenge to test the SSAP algorithm. Our results show that our deployment algorithm outperforms other competing algorithms, such as the classic secretary approach and baseline partitioning methods, and is comparable to an offline oracle algorithm.


title: Stochastic Assignment for Deploying Multiple Marsupial Robots

abstract: Marsupial robot teams consist of carrier robots that transport and deploy multiple passenger robots, such as a team of ground robots that carry and deploy multiple aerial robots, to rapidly explore complex environments. We specifically address the problem of planning the deployment times and locations of the carrier robots to best meet the objectives of a mission while reasoning over uncertain future observations and rewards. While prior work proposed optimal, polynomial-time solutions to single-carrier robot systems, the multiple-carrier robot deployment problem is fundamentally harder as it requires addressing conflicts and dependencies between deployments of multiple passenger robots. We propose a centralized heuristic search algorithm for the multiple-carrier robot deployment problem that combines Monte Carlo Tree Search with a dynamic programming-based solution to the Sequential Stochastic Assignment Problem as a rollout action-selection policy. Our results with both procedurally-generated data and data drawn from the DARPA Subterranean Challenge Urban Circuit show the viability of our approach and substantial exploration performance improvements over alternative algorithms.


title: Native defects in monolayer GaS and GaSe: electrical properties and thermodynamic stability

abstract: Structural, electronic and thermodynamic properties of native defects in GaS and GaSe monolayers are investigated by means of accurate ab-initio calculations. Based on their charge transition levels we assess the influence of the studied defects on the electrical properties of the monolayers. Specifically, we show that native defects do not behave as shallow dopants and their presence cannot account for the experimentally observed intrinsic doping. In addition, we predict that native defects are efficient compensation and recombination centers. Besides pointing out their detrimental nature, we also calculate the corresponding finite temperature formation energies and provide a window of growth conditions able to reduce the concentration of all relevant native defects.


title: 1st Place Solution for the UVO Challenge on Image-based Open-World Segmentation 2021

abstract: We describe our two-stage instance segmentation framework we use to compete in the challenge. The first stage of our framework consists of an object detector, which generates object proposals in the format of bounding boxes. Then, the images and the detected bounding boxes are fed to the second stage, where a segmentation network is applied to segment the objects in the bounding boxes. We train all our networks in a class-agnostic way. Our approach achieves the first place in the UVO 2021 Image-based Open-World Segmentation Challenge.


title: Densities probabilities of a Bose-Fermi Mixture in 1D Double well potential

abstract: We use the two mode approximation for a interacting one-dimensional spinless soft core bosons and one half spin fermions in a double-well potential with a large central barrier. We include all the on-site boson-boson, fermion-fermion and boson-fermion repulsive contact potential represented by delta-function and considered bosonic and fermionic isotopes of ytterbium(Yb) $^{170}Yb$ and $^{171}Yb$ respectively. By means of the approximation, we find that in the regime $U_{BF}>U_{BB}$ give rise to a immiscible phase and in the regime $U_{BB}\geq U_{BF}$ give rise to a miscible phase, that is characterized by a temporal overlap of the bosonic and fermionic probability densities. We also report that due to the Bose-Fermi interaction, the system presents an apparent destruction of the collapse-revival oscillation of boson density probability at least in the ranges investigated.


title: Structural effect on phonon attenuation in metallic liquids and glasses

abstract: The attenuation rate of vibrational excitations in various metallic liquids and glasses has been reported to change from the quadratic dependence on wavevector at low wavevectors to the linear dependence at high wavevectors. However, the origin of this behavior is not clear. Here, the analysis of this phenomenon through molecular dynamics is presented for prototypical metallic liquids, Cu56Zr44 and Fe. It is shown that the crossover wavevector is strongly correlated with the structural coherence length characterizing coarse-grained density correlations. We suggest that the linear dependence is caused by scattering of vibrational excitations by structural activation processes with low activation energies which are distinctively observed in metallic systems.


title: A New Automatic Change Detection Frame-work Based on Region Growing and Weighted Local Mutual Information: Analysis of Breast Tumor Response to Chemotherapy in Serial MR Images

abstract: The automatic analysis of subtle changes between longitudinal MR images is an important task as it is still a challenging issue in scope of the breast medical image processing. In this paper we propose an effective automatic change detection framework composed of two phases since previously used methods have features with low distinctive power. First, in the preprocessing phase an intensity normalization method is suggested based on Hierarchical Histogram Matching (HHM) that is more robust to noise than previous methods. To eliminate undesirable changes and extract the regions containing significant changes the proposed Extraction Region of Changes (EROC) method is applied based on intensity distribution and Hill-Climbing algorithm. Second, in the detection phase a region growing-based approach is suggested to differentiate significant changes from unreal ones. Due to using proposed Weighted Local Mutual Information (WLMI) method to extract high level features and also utilizing the principle of the local consistency of changes, the proposed approach enjoys reasonable performance. The experimental results on both simulated and real longitudinal Breast MR Images confirm the effectiveness of the proposed framework. Also, this framework outperforms the human expert in some cases which can detect many lesion evolutions that are missed by expert.


title: RUNMON-RIFT: Adaptive Configuration and Healing for Large-Scale Parameter Inference

abstract: Gravitational wave parameter inference pipelines operate on data containing unknown sources on distributed hardware with unreliable performance. For one specific analysis pipeline (RIFT), we have developed a flexible tool (RUNMON-RIFT) to mitigate the most common challenges introduced by these two uncertainties. On the one hand, RUNMON provides several mechanisms to identify and redress unreliable computing environments. On the other hand, RUNMON provides mechanisms to adjust pipeline-specific run settings, including prior ranges, to ensure the analysis completes and encompasses the physical source parameters. We demonstrate both general features with two controlled examples.


title: Schmidt rank of quartics over perfect fields

abstract: Let $k$ be a perfect field of characteristic $\neq 2$. We prove that the Schmidt rank (also known as strength) of a quartic polynomial $f$ over $k$ is bounded above in terms of only the Schmidt rank of $f$ over $\overline{k}$, an algebraic closure of $k$.


title: Regret Minimization in Isotonic, Heavy-Tailed Contextual Bandits via Adaptive Confidence Bands

abstract: In this paper we initiate a study of non parametric contextual bandits under shape constraints on the mean reward function. Specifically, we study a setting where the context is one dimensional, and the mean reward function is isotonic with respect to this context. We propose a policy for this problem and show that it attains minimax rate optimal regret. Moreover, we show that the same policy enjoys automatic adaptation; that is, for subclasses of the parameter space where the true mean reward functions are also piecewise constant with $k$ pieces, this policy remains minimax rate optimal simultaneously for all $k \geq 1.$ Automatic adaptation phenomena are well-known for shape constrained problems in the offline setting; %The phenomenon of automatic adaptation of shape constrained methods is known to occur in offline problems; we show that such phenomena carry over to the online setting. The main technical ingredient underlying our policy is a procedure to derive confidence bands for an underlying isotonic function using the isotonic quantile estimator. The confidence band we propose is valid under heavy tailed noise, and its average width goes to $0$ at an adaptively optimal rate. We consider this to be an independent contribution to the isotonic regression literature.


title: 2020 State of the Octoverse: Securing the World's Software

abstract: Open source is the connective tissue for much of the information economy. You would be hard-pressed to find a scenario where your data does not pass through at least one open source component. Many of the services and technology we all rely on, from banking to healthcare, also rely on open source software. The artifacts of open source code serve as critical i infrastructure for much of the global economy, making the security of open source software mission-critical to the world.


title: Laser Communication with Proxima and Alpha Centauri using the Solar Gravitational Lens

abstract: A search was conducted for laser signals, both sub-second pulses and continuous emission, from the regions of the sky opposite Proxima and Alpha Centauri. These regions are located at the foci of the gravitational lensing caused by the Sun, ideal for amplifying transmissions between our Solar System and those two nearest stellar neighbors. During six months in 2020 and 2021, 88000 exposures for Proxima Cen and 47000 exposures for Alpha Cen were obtained. No evidence was detected of light pulses or continuous laser emission in the wavelength range of 380 to 950 nm. We would have detected a laser having a power of just 100 Watts.


title: 2020 State of the Octoverse: Finding Balance Between Work and Play

abstract: Over the past year, many developers and other technology professionals have transitioned to a remote-first world, as COVID-19 pressed organizations to support working from home whenever possible. This shift quickly changed the routines and environments where we work and learn, redrawing the lines between personal and professional lives. How does this affect the ways we develop and deliver software, both at work and in our open source projects?


title: Neural Stochastic PDEs: Resolution-Invariant Learning of Continuous Spatiotemporal Dynamics

abstract: Stochastic partial differential equations (SPDEs) are the mathematical tool of choice for modelling spatiotemporal PDE-dynamics under the influence of randomness. Based on the notion of mild solution of an SPDE, we introduce a novel neural architecture to learn solution operators of PDEs with (possibly stochastic) forcing from partially observed data. The proposed Neural SPDE model provides an extension to two popular classes of physics-inspired architectures. On the one hand, it extends Neural CDEs and variants -- continuous-time analogues of RNNs -- in that it is capable of processing incoming sequential information arriving at arbitrary spatial resolutions. On the other hand, it extends Neural Operators -- generalizations of neural networks to model mappings between spaces of functions -- in that it can parameterize solution operators of SPDEs depending simultaneously on the initial condition and a realization of the driving noise. By performing operations in the spectral domain, we show how a Neural SPDE can be evaluated in two ways, either by calling an ODE solver (emulating a spectral Galerkin scheme), or by solving a fixed point problem. Experiments on various semilinear SPDEs, including the stochastic Navier-Stokes equations, demonstrate how the Neural SPDE model is capable of learning complex spatiotemporal dynamics in a resolution-invariant way, with better accuracy and lighter training data requirements compared to alternative models, and up to 3 orders of magnitude faster than traditional solvers.


title: Bayesian spatial voting model to characterize the legislative behavior of the Colombian Senate 2010-2014

abstract: This paper applies Bayesian methodologies to characterize the legislative behavior of the Colombian Senate during the period 2010-2014. The analysis is carried out through the plenary roll call votes of this legislative chamber. In addition, parliamentary electoral behavior is operationalized by implementing the one-dimensional standard Bayesian ideal point estimator via the Markov chain Monte Carlo algorithms. The results contribute mainly to two points: political space dimensionality and the identification of pivot legislators. The pattern revealed by the estimated ideal points suggests a latent non-ideological trait (opposition - non-opposition) underlying the vote of deputies in the Senate. Thus, in addition to providing empirical evidence for a better understanding of legislative policy in Colombia during the period under analysis, this work also offers methodological and theoretical tools to guide the analysis of roll call vote data in contexts of unbalanced parliaments (as opposed to the U.S. parliament), taking the particular case of the Colombian Senate as a reference.


title: Higher Coleman Theory

abstract: We develop local cohomology techniques to study the finite slope part of the coherent cohomology of Shimura varieties. The local cohomology groups we consider are a generalization of overconvergent modular forms, and they are defined by using a stratification on the Shimura variety obtained from the Bruhat stratification on a flag variety via the Hodge-Tate period map. We construct a spectral sequence from the local cohomologies to the classical cohomology and use it to obtain classicality and vanishing results. We also develop a theory of p-adic families and construct eigenvarieties. As an application, we prove some new properties of Galois representations arising from certain non-regular algebraic cuspidal automorphic representations.


title: An inductive approach to generalized abundance using nef reduction

abstract: We use the canonical bundle formula for parabolic fibrations to give an inductive approach to the generalized abundance conjecture using nef reduction. In particular, we observe that generalized abundance holds for a klt pair $(X,B)$ if the nef dimension $n(K_X+B+L)=2$ and $K_X+B \geq 0$ or $n(K_X+B+L)=3$ and $\kappa(K_X+B )>0$.


title: First global QCD analysis of the TMD $\boldsymbol{g_{1T}}$ from semi-inclusive DIS data

abstract: The worm-gear transverse momentum dependent (TMD) function $g_{1T}$ is one of the eight leading-twist TMDs and has the probabilistic interpretation of finding a longitudinally polarized quark inside a transversely polarized hadron. In this work, we present the first simultaneous extraction of $g_{1T}$ from all the available experimental measurements. The study involves the analysis of COMPASS, HERMES, and Jefferson Lab data on semi-inclusive deep-inelastic scattering using Monte Carlo techniques. We also compare $g_{1T}$ obtained from this experimental data with different theoretical approaches, such as the large-$N_c$ approximation, the Wandzura-Wilczek-type approximation, and lattice QCD.


title: Approaching the Chiral and Continuum Limit of Large-N QCD

abstract: We present preliminary results from our calculation of the low energy constants (LECs) of the chiral effective theory for 3, 4 and 5 color QCD with $N_f=2$ dynamical fermion flavors. We simulate with clover fermions over a range of lattice couplings and quark masses. We observe the expected $N_c$ scaling for the LECs appropriate to the condensate $B$ and pseudoscalar decay constant $F$. The range of quark masses over which leading order chiral perturbation theory describes the data grows as $N_c$ rises.


title: A Simple Approach to Continual Learning by Transferring Skill Parameters

abstract: In order to be effective general purpose machines in real world environments, robots not only will need to adapt their existing manipulation skills to new circumstances, they will need to acquire entirely new skills on-the-fly. A great promise of continual learning is to endow robots with this ability, by using their accumulated knowledge and experience from prior skills. We take a fresh look at this problem, by considering a setting in which the robot is limited to storing that knowledge and experience only in the form of learned skill policies. We show that storing skill policies, careful pre-training, and appropriately choosing when to transfer those skill policies is sufficient to build a continual learner in the context of robotic manipulation. We analyze which conditions are needed to transfer skills in the challenging Meta-World simulation benchmark. Using this analysis, we introduce a pair-wise metric relating skills that allows us to predict the effectiveness of skill transfer between tasks, and use it to reduce the problem of continual learning to curriculum selection. Given an appropriate curriculum, we show how to continually acquire robotic manipulation skills without forgetting, and using far fewer samples than needed to train them from scratch.


title: Estimating energy levels of a three-level atom in single and multi-parameter metrological schemes

abstract: Determining the energy levels of a quantum system is a significant task, for instance, to analyze reaction rates in drug discovery and catalysis or characterize the compatibility of materials. In this paper we exploit quantum metrology, the research field focusing on the estimation of unknown parameters exploiting quantum resources, to address this problem for a three-level system interacting with laser fields. The performance of simultaneous estimation of the levels compared to independent one is also investigated in various scenarios. Moreover, we introduce, the Hilbert-Schmidt speed (HSS), a special type of quantum statistical speed, as a powerful figure of merit for enhancing estimation of energy spectrum. This measure is easily computable, because it does not require diagonalization of the system state, verifying its efficiency in high-dimensional systems.


title: Koszul duality in quantum field theory

abstract: In this article, we introduce basic aspects of the algebraic notion of Koszul duality for a physics audience. We then review its appearance in the physical problem of coupling QFTs to topological line defects, and illustrate the concept with some examples drawn from twists of various simple supersymmetric theories. Though much of the content of this article is well-known to experts, the presentation and examples have not, to our knowledge, appeared in the literature before. Our aim is to provide an elementary introduction for those interested in the appearance of Koszul duality in supersymmetric gauge theories with line defects and, ultimately, its generalizations to higher-dimensional defects and twisted holography.


title: Scalar field dark matter with two components: combined approach from particle physics and cosmology

abstract: In this work we explore the possibility of incorporating particle physics motivated scalar fields to the dark matter cosmological model. In this landscape, we consider the classical complex scalar field in a certain region in the parameter space of the model that increases the number of neutrino species $N_{\mathrm{eff}}$, in order to be consistent with the observed abundance of light elements produced at Big Bang Nucleosynthesis (BBN). We perform analyses using one and two scalar fields. We examine the difference between these models and the priors considered at the edges of the cosmic ladder, this with the purpose of studying the impact of such models on the Hubble cosmic flow. In the two scalar field models we explore the possibility of combining an axion and a Higgs-like field as well as a Higgs-like field and the classical field, we show that in the first case there is no set of parameters that allows us to be consistent with $N_\mathrm{eff}$, while in the second case a strong restriction to the set of parameters is obtained. This last restriction is given in terms of a maximum bound of the fraction of Higgs-like field that can be incorporated together with the classical field. Our results could be relevant in the direct dark matter detection programs.


title: A Neural-Network-defined Gaussian Mixture Model for particle identification applied to the LHCb fixed-target programme

abstract: Particle identification in large high-energy physics experiments typically relies on classifiers obtained by combining many experimental observables. Predicting the probability density function (pdf) of such classifiers in the multivariate space covering the relevant experimental features is usually challenging. The detailed simulation of the detector response from first principles cannot provide the reliability needed for the most precise physics measurements. Data-driven modelling is usually preferred, though sometimes limited by the available data size and different coverage of the feature space by the control channels. In this paper, we discuss a novel approach to the modelling of particle identification classifiers using machine-learning techniques. The marginal pdf of the classifiers is described with a Gaussian Mixture Model, whose parameters are predicted by Multi Layer Perceptrons trained on calibration data. As a proof of principle, the method is applied to the data acquired by the LHCb experiment in its fixed-target configuration. The model is trained on a data sample of proton-neon collisions and applied to smaller data samples of proton-helium and proton-argon collisions collected at different centre-of-mass energies. The method is shown to perform better than a detailed simulation-based approach, to be fast and suitable to be applied to a large variety of use cases.


title: Collective gradient sensing with limited positional information

abstract: Eukaryotic cells sense chemical gradients to decide where and when to move. Clusters of cells can sense gradients more accurately than individual cells by integrating measurements of the concentration made across the cluster. Is this gradient sensing accuracy impeded when cells have limited knowledge of their position within the cluster, i.e. limited positional information? We apply maximum likelihood estimation to study gradient sensing accuracy of a cluster of cells with finite positional information. If cells must estimate their location within the cluster, this lowers the accuracy of collective gradient sensing. We compare our results with a tug-of-war model where cells respond to the gradient by a mechanism of collective guidance without relying on their positional information. As the cell positional uncertainty increases, there is a trade-off where the tug-of-war model responds more accurately to the chemical gradient. However, for sufficiently large cell clusters or shallow chemical gradients, the tug-of-war model will always be suboptimal to one that integrates information from all cells, even if positional uncertainty is high.


title: Learning Domain Specific Language Models for Automatic Speech Recognition through Machine Translation

abstract: Automatic Speech Recognition (ASR) systems have been gaining popularity in the recent years for their widespread usage in smart phones and speakers. Building ASR systems for task-specific scenarios is subject to the availability of utterances that adhere to the style of the task as well as the language in question. In our work, we target such a scenario wherein task-specific text data is available in a language that is different from the target language in which an ASR Language Model (LM) is expected. We use Neural Machine Translation (NMT) as an intermediate step to first obtain translations of the task-specific text data. We then train LMs on the 1-best and N-best translations and study ways to improve on such a baseline LM. We develop a procedure to derive word confusion networks from NMT beam search graphs and evaluate LMs trained on these confusion networks. With experiments on the WMT20 chat translation task dataset, we demonstrate that NMT confusion networks can help to reduce the perplexity of both n-gram and recurrent neural network LMs compared to those trained only on N-best translations.


title: Dark Matter Axion Search Using a Josephson Traveling Wave Parametric Amplifier

abstract: We present a new exclusion bound of axion-like particle dark matter with axion-photon couplings above $\mathrm{10^{-13}}$ $\mathrm{GeV^{-1}}$ over the frequency range 4796.7--4799.5 MHz, corresponding to a narrow range of axion masses centered around 19.84 $\mu$eV. This measurement represents the first implementation of a Josephson Traveling Wave Parametric Amplifier (JTWPA) in a dark matter search. The JTWPA was operated in the insert of the Axion Dark Matter eXperiment (ADMX) as part of an independent receiver chain that was attached to a 0.588-liter cavity. The ability of the JTWPA to deliver high gain over a wide (3 GHz) bandwidth has engendered interest from those aiming to perform broadband axion searches, a longstanding goal in this field.


title: Current status of PAPYRUS : the pyramid based adaptive optics system at LAM/OHP

abstract: The Provence Adaptive optics Pyramid Run System (PAPYRUS) is a pyramid-based Adaptive Optics (AO) system that will be installed at the Coude focus of the 1.52m telescope (T152) at the Observatoire de Haute Provence (OHP). The project is being developed by PhD students and Postdocs across France with support from staff members consolidating the existing expertise and hardware into an R&D testbed. This testbed allows us to run various pyramid wavefront sensing (WFS) control algorithms on-sky and experiment on new concepts for wavefront control with additional benefit from the high number of available nights at this telescope. It will also function as a teaching tool for students during the planned AO summer school at OHP. To our knowledge, this is one of the first pedagogic pyramid-based AO systems on-sky. The key components of PAPYRUS are a 17x17 actuators Alpao deformable mirror with a Alpao RTC, a very low noise camera OCAM2k, and a 4-faces glass pyramid. PAPYRUS is designed in order to be a simple and modular system to explore wavefront control with a pyramid WFS on sky. We present an overview of PAPYRUS, a description of the opto-mechanical design and the current status of the project.


title: $r$-Hyperideals and Generalizations of $r$-Hyperideals in Krasner Hyperrings

abstract: In this study, we examine some properties of $r$-hyperideals in the commutative Krasner hyperrings. Some properties of $pr$-hyperideals are also studied. The relation between prime hyperideals and $r$-hyperideals is investigated. We show that the image and the inverse image of an $r$-hyperideal is also an $r$-hyperideal. We also introduce a generalization of r-hyperideals and we prove some properties of them.


title: H\"older continuity of the Lyapunov exponents of linear cocycles over hyperbolic maps

abstract: Given a hyperbolic homeomorphism on a compact metric space, consider the space of linear cocycles over this base dynamics which are H\"older continuous and whose projective actions are partially hyperbolic dynamical systems. We prove that locally near any typical cocycle, the Lyapunov exponents are H\"older continuous functions relative to the uniform topology. This result is obtained as a consequence of a uniform large deviations type estimate in the space of cocycles. As a byproduct of our approach, we also establish other statistical properties for the iterates of such cocycles, namely a central limit theorem and a large deviations principle.


title: Addressing Positivity Violations in Causal Effect Estimation using Gaussian Process Priors

abstract: In observational studies, causal inference relies on several key identifying assumptions. One identifiability condition is the positivity assumption, which requires the probability of treatment be bounded away from 0 and 1. That is, for every covariate combination, it should be possible to observe both treated and control subjects, i.e., the covariate distributions should overlap between treatment arms. If the positivity assumption is violated, population-level causal inference necessarily involves some extrapolation. Ideally, a greater amount of uncertainty about the causal effect estimate should be reflected in such situations. With that goal in mind, we construct a Gaussian process model for estimating treatment effects in the presence of practical violations of positivity. Advantages of our method include minimal distributional assumptions, a cohesive model for estimating treatment effects, and more uncertainty associated with areas in the covariate space where there is less overlap. We assess the performance of our approach with respect to bias and efficiency using simulation studies. The method is then applied to a study of critically ill female patients to examine the effect of undergoing right heart catheterization.


title: Recognizability of morphisms

abstract: We investigate several questions related to the notion of recognizable morphism. The main result is a new proof of Moss\'e's theorem and actually of a generalization to non primitive morphisms due to Berth\'e et al. We actually prove the result of Berth\'e et al. for the most general class of morphisms, including ones with erasable letters. It is derived from a result concerning elementary morphisms for which we also provide a new proof. We also show how to decide whether an injective morphism is recognizable on the full shift for aperiodic points.


title: Kinematic distributions of the $\eta_c$ leptoproduction in association with light hadrons

abstract: The $\eta_c$ meson leptoproduction is calculated within the nonrelativistic QCD framework for the first time. It is found that the colour-singlet channel, although suppressed by a factor of $\alpha_s$ relative to the colour-octet ones, provides the most important contribution for almost all the experimental conditions, which disagrees with some of the expectations before computation. We present the differential cross sections with respect to $p_t^2$, $p_t^{\star2}$, $Q^2$, $W$, and $z$, for both HERA and EIC experimental conditions as a reference for future studies. The scale dependence and long-distance-matrix-element dependence are also investigated in this paper.


title: Performance Bounds for PDE-Constrained Optimization under Uncertainty

abstract: Computational approaches to PDE-constrained optimization under uncertainty may involve finite-dimensional approximations of control and state spaces, sample average approximations of measures of risk and reliability, smooth approximations of nonsmooth functions, penalty approximations of constraints as well as many other kinds of inaccuracies. In this paper, we analyze the performance of controls obtained by an approximation-based algorithm and in the process develop estimates of optimality gaps for general optimization problems defined on metric spaces. Under mild assumptions, we establish that limiting controls have arbitrarily small optimality gaps provided that the inaccuracies in the various approximations vanish. We carry out the analysis for a broad class of problems with multiple expectation, risk, and reliability functions involving PDE solutions and appearing in objective as well as constraint expressions. In particular, we address problems with buffered failure probability constraints approximated via an augmented Lagrangian. We demonstrate the framework on an elliptic PDE with a random coefficient field and a distributed control function.


title: Thermal effects on tidal deformability of a coalescing binary neutron star system

abstract: The study of neutron star mergers by the detection of the emitted gravitational waves is one of the most promised tools to study the properties of dense nuclear matter at high densities. It is worth claiming that, at the moment, strong evidence that the temperature of the stars is zero during the last orbits before coalescing does not exist. Contrariwise, there are some theoretical predictions suggesting that the star's temperature might even be a few MeV. According to the main theory, the tides transfer mechanical energy and angular momentum to the star at the expense of the orbit, where friction within the star converts the mechanical energy into heat. During the inspiral these effects are potentially detectable. Different treatments have been used to estimate the transfer of the mechanical energy and the size of the tidal friction, leading to different conclusions about the importance of pre-merger tidal effects. The present work is dedicated to the study of the effect of temperature on the tidal deformability of neutron stars during the inspiral of a neutron star system just before the merger. We applied a class of hot equations of state originated from various nuclear models and found that even for low values of temperature ($T<1$ MeV) the effects on the basic ingredients of tidal deformability are not negligible. However, according to the main finding, the effect of the temperature on the tidal deformability is indistinguishable. The consequences of this unexpected result are discussed and analyzed.


title: The time-domain Landau-Zener-Stuckelberg-Majorana interference in optical lattice clock

abstract: The interference between a sequence of Landau-Zener (LZ) transitions can produce Rabi oscillations(LZROs). This phenomenon is a kind of time-domain Landau-Zener-Stuckelberg-Majorana (LZSM) interference. However, experimental demonstrations of this LZSM interference induced Rabi oscillation are extremely hard due to the short coherence time of the driven quantum system. Here, we study theoretically the time-domain LZSM interference between the clock transition in one-dimensional (1D) Sr-87 optical lattice clock (OLC) system. With the help of both the adiabatic-impulse model and Floquet numerical simulation method, the LZROs with special step-like structure are clearly found both in fast- and slow-passage limit in the real experiment parameter regions. In addition, the dephasing effect caused by the system temperature can be suppressed with destructive interference in the slow-passage limit. Finally, we discuss the possible Bloch-Siegert shift while the pulse time is away from the half-integer and integer periods. Our work provides a clear roadmap to observe the LZROs on the OLC platform.


title: An application of a small area procedure with correlation between measurement error and sampling error to the Conservation Effects Assessment Project

abstract: County level estimates of mean sheet and rill erosion from the Conservation Effects Assessment Project (CEAP) are useful for program development and evaluation. Since county sample sizes in the CEAP survey are insufficient to support reliable direct estimators, small area estimation procedures are needed. The quantity of water runoff is a useful covariate but is unavailable for the full population. We use an estimate of mean runoff from the CEAP survey as a covariate in a small area model with sheet and rill erosion as the response. As the runoff and sheet and rill erosion are estimators from the same survey, the measurement error in the covariate is important as is the correlation between the measurement error and the sampling error. We conduct a detailed investigation of small area estimation in the presence of a correlation between the measurement error in the covariate and the sampling error in the response. In simulations, the proposed predictor is superior to small area predictors that assume the response and covariate are uncorrelated or that ignore the measurement error entirely.


title: A Vertex Model for Supersymmetric LLT Polynomials

abstract: We describe a Yang-Baxter integrable vertex model, which can be realized as a degeneration of a vertex model introduced by Aggarwal, Borodin, and Wheeler. From this vertex model, we construct a certain class of partition functions that we show are essentially equal to the super ribbon functions of Lam. Using the vertex model formalism, we give proofs of many properties of these polynomials, namely a Cauchy identity and generalizations of known identities for supersymmetric Schur polynomials.


title: Effect of the electronic charge gap on LO bond-stretching phonons in undoped La$_2$CuO$_4$ calculated using LDA+U

abstract: Typical density-functional theory calculations that wrongly predict undoped cuprates to be metallic also predict Cu-O half- and full-breathing phonon energies that are significantly softer than observed, presumably because of weak on-site Coulomb repulsion on the Cu 3d orbitals. We used DFT+U calculations with antiferromagnetic supercells of La$_2$CuO$_4$ to establish correlation between the on-site repulsion strength, tuned via adjusting the value of U, and phonon dispersions. We find that breathing and half-breathing phonons reach experimental values when U is tuned to obtain the correct optical gap and magnetic moments. We demonstrate that using distorted supercells within DFT+U is a promising framework to model phonons in undoped cuprates and other perovskite oxides with complex, interrelated structural and electronic degrees of freedom.