title: Successive deuteration in low-mass star-forming regions: the case of D$_{2}$-methanol (CHD$_{2}$OH) in IRAS 16293-2422

abstract: Accurate quantification of the column density of di-deuterated methanol is a key missing puzzle piece in the otherwise thoroughly constrained family of D-bearing methanol in the deeply embedded low-mass protostellar system and astrochemical template source IRAS16293-2422. A spectroscopic dataset for astrophysical purposes is built for CHD$_{2}$OH and made publicly available to facilitate accurate characterization of this species in astrochemical surveys. The newly computed line list and partition function are used to search for CHD$_{2}$OH towards IRAS16293-2422 A and B in data from ALMA-PILS. Only non-blended, optically thin lines of CHD$_{2}$OH are used for the synthetic spectral fitting. The constructed spectroscopic database contains line frequencies and strengths for 7417 transitions in the 0 to 500 GHz frequency range. ALMA-PILS observations in the 329-363 GHz range are used to identify 105 unique, non-blended, optically thin line frequencies of CHD$_{2}$OH for synthetic spectral fitting. The derived excitation temperatures and column densities yield high D/H ratios of CHD$_{2}$OH in IRAS 16293-2422 A and B of 7.5$\pm$1.1% and 7.7$\pm$1.2%, respectively. Deuteration in IRAS 16293-2422 is not higher than in other low-mass star-forming regions. Di-deuterated molecules consistently have higher D/H ratios than their monodeuterated counterparts in all low-mass protostars, which may be a natural consequence of H-D substitution reactions as seen in laboratory experiments. The Solar System's natal cloud, as traced by comet 67P/Churyumov-Gerasimenko, may have had a lower initial abundance of D, been warmer than the cloud of IRAS16293-2422, or been partially reprocessed. In combination with accurate spectroscopy, a careful spectral analysis, and a consideration of the underlying assumptions, successive deuteration is a robust window on the physicochemical provenance of star-forming systems.


title: Spin generalizations of the Benjamin-Ono equation

abstract: We present new soliton equations related to the $A$-type spin Calogero-Moser (CM) systems introduced by Gibbons and Hermsen. These equations are spin generalizations of the Benjamin-Ono (BO) equation and the recently introduced non-chiral intermediate long-wave (ncILW) equation. We obtain multi-soliton solutions of these spin generalizations of the BO equation and the ncILW equation via a spin-pole ansatz where the spin-pole dynamics is governed by the spin CM system in the rational and hyperbolic cases, respectively. We also propose physics applications of the new equations, and we introduce a spin generalization of the standard intermediate long-wave equation which interpolates between the matrix Korteweg-de Vries equation, the Heisenberg ferromagnet equation, and the spin BO equation.


title: Random walk on nonnegative integers in beta distributed random environment

abstract: We consider random walks on the nonnegative integers in a space-time dependent random environment. We assume that transition probabilities are given by independent $\mathrm{Beta}(\mu,\mu)$ distributed random variables, with a specific behaviour at the boundary, controlled by an extra parameter $\eta$. We show that this model is exactly solvable and prove a formula for the mixed moments of the random heat kernel. We then provide two formulas that allow us to study the large-scale behaviour. The first involves a Fredholm Pfaffian, which we use to prove a local limit theorem describing how the boundary parameter $\eta$ affects the return probabilities. The second is an explicit series of integrals, and we show that non-rigorous critical point asymptotics suggest that the large deviation behaviour of this half-space random walk in random environment is the same as for the analogous random walk on $\mathbb{Z}$.

categories: math.PR cond-mat.stat-mech math-ph math.MP


title: A fast radio burst progenitor born in a galaxy merger

abstract: We report a Giant Metrewave Radio Telescope HI 21cm mapping study of the neutral atomic hydrogen (HI) in the host galaxy of the fast radio burst (FRB) FRB20180916B at $z \approx 0.03399$. We find that the FRB host has an HI mass of $\rm M_{HI} = (2.74 \pm 0.33) \times 10^9 \ M_\odot$ and a high HI-to-stellar mass ratio, $\approx 1.3$. The FRB host is thus a gas-rich but near-quiescent galaxy, that is likely to have acquired a significant mass of HI in the recent past. The HI distribution is disturbed, with extended HI 21cm emission detected in a north-eastern tail, a counter-tail towards the south, an HI hole between the galaxy centre and the FRB location, and a high HI column density measured close to the FRB position. The FRB host is part of a group with four companions detected in their HI 21cm emission, the nearest of which is only 22~kpc from the FRB location. The gas-richness and disturbed HI distribution indicate that the FRB host has recently undergone a minor merger, which increased its HI mass, disturbed the HI in the galaxy disk, and compressed the HI near the FRB location to increase its surface density. We propose that this merger caused the burst of star-formation in the outskirts of the galaxy that gave rise to the FRB progenitor. The evidence for a minor merger is consistent with scenarios in which the FRB progenitor is a massive star, formed due to the merger event.


title: Lambda the Ultimate SSA: Optimizing Functional Programs in SSA

abstract: Static Single Assignment (SSA) is the workhorse of modern optimizing compilers for imperative programming languages. However, functional languages have been slow to adopt SSA and prefer to use intermediate representations based on minimal lambda calculi due to SSA's inability to express higher order constructs. We exploit a new SSA construct -- regions -- in order to express functional optimizations via classical SSA based reasoning. Region optimization currently relies on ad-hoc analyses and transformations on imperative programs. These ad-hoc transformations are sufficient for imperative languages as regions are used in a limited fashion. In contrast, we use regions pervasively to model sub-expressions in our functional IR. This motivates us to systematize region optimizations. We extend classical SSA reasoning to regions for functional-style analyses and transformations. We implement a new SSA+regions based backend for LEAN4, a theorem prover that implements a purely functional, dependently typed programming language. Our backend is feature-complete and handles all constructs of LEAN4's functional intermediate representation {\lambda}rc within the SSA framework. We evaluate our proposed region optimizations by optimizing {\lambda}rc within an SSA+regions based framework implemented in MLIR and demonstrating performance parity with the current LEAN4 backend. We believe our work will pave the way for a unified optimization framework capable of representing, analyzing, and optimizing both functional and imperative languages.


title: Exotic hadrons in heavy ion collisions

abstract: Nucleus-nucleus collisions offer a great opportunity for analyzing and determining the intrinsic nature of heavy and exotic hadrons. In this sense, here we discuss how the production and dissociation of hadron states are affected by reactions during the expansion of hadronic matter in a heavy-ion collision environment. We give emphasis to recent works on exotic states, revisiting as a case study the time evolution of the abundances of the $X_{0,1}(2900)$ states.


title: Online Generation of Proofs Without Words

abstract: Understanding geometric relationships with little mathematical knowledge can be challenging for today's students and teachers. A new toolset is introduced that is able to create a proof without words by combining the benefits of the Geometric Deduction Database method (to obtain a readable proof of a geometric statement) and the GeoGebra framework (that makes it possible to export these data as an online applet in a simple way).


title: Automated Theorem Proving in the Classroom

abstract: We report on several scenarios of using automated theorem proving software in university education. In particular, we focus on using the Theorema system in a software-enhanced logic-course for students in computer science or artificial intelligence. The purpose of using logic-software in our teaching is not to teach students the proper use of a particular piece of software. In contrast, we try to employ certain software in order to spark students' motivation and to support their understanding of logic principles they are supposed to understand after having passed the course. In a sense, we try to let the software act as a logic-tutor, the software is not an additional subject we teach.


title: Large deviation principle for geometric and topological functionals and associated point processes

abstract: We prove a large deviation principle for the point process associated to $k$-element connected components in $\mathbb R^d$ with respect to the connectivity radii $r_n\to\infty$. The random points are generated from a homogeneous Poisson point process, so that $(r_n)_{n\ge1}$ satisfies $n^kr_n^{d(k-1)}\to\infty$ and $nr_n^d\to0$ as $n\to\infty$ (i.e., sparse regime). The rate function for the obtained large deviation principle can be represented as relative entropy. As an application, we deduce large deviation principles for various functionals and point processes appearing in stochastic geometry and topology. As concrete examples of topological invariants, we consider persistent Betti numbers of geometric complexes and the number of Morse critical points of the min-type distance function.


title: Measuring correlations from the collective spin fluctuations of a large ensemble of lattice-trapped dipolar spin-3 atoms

abstract: We perform collective spin measurements to study the buildup of two-body correlations between $\approx10^4$ spin $s=3$ chromium atoms pinned in a 3D optical lattice. The spins interact via long range and anisotropic dipolar interactions. From the fluctuations of total magnetization, measured at the standard quantum limit, we estimate the dynamical growth of the connected pairwise correlations associated with magnetization. The quantum nature of the correlations is assessed by comparisons with short and long time expansions, and numerical simulations. Our work shows that measuring fluctuations of spin populations provides new ways to characterize correlations in quantum many-body systems, for $s>1/2$ spins


title: Functons of perturbed pairs of dissipative operators

abstract: Let $f$ be a function in the inhomogeneous analytic Besov space $B_{\infty,1}^1$. For a pair $(L,M)$ of not necessarily commuting maximal dissipative operators, we define the function $f(L,M)$ of $L$ and $M$ as a densely defined linear operator. We prove for $p\in[1,2]$ that if $(L_1,M_1)$ and $(L_2,M_2)$ are pairs of not necessarily commuting maximal dissipative operators such that both differences $L_1-L_2$ and $M_1-M_2$ belong to the Schatten--von Neumann class $\boldsymbol{S}_p$ than for an arbitrary function $f$ in the inhomogeneous analytic Besov space $B_{\infty,1}^1$, the operator difference $f(L_1,M_1)-f(L_2,M_2)$ belongs to $\boldsymbol{S}_p$ and the following Lipschitz type estimate holds: $$ \|f(L_1,M_1)-f(L_2,M_2)\|_{\boldsymbol{S}_p} \le\operatorname{const}\|f\|_{B_{\infty,1}^1}\max\big\{\|L_1-L_2\|_{\boldsymbol{S}_p},\|M_1-M_2\|_{\boldsymbol{S}_p}\big\}. $$


