title: Bounds for the degree and Betti sequences along a graded resolution

abstract: The main goal of this paper is to size up the minimal graded free resolution of a homogeneous ideal in terms of its generating degrees. By and large, this is too ambitious an objective. As understood, sizing up means looking closely at the two available parameters: the shifts and the Betti numbers. Since, in general, bounds for the shifts can behave quite steeply, we filter the difficulty by the subadditivity of the syzygies. The method we applied is hopefully new and sheds additional light on the structure of the minimal free resolution. We use the Boij-S\"oderberg techniques for the Betti numbers to get polynomial upper bounds for them. It is expected that the landscape of hypersurface singularities already contains most of the difficult corners of the arbitrary case. In this regard, we treat some facets of this case, including an improvement of one of the basic results of a recent work by Bus\'e--Dimca--Schenck--Sticlaru.


title: Far and extreme UV radiation feedback in molecular clouds and its influence on the mass and size of star clusters

abstract: We study the formation of star clusters in molecular clouds by performing three-dimensional radiation hydrodynamics simulations with far ultraviolet (FUV; $6 ~{\rm eV} \leqq h \nu \leqq 13.6 ~{\rm eV}$) and extreme ultraviolet (EUV; $h\nu \geqq 13.6~{\rm eV}$) radiative feedback. We find that the FUV feedback significantly suppresses the star formation in diffuse clouds with the initial surface densities of $\Sigma_{\rm cl} \lesssim \rm 50~M_{\odot} \; pc^{-2}$. In the cases of clouds with $\Sigma_{\rm cl} \sim \rm 100-200~M_{\odot} \; pc^{-2}$, the EUV feedback plays a main role and decrease the star formation efficiencies less than $0.3$. We show that thermal pressure from PDRs or H{\sc ii} regions disrupts the clouds and makes the size of the star clusters larger. Consequently, the clouds with the mass $M_{\rm cl} \lesssim 10^{5}~\rm M_{\odot}$ and the surface density $\Sigma_{\rm cl} \lesssim 200~\rm M_{\odot}\; pc^{-2}$ remain the star clusters with the stellar densities of $\sim 100~\rm M_{\odot}\; pc^{-3}$ that nicely match the observed open clusters in the Milky Way. If the molecular clouds are massive ($M_{\rm cl} \gtrsim 10^{5}~\rm M_{\odot}$) and compact ($\Sigma \gtrsim 400~\rm M_{\odot}\; pc^{-2}$), the radiative feedback is not effective and they form massive dense cluster with the stellar densities of $\sim 10^{4}~\rm M_{\odot}\; pc^{-3}$ like observed globular clusters or young massive star clusters. Thus, we suggest that the radiative feedback and the initial conditions of molecular clouds are key factors inducing the variety of the observed star clusters.


title: Patapasco: A Python Framework for Cross-Language Information Retrieval Experiments

abstract: While there are high-quality software frameworks for information retrieval experimentation, they do not explicitly support cross-language information retrieval (CLIR). To fill this gap, we have created Patapsco, a Python CLIR framework. This framework specifically addresses the complexity that comes with running experiments in multiple languages. Patapsco is designed to be extensible to many language pairs, to be scalable to large document collections, and to support reproducible experiments driven by a configuration file. We include Patapsco results on standard CLIR collections using multiple settings.


title: Razmecheno: Named Entity Recognition from Digital Archive of Diaries "Prozhito"

abstract: The vast majority of existing datasets for Named Entity Recognition (NER) are built primarily on news, research papers and Wikipedia with a few exceptions, created from historical and literary texts. What is more, English is the main source for data for further labelling. This paper aims to fill in multiple gaps by creating a novel dataset "Razmecheno", gathered from the diary texts of the project "Prozhito" in Russian. Our dataset is of interest for multiple research lines: literary studies of diary texts, transfer learning from other domains, low-resource or cross-lingual named entity recognition. Razmecheno comprises 1331 sentences and 14119 tokens, sampled from diaries, written during the Perestroika. The annotation schema consists of five commonly used entity tags: person, characteristics, location, organisation, and facility. The labelling is carried out on the crowdsourcing platfrom Yandex.Toloka in two stages. First, workers selected sentences, which contain an entity of particular type. Second, they marked up entity spans. As a result 1113 entities were obtained. Empirical evaluation of Razmecheno is carried out with off-the-shelf NER tools and by fine-tuning pre-trained contextualized encoders. We release the annotated dataset for open access.


title: On module categories related to Sp(N-1) \subset Sl(N)

abstract: Let $V=\C^N$ with $N$ odd. We construct a $q$-deformation of $\End_{Sp(N-1)}(V^{\otimes n})$ which contains $\End_{U_q\sl_N}(V^{\otimes n})$. It is a quotient of an abstract two-variable algebra which is defined by adding one more generator to the generators of the Hecke algebras $H_n$. These results suggest the existence of module categories of $Rep(U_q\sl_N)$ which may not come from already known coideal subalgebras of $U_q\sl_N$. We moreover indicate how this can be used to construct module categories of the associated fusion tensor categories as well as subfactors, along the lines of previous work for inclusions $Sp(N)\subset SL(N)$.


title: Quantum Limits to Classically Spoofing an Electromagnetic Signal

abstract: Spoofing an electromagnetic signal involves measuring its properties and preparing a spoof signal that is a close enough copy to fool a receiver. A classic application of spoofing is in radar where an airborne target attempts to avoid being tracked by a ground-based radar by emitting pulses indicating a false range or velocity. In certain scenarios it has been shown that a sensor can exploit quantum mechanics to detect spoofing at the single-photon level. Here we analyze an idealized spoofing scenario where a transmitter-receiver pair, seeking to detect spoofing, utilizes a signal chosen randomly from a set of non-orthogonal, coherent states. We show that a spoofer optimally employing classical information on the state of the transmitted signal (i.e. the best measure-and-prepare strategy allowed by quantum mechanics) inevitably emits imperfect spoofs that can be exploited by the receiver to reveal the presence of the spoofer, or to discriminate between true reflections and spoofs. Importantly, we show that the quantum limitations on classical spoofing remain significant even in the large mean-photon-number regime.


abstract: Given a union of non-linear manifolds, non-linear subspace clustering or manifold clustering aims to cluster data points based on manifold structures and also learn to parameterize each manifold as a linear subspace in a feature space. Deep neural networks have the potential to achieve this goal under highly non-linear settings given their large capacity and flexibility. We argue that achieving manifold clustering with neural networks requires two essential ingredients: a domain-specific constraint that ensures the identification of the manifolds, and a learning algorithm for embedding each manifold to a linear subspace in the feature space. This work shows that many constraints can be implemented by data augmentation. For subspace feature learning, Maximum Coding Rate Reduction (MCR$^2$) objective can be used. Putting them together yields {\em Neural Manifold Clustering and Embedding} (NMCE), a novel method for general purpose manifold clustering, which significantly outperforms autoencoder-based deep subspace clustering. Further, on more challenging natural image datasets, NMCE can also outperform other algorithms specifically designed for clustering. Qualitatively, we demonstrate that NMCE learns a meaningful and interpretable feature space. As the formulation of NMCE is closely related to several important Self-supervised learning (SSL) methods, we believe this work can help us build a deeper understanding on SSL representation learning.


title: E-ADDA: Unsupervised Adversarial Domain Adaptation Enhanced by a New Mahalanobis Distance Loss for Smart Computing

abstract: In smart computing, the labels of training samples for a specific task are not always abundant. However, the labels of samples in a relevant but different dataset are available. As a result, researchers have relied on unsupervised domain adaptation to leverage the labels in a dataset (the source domain) to perform better classification in a different, unlabeled dataset (target domain). Existing non-generative adversarial solutions for UDA aim at achieving domain confusion through adversarial training. The ideal scenario is that perfect domain confusion is achieved, but this is not guaranteed to be true. To further enforce domain confusion on top of the adversarial training, we propose a novel UDA algorithm, \textit{E-ADDA}, which uses both a novel variation of the Mahalanobis distance loss and an out-of-distribution detection subroutine. The Mahalanobis distance loss minimizes the distribution-wise distance between the encoded target samples and the distribution of the source domain, thus enforcing additional domain confusion on top of adversarial training. Then, the OOD subroutine further eliminates samples on which the domain confusion is unsuccessful. We have performed extensive and comprehensive evaluations of E-ADDA in the acoustic and computer vision modalities. In the acoustic modality, E-ADDA outperforms several state-of-the-art UDA algorithms by up to 29.8%, measured in the f1 score. In the computer vision modality, the evaluation results suggest that we achieve new state-of-the-art performance on popular UDA benchmarks such as Office-31 and Office-Home, outperforming the second best-performing algorithms by up to 17.9%.


title: Heat Conduction Plate Layout Optimization using Physics-driven Convolutional Neural Networks

abstract: The layout optimization of the heat conduction is essential during design in engineering, especially for thermal sensible products. When the optimization algorithm iteratively evaluates different loading cases, the traditional numerical simulation methods used usually lead to a substantial computational cost. To effectively reduce the computational effort, data-driven approaches are used to train a surrogate model as a mapping between the prescribed external loads and various geometry. However, the existing model are trained by data-driven methods which requires intensive training samples that from numerical simulations and not really effectively solve the problem. Choosing the steady heat conduction problems as examples, this paper proposes a Physics-driven Convolutional Neural Networks (PD-CNN) method to infer the physical field solutions for random varied loading cases. After that, the Particle Swarm Optimization (PSO) algorithm is used to optimize the sizes and the positions of the hole masks in the prescribed design domain, and the average temperature value of the entire heat conduction field is minimized, and the goal of minimizing heat transfer is achieved. Compared with the existing data-driven approaches, the proposed PD-CNN optimization framework not only predict field solutions that are highly consistent with conventional simulation results, but also generate the solution space with without any pre-obtained training data.


title: Reinforcement Learning Your Way: Agent Characterization through Policy Regularization