title: Gravitational entropy and the flatness, homogeneity and isotropy puzzles

abstract: We suggest a new explanation for the observed large scale flatness, homogeneity and isotropy of the universe. The basic ingredients are elementary and well-known, namely Einstein's theory of gravity and Hawking's method of computing gravitational entropy. The new twist is provided by the boundary conditions we recently proposed for "big bang" type singularities dominated by conformal matter, enforcing $CPT$ symmetry and analyticity. Here, we show that, besides allowing us to describe the big bang, these boundary conditions allow new gravitational instantons, enabling us to calculate the gravitational entropy of cosmologies which include radiation, dark energy and space curvature of either sign. We find the gravitational entropy of these universes, $S_g \sim S_\Lambda^{1/ 4} S_r$, where $S_\Lambda$ is the famous de Sitter entropy and $S_r$ is the total entropy in radiation. To the extent that $S_g$ exceeds $S_\Lambda$, the most probable universe is flat. By analysing the perturbations about our new instantons, we argue it is also homogeneous and isotropic on large scales.


title: Causality in Configurable Software Systems

abstract: Detecting and understanding reasons for defects and inadvertent behavior in software is challenging due to their increasing complexity. In configurable software systems, the combinatorics that arises from the multitude of features a user might select from adds a further layer of complexity. We introduce the notion of feature causality, which is based on counterfactual reasoning and inspired by the seminal definition of actual causality by Halpern and Pearl. Feature causality operates at the level of system configurations and is capable of identifying features and their interactions that are the reason for emerging functional and non-functional properties. We present various methods to explicate these reasons, in particular well-established notions of responsibility and blame that we extend to the feature-oriented setting. Establishing a close connection of feature causality to prime implicants, we provide algorithms to effectively compute feature causes and causal explications. By means of an evaluation on a wide range of configurable software systems, including community benchmarks and real-world systems, we demonstrate the feasibility of our approach: We illustrate how our notion of causality facilitates to identify root causes, estimate the effects of features, and detect feature interactions.


title: Annotating the Tweebank Corpus on Named Entity Recognition and Building NLP Models for Social Media Analysis

abstract: Social media data such as Twitter messages ("tweets") pose a particular challenge to NLP systems because of their short, noisy, and colloquial nature. Tasks such as Named Entity Recognition (NER) and syntactic parsing require highly domain-matched training data for good performance. To date, there is no complete training corpus for both NER and syntactic analysis (e.g., part of speech tagging, dependency parsing) of tweets. While there are some publicly available annotated NLP datasets of tweets, they are only designed for individual tasks. In this study, we aim to create Tweebank-NER, an English NER corpus based on Tweebank V2 (TB2), train state-of-the-art (SOTA) Tweet NLP models on TB2, and release an NLP pipeline called Twitter-Stanza. We annotate named entities in TB2 using Amazon Mechanical Turk and measure the quality of our annotations. We train the Stanza pipeline on TB2 and compare with alternative NLP frameworks (e.g., FLAIR, spaCy) and transformer-based models. The Stanza tokenizer and lemmatizer achieve SOTA performance on TB2, while the Stanza NER tagger, part-of-speech (POS) tagger, and dependency parser achieve competitive performance against non-transformer models. The transformer-based models establish a strong baseline in Tweebank-NER and achieve the new SOTA performance in POS tagging and dependency parsing on TB2. We release the dataset and make both the Stanza pipeline and BERTweet-based models available "off-the-shelf" for use in future Tweet NLP research. Our source code, data, and pre-trained models are available at: \url{https://github.com/social-machines/TweebankNLP}.


title: Galactic Mass-to-Light Ratios With Superfluid Dark Matter

abstract: We make rotation curve fits to test the superfluid dark matter model. In addition to verifying that the resulting fits match the rotation curve data reasonably well, we aim to evaluate how satisfactory they are with respect to two criteria, namely, how reasonable the resulting stellar mass-to-light ratios are and whether the fits end up in the regime of superfluid dark matter where the model resembles modified Newtonian dynamics (MOND). We fitted the superfluid dark matter model to the rotation curves of 169 galaxies in the SPARC sample. We found that the mass-to-light ratios obtained with superfluid dark matter are generally acceptable in terms of stellar populations. However, the best-fit mass-to-light ratios have an unnatural dependence on the size of the galaxy in that giant galaxies have systematically lower mass-to-light ratios than dwarf galaxies. A second finding is that the superfluid often fits the rotation curves best in the regime where the superfluid's force cannot resemble that of MOND without adjusting a boundary condition separately for each galaxy. In that case, we can no longer expect superfluid dark matter to reproduce the phenomenologically observed scaling relations that make MOND appealing. If, on the other hand, we consider only solutions whose force approximates MOND well, then the total mass of the superfluid is in tension with gravitational lensing data. We conclude that even the best fits with superfluid dark matter are still unsatisfactory for two reasons. First, the resulting stellar mass-to-light ratios show an unnatural trend with galaxy size. Second, the fits do not end up in the regime that automatically resembles MOND, and if we force the fits to do so, the total dark matter mass is in tension with strong lensing data.


title: Biophysical and Biochemical mechanisms underlying Collective Cell Migration in Cancer Metastasis

abstract: Multicellular collective migration is a ubiquitous strategy of cells to translocate spatially in diverse tissue environments to accomplish a wide variety of biological phenomena, viz. embryonic development, wound healing, and tumor progression. Diverse cellular functions and behaviors, for instance, cell protrusions, active contractions, cell-cell adhesion, biochemical signaling, remodeling of tissue micro-environment, etc., play their own role concomitantly to have a single concerted consequence of multicellular migration. Thus unveiling the driving principles, both biochemical and biophysical, of the inherently complex process of collective cell migration is an insurmountable task. Mathematical and computational models, in tandem with experimental data, help in shedding some light on it. Here we review different factors influencing Collective Cell Migration and then focus on different mathematical and computational models - discrete, hybrid, and continuum - which helps in revealing different aspects of multicellular migration. Finally, we discuss the applications of these modeling frameworks specific to cancer


title: TranAD: Deep Transformer Networks for Anomaly Detection in Multivariate Time Series Data

abstract: Efficient anomaly detection and diagnosis in multivariate time-series data is of great importance for modern industrial applications. However, building a system that is able to quickly and accurately pinpoint anomalous observations is a challenging problem. This is due to the lack of anomaly labels, high data volatility and the demands of ultra-low inference times in modern applications. Despite the recent developments of deep learning approaches for anomaly detection, only a few of them can address all of these challenges. In this paper, we propose TranAD, a deep transformer network based anomaly detection and diagnosis model which uses attention-based sequence encoders to swiftly perform inference with the knowledge of the broader temporal trends in the data. TranAD uses focus score-based self-conditioning to enable robust multi-modal feature extraction and adversarial training to gain stability. Additionally, model-agnostic meta learning (MAML) allows us to train the model using limited data. Extensive empirical studies on six publicly available datasets demonstrate that TranAD can outperform state-of-the-art baseline methods in detection and diagnosis performance with data and time-efficient training. Specifically, TranAD increases F1 scores by up to 17%, reducing training times by up to 99% compared to the baselines.


title: A zero-inflated endemic-epidemic model with an application to measles time series in Germany

abstract: Count data with excessive zeros are often encountered when modelling infectious disease occurrence. The degree of zero inflation can vary over time due to non-epidemic periods as well as by age group or region. The existing endemic-epidemic modelling framework (aka HHH) lacks a proper treatment for surveillance data with excessive zeros as it is limited to Poisson and negative binomial distributions. In this paper, we propose a multivariate zero-inflated endemic-epidemic model with random effects to extend HHH. Parameters of the new zero-inflation and the HHH part of the model can be estimated jointly and efficiently via (penalized) maximum likelihood inference using analytical derivatives. A simulation study confirms proper convergence and coverage probabilities of confidence intervals. Applying the model to measles counts in the 16 German states, 2005--2018, shows that the added zero-inflation improves probabilistic forecasts.


title: Conservative Distributional Reinforcement Learning with Safety Constraints

abstract: Safety exploration can be regarded as a constrained Markov decision problem where the expected long-term cost is constrained. Previous off-policy algorithms convert the constrained optimization problem into the corresponding unconstrained dual problem by introducing the Lagrangian relaxation technique. However, the cost function of the above algorithms provides inaccurate estimations and causes the instability of the Lagrange multiplier learning. In this paper, we present a novel off-policy reinforcement learning algorithm called Conservative Distributional Maximum a Posteriori Policy Optimization (CDMPO). At first, to accurately judge whether the current situation satisfies the constraints, CDMPO adapts distributional reinforcement learning method to estimate the Q-function and C-function. Then, CDMPO uses a conservative value function loss to reduce the number of violations of constraints during the exploration process. In addition, we utilize Weighted Average Proportional Integral Derivative (WAPID) to update the Lagrange multiplier stably. Empirical results show that the proposed method has fewer violations of constraints in the early exploration process. The final test results also illustrate that our method has better risk control.


title: Polar Coded Merkle Tree: Improved Detection of Data Availability Attacks in Blockchain Systems

abstract: Light nodes in blockchain systems are known to be vulnerable to data availability (DA) attacks where they accept an invalid block with unavailable portions. Previous works have used LDPC and 2-D Reed Solomon (2D-RS) codes with Merkle Trees to mitigate DA attacks. While these codes have demonstrated improved performance across a variety of metrics such as DA detection probability, they are difficult to apply to blockchains with large blocks due to generally intractable code guarantees for large codelengths (LDPC), large decoding complexity (2D-RS), or large coding fraud proof sizes (2D-RS). We address these issues by proposing the novel Polar Coded Merkle Tree (PCMT) which is a Merkle Tree built from the encoding graphs of polar codes and a specialized polar code construction called Sampling-Efficient Freezing (SEF). We demonstrate that the PCMT with SEF polar codes performs well in detecting DA attacks for large block sizes.


title: Extending the Vocabulary of Fictional Languages using Neural Networks

abstract: Fictional languages have become increasingly popular over the recent years appearing in novels, movies, TV shows, comics, and video games. While some of these fictional languages have a complete vocabulary, most do not. We propose a deep learning solution to the problem. Using style transfer and machine translation tools, we generate new words for a given target fictional language, while maintaining the style of its creator, hence extending this language vocabulary.


title: Sparsification of Decomposable Submodular Functions

abstract: Submodular functions are at the core of many machine learning and data mining tasks. The underlying submodular functions for many of these tasks are decomposable, i.e., they are sum of several simple submodular functions. In many data intensive applications, however, the number of underlying submodular functions in the original function is so large that we need prohibitively large amount of time to process it and/or it does not even fit in the main memory. To overcome this issue, we introduce the notion of sparsification for decomposable submodular functions whose objective is to obtain an accurate approximation of the original function that is a (weighted) sum of only a few submodular functions. Our main result is a polynomial-time randomized sparsification algorithm such that the expected number of functions used in the output is independent of the number of underlying submodular functions in the original function. We also study the effectiveness of our algorithm under various constraints such as matroid and cardinality constraints. We complement our theoretical analysis with an empirical study of the performance of our algorithm.


title: The Center-to-Limb Variation of Non-Thermal Velocities using IRIS Si IV

abstract: We study the non-thermal velocities in the quiet-sun using various high spatial, temporal, and spectral resolution observations from the Interface Region Imaging Spectrograph (IRIS). We focus our analysis on the transition region using the optically thin line (Si IV 1393.7 \AA), and select line profiles that are nearly Gaussian. We find evidence of a centre-to-limb variation using different observations having different exposure times, ranging from 5-30 s. The distribution of non-thermal velocities close to the limb are observed to peak around 20 km s$^{-1}$ while the disc observations show a peak around 15 km s$^{-1}$. The distributions are also different. The overall variation in the non-thermal velocities are correlated with the intensity of the line, as found previously. The on-disc velocities are smaller than most previous observations. In general, we find that the non-thermal velocities are independent of the selected exposure times. The Si IV lines didn't seem to exhibit any significant opacity effects. We conclude that these Doppler motions are mostly transverse to the radial direction. The possibility of swaying/torsional motions leading to such variations are validated from these IRIS observations.


title: Accelerating Bayesian inference of dependency between complex biological traits

abstract: Inferring dependencies between complex biological traits while accounting for evolutionary relationships between specimens is of great scientific interest yet remains infeasible when trait and specimen counts grow large. The state-of-the-art approach uses a phylogenetic multivariate probit model to accommodate binary and continuous traits via a latent variable framework, and utilizes an efficient bouncy particle sampler (BPS) to tackle the computational bottleneck -- integrating many latent variables from a high-dimensional truncated normal distribution. This approach breaks down as the number of specimens grows and fails to reliably characterize conditional dependencies between traits. Here, we propose an inference pipeline for phylogenetic probit models that greatly outperforms BPS. The novelty lies in 1) a combination of the recent Zigzag Hamiltonian Monte Carlo (Zigzag-HMC) with linear-time gradient evaluations and 2) a joint sampling scheme for highly correlated latent variables and correlation matrix elements. In an application exploring HIV-1 evolution from 535 viruses, the inference requires joint sampling from an 11,235-dimensional truncated normal and a 24-dimensional covariance matrix. Our method yields a 5-fold speedup compared to BPS and makes it possible to learn partial correlations between candidate viral mutations and virulence. Computational speedup now enables us to tackle even larger problems: we study the evolution of influenza H1N1 glycosylations on around 900 viruses. For broader applicability, we extend the phylogenetic probit model to incorporate categorical traits, and demonstrate its use to study Aquilegia flower and pollinator co-evolution.