abstract: The increased complexity of state-of-the-art reinforcement learning (RL) algorithms have resulted in an opacity that inhibits explainability and understanding. This has led to the development of several post-hoc explainability methods that aim to extract information from learned policies thus aiding explainability. These methods rely on empirical observations of the policy and thus aim to generalize a characterization of agents' behaviour. In this study, we have instead developed a method to imbue a characteristic behaviour into agents' policies through regularization of their objective functions. Our method guides the agents' behaviour during learning which results in an intrinsic characterization; it connects the learning process with model explanation. We provide a formal argument and empirical evidence for the viability of our method. In future work, we intend to employ it to develop agents that optimize individual financial customers' investment portfolios based on their spending personalities.


title: Thixotropic spectra and Ashby-style charts for thixotropy

abstract: There is no universal model for thixotropy, and comparing thixotropic effects between different fluids is a subtle yet challenging problem. We introduce a generalized (model-insensitive) framework for comparing thixotropic properties based on thixotropic spectra. A superposition of exponential stress modes distributed over thixotropic timescales is used to quantify buildup and breakdown times and mode strengths in response to step-change input. This mathematical framework is tested with several experimental step-shear rate data on colloidal suspensions. Low-dimensional metrics based on moments of the distribution reveal characteristic average thixotropic properties which are visualized on Ashby-style diagrams. This method outlines a framework for describing thixotropy across a diverse range of microstructures, supporting scientific studies as well as material selection for engineering design applications.


title: Text and Code Embeddings by Contrastive Pre-Training

abstract: Text embeddings are useful features in many applications such as semantic search and computing text similarity. Previous work typically trains models customized for different use cases, varying in dataset choice, training objective and model architecture. In this work, we show that contrastive pre-training on unsupervised data at scale leads to high quality vector representations of text and code. The same unsupervised text embeddings that achieve new state-of-the-art results in linear-probe classification also display impressive semantic search capabilities and sometimes even perform competitively with fine-tuned models. On linear-probe classification accuracy averaging over 7 tasks, our best unsupervised model achieves a relative improvement of 4% and 1.8% over previous best unsupervised and supervised text embedding models respectively. The same text embeddings when evaluated on large-scale semantic search attains a relative improvement of 23.4%, 14.7%, and 10.6% over previous best unsupervised methods on MSMARCO, Natural Questions and TriviaQA benchmarks, respectively. Similarly to text embeddings, we train code embedding models on (text, code) pairs, obtaining a 20.8% relative improvement over prior best work on code search.


title: Quantum density estimation with density matrices: Application to quantum anomaly detection

abstract: Density estimation is a central task in statistics and machine learning. This problem aims to determine the underlying probability density function that best aligns with an observed data set. Some of its applications include statistical inference, unsupervised learning, and anomaly detection. Despite its relevance, few works have explored the application of quantum computing to density estimation. In this article, we present a novel quantum-classical density matrix density estimation model, called Q-DEMDE, based on the expected values of density matrices and a novel quantum embedding called quantum Fourier features. The method uses quantum hardware to build probability distributions of training data via mixed quantum states. As a core subroutine, we propose a new algorithm to estimate the expected value of a mixed density matrix from its spectral decomposition on a quantum computer. In addition, we present an application of the method for quantum-classical anomaly detection. We evaluated the density estimation model with quantum random and quantum adaptive Fourier features on different data sets on a quantum simulator and a real quantum computer. An important result of this work is to show that it is possible to perform density estimation and anomaly detection with high performance on present-day quantum computers.


title: Evolutionary timeline of a modeled cell

abstract: A theoretical study of cell evolution is presented here. By using a toolbox containing an intracellular catalytic reaction network model and a mutation-selection process, four distinct phases of self-organization were unveiled. First, the nutrients prevail as the central substrate of the chemical reactions. Second, the cell becomes a small-world. Third, a highly connected core component emerges, concurrently with the nutrient carriers becoming the central product of reactions. Finally, the cell reaches a steady configuration where the concentrations of the core chemical species are described by Zipf's law.


title: Rheological design of thickened alcohol-based hand rubs

abstract: The handleability and sensory perception of hand sanitisers by consumers affect the hygiene outcome. Spillage may result in under-dosing and poor sensory properties can lead to under-utilisation. We first propose four principles (low run off, spreadability, smoothness and non-stickiness) for designing the rheology of thickened alcohol-based hand rubs with acceptable handleability and hand feel. We then evaluate a commercial hand gel and a variety of simplified formulations thickened with microgels (Carbopol 974P, Carbopol Ultrez 20 and Sepimax Zen), or linear polymers (Jaguar HP 120 COS), and evaluate them against these design criteria. All four additives provide acceptable spreadability by shear thinning to $\eta \approx 10^{-1} Pa \, s$ at $\dot\gamma \sim 10^3 s^{-1}$. Either the finite yield stress conferred by the microgels ($\sigma_y \gtrsim 10 Pa$) or the increase in low-shear viscosity provided by the linear polymer ($\eta \gtrsim 1 Pa \, s$ at $\dot\gamma \lesssim 0.1 s^{-1}$) give rise to acceptably low run-off. However, the formulation using the linear polymer shows a filament breakage time of $\tau_{\rm b} \approx 1 s$ in capillary rheology, which may result in stickiness and therefore a less than optimal hand feel.


title: Photoinduced structural dynamics of multiferroic TbMnO$_3$

abstract: We use time-resolved hard x-ray diffraction to investigate the structural dynamics of the multiferroic insulator TbMnO$_3$ in the low temperature antiferromagnetic and ferroelectrically ordered phase. The lattice response following photoexcitation at 1.55 eV is detected by measuring the (0 2 4) and (1 3 -5) Bragg reflections. A 0.02% tensile strain, normal to the surface, is seen to arise within 20 - 30 ps. The magnitude of this transient strain is over an order of magnitude lower than that predicted from laser-induced heating, which we attribute to a bottleneck in the energy transfer between the electronic and lattice subsystems. The timescale for the transient expansion is consistent with that of previously reported demagnetization dynamics. We discuss a possible relationship between structural and demagnetization dynamics in TbMnO$_3$, in which photoinduced atomic motion modulates the exchange interaction, leading to a destruction of the magnetic order in the system.


title: Low-energy physics for an iron phthalocyanine molecule on Au(111)

abstract: The system of an iron phthalocyanine molecule on the Au(111) surface, has been studied recently due to its peculiar properties. In particular, several surprising results of scanning tunneling spectroscopy changing the position of the molecule and applying magnetic field can be explained by the {\it non-Landau} Fermi liquid state of a 2-channel spin-1 Kondo model with anisotropy. The localized orbitals near the Fermi level are three, one of symmetry $z^2$ and two (nearly) degenerate $\pi$ orbitals of symmetry $xz$ and $yz$. Previous studies using the numerical renormalization group neglected one of these orbitals to render the problem tractable. Here we investigate, using a slave-boson mean-field approximation, if the splitting $S$ between $\pi$ orbitals caused by spin-orbit coupling (SOC) justifies this approximation. We obtain an abrupt transition from a 3-band regime to a 2-band one at a value of $S$ which is about 1/3 of the atomic SOC for Fe, justifying the 2-band model for the system.


title: Constant Inapproximability for PPA

abstract: In the $\varepsilon$-Consensus-Halving problem, we are given $n$ probability measures $v_1, \dots, v_n$ on the interval $R = [0,1]$, and the goal is to partition $R$ into two parts $R^+$ and $R^-$ using at most $n$ cuts, so that $|v_i(R^+) - v_i(R^-)| \leq \varepsilon$ for all $i$. This fundamental fair division problem was the first natural problem shown to be complete for the class PPA, and all subsequent PPA-completeness results for other natural problems have been obtained by reducing from it. We show that $\varepsilon$-Consensus-Halving is PPA-complete even when the parameter $\varepsilon$ is a constant. In fact, we prove that this holds for any constant $\varepsilon < 1/5$. As a result, we obtain constant inapproximability results for all known natural PPA-complete problems, including Necklace-Splitting, the Discrete-Ham-Sandwich problem, two variants of the pizza sharing problem, and for finding fair independent sets in cycles and paths.


title: First-Order Game Logic and Modal Mu-Calculus

abstract: This paper investigates first-order game logic and first-order modal mu-calculus, which extend their propositional modal logic counterparts with first-order modalities of interpreted effects such as variable assignments. Unlike in the propositional case, both logics are shown to have the same expressive power and their proof calculi to have the same deductive power. Both calculi are also mutually relatively complete. In the presence of differential equations, corollaries obtain usable and complete translations between differential game logic, a logic for the deductive verification of hybrid games, and the differential mu-calculus, the modal mu-calculus for hybrid systems. The differential mu-calculus is complete with respect to first-order fixpoint logic and differential game logic is complete with respect to its ODE-free fragment.


title: Effects of operator nonlocality on closures for multicomponent reactive flows based on dispersion analysis

abstract: Algebraic closure models with spatially nonlocal operators that are associated with both unresolved advective transport and nonlinear reaction terms in a Reynolds-averaged Navier-Stokes context are presented in this work. In particular, a system of two species subject to binary reaction and transport by advection and diffusion are examined by expanding upon analysis originally developed for binary reactions in the context of Taylor dispersion of scalars. This work extends model forms from weakly-nonlinear extensions of that dispersion theory and the role of nonlocality in the presence of reactions is studied and captured by analytic expressions. These expressions can be incorporated into an eddy diffusivity matrix that explicitly capture the influence of chemical kinetics and flow conditions on the closure operators and we demonstrate that the model form derived in a laminar context can be directly translated to an analogous setup in homogeneous isotropic turbulence, which has implications as a subgrid scale model. We show that this framework can improve prediction of mean quantities compared to previous purely local results, but does not fully close unresolved terms.


title: Zero-Truncated Poisson Regression for Sparse Multiway Count Data Corrupted by False Zeros

abstract: We propose a novel statistical inference methodology for multiway count data that is corrupted by false zeros that are indistinguishable from true zero counts. Our approach consists of zero-truncating the Poisson distribution to neglect all zero values. This simple truncated approach dispenses with the need to distinguish between true and false zero counts and reduces the amount of data to be processed. Inference is accomplished via tensor completion that imposes low-rank tensor structure on the Poisson parameter space. Our main result shows that an $N$-way rank-$R$ parametric tensor $\boldsymbol{\mathscr{M}}\in(0,\infty)^{I\times \cdots\times I}$ generating Poisson observations can be accurately estimated by zero-truncated Poisson regression from approximately $IR^2\log_2^2(I)$ non-zero counts under the nonnegative canonical polyadic decomposition. Our result also quantifies the error made by zero-truncating the Poisson distribution when the parameter is uniformly bounded from below. Therefore, under a low-rank multiparameter model, we propose an implementable approach guaranteed to achieve accurate regression in under-determined scenarios with substantial corruption by false zeros. Several numerical experiments are presented to explore the theoretical results.

categories: stat.ME cs.MS cs.NA math.NA math.ST stat.ML stat.TH


title: Automatic Recognition and Digital Documentation of Cultural Heritage Hemispherical Domes using Images

abstract: Advancements in optical metrology has enabled documentation of dense 3D point clouds of cultural heritage sites. For large scale and continuous digital documentation, processing of dense 3D point clouds becomes computationally cumbersome, and often requires additional hardware for data management, increasing the time cost, and complexity of projects. To this end, this manuscript presents an original approach to generate fast and reliable semantic digital models of heritage hemispherical domes using only two images. New closed formulations were derived to establish the relationships between spheres and their projected ellipses onto images, which fostered the development of a new automatic framework for as-built generation of spheres. The effectiveness of the proposed method was evaluated under both laboratory and real-world datasets. The results revealed that the proposed method achieved as-built modeling accuracy of around 6mm, while improving the computation time by a factor of 7, when compared to established point cloud processing methods.


title: First direct measurement of the $^{13}$N($\alpha$,$p$)$^{16}$O reaction relevant for core-collapse supernovae nucleosynthesis

abstract: Understanding the explosion mechanism of a core-collapse supernova (CCSN) is important to accurately model CCSN scenarios for different progenitor stars using model-observation comparisons. The uncertainties of various nuclear reaction rates relevant for CCSN scenarios strongly affect the accuracy of these stellar models. Out of these reactions, the $^{13}$N($\alpha$,$p$)$^{16}$O reaction has been found to affect various stages of a CCSN at varying temperatures. This work presents the first direct measurement of the $^{13}$N($\alpha$,$p$)$^{16}$O reaction performed using a 34.6 MeV beam of radioactive $^{13}$N ions and the active-target detector MUSIC (MUlti-Sampling Ionization Chamber) at Argonne National Laboratory. The resulting total $^{13}$N($\alpha$,$p$)$^{16}$O reaction cross sections from this measurement in the center-of-mass energy range of 3.26 - 6.02 MeV are presented and compared with calculations using the Hauser-Feshbach formalism. Uncertainties in the reaction rate have been dramatically reduced at CCSN temperatures.


title: Online Convex Optimization Using Coordinate Descent Algorithms

abstract: This paper considers the problem of online optimization where the objective function is time-varying. In particular, we extend coordinate descent type algorithms to the online case, where the objective function varies after a finite number of iterations of the algorithm. Instead of solving the problem exactly at each time step, we only apply a finite number of iterations at each time step. Commonly used notions of regret are used to measure the performance of the online algorithm. Moreover, coordinate descent algorithms with different updating rules are considered, including both deterministic and stochastic rules that are developed in the literature of classical offline optimization. A thorough regret analysis is given for each case. Finally, numerical simulations are provided to illustrate the theoretical results.


title: The VMC survey -- XLV. Proper motion of the outer LMC and the impact of the SMC