title: The Pseudo-Lindley Alpha Power transformed distribution, mathematical characterizations and asymptotic properties

abstract: We introduce a new generalization of the Pseudo-Lindley distribution by applying alpha power transformation. The obtained distribution is referred as the Pseudo-Lindley alpha power transformed distribution (\textit{PL-APT}). Some tractable mathematical properties of the \textit{PL-APT} distribution as reliability, hazard rate, order statistics and entropies are provided. The maximum likelihood method is used to obtain the parameters' estimation of the \textit{PL-APT} distribution. The asymptotic properties of the proposed distribution are discussed. Also, a simulation study is performed to compare the modeling capability and flexibility of \textit{PL-APT} with Lindley and Pseudo-Lindley distributions. The \textit{PL-APT} provides a good fit as the Lindley and the Pseudo-Lindley distribution. The extremal domain of attraction of \textit{PL-APT} is found and its quantile and extremal quantile functions studied. Finally, the extremal value index is estimated by the double-indexed Hill's estimator (Ngom and Lo, 2016) and related asymptotic statistical tests are provided and characterized.


title: Landscape of heavy baryons from the perspective of the chiral quark-soliton model

abstract: We employ the chiral quark-soliton model and the heavy quark symmetry to describe spectra of charm and beauty baryons. Heavy baryons can be classified according to the SU(3) representations of the light sector. We argue that recently discovered $\Xi_b$ states can be interpreted as negative parity excited anti-triplets or sextets, and the $\Sigma_b$ states as negative parity sextets. Consequences of such assignments for the decay patterns are discussed and also predictions of masses of the yet unmeasured sextet members are given.


title: Daily-Resolved Lightning Climatology of the Eastern Alpine Region at the Kilometer Scale

abstract: Lightning flashes are rare albeit hazardous events. Despite this scarcity, generalized additive models (GAMs) succeed in producing a climatology of lightning occurrence for the eastern Alps and surrounding lowlands at an unprecedented resolution of 1\,km$^2$ for each day of April through September with data from the ALDIS lightning location system. The GAM adds the effects of seasonality, jaggedness of the terrain, and seasonally varying effects of elevation and region, thus combining information from analysis cells sharing similar characteristics. The probability of a cloud-to-ground discharge over 1\,km$^2$ on a given day is typically less than 1\,\% with a rapid increase in spring, followed by a plateau and a gentler tapering-off in fall. Probabilities early are lower at high elevations but increase once their snow cover is gone. Regional patterns of lightning also vary with season with an overall southward shift later in the year but more complex details. Grid cells with jagged topography have a higher probability of lightning.


title: A detailed look at the Calabi-Eguchi-Hanson spaces

abstract: This article takes a detailed look at the Ricci-flat metrics introduced by Eguchi-Hanson and Calabi on the canonical line bundle of complex projective space. We give a description of these spaces as resolutions of certain orbifold singularities. We then compute the curvature explicitly and show that all compact, minimal submanifolds are contained in the zero section. This extends a result by Tsai and Wang.


title: Convergence of Policy Gradient for Entropy Regularized MDPs with Neural Network Approximation in the Mean-Field Regime

abstract: We study the global convergence of policy gradient for infinite-horizon, continuous state and action space, and entropy-regularized Markov decision processes (MDPs). We consider a softmax policy with (one-hidden layer) neural network approximation in a mean-field regime. Additional entropic regularization in the associated mean-field probability measure is added, and the corresponding gradient flow is studied in the 2-Wasserstein metric. We show that the objective function is increasing along the gradient flow. Further, we prove that if the regularization in terms of the mean-field measure is sufficient, the gradient flow converges exponentially fast to the unique stationary solution, which is the unique maximizer of the regularized MDP objective. Lastly, we study the sensitivity of the value function along the gradient flow with respect to regularization parameters and the initial condition. Our results rely on the careful analysis of the non-linear Fokker-Planck-Kolmogorov equation and extend the pioneering work of Mei et al. 2020 and Agarwal et al. 2020, which quantify the global convergence rate of policy gradient for entropy-regularized MDPs in the tabular setting.


title: Strain-induced orbital energy shift in antiferromagnetic RuO2 revealed by resonant elastic x-ray scattering

abstract: In its ground state, RuO2 was long thought to be an ordinary metallic paramagnet. Recent neutron and x-ray diffraction revealed that bulk RuO2 is an antiferromagnet (AFM) with TN above 300 K. Furthermore, epitaxial strain induces novel superconductivity in thin films of RuO2 below 2 K. Here, we present a resonant elastic x-ray scattering (REXS) study at the Ru L2 edge of the strained RuO2 films exhibiting the strain-induced superconductivity. We observe an azimuthal modulation of the 100 Bragg peak consistent with canted AFM found in bulk. Most notably, in the strained films displaying novel superconductivity, we observe a ~1 eV shift of the Ru eg orbitals to a higher energy. The energy shift is smaller in thicker, relaxed films and films with a different strain direction. Our results provide further evidence of the utility of epitaxial strain as a tuning parameter in complex oxides.

categories: cond-mat.str-el cond-mat.mtrl-sci cond-mat.other cond-mat.supr-con


title: Angular-momentum projection in coupled-cluster theory: structure of $^{34}$Mg

abstract: Single-reference coupled-cluster theory is an accurate and affordable computational method for the nuclear many-body problem. For open-shell nuclei, the reference state typically breaks rotational invariance and angular momentum must be restored as a good quantum number. We perform angular-momentum projection after variation and employ the disentangled coupled-cluster formalism and a Hermitian approach. We compare our results with benchmarks for $^8$Be and $^{20}$Ne using a two-nucleon interaction from chiral effective field theory and for $pf$-shell nuclei within the traditional shell model. We compute the rotational band in the exotic nucleus $^{34}$Mg and find agreement with data.


title: Knot Floer homology and surgery on equivariant knots

abstract: Given an equivariant knot $K$ of order $2$, we study the induced action of the symmetry on the knot Floer homology. We relate this action with the induced action of the symmetry on the Heegaard Floer homology of large surgeries on $K$. This surgery formula can be thought of as an equivariant analog of the involutive large surgery formula proved by Hendricks and Manolescu. As a consequence, we obtain that for certain double branched covers of $S^{3}$ and corks, the induced action of the involution on Heegaard Floer homology can be identified with an action on the knot Floer homology. As an application, we calculate equivariant correction terms which are invariants of a generalized version of the spin rational homology cobordism group, and define two knot concordance invariants. We also compute the action of the symmetry on the knot Floer complex of $K$ for several equivariant knots.


title: 3-dimensional mixed BF theory and Hitchin's integrable system

abstract: The affine Gaudin model, associated with an untwisted affine Kac-Moody algebra, is known to arise from a certain gauge fixing of 4-dimensional mixed topological-holomorphic Chern-Simons theory in the Hamiltonian framework. We show that the finite Gaudin model, associated with a finite-dimensional semisimple Lie algebra, or more generally the tamely ramified Hitchin system on an arbitrary Riemann surface, can likewise be obtained from a similar gauge fixing of 3-dimensional mixed BF theory in the Hamiltonian framework.


title: Measurement of the inclusive and differential $\mathrm{t\bar{t}}\gamma$ cross sections in the dilepton channel and effective field theory interpretation in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: The production cross section of a top quark pair in association with a photon is measured in proton-proton collisions in the decay channel with two oppositely charged leptons (e$^\pm\mu^\mp$, e$^+$e$^-$, or $\mu^+\mu^-$). The measurement is performed using 138 fb$^{-1}$ of proton-proton collision data recorded by the CMS experiment at $\sqrt{s} =$ 13 TeV during the 2016-2018 data-taking period of the CERN LHC. A fiducial phase space is defined such that photons radiated by initial-state particles, top quarks, or any of their decay products are included. An inclusive cross section of 175.2 $\pm$ 2.5 (stat) $\pm$ 6.3 (syst) fb is measured in a signal region with at least one jet coming from the hadronization of a bottom quark and exactly one photon with transverse momentum above 20 GeV. Differential cross sections are measured as functions of several kinematic observables of the photon, leptons, and jets, and compared to standard model predictions. The measurements are also interpreted in the standard model effective field theory framework, and limits are found on the relevant Wilson coefficients from these results alone and in combination with a previous CMS measurement of the $\mathrm{t\bar{t}}\gamma$ production process using the lepton+jets final state.


title: Efficient emulation of relativistic heavy ion collisions with transfer learning

abstract: Measurements from the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC) can be used to study the properties of quark-gluon plasma. Systematic constraints on these properties must combine measurements from different collision systems and methodically account for experimental and theoretical uncertainties. Such studies require a vast number of costly numerical simulations. While computationally inexpensive surrogate models ("emulators") can be used to efficiently approximate the predictions of heavy ion simulations across a broad range of model parameters, training a reliable emulator remains a computationally expensive task. We use transfer learning to map the parameter dependencies of one model emulator onto another, leveraging similarities between different simulations of heavy ion collisions. By limiting the need for large numbers of simulations to only one of the emulators, this technique reduces the numerical cost of comprehensive uncertainty quantification when studying multiple collision systems and exploring different models.


title: Large Hybrid Time-Varying Parameter VARs

abstract: Time-varying parameter VARs with stochastic volatility are routinely used for structural analysis and forecasting in settings involving a few endogenous variables. Applying these models to high-dimensional datasets has proved to be challenging due to intensive computations and over-parameterization concerns. We develop an efficient Bayesian sparsification method for a class of models we call hybrid TVP-VARs--VARs with time-varying parameters in some equations but constant coefficients in others. Specifically, for each equation, the new method automatically decides whether the VAR coefficients and contemporaneous relations among variables are constant or time-varying. Using US datasets of various dimensions, we find evidence that the parameters in some, but not all, equations are time varying. The large hybrid TVP-VAR also forecasts better than many standard benchmarks.


title: Detection of Chemo-Kinematical Structures in Leo I

abstract: Context: A variety of formation models for dwarf spheroidal (dSph) galaxies have been proposed in the literature, but generally they have not been quantitatively compared with observations. Aims: We search for chemodynamical patterns in our observational data set and compare the results with mock galaxies consisting of pure random motions, and simulated dwarfs formed via the dissolving star cluster and tidal stirring models. Methods: We made use of a new spectroscopic data set for the Milky Way dSph Leo I, combining 288 stars observed with Magellan/IMACS and existing Keck/DEIMOS data, to provide velocity and metallicity measurements for 953 Leo I member stars. We used a specially developed algorithm called Beacon to detect chemo-kinematical patterns in the observed and simulated data. Results: After analysing the Leo I data, we report the detection of 14 candidate streams of stars that may have originated in disrupted star clusters. The angular momentum vectors of these streams are randomly oriented, consistent with the lack of rotation in Leo I. These results are consistent with the predictions of the dissolving cluster model. In contrast, we find fewer candidate stream signals in mock data sets that lack coherent motions ~99% of the time. The chemodynamical analysis of the tidal stirring simulation produces streams that share a common orientation of their angular momenta, which is inconsistent with the Leo I data. Conclusions: Even though it is very difficult to distinguish which of the detected streams are real and which are only noise, we can be certain that there are more streams detected in the observational data of Leo I than expected in pure random data.


title: Proportional Ranking in Primary Elections: A Case Study