abstract: The Large Magellanic Cloud (LMC) is the most luminous satellite galaxy of the Milky Way and owing to its companion, the Small Magellanic Cloud (SMC), represents an excellent laboratory to study the interaction of dwarf galaxies. The aim of this study is to investigate the kinematics of the outer regions of the LMC by using stellar proper motions to understand the impact of interactions, e.g. with the SMC about 250 Myr ago. {We calculate proper motions using multi-epoch $K_\mathrm{s}$-band images from the VISTA survey of the Magellanic Clouds system (VMC). Observations span a time baseline of 2$-$5 yr. We combine the VMC data with data from the Gaia early Data Release 3 and introduce a new method to distinguish between Magellanic and Milky Way stars based on a machine learning algorithm. This new technique enables a larger and cleaner sample selection of fainter sources as it reaches below the red clump of the LMC. We investigate the impact of the SMC on the rotational field of the LMC and find hints of stripped SMC debris. The south east region of the LMC shows a slow rotational speed compared to the overall rotation. $N$-body simulations suggest that this could be caused by a fraction of stripped SMC stars, located in that particular region, that move opposite to the expected rotation.


title: Hydrodynamic effects on the liquid-hexatic transition of active colloids

abstract: We study numerically the role of hydrodynamics in the liquid-hexatic transition of active colloids at intermediate activity, where motility induced phase separation (MIPS) does not occur. We show that in the case of active Brownian particles (ABP), the critical density of the transition decreases upon increasing the particle's mass, enhancing ordering, while self-propulsion has the opposite effect in the activity regime considered. Active hydrodynamic particles (AHP), instead, undergo the liquid-hexatic transition at higher values of packing fraction $\phi$ than the corresponding ABP, suggesting that hydrodynamics have the net effect of disordering the system. At increasing densities, close to the hexatic-liquid transition, we found in the case of AHP the appearance of self-sustained organized motion with clusters of particles moving coherently.


title: The California-Kepler Survey. X. The Radius Gap as a Function of Stellar Mass, Metallicity, and Age

abstract: In 2017, the California-Kepler Survey (CKS) published its first data release (DR1) of high-resolution optical spectra of 1305 planet hosts. Refined CKS planet radii revealed that small planets are bifurcated into two distinct populations: super-Earths (smaller than 1.5 $R_E$) and sub-Neptunes (between 2.0 and 4.0 $R_E$), with few planets in between (the "Radius Gap.") Several theoretical models of the Radius Gap predict variation with stellar mass, but testing these predictions are challenging with CKS DR1 due to its limited $M_\star$ range of 0.8-1.4 $M_\odot$. Here, we present CKS DR2 with 411 additional spectra and derived properties focusing on stars of 0.5-0.8 $M_\odot$. We found the Radius Gap follows $R_p \propto P^m$ with $m = -0.10 \pm 0.03$, consistent with predictions of XUV- and core-powered mass-loss mechanisms. We found no evidence that $m$ varies with $M_\star$. We observed a correlation between the average sub-Neptune size and $M_\star$. Over 0.5 to 1.4 $M_\odot$, the average sub-Neptune grows from 2.1 to 2.6 $R_E$, following $R_p \propto M_\star^\alpha$ with $\alpha = 0.25 \pm 0.03$. In contrast, there is no detectable change for super-Earths. These $M_\star$-$R_p$ trends suggests that protoplanetary disks can efficiently produce cores up to a threshold mass of $M_c$, which grows linearly with stellar mass according to $M_c \approx 10 M_E~(M_\star / M_\odot)$. There is no significant correlation between sub-Neptune size and stellar metallicity (over $-$0.5 to $+$0.5 dex) suggesting a weak relationship between planet envelope opacity and stellar metallicity. Finally, there is no significant variation in sub-Neptune size with stellar age (over 1 to 10 Gyr), which suggests that the majority of envelope contraction concludes after $\sim$1 Gyr.


title: Trap Frequency Measurement with a Pulsed Atom Laser

abstract: We describe a novel method of single-shot trap frequency measurement for a confined Bose-Einstein Condensate, which uses an atom laser to repeatedly sample the mean velocity of trap oscillations as a function of time. The method is able to determine the trap frequency to an accuracy of 39~ppm (16~mHz) in a single experimental realization, improving on the literature by a factor of three. Further, we show that by employing a reconstructive aliasing approach our method can be applied to trap frequencies more than a factor of 3 greater than the sampling frequency.


title: Affine Body Dynamics: Fast, Stable & Intersection-free Simulation of Stiff Materials

abstract: Simulating stiff materials in applications where deformations are either not significant or can safely be ignored is a pivotal task across fields. Rigid body modeling has thus long remained a fundamental tool and is, by far, the most popular simulation strategy currently employed for modeling stiff solids. At the same time, numerical models of a rigid body continue to pose a number of known challenges and trade-offs including intersections, instabilities, inaccuracies, and/or slow performances that grow with contact-problem complexity. In this paper we revisit this problem and present ABD, a simple and highly effective affine body dynamics framework, which significantly improves state-of-the-art stiff simulations. We trace the challenges in the rigid-body IPC (incremental potential contact) method to the necessity of linearizing piecewise-rigid (SE(3)) trajectories and subsequent constraints. ABD instead relaxes the unnecessary (and unrealistic) constraint that each body's motion be exactly rigid with a stiff orthogonality potential, while preserving the rigid body model's key feature of a small coordinate representation. In doing so ABD replaces piecewise linearization with piecewise linear trajectories. This, in turn, combines the best from both parties: compact coordinates ensure small, sparse system solves, while piecewise-linear trajectories enable efficient and accurate constraint (contact and joint) evaluations. Beginning with this simple foundation, ABD preserves all guarantees of the underlying IPC model e.g., solution convergence, guaranteed non-intersection, and accurate frictional contact. Over a wide range and scale of simulation problems we demonstrate that ABD brings orders of magnitude performance gains (two- to three-order on the CPU and an order more utilizing the GPU, which is 10,000x speedups) over prior IPC-based methods with a similar or higher simulation quality.


title: On the location of the excess wing relative to the a-loss peak in the susceptibility spectra from simulations with the swap Monte Carlo algorithm

abstract: An advance was made by Guiselin et al. [arXiv:2103.01569 (2021)] in molecular dynamics simulations of the equilibrium dynamics of supercooled liquids near the experimental glass transition by utilizing the giant equilibration speedup provided by the swap Monte Carlo algorithm. The found emergence of a power law in relaxation spectra at lower temperatures on the high frequency flank of the alpha-loss peak in analogy to the excess wings observed experimentally in molecular glass-formers. Their remarkable finding leads to the question of where the excess wing is located relative to the alpha-loss peak in the susceptibility spectrum. I provide an answer by identifying the excess wing as the unresolved Johari-Goldstein (JG) b-relaxation and using the reciprocal of its relaxation time tau_JG to assess the location of the excess wing. The Coupling Model (CM) has a history of being successful in determining approximately the values of tau_JG(T) in molecular liquids whether the JG beta-relaxation is resolved or not (i.e., excess wing). It is applied to the simulation data and the results of tau_JG(T) successfully account for the locations of the excess wings in the simulation spectra at different temperatures. The time evolution of the dynamics of the distribution of processes composing the JG b-relaxation suggested by the CM based on experimental data in molecular glass-formers are in agreement with that given by Guiselin et al. from simulations of a size-polydisperse mixture of N soft spheres interacting with a repulsive power law pair potential. The agreement brings their results closer to experimental data of real molecular glass-formers.


title: Proteome-scale Deployment of Protein Structure Prediction Workflows on the Summit Supercomputer

abstract: Deep learning has contributed to major advances in the prediction of protein structure from sequence, a fundamental problem in structural bioinformatics. With predictions now approaching the accuracy of crystallographic resolution in some cases, and with accelerators like GPUs and TPUs making inference using large models rapid, fast genome-level structure prediction becomes an obvious aim. Leadership-class computing resources can be used to perform genome-scale protein structure prediction using state-of-the-art deep learning models, providing a wealth of new data for systems biology applications. Here we describe our efforts to efficiently deploy the AlphaFold2 program, for full-proteome structure prediction, at scale on the Oak Ridge Leadership Computing Facility's resources, including the Summit supercomputer. We performed inference to produce the predicted structures for 35,634 protein sequences, corresponding to three prokaryotic proteomes and one plant proteome, using under 4,000 total Summit node hours, equivalent to using the majority of the supercomputer for one hour. We also designed an optimized structure refinement that reduced the time for the relaxation stage of the AlphaFold pipeline by over 10X for longer sequences. We demonstrate the types of analyses that can be performed on proteome-scale collections of sequences, including a search for novel quaternary structures and implications for functional annotation.


title: On the seminormal bases and dual seminormal bases of the cyclotomic Hecke algebras of type $G(\ell,1,n)$

abstract: This paper studies the seminormal bases $\{f_{\mathfrak{s}\mathfrak{t}}\}$ and the dual seminormal bases $\{g_{\mathfrak{s}\mathfrak{t}}\}$ of the non-degenerate and the degenerate cyclotomic Hecke algebras ${H}_{\ell,n}$ of type $G(\ell,1,n)$. We present some explicit formulae for the constants $\alpha_{\mathfrak{s}\mathfrak{t}}:=g_{\mathfrak{s}\mathfrak{t}}/f_{\mathfrak{s}\mathfrak{t}}\in K^\times$ in terms of the $\gamma$-coefficients $\{\gamma_{\mathfrak{u}}, \gamma'_{\mathfrak{u}}\}$ of $H_{\ell,n}$. In particular, we answer a question of Mathas on the rationality of square roots of some quotients of products of $\gamma$-coefficients. We obtain some explicit formulae for the expansion of each seminormal bases of $H_{\ell,n-1}$ as a linear combination of the seminormal bases of $H_{\ell,n}$ under the natural inclusion $H_{\ell,n-1}\hookrightarrow H_{\ell,n}$.


title: Mysterious Odd Radio Circle near the Large Magellanic Cloud -- An Intergalactic Supernova Remnant?

abstract: We report the discovery of J0624-6948, a low-surface brightness radio ring, lying between the Galactic Plane and the Large Magellanic Cloud (LMC). It was first detected at 888 MHz with the Australian Square Kilometre Array Pathfinder (ASKAP), and with a diameter of ~196 arcsec. This source has phenomenological similarities to Odd Radio Circles (ORCs). Significant differences to the known ORCs - a flatter radio spectral index, the lack of a prominent central galaxy as a possible host, and larger apparent size - suggest that J0624-6948 may be a different type of object. We argue that the most plausible explanation for J0624-6948 is an intergalactic supernova remnant due to a star that resided in the LMC outskirts that had undergone a single-degenerate type Ia supernova, and we are seeing its remnant expand into a rarefied, intergalactic environment. We also examine if a massive star or a white dwarf binary ejected from either galaxy could be the supernova progenitor. Finally, we consider several other hypotheses for the nature of the object, including the jets of an active galactic nucleus (AGN) or the remnant of a nearby stellar super-flare.


title: A solution for the quasi-one-dimensional linearised Euler equations with heat transfer

abstract: The unsteady response of nozzles with steady heat transfer forced by acoustic and/or entropy waves is modelled. The approach is based on the quasi-one-dimensional linearised Euler equations. The equations are cast in terms of three variables, namely the dimensionless mass, stagnation temperature and entropy fluctuations, which are invariants of the system at zero frequency and with no heat transfer. The resulting first-order system of differential equations is then solved using the Magnus expansion method, where the perturbation parameters are the normalised frequency and the volumetric heat transfer. In this work, a measure of the flow non-isentropicity (in this case the steady heat transfer) is used for the first time as an expansion parameter. The solution method was applied to a converging-diverging nozzle with constant heat transfer for both sub-critical and super-critical flow cases, showing good agreement with numerical predictions. It was observed that the acoustic and entropy transfer functions of the nozzle strongly depend on the frequency and heat transfer.


title: 4.4 kV $\beta$-Ga$_2$O$_3$ Power MESFETs with Lateral Figure of Merit exceeding 100 MW/cm$^2$

abstract: Field-plated (FP) depletion-mode MOVPE-grown $\beta$-Ga$_2$O$_3$ lateral MESFETs are realized with superior reverse breakdown voltages and ON currents. A sandwiched SiN$_x$ dielectric field plate design was utilized that prevents etching-related damage in the active region and a deep mesa-etching was used to reduce reverse leakage. The device with L$_{GD}$ = 34.5 $\mu$m exhibits an ON current (I$_{DMAX}$) of 56 mA/mm, a high I$_{ON}$/I$_{OFF}$ ratio $>$ 10$^8$ and a very low reverse leakage until catastrophic breakdown at $\sim$ 4.4kV. The highest measurable V$_{BR}$ recorded was 4.57 kV (L$_{GD}$ = 44.5 $\mu$m). An LFOM of 132 MW/cm$^2$ was calculated for a V$_{BR}$ of $\sim$ 4.4 kV. The reported results are the first $>$ 4kV-class Ga$_2$O$_3$ transistors to surpass the theoretical FOM of Silicon. These are also the highest I$_{DMAX}$ and lowest R$_{ON}$ values achieved simultaneously for any $\beta$-Ga$_2$O$_3$ device with V$_{BR}$ $>$ 4kV to date. This work highlights that high breakdown voltages (V$_{BR}$), high lateral figure of merit (LFOM) and high ON currents can be achieved simultaneously in $\beta$-Ga$_2$O$_3$ lateral transistors.


title: PONI: Potential Functions for ObjectGoal Navigation with Interaction-free Learning

abstract: State-of-the-art approaches to ObjectGoal navigation rely on reinforcement learning and typically require significant computational resources and time for learning. We propose Potential functions for ObjectGoal Navigation with Interaction-free learning (PONI), a modular approach that disentangles the skills of `where to look?' for an object and `how to navigate to (x, y)?'. Our key insight is that `where to look?' can be treated purely as a perception problem, and learned without environment interactions. To address this, we propose a network that predicts two complementary potential functions conditioned on a semantic map and uses them to decide where to look for an unseen object. We train the potential function network using supervised learning on a passive dataset of top-down semantic maps, and integrate it into a modular framework to perform ObjectGoal navigation. Experiments on Gibson and Matterport3D demonstrate that our method achieves the state-of-the-art for ObjectGoal navigation while incurring up to 1,600x less computational cost for training. Code and pre-trained models are available: https://vision.cs.utexas.edu/projects/poni/


title: The Pop-stack-sorting Operator on Tamari Lattices

abstract: Motivated by the pop-stack-sorting map on the symmetric groups, Defant defined an operator $\mathsf{Pop}_M : M \to M$ for each complete meet-semilattice $M$ by $$\mathsf{Pop}_M(x)=\bigwedge(\{y\in M: y\lessdot x\}\cup \{x\}).$$ This paper concerns the dynamics of $\mathsf{Pop}_{\mathrm{Tam}_n}$, where $\mathrm{Tam}_n$ is the $n$-th Tamari lattice. We say an element $x\in \mathrm{Tam}_n$ is $t$-$\mathsf{Pop}$-sortable if $\mathsf{Pop}_M^t (x)$ is the minimal element and we let $h_t(n)$ denote the number of $t$-$\mathsf{Pop}$-sortable elements in $\mathrm{Tam}_n$. We find an explicit formula for the generating function $\sum_{n\ge 1}h_t(n)z^n$ and verify Defant's conjecture that it is rational. We furthermore prove that the size of the image of $\mathsf{Pop}_{\mathrm{Tam}_n}$ is the Motzkin number $M_n$, settling a conjecture of Defant and Williams.


title: On the Crawford number attaining operators

abstract: We study the denseness of Crawford number attaining operators on Banach spaces. Mainly, we prove that if a Banach space has the RNP, then the set of Crawford number attaining operators is dense in the space of bounded linear operators. We also see among others that the set of Crawford number attaining operators may be dense in the space of all bounded linear operators while they do not coincide, by observing the case of compact operators when the Banach space has a 1-unconditional basis. Furthermore, we show a Bishop-Phelps-Bollob\'as type property for the Crawford number for certain Banach spaces, and we finally discuss some difficulties and possible problems on the topic.


title: Variational Autoencoders for Reliability Optimization in Multi-Access Edge Computing Networks

abstract: Multi-access edge computing (MEC) is viewed as an integral part of future wireless networks to support new applications with stringent service reliability and latency requirements. However, guaranteeing ultra-reliable and low-latency MEC (URLL MEC) is very challenging due to uncertainties of wireless links, limited communications and computing resources, as well as dynamic network traffic. Enabling URLL MEC mandates taking into account the statistics of the end-to-end (E2E) latency and reliability across the wireless and edge computing systems. In this paper, a novel framework is proposed to optimize the reliability of MEC networks by considering the distribution of E2E service delay, encompassing over-the-air transmission and edge computing latency. The proposed framework builds on correlated variational autoencoders (VAEs) to estimate the full distribution of the E2E service delay. Using this result, a new optimization problem based on risk theory is formulated to maximize the network reliability by minimizing the Conditional Value at Risk (CVaR) as a risk measure of the E2E service delay. To solve this problem, a new algorithm is developed to efficiently allocate users' processing tasks to edge computing servers across the MEC network, while considering the statistics of the E2E service delay learned by VAEs. The simulation results show that the proposed scheme outperforms several baselines that do not account for the risk analyses or statistics of the E2E service delay.


title: Derived category and ACM bundles of moduli space of vector bundles on a curve

abstract: We show that the derived category of a curve is embedded into the derived category of the moduli space of vector bundles on the curve of coprime rank and degree. We also generalize the semiorthogonal decomposition constructed by Narasimhan and Belmans-Mukhopadhyay. Finally, we produce a one-dimensional family of ACM bundles over the moduli space.


title: Self-Supervised Point Cloud Registration with Deep Versatile Descriptors

abstract: As a fundamental yet challenging problem in intelligent transportation systems, point cloud registration attracts vast attention and has been attained with various deep learning-based algorithms. The unsupervised registration algorithms take advantage of deep neural network-enabled novel representation learning while requiring no human annotations, making them applicable to industrial applications. However, unsupervised methods mainly depend on global descriptors, which ignore the high-level representations of local geometries. In this paper, we propose to jointly use both global and local descriptors to register point clouds in a self-supervised manner, which is motivated by a critical observation that all local geometries of point clouds are transformed consistently under the same transformation. Therefore, local geometries can be employed to enhance the representation ability of the feature extraction module. Moreover, the proposed local descriptor is flexible and can be integrated into most existing registration methods and improve their performance. Besides, we also utilize point cloud reconstruction and normal estimation to enhance the transformation awareness of global and local descriptors. Lastly, extensive experimental results on one synthetic and three real-world datasets demonstrate that our method outperforms existing state-of-art unsupervised registration methods and even surpasses supervised ones in some cases. Robustness and computational efficiency evaluations also indicate that the proposed method applies to intelligent vehicles.


title: Maximizing information from chemical engineering data sets: Applications to machine learning

abstract: It is well-documented how artificial intelligence can have (and already is having) a big impact on chemical engineering. But classical machine learning approaches may be weak for many chemical engineering applications. This review discusses how challenging data characteristics arise in chemical engineering applications. We identify four characteristics of data arising in chemical engineering applications that make applying classical artificial intelligence approaches difficult: (1) high variance, low volume data, (2) low variance, high volume data, (3) noisy/corrupt/missing data, and (4) restricted data with physics-based limitations. For each of these four data characteristics, we discuss applications where these data characteristics arise and show how current chemical engineering research is extending the fields of data science and machine learning to incorporate these challenges. Finally, we identify several challenges for future research.


title: Carbon Explorer: A Holistic Approach for Designing Carbon Aware Datacenters

abstract: Technology companies have been leading the way to a renewable energy transformation, by investing in renewable energy sources to reduce the carbon footprint of their datacenters. In addition to helping build new solar and wind farms, companies make power purchase agreements or purchase carbon offsets, rather than relying on renewable energy every hour of the day, every day of the week (24/7). Relying on renewable energy 24/7 is challenging due to the intermittent nature of wind and solar energy. Inherent variations in solar and wind energy production causes excess or lack of supply at different times. To cope with the fluctuations of renewable energy generation, multiple solutions must be applied. These include: capacity sizing with a mix of solar and wind power, energy storage options, and carbon aware workload scheduling. However, depending on the region and datacenter workload characteristics, the carbon-optimal solution varies. Existing work in this space does not give a holistic view of the trade-offs of each solution and often ignore the embodied carbon cost of the solutions. In this work, we provide a framework, Carbon Explorer, to analyze the multi-dimensional solution space by taking into account operational and embodided footprint of the solutions to help make datacenters operate on renewable energy 24/7. The solutions we analyze include capacity sizing with a mix of solar and wind power, battery storage, and carbon aware workload scheduling, which entails shifting the workloads from times when there is lack of renewable supply to times with abundant supply.


title: Diversity Enhancement via Magnitude

abstract: Promoting and maintaining diversity of candidate solutions is a key requirement of evolutionary algorithms in general and multi-objective evolutionary algorithms in particular. In this paper, we use the recently developed theory of magnitude to construct a gradient flow and similar notions that systematically manipulate finite subsets of Euclidean space to enhance their diversity, and apply the ideas in service of multi-objective evolutionary algorithms. We demonstrate diversity enhancement on benchmark problems using leading algorithms, and discuss extensions of the framework.


title: Pairing symmetry in infinite-layer nickelate superconductor

abstract: The superconducting infinite-layer nickelate family has risen as a promising platform for revealing the mechanism of high-temperature superconductivity. However, its challenging material synthesis has obscured effort in understanding the nature of its ground state and low-lying excitations, which is a prerequisite for identifying the origin of the Cooper pairing in high-temperature superconductors. In particular, the superconducting gap symmetry of nickelates has hardly been investigated and remains controversial. Here, we report the pairing symmetry of the infinite-layer nickelates determined by London penetration depth measurements in neodymium-based (Nd,Sr)NiO$_2$ and lanthanide-based (La,Ca)NiO$_2$ thin films of high crystallinity. A rare-earth-specific order parameter is observed. While the lanthanide nickelates follow dirty line-node behaviour, the neodymium-counterpart exhibits nodeless order parameters such as the $(d+is)$ wave. In contrast to the cuprates, our results suggest that the superconducting order parameter in nickelates is beyond a single $d_(x^2-y^2 )$-wave gap. Furthermore, the superfluid density shows a long tail near the superconducting transition temperature which is consistent with the emergence of a two-dimensional to three-dimensional crossover in the superconducting state. These observations challenge the early theoretical framework and propel further experimental and theoretical interests in the pairing nature of the infinite-layer nickelate family.


title: p-orbital disclination states in non-Euclidean geometries

abstract: Disclinations are ubiquitous lattice defects existing in almost all crystalline materials. In two-dimensional nanomaterials, disclinations lead to the warping and deformation of the hosting material, yielding non-Euclidean geometries. However, such geometries have never been investigated experimentally in the context of topological phenomena. Here, by creating the physical realization of disclinations in conical and saddle-shaped acoustic systems, we demonstrate that disclinations can lead to topologically protected bound modes in non-Euclidean surfaces. In the designed honeycomb sonic crystal for p-orbital acoustic waves, non-Euclidean geometry interplay with the p-orbital physics and the band topology, showing intriguing emergent features as confirmed by consistent experiments and simulations. Our study opens a pathway towards topological phenomena in non-Euclidean geometries that may inspire future studies on, e.g., electrons and phonons in nanomaterials with curved surfaces.


title: IsoDAR@Yemilab: A Report on the Technology, Capabilities, and Deployment

abstract: IsoDAR@Yemilab is a novel isotope-decay-at-rest experiment that has preliminary approval to run at the Yemi underground laboratory (Yemilab) in Jeongseon-gun, South Korea. In this technical report, we describe in detail the considerations for installing this compact particle accelerator and neutrino target system at the Yemilab underground facility. Specifically, we describe the caverns being prepared for IsoDAR, and address installation, hielding, and utilities requirements. To give context and for completeness, we also briefly describe the physics opportunities of the IsoDAR neutrino source when paired with the Liquid Scintillator Counter (LSC) at Yemilab, and review the technical design of the neutrino source.


title: Migration of self-propelling agent in a turbulent environment with minimal energy consumption

abstract: We present a numerical study of training a self-propelling agent to migrate in the unsteady flow environment. We control the agent to utilize the background flow structure by adopting the reinforcement learning algorithm to minimize energy consumption. We considered the agent migrating in two types of flows: one is simple periodical double-gyre flow as a proof-of-concept example, while the other is complex turbulent Rayleigh-B\'enard convection as a paradigm for migrating in the convective atmosphere or the ocean. The results show that the smart agent in both flows can learn to migrate from one position to another while utilizing background flow currents as much as possible to minimize the energy consumption, which is evident by comparing the smart agent with a naive agent that moves straight from the origin to the destination. In addition, we found that compared to the double-gyre flow, the flow field in the turbulent Rayleigh-B\'enard convection exhibits more substantial fluctuations, and the training agent is more likely to explore different migration strategies; thus, the training process is more difficult to converge. Nevertheless, we can still identify an energy-efficient trajectory that corresponds to the strategy with the highest reward received by the agent. These results have important implications for many migration problems such as unmanned aerial vehicles flying in a turbulent convective environment, where planning energy-efficient trajectories are often involved.


title: Performance Analysis of Multiple-Antenna Ambient Backscatter Systems at Finite Blocklengths

abstract: This paper analyzes the maximal achievable rate for a given blocklength and error probability over a multiple-antenna ambient backscatter channel with perfect channel state information at the receiver. The result consists of a finite blocklength channel coding achievability bound and a converse bound based on the Neyman-Pearson test and the normal approximation based on the Berry- Esseen Theorem. Numerical evaluation of these bounds shows fast convergence to the channel capacity as the blocklength increases and also proves that the channel dispersion is an accurate measure of the backoff from capacity due to finite blocklength.


title: NAPA: Neighborhood-Assisted and Posterior-Adjusted Two-sample Inference

abstract: Two-sample multiple testing problems of sparse spatial data are frequently arising in a variety of scientific applications. In this article, we develop a novel neighborhood-assisted and posterior-adjusted (NAPA) approach to incorporate both the spatial smoothness and sparsity type side information to improve the power of the test while controlling the false discovery of multiple testing. We translate the side information into a set of weights to adjust the $p$-values, where the spatial pattern is encoded by the ordering of the locations, and the sparsity structure is encoded by a set of auxiliary covariates. We establish the theoretical properties of the proposed test, including the guaranteed power improvement over some state-of-the-art alternative tests, and the asymptotic false discovery control. We demonstrate the efficacy of the test through intensive simulations and two neuroimaging applications.


title: ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- VIII. A search for hot cores by using C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH lines

abstract: Hot cores characterized by rich lines of complex organic molecules are considered as ideal sites for investigating the physical and chemical environments of massive star formation. We present a search for hot cores by using typical nitrogen- and oxygen-bearing complex organic molecules (C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH), based on ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS). The angular resolutions and line sensitivities of the ALMA observations are better than 2 arcsec and 10 mJy/beam, respectively. A total of 60 hot cores are identified with 45 being newly detected, in which the complex organic molecules have high gas temperatures ($>$ 100 K) and small source sizes ($<$ 0.1 pc). So far this is the largest sample of hot cores observed with similar angular resolution and spectral coverage. The observations have also shown nitrogen and oxygen differentiation in both line emission and gas distribution in 29 hot cores. Column densities of CH$_3$OH and CH$_3$OCHO increase as rotation temperatures rise. The column density of CH$_3$OCHO correlates tightly with that of CH$_3$OH. The pathways for production of different species are discussed. Based on the spatial position difference between hot cores and UC~H{\sc ii} regions, we conclude that 24 hot cores are externally heated while the other hot cores are internally heated. The observations presented here will potentially help establish a hot core template for studying massive star formation and astrochemistry.


title: Thin-gap averaging of variable-viscosity flows: application to thermoviscous fingering

abstract: A consistent averaging technique, using the weighted residual integral boundary layer (WRIBL) method, is presented for flow through a thin-gap geometry wherein the fluid's viscosity varies across the gap. In such situations, the flow has a non-parabolic cross-gap velocity profile -- an effect that is ignored by Darcy models conventionally used for such Hele-Shaw flows. The WRIBL technique systematically accounts for the cross-gap variation of viscosity and yields reduced-order equations for the gap-averaged fluid flow rate. As a test case, we consider a fluid with a temperature-dependent viscosity and analyse the previously-studied problem of thermoviscous fingering: a hot fluid flowing through a Hele-Shaw geometry with cold walls spontaneously forms channels of low-viscosity, hot fluid, separated by regions of high-viscosity, cold fluid. The temperature of the cold walls is assumed either to be constant, a scenario that mimics the upward flow of magma through fissures in the Earth's crust, or to vary linearly along the direction of flow. In both cases, the predictions of the WRIBL model, regarding the multiplicity of uniform steady flow states and their linear stability, are compared with that of the hitherto used, ad hoc, Darcy model (Helfrich, J. Fluid Mech., vol. 305, 1995, pp. 219-238), as well as with calculations of the full three-dimensional governing equations (Wylie and Lister, J. Fluid Mech., vol. 305, 1995, pp. 239-261). Though the results are qualitatively similar, the WRIBL model is found to be much more accurate than the Darcy model. The averaging method is presented in a general manner to facilitate its application to other physical situations where, for example, the viscosity depends on solute/particle concentration in a solution/suspension.


title: On evolving natural curvature for an inextensible, unshearable, viscoelastic rod

abstract: We formulate and consider the problem of an inextensible, unshearable, viscoelastic rod, with evolving natural configuration, moving on a plane. We prove that the dynamic equations describing quasistatic motion of an Eulerian strut, an infinite dimensional dynamical system, are globally well-posed. For every value of the terminal thrust, these equations contain a smooth embedded curve of static solutions (equilibrium points). We characterize the spectrum of the linearized equations about an arbitrary equilibrium point, and using this information and a convergence result for dynamical systems due to Brunovsk\'y and Pol\'acik, we prove that every solution to the quasistatic equations of motion converges to an equilibrium point as time goes to infinity.


title: Are Commercial Face Detection Models as Biased as Academic Models?

abstract: As facial recognition systems are deployed more widely, scholars and activists have studied their biases and harms. Audits are commonly used to accomplish this and compare the algorithmic facial recognition systems' performance against datasets with various metadata labels about the subjects of the images. Seminal works have found discrepancies in performance by gender expression, age, perceived race, skin type, etc. These studies and audits often examine algorithms which fall into two categories: academic models or commercial models. We present a detailed comparison between academic and commercial face detection systems, specifically examining robustness to noise. We find that state-of-the-art academic face detection models exhibit demographic disparities in their noise robustness, specifically by having statistically significant decreased performance on older individuals and those who present their gender in a masculine manner. When we compare the size of these disparities to that of commercial models, we conclude that commercial models - in contrast to their relatively larger development budget and industry-level fairness commitments - are always as biased or more biased than an academic model.


title: A First Look at First-Passage Processes

abstract: These notes are based on the lectures that I gave (virtually) at the Bruneck Summer School in 2021 on first-passage processes and some applications of the basic theory. I begin by defining what is a first-passage process and presenting the connection between the first-passage probability and the familiar occupation probability. Some basic features of first passage on the semi-infinite line and a finite interval are then discussed, such as splitting probabilities and first-passage times. I also treat the fundamental connection between first passage and electrostatics. A number of applications of first-passage processes are then presented, including the hitting probability for a sphere in greater than two dimensions, reaction rate theory and its extension to receptors on a cell surface, first-passage inside an infinite absorbing wedge in two dimensions, stochastic hunting processes in one dimension, the survival of a diffusing particle in an expanding interval, and finally the dynamics of the classic birth-death process.


title: An Efficient Algorithm for the Partitioning Min-Max Weighted Matching Problem

abstract: The Partitioning Min-Max Weighted Matching (PMMWM) problem is an NP-hard problem that combines the problem of partitioning a group of vertices of a bipartite graph into disjoint subsets with limited size and the classical Min-Max Weighted Matching (MMWM) problem. Kress et al. proposed this problem in 2015 and they also provided several algorithms, among which MP$_{\text{LS}}$ is the state-of-the-art. In this work, we observe there is a time bottleneck in the matching phase of MP$_{\text{LS}}$. Hence, we optimize the redundant operations during the matching iterations, and propose an efficient algorithm called the MP$_{\text{KM-M}}$ that greatly speeds up MP$_{\text{LS}}$. The bottleneck time complexity is optimized from $O(n^3)$ to $O(n^2)$. We also prove the correctness of MP$_{\text{KM-M}}$ by the primal-dual method. To test the performance on diverse instances, we generate various types and sizes of benchmarks, and carried out an extensive computational study on the performance of MP$_{\text{KM-M}}$ and MP$_{\text{LS}}$. The evaluation results show that our MP$_{\text{KM-M}}$ greatly shortens the runtime as compared with MP$_{\text{LS}}$ while yielding the same solution quality.


title: Positive solutions to multi-critical elliptic problems

abstract: In this paper, we investigate the existence of multiple solutions to the following multi-critical elliptic problem \begin{equation}\label{eq:0.1} \left\{\begin{aligned} -\Delta u & =\lambda |u|^{p-2}u +\sum_{i=1}^k(|x|^{-(N-\alpha_i)}*|u|^{2^*_i})|u|^{2^*_i-2}u\quad {\rm in}\quad \Omega,\\ &u\in H^1_0(\Omega)\\ \end{aligned}\right. \end{equation} in connection with the topology of the bounded domain $\Omega\subset \mathbb{R}^N, \,N\geq 4$, where $\lambda>0$, $2^*_i=\frac{N+\alpha_i}{N-2}$ with $N-4<\alpha_i<N,\ \ i=1,2,\cdot\cdot\cdot, k$ are critical Hardy-Littlewood-Sobolev exponents and $2<p<22^*_{min}$ with $2^*_{min}=\min\{2^*_i, \ i=1,2,\cdot\cdot\cdot, k\}$. We show that there is $\lambda^*>0$ such that if $0<\lambda<\lambda^*$ problem \eqref{eq:0.1} possesses at least $cat_\Omega(\Omega)$ positive solutions. We also study the existence and uniqueness of solutions for the limit problem of \eqref{eq:0.1}.


title: MonLAD: Money Laundering Agents Detection in Transaction Streams

abstract: Given a stream of money transactions between accounts in a bank, how can we accurately detect money laundering agent accounts and suspected behaviors in real-time? Money laundering agents try to hide the origin of illegally obtained money by dispersive multiple small transactions and evade detection by smart strategies. Therefore, it is challenging to accurately catch such fraudsters in an unsupervised manner. Existing approaches do not consider the characteristics of those agent accounts and are not suitable to the streaming settings. Therefore, we propose MonLAD and MonLAD-W to detect money laundering agent accounts in a transaction stream by keeping track of their residuals and other features; we devise AnoScore algorithm to find anomalies based on the robust measure of statistical deviation. Experimental results show that MonLAD outperforms the state-of-the-art baselines on real-world data and finds various suspicious behavior patterns of money laundering. Additionally, several detected suspected accounts have been manually-verified as agents in real money laundering scenario.


title: Direct measurement of non-thermal electron acceleration from magnetically driven reconnection in a laboratory plasma

abstract: Magnetic reconnection is a ubiquitous astrophysical process that rapidly converts magnetic energy into some combination of plasma flow energy, thermal energy, and non-thermal energetic particles, including energetic electrons. Various reconnection acceleration mechanisms in different low-$\beta$ (plasma-to-magnetic pressure ratio) and collisionless environments have been proposed theoretically and studied numerically, including first- and second-order Fermi acceleration, betatron acceleration, parallel electric field acceleration along magnetic fields, and direct acceleration by the reconnection electric field. However, none of them have been heretofore confirmed experimentally, as the direct observation of non-thermal particle acceleration in laboratory experiments has been difficult due to short Debye lengths for \textit{in-situ} measurements and short mean free paths for \textit{ex-situ} measurements. Here we report the direct measurement of accelerated non-thermal electrons from low-$\beta$ magnetically driven reconnection in experiments using a laser-powered capacitor coil platform. We use kiloJoule lasers to drive parallel currents to reconnect MegaGauss-level magnetic fields in a quasi-axisymmetric geometry. The angular dependence of the measured electron energy spectrum and the resulting accelerated energies, supported by particle-in-cell simulations, indicate that the mechanism of direct electric field acceleration by the out-of-plane reconnection electric field is at work. Scaled energies using this mechanism show direct relevance to astrophysical observations. Our results therefore validate one of the proposed acceleration mechanisms by reconnection, and establish a new approach to study reconnection particle acceleration with laboratory experiments in relevant regimes.

categories: physics.plasm-ph astro-ph.HE astro-ph.IM astro-ph.SR physics.space-ph


title: Learning Resource Allocation Policies from Observational Data with an Application to Homeless Services Delivery

abstract: We study the problem of learning, from observational data, fair and interpretable policies that effectively match heterogeneous individuals to scarce resources of different types. We model this problem as a multi-class multi-server queuing system where both individuals and resources arrive stochastically over time. Each individual, upon arrival, is assigned to a queue where they wait to be matched to a resource. The resources are assigned in a first come first served (FCFS) fashion according to an eligibility structure that encodes the resource types that serve each queue. We propose a methodology based on techniques in modern causal inference to construct the individual queues as well as learn the matching outcomes and provide a mixed-integer optimization (MIO) formulation to optimize the eligibility structure. The MIO problem maximizes policy outcome subject to wait time and fairness constraints. It is very flexible, allowing for additional linear domain constraints. We conduct extensive analyses using synthetic and real-world data. In particular, we evaluate our framework using data from the U.S. Homeless Management Information System (HMIS). We obtain wait times as low as an FCFS policy while improving the rate of exit from homelessness for underserved or vulnerable groups (7% higher for the Black individuals and 15% higher for those below 17 years old) and overall.


title: Estimates in $L_{p}$ for solutions of SPDEs with coefficients in Morrey classes

abstract: For solutions of a certain class of SPDEs in divergence form we present some estimates of their $L_{p}$-norms and the $L_{p}$-norms of their first-order derivatives. The main novelty is that the low-order coefficients are supposed to belong to certain Morrey classes instead of $L_{p}$-spaces. Our results are new even if there are no stochastic terms in the equation.


title: Identifying a Training-Set Attack's Target Using Renormalized Influence Estimation

abstract: Targeted training-set attacks inject malicious instances into the training set to cause a trained model to mislabel one or more specific test instances. This work proposes the task of target identification, which determines whether a specific test instance is the target of a training-set attack. Target identification can be combined with adversarial-instance identification to find (and remove) the attack instances, mitigating the attack with minimal impact on other predictions. Rather than focusing on a single attack method or data modality, we build on influence estimation, which quantifies each training instance's contribution to a model's prediction. We show that existing influence estimators' poor practical performance often derives from their over-reliance on training instances and iterations with large losses. Our renormalized influence estimators fix this weakness; they far outperform the original estimators at identifying influential groups of training examples in both adversarial and non-adversarial settings, even finding up to 100% of adversarial training instances with no clean-data false positives. Target identification then simplifies to detecting test instances with anomalous influence values. We demonstrate our method's effectiveness on backdoor and poisoning attacks across various data domains, including text, vision, and speech, as well as against a gray-box, adaptive attacker that specifically optimizes the adversarial instances to evade our method. Our source code is available at https://github.com/ZaydH/target_identification.


title: Underwater Acoustic Communication Channel Modeling using Deep Learning

abstract: With the recent increase in the number of underwater activities, having effective underwater communication systems has become increasingly important. Underwater acoustic communication has been widely used but greatly impaired due to the complicated nature of the underwater environment. In a bid to better understand the underwater acoustic channel so as to help in the design and improvement of underwater communication systems, attempts have been made to model the underwater acoustic channel using mathematical equations and approximations under some assumptions. In this paper, we explore the capability of machine learning and deep learning methods to learn and accurately model the underwater acoustic channel using real underwater data collected from a water tank with disturbance and from lake Tahoe. Specifically, Deep Neural Network (DNN) and Long Short Term Memory (LSTM) are applied to model the underwater acoustic channel. Experimental results show that these models are able to model the underwater acoustic communication channel well and that deep learning models, especially LSTM are better models in terms of mean absolute percentage error.


title: On the Optimality of Linear Index Coding over the Fields with Characteristic Three

abstract: It has been known that the insufficiency of linear coding in achieving the optimal rate of the general index coding problem is rooted in its rate's dependency on the field size. However, this dependency has been described only through the two well-known matroid instances, namely the Fano and non-Fano matroids, which, in turn, limits its scope only to the fields with characteristic two. In this paper, we extend this scope to demonstrate the reliance of linear index coding rate on fields with characteristic three. By constructing two index coding instances of size 29, we prove that for the first instance, linear coding is optimal only over the fields with characteristic three, and for the second instance, linear coding over any field with characteristic three can never be optimal. Then, a variation of the second instance is designed as the third index coding instance of size 58. For this instance, it is proved that while linear coding over any field with characteristic three cannot be optimal, there exists a nonlinear code over the fields with characteristic three, which achieves its optimal rate. Connecting the first and third index coding instances in two specific ways, called no-way and two-way connections, will lead to two new index coding instances of size 87 and 91, for which linear coding is outperformed by nonlinear codes. Another main contribution of this paper is the reduction of the key constraints on the space of the linear coding for the first and second index coding instances, each of size 29, into a matroid instance with the ground set of size 9, whose linear representability is dependent on the fields with characteristic three. The proofs and discussions provided in this paper through using these two relatively small matroid instances will shed light on the underlying reason causing the linear coding to become insufficient for the general index coding problem.


title: GENGA II: GPU planetary N-body simulations with non-Newtonian forces and high number of particles

abstract: We present recent updates and improvements of the graphical processing unit (GPU) N-body code GENGA. Modern state-of-the-art simulations of planet formation require the use of a very high number of particles to accurately resolve planetary growth and to quantify the effect of dynamical friction. At present the practical upper limit is in the range of 30,000 - 60,000 fully interactive particles; possibly a little more on the latest GPU devices. While the original hybrid symplectic integration method has difficulties to scale up to these numbers, we have improved the integration method by i) introducing higher level changeover functions and ii) code improvements to better use the most recent GPU hardware efficiently for such large simulations. We added treatments of non-Newtonian forces such as general relativity, tidal interaction, rotational deformation, the Yarkovsky effect, and Poynting-Robertson drag, as well as a new model to treat virtual collisions of small bodies in the solar system. We added new tools to GENGA, such as semi-active test particles that feel more massive bodies but not each other, a more accurate collision handling and a real-time openGL visualization. We present example simulations, including a 1.5 billion year terrestrial planet formation simulation that initially started with 65,536 particles, a 3.5 billion year simulation without gas giants starting with 32,768 particles, the evolution of asteroid fragments in the solar system, and the planetesimal accretion of a growing Jupiter simulation. GENGA runs on modern NVIDIA and AMD GPUs.