abstract: Many democratic political parties hold primary elections, which nicely reflects their democratic nature and promote, among other things, the democratic value of inclusiveness. However, the methods currently used for holding such primary elections may not be the most suitable, especially if some form of proportional ranking is desired. In this paper, we compare different algorithmic methods for holding primaries (i.e., different aggregation methods for voters' ballots), by evaluating the degree of proportional ranking that is achieved by each of them using real-world data. In particular, we compare six different algorithms by analyzing real-world data from a recent primary election conducted by the Israeli Democratit party. Technically, we analyze unique voter data and evaluate the proportionality achieved by means of cluster analysis, aiming at pinpointing the representation that is granted to different voter groups under each of the algorithmic methods considered. Our finding suggest that, contrary to the most-prominent primaries algorithm used (i.e., Approval), other methods such as Sequential Proportional Approval or Phragmen can bring about better proportional ranking and thus may be better suited for primary elections in practice.


title: Bregman Deviations of Generic Exponential Families

abstract: We revisit the method of mixture technique, also known as the Laplace method, to study the concentration phenomenon in generic exponential families. Combining the properties of Bregman divergence associated with log-partition function of the family with the method of mixtures for super-martingales, we establish a generic bound controlling the Bregman divergence between the parameter of the family and a finite sample estimate of the parameter. Our bound is time-uniform and makes appear a quantity extending the classical information gain to exponential families, which we call the Bregman information gain. For the practitioner, we instantiate this novel bound to several classical families, e.g., Gaussian, Bernoulli, Exponential, Weibull, Pareto, Poisson and Chi-square yielding explicit forms of the confidence sets and the Bregman information gain. We further numerically compare the resulting confidence bounds to state-of-the-art alternatives for time-uniform concentration and show that this novel method yields competitive results. Finally, we highlight the benefit of our concentration bounds on some illustrative applications.


title: Visualizing fluid flows via regularized optimal mass transport with applications to neuroscience

abstract: Regularized optimal mass transport (rOMT) problem adds a diffusion term to the continuity equation in the original dynamic formulation of the optimal mass transport (OMT) problem proposed by Benamou and Brenier. We show that the rOMT model serves as a powerful tool in computational fluid dynamics (CFD) for visualizing fluid flows in the glymphatic system. In the present work, we describe how to modify the previous numerical method for efficient implementation, resulting in a significant reduction in computational runtime. Numerical results applied to synthetic and real-data are provided.


title: Enabling Deep Reinforcement Learning on Energy Constrained Devices at the Edge of the Network

abstract: Deep Reinforcement Learning (DRL) solutions are becoming pervasive at the edge of the network as they enable autonomous decision-making in a dynamic environment. However, to be able to adapt to the ever-changing environment, the DRL solution implemented on an embedded device has to continue to occasionally take exploratory actions even after initial convergence. In other words, the device has to occasionally take random actions and update the value function, i.e., re-train the Artificial Neural Network (ANN), to ensure its performance remains optimal. Unfortunately, embedded devices often lack processing power and energy required to train the ANN. The energy aspect is particularly challenging when the edge device is powered only by a means of Energy Harvesting (EH). To overcome this problem, we propose a two-part algorithm in which the DRL process is trained at the sink. Then the weights of the fully trained underlying ANN are periodically transferred to the EH-powered embedded device taking actions. Using an EH-powered sensor, real-world measurements dataset, and optimizing for Age of Information (AoI) metric, we demonstrate that such a DRL solution can operate without any degradation in the performance, with only a few ANN updates per day.


title: OSSID: Online Self-Supervised Instance Detection by (and for) Pose Estimation

abstract: Real-time object pose estimation is necessary for many robot manipulation algorithms. However, state-of-the-art methods for object pose estimation are trained for a specific set of objects; these methods thus need to be retrained to estimate the pose of each new object, often requiring tens of GPU-days of training for optimal performance. In this paper, we propose the OSSID framework, leveraging a slow zero-shot pose estimator to self-supervise the training of a fast detection algorithm. This fast detector can then be used to filter the input to the pose estimator, drastically improving its inference speed. We show that this self-supervised training exceeds the performance of existing zero-shot detection methods on two widely used object pose estimation and detection datasets, without requiring any human annotations. Further, we show that the resulting method for pose estimation has a significantly faster inference speed, due to the ability to filter out large parts of the image. Thus, our method for self-supervised online learning of a detector (trained using pseudo-labels from a slow pose estimator) leads to accurate pose estimation at real-time speeds, without requiring human annotations. Supplementary materials and code can be found at https://georgegu1997.github.io/OSSID/


title: Subfactors from regular graphs induced by association schemes

abstract: We clarify the relations between the mathematical structures that enable fashioning quantum walks on regular graphs and their realizations in anyonic systems. Our protagonist is association schemes that may be synthesized from type-II matrices which have a canonical construction of subfactors. This way we set up quantum walks on growing distance-regular graphs induced by association schemes via interacting Fock spaces and relate them to anyon systems described by subfactors. We discuss in detail a large family of graphs that may be treated within this approach. Classification of association schemes and realizable anyon systems are complex combinatorial problems and we tackle a part of it with a quantum walk application based approach.


title: Datasheet for the Pile

abstract: This datasheet describes the Pile, a 825 GiB dataset of human-authored text compiled by EleutherAI for use in large-scale language modeling. The Pile is comprised of 22 different text sources, ranging from original scrapes done for this project, to text data made available by the data owners, to third-party scrapes available online.


title: Model-driven Cluster Resource Management for AI Workloads in Edge Clouds

abstract: Since emerging edge applications such as Internet of Things (IoT) analytics and augmented reality have tight latency constraints, hardware AI accelerators have been recently proposed to speed up deep neural network (DNN) inference run by these applications. Resource-constrained edge servers and accelerators tend to be multiplexed across multiple IoT applications, introducing the potential for performance interference between latency-sensitive workloads. In this paper, we design analytic models to capture the performance of DNN inference workloads on shared edge accelerators, such as GPU and edgeTPU, under different multiplexing and concurrency behaviors. After validating our models using extensive experiments, we use them to design various cluster resource management algorithms to intelligently manage multiple applications on edge accelerators while respecting their latency constraints. We implement a prototype of our system in Kubernetes and show that our system can host 2.3X more DNN applications in heterogeneous multi-tenant edge clusters with no latency violations when compared to traditional knapsack hosting algorithms.


title: Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

abstract: For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.


title: Momentum-resolved time evolution with matrix product states

abstract: We introduce a method based on matrix product states (MPS) for computing spectral functions of (quasi) one-dimensional spin chains, working directly in momentum space in the thermodynamic limit. We simulate the time evolution after applying a momentum operator to an MPS ground state by working with the momentum superposition of a window MPS. We show explicitly for the spin-1 Heisenberg chain that the growth of entanglement is smaller in momentum space, even inside a two-particle continuum, such that we can attain very accurate spectral functions with relatively small bond dimension. We apply our method to compute spectral lineshapes of the gapless XXZ chain and the square-lattice J1-J2 Heisenberg model on a six-leg cylinder.


title: Deuterium chemistry and D/H ratios in Class 0/I proto-brown dwarfs

abstract: We have conducted the first extensive observational survey of several deuterated species in 16 Class 0/I proto-brown dwarfs (proto-BDs) and 4 Class Flat/Class II brown dwarfs. Observations were obtained with the IRAM 30m telescope in the DCO$^{+}$ (3-2), DCN (3-2), DNC (3-2), and N$_{2}$D$^{+}$ (3-2) lines. The DCO$^{+}$/H$^{13}$CO$^{+}$, DCN/H$^{13}$CN, and DNC/HN$^{13}$C ratios are comparatively higher and show a narrower range than the DCO$^{+}$/HCO$^{+}$, DCN/HCN, and DNC/HNC ratios, respectively. The mean D/H ratios for the proto-BDs derived from these molecules are in the range of $\sim$0.02-3. Both low-temperature gas-phase ion-molecule deuteron transfer and grain surface reactions are required to explain the enhanced deuterium fractionation. The very dense and cold ($n_{H_{2}}$ $\geq$10$^{6}$ cm$^{-3}$, T $\leq$10 K) interior of the proto-BDs provide the suitable conditions for efficient deuterium fractionation in these cores. There is no correlation between the D/H ratios and the CO depletion factor, with the exception of the DCN/HCN ratios that show a strong anti-correlation possibly due to the difference in the peak emitting regions of the DCN and HCN molecules. Over a wide range in the bolometric luminosities spanning $\sim$0.002--40 L$_{\odot}$, we find a trend of higher DCO$^{+}$/HCO$^{+}$ (r = -0.7) and DCN/HCN (r = -0.6) ratios, nearly constant DNC/HNC (r = -0.4) and DNC/HN$^{13}$C (r = -0.1) ratios, lower N$_{2}$D$^{+}$/N$_{2}$H$^{+}$ ratios (r = 0.6) in the proto-BDs compared to protostars. Only one Class II brown dwarf shows emission in the DCO$^{+}$ (3-2) line. No correlation is seen between the D/H ratios and the evolutionary stage.


title: Towards Federated Clustering: A Federated Fuzzy $c$-Means Algorithm (FFCM)

abstract: Federated Learning (FL) is a setting where multiple parties with distributed data collaborate in training a joint Machine Learning (ML) model while keeping all data local at the parties. Federated clustering is an area of research within FL that is concerned with grouping together data that is globally similar while keeping all data local. We describe how this area of research can be of interest in itself, or how it helps addressing issues like non-independently-identically-distributed (i.i.d.) data in supervised FL frameworks. The focus of this work, however, is an extension of the federated fuzzy $c$-means algorithm to the FL setting (FFCM) as a contribution towards federated clustering. We propose two methods to calculate global cluster centers and evaluate their behaviour through challenging numerical experiments. We observe that one of the methods is able to identify good global clusters even in challenging scenarios, but also acknowledge that many challenges remain open.


title: A Privacy-Preserving Unsupervised Domain Adaptation Framework for Clinical Text Analysis

abstract: Unsupervised domain adaptation (UDA) generally aligns the unlabeled target domain data to the distribution of the source domain to mitigate the distribution shift problem. The standard UDA requires sharing the source data with the target, having potential data privacy leaking risks. To protect the source data's privacy, we first propose to share the source feature distribution instead of the source data. However, sharing only the source feature distribution may still suffer from the membership inference attack who can infer an individual's membership by the black-box access to the source model. To resolve this privacy issue, we further study the under-explored problem of privacy-preserving domain adaptation and propose a method with a novel differential privacy training strategy to protect the source data privacy. We model the source feature distribution by Gaussian Mixture Models (GMMs) under the differential privacy setting and send it to the target client for adaptation. The target client resamples differentially private source features from GMMs and adapts on target data with several state-of-art UDA backbones. With our proposed method, the source data provider could avoid leaking source data privacy during domain adaptation as well as reserve the utility. To evaluate our proposed method's utility and privacy loss, we apply our model on a medical report disease label classification task using two noisy challenging clinical text datasets. The results show that our proposed method can preserve source data's privacy with a minor performance influence on the text classification task.


title: On the Molecular Picture and Interfacial Temperature Discontinuity During Evaporation and Condensation

abstract: Although it has been shown experimentally that a temperature discontinuity exists at the liquid-vapor interface during evaporation and condensation, quantitatively modeling this temperature jump has been difficult. The classical Schrage equation does not give enough information to determine the interfacial temperature jump. Starting from the Boltzmann transport equation, this paper establishes three interfacial boundary conditions to connect the temperature, density, and pressure jumps at the liquid-vapor interface to the interfacial mass and heat fluxes: one for the mass flux (the Schrage equation), one for the heat flux, and the third for the density discontinuities. These expressions can be readily coupled to heat and mass transport equations in the continuum of the liquid and the vapor phases, enabling one to determine the values of the interfacial temperature, density, and pressure jumps. Comparison with past experiments is favorable. A thermomolecular emission model, mimicking thermionic emission of electrons, is also presented to gain more molecular-level insights on the thermal evaporation processes.


title: Asymptotic properties of Bayesian inference in linear regression with a structural break

abstract: This paper studies large sample properties of a Bayesian approach to inference about slope parameters $\gamma$ in linear regression models with a structural break. In contrast to the conventional approach to inference about $\gamma$ that does not take into account the uncertainty of the unknown break location $\tau$, the Bayesian approach that we consider incorporates such uncertainty. Our main theoretical contribution is a Bernstein-von Mises type theorem (Bayesian asymptotic normality) for $\gamma$ under a wide class of priors, which essentially indicates an asymptotic equivalence between the conventional frequentist and Bayesian inference. Consequently, a frequentist researcher could look at credible intervals of $\gamma$ to check robustness with respect to the uncertainty of $\tau$. Simulation studies show that the conventional confidence intervals of $\gamma$ tend to undercover in finite samples whereas the credible intervals offer more reasonable coverages in general. As the sample size increases, the two methods coincide, as predicted from our theoretical conclusion. Using data from Paye and Timmermann (2006) on stock return prediction, we illustrate that the traditional confidence intervals on $\gamma$ might underrepresent the true sampling uncertainty.


title: Critical periods and Autism Spectrum Disorders, a role for sleep

abstract: Brain development relies on both experience and genetically defined programs. Time windows where certain brain circuits are particularly receptive to external stimuli, resulting in heightened plasticity, are referred to as critical periods. Sleep is thought to be essential for normal brain development. Importantly, studies have shown that sleep enhances critical period plasticity and promotes experience-dependent synaptic pruning in the developing mammalian brain. Therefore, normal plasticity during critical periods depends on proper sleep. Problems falling and staying asleep occur at a higher rate in Autism Spectrum Disorder (ASD) relative to typical development. In this review, we explore the potential link between sleep, critical period plasticity, and ASD. First, we review the importance of critical period plasticity in typical development and the role of sleep in this process. Next, we summarize the evidence linking ASD with deficits in synaptic plasticity in rodent models of high-confident ASD gene candidates. We then show that almost all the high-confidence rodent models of ASD that show sleep deficits also display plasticity deficits. Given how important sleep is for critical period plasticity, it is essential to understand the connections between synaptic plasticity, sleep, and brain development in ASD. However, studies investigating sleep or plasticity during critical periods in ASD mouse models are lacking. Therefore, we highlight an urgent need to consider developmental trajectory in studies of sleep and plasticity in neurodevelopmental disorders.


title: VaxEquity: A Data-Driven Risk Assessment and Optimization Framework for Equitable Vaccine Distribution

abstract: With the continuous rise of the COVID-19 cases worldwide, it is imperative to ensure that all those vulnerable countries lacking vaccine resources can receive sufficient support to contain the risks. COVAX is such an initiative operated by the WHO to supply vaccines to the most needed countries. One critical problem faced by the COVAX is how to distribute the limited amount of vaccines to these countries in the most efficient and equitable manner. This paper aims to address this challenge by first proposing a data-driven risk assessment and prediction model and then developing a decision-making framework to support the strategic vaccine distribution. The machine learning-based risk prediction model characterizes how the risk is influenced by the underlying essential factors, e.g., the vaccination level among the population in each COVAX country. This predictive model is then leveraged to design the optimal vaccine distribution strategy that simultaneously minimizes the resulting risks while maximizing the vaccination coverage in these countries targeted by COVAX. Finally, we corroborate the proposed framework using case studies with real-world data.


title: Transparent Single-Cell Set Classification with Kernel Mean Embeddings

abstract: Modern single-cell flow and mass cytometry technologies measure the expression of several proteins of the individual cells within a blood or tissue sample. Each profiled biological sample is thus represented by a set of hundreds of thousands of multidimensional cell feature vectors, which incurs a high computational cost to predict each biological sample's associated phenotype with machine learning models. Such a large set cardinality also limits the interpretability of machine learning models due to the difficulty in tracking how each individual cell influences the ultimate prediction. We propose using Kernel Mean Embedding to encode the cellular landscape of each profiled biological sample. Although our foremost goal is to make a more transparent model, we find that our method achieves comparable or better accuracies than the state-of-the-art gating-free methods through a simple linear classifier. As a result, our model contains few parameters but still performs similarly to deep learning models with millions of parameters. In contrast with deep learning approaches, the linearity and sub-selection step of our model makes it easy to interpret classification results. Analysis further shows that our method admits rich biological interpretability for linking cellular heterogeneity to clinical phenotype.


title: A Multi-factor Multi-level and Interaction based (M2I) Authentication Framework for Internet of Things (IoT) Applications

abstract: Existing authentication solutions proposed for Internet of Things (IoT) provide a single Level of Assurance (LoA) regardless of the sensitivity levels of the resources or interactions between IoT devices being protected. For effective (with adequate level of protection) and efficient (with as low overhead costs as possible) protections, it may be desirable to tailor the protection level in response to the sensitivity level of the resources, as a stronger protection level typically imposes a higher level of overheads costs. In this paper, we investigate how to facilitate multi-LoA authentication for IoT by proposing a multi-factor multi-level and interaction based (M2I) authentication framework. The framework implements LoA linked and interaction based authentication. Two interaction modes are investigated, P2P (Peer-to-Peer) and O2M (One-to-Many) via the design of two corresponding protocols. Evaluation results show that adopting the O2M interaction mode in authentication can cut communication cost significantly; compared with that of the Kerberos protocol, the O2M protocol reduces the communication cost by 42% ~ 45%. The protocols also introduce less computational cost. The P2P and O2M protocol, respectively, reduce the computational cost by 70% ~ 72% and 81% ~ 82% in comparison with that of Kerberos. Evaluation results also show that the two factor authentication option costs twice as much as that of the one-factor option.


title: Positron emission in elastic collisions of fully ionized high-Z heavy ions

abstract: Recent theoretical investigations of pair-production probabilities for given sets of nuclear trajectories predict the observability of spontaneous positron emission in elastic collisions of heavy ions. These calculations are extended to the investigation of ratios $P_{y}(\eta)/P_{y}(1)$ for partial probabilities of positron emission at small positron energies. The results suggest a possible observability of spontaneous positron emission in elastic collisions of bare heavy ions by analyzing the slope of $P_{y}(\eta)/P_{y}(1)$ for values of $R_{min}$ $~_{\sim}\!\!\!\!^<$ 0.75 $R_{cr}$. In addition calculations of cross sections for positron emission in U+Cm and Cm+Cm collisions of fully ionized heavy ions are compared with former results with partly ionized projectiles and show an increase by factors of 57 to 72.


title: FMM-LU: A fast direct solver for multiscale boundary integral equations in three dimensions

abstract: We present a fast direct solver for boundary integral equations on complex surfaces in three dimensions using an extension of the recently introduced recursive strong skeletonization scheme. For problems that are not highly oscillatory, our algorithm computes an ${LU}$-like hierarchical factorization of the dense system matrix, permitting application of the inverse in $\mathcal O(n)$ time, where $n$ is the number of unknowns on the surface. The factorization itself also scales linearly with the system size, albeit with a somewhat larger constant. The scheme is built on a level-restricted adaptive octree data structure, and therefore it is compatible with highly nonuniform discretizations. Furthermore, the scheme is coupled with high-order accurate locally-corrected Nystr\"om quadrature methods to integrate the singular and weakly-singular Green's functions used in the integral representations. Our method has immediate applications to a variety of problems in computational physics. We concentrate here on studying its performance in acoustic scattering (governed by the Helmholtz equation) at low to moderate frequencies, and provide rigorous justification for compression of submatrices via proxy surfaces.


title: Automated Counting and Statistical Analysis of Labeled Trees with Degree Restrictions

abstract: Arthur Cayley famously proved that there are n to the power n-2 labeled trees on n vertices. Here we go much further and show how to enumerate, fully automatically, labeled trees such that every vertex has a number of neighbors that belongs to a specified finite set, and also count trees where the number of neighbors is not allowed to be in a given finite set. We also give detailed statistical analysis, and show that in the sample space of labeled trees with n vertices, the random variable "number of vertices with d neighbors" is asymptotically normal, and for any different degrees, are jointly asymptotically normal, but of course, not independently so (except for the pair (1,3), i.e. the number of leaves and the number of degree-3 vertices, where there are asymptotically independent). We also give new proofs to Amram Meir and John Noon's expressions for the limiting expectation and variance for these, and derive an explicit expression for the covariance.


title: Hermite-based, One-step, Variational and Galerkin Time Integrators for Mechanical Systems

abstract: In this paper, we present two Hermite polynomial based approaches to derive one-step numerical integrators for mechanical systems. These methods are based on discretizing the configuration using Hermite polynomials which leads to numerical trajectories continuous in both configuration and velocity. First, we incorporate Hermite polynomials for time-discretization and derive one-step variational methods by discretizing the Lagrange-d'Alembert principle over a single time step. Second, we present the Galerkin approach to derive one-step numerical integrators by setting the weighted average of the residual of the equations of motion over a time step to zero. We consider three numerical examples to understand the numerical performance of the one-step variational and Galerkin methods. We first study a particle in a double-well potential and compare the variational approach results with the corresponding results for the Galerkin approach. We then study the Duffing oscillator to understand the numerical behavior in presence of dissipative forces. Finally, we apply the proposed methods to a nonlinear aeroelastic system with two degrees of freedom. Both variational and Galerkin one-step methods capture conservative and nonconservative dynamics accurately with excellent energy behavior. The one-step Galerkin methods exhibit better trajectory and energy performance than the one-step variational methods and the variational integrators.


title: Cutting Through the Noise to Infer Autonomous System Topology

abstract: The Border Gateway Protocol (BGP) is a distributed protocol that manages interdomain routing without requiring a centralized record of which autonomous systems (ASes) connect to which others. Many methods have been devised to infer the AS topology from publicly available BGP data, but none provide a general way to handle the fact that the data are notoriously incomplete and subject to error. This paper describes a method for reliably inferring AS-level connectivity in the presence of measurement error using Bayesian statistical inference acting on BGP routing tables from multiple vantage points. We employ a novel approach for counting AS adjacency observations in the AS-PATH attribute data from public route collectors, along with a Bayesian algorithm to generate a statistical estimate of the AS-level network. Our approach also gives us a way to evaluate the accuracy of existing reconstruction methods and to identify advantageous locations for new route collectors or vantage points.

categories: cs.NI cs.SI physics.data-an physics.soc-ph


title: On the minimax rate of the Gaussian sequence model under bounded convex constraints

abstract: We determine the exact minimax rate of a Gaussian sequence model under bounded convex constraints, purely in terms of the local geometry of the given constraint set $K$. Our main result shows that the minimax risk (up to constant factors) under the squared $\ell_2$ loss is given by $\epsilon^{*2} \wedge \operatorname{diam}(K)^2$ with \begin{align*} \epsilon^* = \sup \bigg\{\epsilon : \frac{\epsilon^2}{\sigma^2} \leq \log M^{\operatorname{loc}}(\epsilon)\bigg\}, \end{align*} where $\log M^{\operatorname{loc}}(\epsilon)$ denotes the local entropy of the set $K$, and $\sigma^2$ is the variance of the noise. We utilize our abstract result to re-derive known minimax rates for some special sets $K$ such as hyperrectangles, ellipses, and more generally quadratically convex orthosymmetric sets. Finally, we extend our results to the unbounded case with known $\sigma^2$ to show that the minimax rate in that case is $\epsilon^{*2}$.


title: Radium Ion Optical Clock

abstract: We report the first operation of a Ra$^{+}$ optical clock, a promising high-performance clock candidate. The clock uses a single trapped $^{226}$Ra$^{+}$ ion and operates on the $7s\ ^2S_{1/2}\rightarrow$ $6d\ ^2D_{5/2}$ electric quadrupole transition. By self-referencing three pairs of symmetric Zeeman transitions, we demonstrate a frequency instability of 1.1$\times10^{-13}$/$\sqrt{\tau}$, where $\tau$ is the averaging time in seconds. The total systematic uncertainty is evaluated to be ${\Delta \nu / \nu = 9 \times 10^{-16}}$. Using the clock, we realize the first measurement of the ratio of the $D_{5/2}$ state to the $S_{1/2}$ state Land\'{e} $g$-factors: $g_{D}/g_{S}$ = 0.5988053(11). A Ra$^{+}$ optical clock could improve limits on the time variation of the fine structure constant, $\dot \alpha / \alpha$, in an optical frequency comparison. The ion also has several features that make it a suitable system for a transportable optical clock.


title: Bubble Nucleation to All Orders

abstract: This paper extends classical results by Langer and Kramers and combines them with modern methods from high-temperature field theory. Assuming Langevin dynamics, the end-product is an all-orders description of bubble-nucleation at high temperatures. Specifically, it's shown that equilibrium and non-equilibrium effects factorize to all orders, and that the nucleation rate splits into a statistical and a dynamical prefactor. The derivation clarifies, and incorporates, higher-order corrections from zero-modes. The rate is also shown to be real to all orders in perturbation theory. The methods are applied to several models. As such, Feynman rules are given; the relevant power-counting is introduced; RG invariance is shown; the connection with the effective action is discussed, and an explicit construction of propagators in an inhomogeneous background is given. The formalism applies to both phase and Sphaleron transitions. While mainly focused on field theory, the methods are applicable to finite-dimensional systems. Finally, as this paper assumes an effective Langevin description, all results only hold within this framework.


title: Fragmentation in turbulence by small eddies

abstract: From air-sea gas exchange, oil pollution, to bioreactors, the ubiquitous fragmentation of bubbles/drops in turbulence has been modelled by relying on the classical Kolmogorov-Hinze paradigm since the 1950s. This framework hypothesizes that bubbles/drops are broken solely by eddies of the same size, even though turbulence is well known for its wide spectrum of scales. Here, by designing an experiment that can physically and cleanly disentangle eddies of various sizes, we report the experimental evidence to challenge this hypothesis and show bubbles are preferentially broken by the sub-bubble-scale eddies. Our work also highlights that fragmentation cannot be quantified solely by the stress criterion or the Weber number; The competition between different time scales is equally important. Instead of being elongated slowly and persistently by flows at their own scales, bubbles are fragmented in turbulence by small eddies via a burst of intense local deformation within a short time.


title: The Newton polytope and Lorentzian property of chromatic symmetric functions

abstract: Chromatic symmetric functions are well-studied symmetric functions in algebraic combinatorics that generalize the chromatic polynomial and are related to Hessenberg varieties and diagonal harmonics. Motivated by the Stanley--Stembridge conjecture, we show that the allowable coloring weights for indifference graphs of Dyck paths are the lattice points of a permutahedron $\mathcal{P}_\lambda$, and we give a formula for the dominant weight $\lambda$. Furthermore, we conjecture that such chromatic symmetric functions are Lorentzian, a property introduced by Br\"and\'en and Huh as a bridge between discrete convex analysis and concavity properties in combinatorics, and we prove this conjecture for abelian Dyck paths. We extend our results on the Newton polytope to incomparability graphs of (3+1)-free posets, and we give a number of conjectures and results stemming from our work, including results on the complexity of computing the coefficients and relations with the $\zeta$ map from diagonal harmonics.


title: FAUST III. Misaligned rotations of the envelope, outflow, and disks in the multiple protostellar system of VLA 1623$-$2417

abstract: We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed.


title: Local Lagrangian reduced-order modeling for Rayleigh-Taylor instability by solution manifold decomposition

abstract: Rayleigh-Taylor instability is a classical hydrodynamic instability of great interest in various disciplines of science and engineering, including astrophyics, atmospheric sciences and climate, geophysics, and fusion energy. Analytical methods cannot be applied to explain the long-time behavior of Rayleigh-Taylor instability, and therefore numerical simulation of the full problem is required. However, in order to capture the growth of amplitude of perturbations accurately, both the spatial and temporal discretization need to be extremely fine for traditional numerical methods, and the long-time simulation may become prohibitively expensive. In this paper, we propose efficient reduced order model techniques to accelerate the simulation of Rayleigh-Taylor instability in compressible gas dynamics. We introduce a general framework for decomposing the solution manifold to construct the temporal domain partition and temporally-local reduced order model construction with varying Atwood number. We propose two practical approaches in this framework, namely decomposition by using physical time and penetration distance respectively. Numerical results are presented to examine the performance of the proposed approaches.


title: Two applications of Grunsky coefficients in the theory of univalent functions

abstract: Let $\mathcal{S}$ denote the class of functions $f$ which are analytic and univalent in the unit disk ${\mathbb D}=\{z:|z|<1\}$ and normalized with $f(z)=z+\sum_{n=2}^{\infty} a_n z^n$. Using a method based on Grusky coefficients we study two problems over the class $\mathcal{S}$: estimate of the fourth logarithmic coefficient and upper bound of the coefficient difference $|a_5|-|a_4|$.


title: Effects of Charge Exchange on the Evaporative Wind of HD 209458b

abstract: The role of charge exchange in shaping exoplanet photoevaporation remains a topic of contention. Exchange of electrons between stellar wind protons from the exoplanet's host star and neutral hydrogen from the planet's wind has been proposed as a mechanism to create "energetic neutral atoms" (ENAs), which could explain the high absorption line velocities observed in systems where mass loss is occurring. In this paper we present results from 3D hydrodynamic simulations of the mass loss of a planet similar to HD 209458b. We self-consistently launch a planetary wind by calculating the ionization and heating resulting from incident high-energy radiation, inject a stellar wind into the simulation, and allow electron exchange between the stellar and planetary winds. We predict the potential production of ENAs by the wind-wind interaction analytically, then present the results of our simulations, which confirm the analytic limits. Within the limits of our hydrodynamic simulation, we find that charge exchange with the stellar wind properties examined here is unable to explain the absorption observed at high Doppler velocities.


title: Controllable Protein Design with Language Models

abstract: The 21st century is presenting humankind with unprecedented environmental and medical challenges. The ability to design novel proteins tailored for specific purposes could transform our ability to respond timely to these issues. Recent advances in the field of artificial intelligence are now setting the stage to make this goal achievable. Protein sequences are inherently similar to natural languages: Amino acids arrange in a multitude of combinations to form structures that carry function, the same way as letters form words and sentences that carry meaning. Therefore, it is not surprising that throughout the history of Natural Language Processing (NLP), many of its techniques have been applied to protein research problems. In the last few years, we have witnessed revolutionary breakthroughs in the field of NLP. The implementation of Transformer pre-trained models has enabled text generation with human-like capabilities, including texts with specific properties such as style or subject. Motivated by its considerable success in NLP tasks, we expect dedicated Transformers to dominate custom protein sequence generation in the near future. Finetuning pre-trained models on protein families will enable the extension of their repertoires with novel sequences that could be highly divergent but still potentially functional. The combination of control tags such as cellular compartment or function will further enable the controllable design of novel protein functions. Moreover, recent model interpretability methods will allow us to open the 'black box' and thus enhance our understanding of folding principles. While early initiatives show the enormous potential of generative language models to design functional sequences, the field is still in its infancy. We believe that protein language models are a promising and largely unexplored field and discuss their foreseeable impact on protein design.


title: Determining the Timescale over Which Stellar Feedback Drives Turbulence in the ISM: A Study of four Nearby Dwarf Irregular Galaxies

abstract: Stellar feedback is fundamental to the modeling of galaxy evolution as it drives turbulence and outflows in galaxies. Understanding the timescales involved are critical for constraining the impact of stellar feedback on the interstellar medium (ISM). We analyzed the resolved star formation histories along with the spatial distribution and kinematics of the atomic and ionized gas of four nearby star-forming dwarf galaxies (NGC 4068, NGC 4163, NGC 6789, UGC 9128) to determine the timescales over which stellar feedback drives turbulence. The four galaxies are within 5 Mpc and have a range of properties including current star formation rates of 0.0005 to 0.01 M$_{\odot}$ yr$^{-1}$, log(M$_*$/M$_{\odot}$) between 7.2 and 8.2, and log(M$_{HI}$/M$_\odot$) between 7.2 and 8.3. Their Color-Magnitude Diagram (CMD) derived star formation histories over the past 500 Myrs were compared to their atomic and ionized gas velocity dispersion and HI energy surface densities as indicators of turbulence. The Spearman's rank correlation coefficient was used to identify any correlations between their current turbulence and their past star formation activity on local scales ($\sim$400 pc). The strongest correlation found was between the HI turbulence measures and the star formation rate 100-200 Myrs ago. This suggests a coupling between the star formation activity and atomic gas on this timescale. No strong correlation between the ionized gas velocity dispersion and the star formation activity between 5-500 Myrs ago was found. The sample and analysis are the foundation of a larger program aimed at understanding the timescales over which stellar feedback drives turbulence.


title: Measuring High-Order Phonon Correlations in an Optomechanical Resonator

abstract: We use single photon detectors to probe the motional state of a superfluid $^4$He resonator of mass $\sim1$ ng. The arrival times of Stokes and anti-Stokes photons (scattered by the resonator's acoustic mode) are used to measure the resonator's phonon coherences up to the fourth order. By post-selecting on photon detection events, we also measure coherences in the resonator when $\leq3$ phonons have been added or subtracted. These measurements are found to be consistent with predictions that assume the acoustic mode to be in thermal equilibrium with a bath through a Markovian coupling.


title: Learning grammar with a divide-and-concur neural network

abstract: We implement a divide-and-concur iterative projection approach to context-free grammar inference. Unlike most state-of-the-art models of natural language processing, our method requires a relatively small number of discrete parameters, making the inferred grammar directly interpretable -- one can read off from a solution how to construct grammatically valid sentences. Another advantage of our approach is the ability to infer meaningful grammatical rules from just a few sentences, compared to the hundreds of gigabytes of training data many other models employ. We demonstrate several ways of applying our approach: classifying words and inferring a grammar from scratch, taking an existing grammar and refining its categories and rules, and taking an existing grammar and expanding its lexicon as it encounters new words in new data.


title: A Deep Learning Approach for Semantic Segmentation of Unbalanced Data in Electron Tomography of Catalytic Materials

abstract: Heterogeneous catalysts possess complex surface and bulk structures, relatively poor intrinsic contrast, and often a sparse distribution of the catalytic nanoparticles (NPs), posing a significant challenge for image segmentation, including the current state-of-the-art deep learning methods. To tackle this problem, we apply a deep learning-based approach for the multi-class semantic segmentation of a $\gamma$-Alumina/Pt catalytic material in a class imbalance situation. Specifically, we used the weighted focal loss as a loss function and attached it to the U-Net's fully convolutional network architecture. We assessed the accuracy of our results using Dice similarity coefficient (DSC), recall, precision, and Hausdorff distance (HD) metrics on the overlap between the ground-truth and predicted segmentations. Our adopted U-Net model with the weighted focal loss function achieved an average DSC score of 0.96 $\pm$ 0.003 in the $\gamma$-Alumina support material and 0.84 $\pm$ 0.03 in the Pt NPs segmentation tasks. We report an average boundary-overlap error of less than 2 nm at the 90th percentile of HD for $\gamma$-Alumina and Pt NPs segmentations. The complex surface morphology of the $\gamma$-Alumina and its relation to the Pt NPs were visualized in 3D by the deep learning-assisted automatic segmentation of a large data set of high-angle annular dark-field (HAADF) scanning transmission electron microscopy (STEM) tomography reconstructions.


title: Classification of left-invariant Einstein metrics on $\mathrm{SL}(2,\mathbb{R})\times\mathrm{SL}(2,\mathbb{R})$ that are bi-invariant under a one-parameter subgroup

abstract: We classify all left-invariant pseudo-Riemannian Einstein metrics on $\mathrm{SL}(2,\mathbb{R})\times \mathrm{SL}(2,\mathbb{R})$ that are bi-invariant under a one-parameter subgroup. We find that there are precisely two such metrics up to homothety, the Killing form and a nearly pseudo-K\"ahler metric.


title: Lung Swapping Autoencoder: Learning a Disentangled Structure-texture Representation of Chest Radiographs

abstract: Well-labeled datasets of chest radiographs (CXRs) are difficult to acquire due to the high cost of annotation. Thus, it is desirable to learn a robust and transferable representation in an unsupervised manner to benefit tasks that lack labeled data. Unlike natural images, medical images have their own domain prior; e.g., we observe that many pulmonary diseases, such as the COVID-19, manifest as changes in the lung tissue texture rather than the anatomical structure. Therefore, we hypothesize that studying only the texture without the influence of structure variations would be advantageous for downstream prognostic and predictive modeling tasks. In this paper, we propose a generative framework, the Lung Swapping Autoencoder (LSAE), that learns factorized representations of a CXR to disentangle the texture factor from the structure factor. Specifically, by adversarial training, the LSAE is optimized to generate a hybrid image that preserves the lung shape in one image but inherits the lung texture of another. To demonstrate the effectiveness of the disentangled texture representation, we evaluate the texture encoder $Enc^t$ in LSAE on ChestX-ray14 (N=112,120), and our own multi-institutional COVID-19 outcome prediction dataset, COVOC (N=340 (Subset-1) + 53 (Subset-2)). On both datasets, we reach or surpass the state-of-the-art by finetuning $Enc^t$ in LSAE that is 77% smaller than a baseline Inception v3. Additionally, in semi-and-self supervised settings with a similar model budget, $Enc^t$ in LSAE is also competitive with the state-of-the-art MoCo. By "re-mixing" the texture and shape factors, we generate meaningful hybrid images that can augment the training set. This data augmentation method can further improve COVOC prediction performance. The improvement is consistent even when we directly evaluate the Subset-1 trained model on Subset-2 without any fine-tuning.


title: A note on the series' of ordinal numbers

abstract: The aim of this short note is to provide a proof to a statement of Sierpi\'nski concerning the number of possible sums of a series (of type $\lambda<\aleph_1$) of arbitrary ordinal numbers.


title: Problems related to Waring-Goldbach problem involving cubes of primes

abstract: In this note, we try to understand the recent development on the Waring-Goldbach problem involving cubes of primes. Especially, we want to determine whether integers that are either primes, squares of primes, cubes of primes, or a cube of an even number can be written as the sum of four cubes of primes. Meanwhile, we raise some problems that may deepen our understanding of the problem about the sum of four cubes of primes. Moreover, some examples suggest that almost all the cubes of integers can be written as the sum of cubes of four integers.


title: Learning physics-constrained subgrid-scale closures in the small-data regime for stable and accurate LES

abstract: We demonstrate how incorporating physics constraints into convolutional neural networks (CNNs) enables learning subgrid-scale (SGS) closures for stable and accurate large-eddy simulations (LES) in the small-data regime (i.e., when the availability of high-quality training data is limited). Using several setups of forced 2D turbulence as the testbeds, we examine the {\it a priori} and {\it a posteriori} performance of three methods for incorporating physics: 1) data augmentation (DA), 2) CNN with group convolutions (GCNN), and 3) loss functions that enforce a global enstrophy-transfer conservation (EnsCon). While the data-driven closures from physics-agnostic CNNs trained in the big-data regime are accurate and stable, and outperform dynamic Smagorinsky (DSMAG) closures, their performance substantially deteriorate when these CNNs are trained with 40x fewer samples (the small-data regime). We show that CNN with DA and GCNN address this issue and each produce accurate and stable data-driven closures in the small-data regime. Despite its simplicity, DA, which adds appropriately rotated samples to the training set, performs as well or in some cases even better than GCNN, which uses a sophisticated equivariance-preserving architecture. EnsCon, which combines structural modeling with aspect of functional modeling, also produces accurate and stable closures in the small-data regime. Overall, GCNN+EnCon, which combines these two physics constraints, shows the best {\it a posteriori} performance in this regime. These results illustrate the power of physics-constrained learning in the small-data regime for accurate and stable LES.


title: Learning Tensor Representations for Meta-Learning

abstract: We introduce a tensor-based model of shared representation for meta-learning from a diverse set of tasks. Prior works on learning linear representations for meta-learning assume that there is a common shared representation across different tasks, and do not consider the additional task-specific observable side information. In this work, we model the meta-parameter through an order-$3$ tensor, which can adapt to the observed task features of the task. We propose two methods to estimate the underlying tensor. The first method solves a tensor regression problem and works under natural assumptions on the data generating process. The second method uses the method of moments under additional distributional assumptions and has an improved sample complexity in terms of the number of tasks. We also focus on the meta-test phase, and consider estimating task-specific parameters on a new task. Substituting the estimated tensor from the first step allows us estimating the task-specific parameters with very few samples of the new task, thereby showing the benefits of learning tensor representations for meta-learning. Finally, through simulation and several real-world datasets, we evaluate our methods and show that it improves over previous linear models of shared representations for meta-learning.


title: Radio Spectra of Luminous, Heavily Obscured WISE-NVSS Selected Quasars

abstract: We present radio spectra spanning $0.1 - 10$ GHz for the sample of heavily obscured luminous quasars with extremely red mid-infrared-optical colors and compact radio emission. The spectra are constructed from targeted 10 GHz observations and archival radio survey data, which together yield $6-11$ flux density measurements for each object. Our suite of Python tools for modeling the radio spectra is publicly available on Github. Our primary result is that most (61%) of the sample have peaked or curved radio spectra and many (36%) could be classified as Gigahertz Peaked Spectrum (GPS) sources. This indicates compact emission regions likely arising from recently triggered radio jets. Assuming synchrotron self-absorption (SSA) generates the peaks, we infer compact source sizes ($3 - 100$ pc) with strong magnetic fields ($6 - 100$ mG) and young ages ($30 - 10^4$ years). Conversely, free-free absorption (FFA) could also create peaks due to the high column densities associated with the deeply embedded nature of the sample. However, we find no correlations between the existence or frequency of the peaks and any parameters of the MIR emission. The high-frequency spectral indices are steep ($\alpha \approx -1$) and correlate, weakly, with the ratio of MIR photon energy density to magnetic energy density, suggesting that the spectral steepening could arise from inverse Compton scattering off the intense MIR photon field. This study provides a foundation for combining multi-frequency and mixed-resolution radio survey data for understanding the impact of young radio jets on the ISM and star formation rates of their host galaxies.