title: Stability of Schr\"odinger Potentials and Convergence of Sinkhorn's Algorithm

abstract: We study the stability of entropically regularized optimal transport with respect to the marginals. Given marginals converging weakly, we establish a strong convergence for the Schr\"odinger potentials describing the density of the optimal couplings. When the marginals converge in total variation, the optimal couplings also converge in total variation. This is applied to show that Sinkhorn's algorithm converges in total variation when costs are quadratic and marginals are subgaussian, or more generally for all continuous costs satisfying an integrability condition.


title: ViT-HGR: Vision Transformer-based Hand Gesture Recognition from High Density Surface EMG Signals

abstract: Recently, there has been a surge of significant interest on application of Deep Learning (DL) models to autonomously perform hand gesture recognition using surface Electromyogram (sEMG) signals. DL models are, however, mainly designed to be applied on sparse sEMG signals. Furthermore, due to their complex structure, typically, we are faced with memory constraints; require large training times and a large number of training samples, and; there is the need to resort to data augmentation and/or transfer learning. In this paper, for the first time (to the best of our knowledge), we investigate and design a Vision Transformer (ViT) based architecture to perform hand gesture recognition from High Density (HD-sEMG) signals. Intuitively speaking, we capitalize on the recent breakthrough role of the transformer architecture in tackling different complex problems together with its potential for employing more input parallelization via its attention mechanism. The proposed Vision Transformer-based Hand Gesture Recognition (ViT-HGR) framework can overcome the aforementioned training time problems and can accurately classify a large number of hand gestures from scratch without any need for data augmentation and/or transfer learning. The efficiency of the proposed ViT-HGR framework is evaluated using a recently-released HD-sEMG dataset consisting of 65 isometric hand gestures. Our experiments with 64-sample (31.25 ms) window size yield average test accuracy of 84.62 +/- 3.07%, where only 78, 210 number of parameters is utilized. The compact structure of the proposed ViT-based ViT-HGR framework (i.e., having significantly reduced number of trainable parameters) shows great potentials for its practical application for prosthetic control.


title: Negative-ResNet: Noisy Ambulatory Electrocardiogram Signal Classification Scheme

abstract: With recently successful applications of deep learning in computer vision and general signal processing, deep learning has shown many unique advantages in medical signal processing. However, data labelling quality has become one of the most significant issues for AI applications, especially when it requires domain knowledge (e.g. medical image labelling). In addition, there might be noisy labels in practical datasets, which might impair the training process of neural networks. In this work, we propose a semi-supervised algorithm for training data samples with noisy labels by performing selected Positive Learning (PL) and Negative Learning (NL). To verify the effectiveness of the proposed scheme, we designed a portable ECG patch -- iRealCare -- and applied the algorithm on a real-life dataset. Our experimental results show that we can achieve an accuracy of 91.0 %, which is 6.2 % higher than a normal training process with ResNet. There are 65 patients in our dataset and we randomly picked 2 patients to perform validation.


title: A sine transform based preconditioned MINRES method for all-at-once systems from constant and variable-coefficient evolutionary PDEs

abstract: In this work, we propose a simple yet generic preconditioned Krylov subspace method for a large class of nonsymmetric block Toeplitz all-at-once systems arising from discretizing evolutionary partial differential equations. Namely, our main result is to propose two novel symmetric positive definite preconditioners, which can be efficiently diagonalized by the discrete sine transform matrix. More specifically, our approach is to first permute the original linear system to obtain a symmetric one, and subsequently develop desired preconditioners based on the spectral symbol of the modified matrix. Then, we show that the eigenvalues of the preconditioned matrix sequences are clustered around $\pm 1$, which entails rapid convergence when the minimal residual method is devised. Alternatively, when the conjugate gradient method on the normal equations is used, we show that our preconditioner is effective in the sense that the eigenvalues of the preconditioned matrix sequence are clustered around unity. An extension of our proposed preconditioned method is given for high-order backward difference time discretization schemes, which can be applied on a wide range of time-dependent equations. Numerical examples are given, also in the variable-coefficient setting, to demonstrate the effectiveness of our proposed preconditioners, which consistently outperforms an existing block circulant preconditioner discussed in the relevant literature.