title: Bamboo Trimming Revisited: Simple Algorithms Can Do Well Too

abstract: The bamboo trimming problem considers $n$ bamboo with growth rates $h_1, h_2, \ldots, h_n$ satisfying $\sum_i h_i = 1$. During a given unit of time, each bamboo grows by $h_i$, and then the bamboo-trimming algorithm gets to trim one of the bamboo back down to height zero. The goal is to minimize the height of the tallest bamboo, also known as the backlog. The bamboo trimming problem is closely related to many scheduling problems, and can be viewed as a variation of the widely-studied fixed-rate cup game, but with constant-factor resource augmentation. Past work has given sophisticated pinwheel algorithms that achieve the optimal backlog of 2 in the bamboo trimming problem. It remained an open question, however, whether there exists a simple algorithm with the same guarantee -- recent work has devoted considerable theoretical and experimental effort to answering this question. Two algorithms, in particular, have appeared as natural candidates: the Reduce-Max algorithm (which always cuts the tallest bamboo) and the Reduce-Fastest$(x)$ algorithm (which cuts the fastest-growing bamboo out of those that have height at least $x$). Both algorithms are conjectured to achieve backlog 2. This paper improves the bounds for both Reduce-Fastest and Reduce-Max. Among other results, we show that the exact optimal backlog for Reduce-Fastest$(x)$ is $x + 1$ for all $x \ge 2$ (proving a conjecture of D'Emidio, Di Stefano, and Navarra in the case of $x = 2$), and we show that Reduce-Fastest$(1)$ does not achieve backlog 2 (disproving a conjecture of D'Emidio, Di Stefano, and Navarra). Finally, we show that there is a different algorithm, which we call the Deadline-Driven Strategy, that is both very simple and achieves the optimal backlog of 2. This resolves the question as to whether there exists a simple worst-case optimal algorithm for the bamboo trimming problem.