title: Varying Coefficient Model via Adaptive Spline Fitting

abstract: The varying coefficient model has received broad attention from researchers as it is a powerful dimension reduction tool for non-parametric modeling. Most existing varying coefficient models fitted with polynomial spline assume equidistant knots and take the number of knots as the hyperparameter. However, imposing equidistant knots appears to be too rigid, and determining the optimal number of knots systematically is also a challenge. In this article, we deal with this challenge by utilizing polynomial splines with adaptively selected and predictor-specific knots to fit the coefficients in varying coefficient models. An efficient dynamic programming algorithm is proposed to find the optimal solution. Numerical results show that the new method can achieve significantly smaller mean squared errors for coefficients compared with the equidistant spline fitting method.


title: Accounting for the Fraction of Carcasses outside the Searched Area and the Estimation of Bird and Bat Fatalities at Wind Energy Facilities

abstract: In estimating bird and bat mortality at wind turbines, it is essential to account for carcasses that lie outside the searched area. In this manuscript we explore some of the difficulties and nuances involved in the spatial prediction of the number of carcasses that lie outside the searched area and provide extensive guidance and documentation for a new R package (dwp) for modeling carcass density as a function of distance from a turbine and estimation of the fraction of carcasses lying within the searched area, which is a critical parameter in fatality estimation software such as GenEst, eoa, acmeR, and carcass.


title: Experimental Observation of Maximum Density Fluctuation in Liquid Te

abstract: We performed small-angle X-ray scattering measurements of liquid Te using a synchrotron radiation facility and observed the maximum scattering intensity near 620 K in the supercooled region (melting temperature 723 K). This result is an experimental observation of the ridge line of the critical density fluctuation associated with the liquid-liquid phase transition that is argued to exist in the supercooled region and verifies the existence of the transition. Similar results have been reported for supercooled liquid water, demonstrating the universality of the liquid-liquid transition concept (or inhomogeneous model) in explaining the thermodynamics of these "anomalous liquids". Moreover, as a more essential issue, it has been demonstrated that liquids can be inhomogeneous at mesoscopic scales, even in their thermal equilibrium state. This will contribute to a greater understanding of the physics of liquids, including glass-forming liquids.


title: Documenting Geographically and Contextually Diverse Data Sources: The BigScience Catalogue of Language Data and Resources

abstract: In recent years, large-scale data collection efforts have prioritized the amount of data collected in order to improve the modeling capabilities of large language models. This prioritization, however, has resulted in concerns with respect to the rights of data subjects represented in data collections, particularly when considering the difficulty in interrogating these collections due to insufficient documentation and tools for analysis. Mindful of these pitfalls, we present our methodology for a documentation-first, human-centered data collection project as part of the BigScience initiative. We identified a geographically diverse set of target language groups (Arabic, Basque, Chinese, Catalan, English, French, Indic languages, Indonesian, Niger-Congo languages, Portuguese, Spanish, and Vietnamese, as well as programming languages) for which to collect metadata on potential data sources. To structure this effort, we developed our online catalogue as a supporting tool for gathering metadata through organized public hackathons. We present our development process; analyses of the resulting resource metadata, including distributions over languages, regions, and resource types; and our lessons learned in this endeavor.


title: Statistical analysis of intermittency and its association with proton heating in the near Sun environment

abstract: We use data from the first six encounters of Parker Solar Probe and employ the Partial Variance of Increments ($PVI$) method to study the statistical properties of coherent structures in the inner heliosphere with the aim of exploring physical connections between magnetic field intermittency and observable consequences such as plasma heating and turbulence dissipation. Our results support proton heating localized in the vicinity of, and strongly correlated with, magnetic structures characterized by $PVI \geq 1$. We show that on average, such events constitute $\approx 19\%$ of the dataset, though variations may occur depending on the plasma parameters. We show that the waiting time distribution ($WT$) of identified events is consistent across all six encounters following a power-law scaling at lower $WTs$. This result indicates that coherent structures are not evenly distributed in the solar wind but rather tend to be tightly correlated and form clusters. We observe that the strongest magnetic discontinuities, $PVI \geq 6$, usually associated with reconnection exhausts, are sites where magnetic energy is locally dissipated in proton heating and are associated with the most abrupt changes in proton temperature. However, due to the scarcity of such events, their relative contribution to energy dissipation is minor. Taking clustering effects into consideration, we show that smaller scale, more frequent structures with PVI between, $1\lesssim PVI \lesssim 6$, play the major role in magnetic energy dissipation. The number density of such events is strongly associated with the global solar wind temperature, with denser intervals being associated with higher $T_{p}$.


title: Continuous-frequency weak electric field measurement with Rydberg atoms

abstract: We demonstrate a continuous frequency electric field measurement based on the far off-resonant AC stark effect in a Rydberg atomic vapor cell. In this configuration, a strong far off-resonant field, denoted as a local oscillator (LO) field, acts as a gain shifting the Rydberg level to a high sensitivity region. An incident weak signal field with a few hundreds of kHz difference from the LO field is mixed with the LO field in Rydberg system to generate an intermediate frequency (IF) signal, which is read out by the Rydberg electromagnetically induced transparency (Rydberg-EIT) spectroscopy. Not like resonant EIT-AT spectra, we realize the electric field measurement of the signal frequency from 2 GHz to 5 GHz using a single Rydberg state. A minimum detectable filed strength is down to 2.31~$\mu$V/cm and a linear dynamic range is over 65~dB. The minimum detectable filed is comparable with a resonant microwave-dressed Rydberg heterodyne receiver using the same system, which is 1.45~$\mu$V/cm. We also show the system has an inherent polarization selectivity feature. Our method can provide a high sensitivity of electric field measurement and be extended to arbitrary frequency measurements.


title: Link Prediction with Contextualized Self-Supervision

abstract: Link prediction aims to infer the link existence between pairs of nodes in networks/graphs. Despite their wide application, the success of traditional link prediction algorithms is hindered by three major challenges -- link sparsity, node attribute noise and dynamic changes -- that are faced by many real-world networks. To address these challenges, we propose a Contextualized Self-Supervised Learning (CSSL) framework that fully exploits structural context prediction for link prediction. The proposed CSSL framework learns a link encoder to infer the link existence probability from paired node embeddings, which are constructed via a transformation on node attributes. To generate informative node embeddings for link prediction, structural context prediction is leveraged as a self-supervised learning task to boost the link prediction performance. Two types of structural context are investigated, i.e., context nodes collected from random walks vs. context subgraphs. The CSSL framework can be trained in an end-to-end manner, with the learning of model parameters supervised by both the link prediction and self-supervised learning tasks. The proposed CSSL is a generic and flexible framework in the sense that it can handle both attributed and non-attributed networks, and operate under both transductive and inductive link prediction settings. Extensive experiments and ablation studies on seven real-world benchmark networks demonstrate the superior performance of the proposed self-supervision based link prediction algorithm over state-of-the-art baselines, on different types of networks under both transductive and inductive settings. The proposed CSSL also yields competitive performance in terms of its robustness to node attribute noise and scalability over large-scale networks.


title: MOORe: Model-based Offline-to-Online Reinforcement Learning

abstract: With the success of offline reinforcement learning (RL), offline trained RL policies have the potential to be further improved when deployed online. A smooth transfer of the policy matters in safe real-world deployment. Besides, fast adaptation of the policy plays a vital role in practical online performance improvement. To tackle these challenges, we propose a simple yet efficient algorithm, Model-based Offline-to-Online Reinforcement learning (MOORe), which employs a prioritized sampling scheme that can dynamically adjust the offline and online data for smooth and efficient online adaptation of the policy. We provide a theoretical foundation for our algorithms design. Experiment results on the D4RL benchmark show that our algorithm smoothly transfers from offline to online stages while enabling sample-efficient online adaption, and also significantly outperforms existing methods.


title: Inflation-induced aneurysm formation and evolution in graded cylindrical tubes of arbitrary thickness

abstract: We study the initiation and evolution of aneurysmal morphology in a pressurized soft tube where the elastic modulus is non-uniform in the radial direction. The primary deformation prior to instability is characterized within the framework of nonlinear elasticity for a general material constitution and a generic modulus gradient. To unravel the influence of modulus gradient on aneurysm formation, we employ the incompressible Gent model and select three representative modulus gradients, including a linear, an exponential, and a sinusoidal function. In particular, the sinusoidal distribution can be used to model actual artery structure. In addition, two prototypical loading conditions are considered, namely, either the resultant axial force or the axial length can be fixed. Based on an explicit bifurcation condition in terms of the internal pressure and the resultant axial force for aneurysm formation or localized bulging, an exhaustive theoretical analysis on bulge initiation is carried out and the effect of geometric and material parameters and modulus gradient on the critical stretch generating localized bulging is revealed. It turns out that the modulus mismatch, as well as the position of maximum modulus, can dramatically affect the onset of localized bulging. Then we analytically elucidate the influence of modulus gradient on bulge propagation and conduct a finite element analysis of bulge evolution based on a robust finite element model established in Abaqus by UHYPER subroutine coding. Interestingly, it is found that a sinusoidally distributed modulus has negligible influence on the critical stretch of bulge initiation, the deformation process of bugle growth, and the maximum size of a bulge. The current analysis can provide useful insight into the biological evolution of human artery and into localized instabilities in graded structures.


title: A Regularity Theory for Static Schr\"odinger Equations on $\mathbb{R}^d$ in Spectral Barron Spaces

abstract: Spectral Barron spaces have received considerable interest recently as it is the natural function space for approximation theory of two-layer neural networks with a dimension-free convergence rate. In this paper we study the regularity of solutions to the whole-space static Schr\"odinger equation in spectral Barron spaces. We prove that if the source of the equation lies in the spectral Barron space $\mathcal{B}^s(\mathbb{R}^d)$ and the potential function admitting a non-negative lower bound decomposes as a positive constant plus a function in $\mathcal{B}^s(\mathbb{R}^d)$, then the solution lies in the spectral Barron space $\mathcal{B}^{s+2}(\mathbb{R}^d)$.