title: Total consensus under high reproductive-variance conditions

abstract: Star-shaped branching patterns of genealogies are common in marine species. High-fecundity marine populations are characterized by low ratios of effective to actual population size, which reflect high variance in reproductive success among parents in mass spawns. When extreme reproduction events occur, offspring from very few parents dominate the population (whereby multiple mergers, or subsets of lineages with star-like trees, arise) and thus, the loss of genetic diversity is significant. Under high reproductive-variance conditions (assuming that reproduction occurs by sampling from the Pareto distribution), this paper explores the distribution of heterozygosity across generations. The result shows that zero heterozygosity is not achieved, implying that the populations may decline without evident loss of genetic variation. It is also found that there are singularities in the heterozygosity distributions. However, in the case of high reproductive variance, the locations of the singular points subtly deviate from those of the case where reproduction occurs by Wright-Fisher sampling.


title: Observation and origin of non-thermal hard X-rays from Jupiter

abstract: Electrons accelerated on Earth by a rich variety of wave scattering or stochastic processes generate hard non-thermal X-ray bremsstrahlung up to >~ 1 MeV and power Earth's various types of aurorae. Although Jupiter's magnetic field is an order of magnitude larger than Earth's, space-based telescopes have previously detected X-rays only up to ~7 keV. On the basis of theoretical models of the Jovian auroral X-ray production, X-ray emission in the ~2-7 keV band has been interpreted as thermal (arising from electrons characterized by a Maxwell-Boltzmann distribution) bremsstrahlung. Here we report the observation of hard X-rays in the 8-20 keV band from the Jovian aurorae, obtained with the NuSTAR X-ray observatory. The X-rays fit to a flat power-law model with slope 0.60+/-0.22 - a spectral signature of non-thermal, hard X-ray bremsstrahlung. We determine the electron flux and spectral shape in the keV to MeV energy range using coeval in situ measurements by the Juno spacecraft's JADE and JEDI instruments. Jovian electron spectra of the form we observe have previously been interpreted to arise in stochastic acceleration, rather than coherent acceleration by electric fields. We reproduce the X-ray spectral shape and approximate flux observed by NuSTAR, and explain the non-detection of hard X-rays by Ulysses, by simulating the non-thermal population of electrons undergoing precipitating electron energy loss, secondary electron generation and bremsstrahlung emission in a model Jovian atmosphere. The results highlight the similarities between the processes generating hard X-ray auroras on Earth and Jupiter, which may be occurring on Saturn, too.


title: Investigating the Future Potential of an Upgraded ALMA to Image Planet Forming Disks at Sub-au Scales

abstract: In recent years, ALMA has been able to observe large-scale substructures within protoplanetary disks. Comparison with the predictions from models of planet-disk interaction has indicated that most of these disk substructures can be explained by the presence of planets with the mass of Neptune or larger at orbital radii of $\approx 5 - 100$ au. Better resolution is needed to observe structures closer to the star, where terrestrial planets are expected to form, as well as structures opened by planets with masses lower than Neptune. We investigate the capabilities of a possible extension to ALMA that would double the longest baseline lengths in the array to detect and resolve disk substructures opened by Earth-mass and Super Earth planets at orbital radii of $1-5$ au. By simulating observations of a family of disk models using this extended configuration in ALMA Bands 6 and 7, we show that an upgraded ALMA would detect gaps in disks formed by super-Earths as close as 1 au, as well as Earth-mass planets down to $2-3$ au from the young host stars in nearby star forming regions.


title: Approximating Output Probabilities of Shallow Quantum Circuits which are Geometrically-local in any Fixed Dimension

abstract: We present a classical algorithm that, for any $D$-dimensional geometrically-local, quantum circuit $C$ of polylogarithmic-depth, and any bit string $x \in {0,1}^n$, can compute the quantity $|<x|C|0^{\otimes n}>|^2$ to within any inverse-polynomial additive error in quasi-polynomial time, for any fixed dimension $D$. This is an extension of the result [CC21], which originally proved this result for $D = 3$. To see why this is interesting, note that, while the $D = 1$ case of this result follows from standard use of Matrix Product States, known for decades, the $D = 2$ case required novel and interesting techniques introduced in [BGM19]. Extending to the case $D = 3$ was even more laborious and required further new techniques introduced in [CC21]. Our work here shows that, while handling each new dimension has historically required a new insight, and fixed algorithmic primitive, based on known techniques for $D \leq 3$, we can now handle any fixed dimension $D > 3$. Our algorithm uses the Divide-and-Conquer framework of [CC21] to approximate the desired quantity via several instantiations of the same problem type, each involving $D$-dimensional circuits on about half the number of qubits as the original. This division step is then applied recursively, until the width of the recursively decomposed circuits in the $D^{th}$ dimension is so small that they can effectively be regarded as $(D-1)$-dimensional problems by absorbing the small width in the $D^{th}$ dimension into the qudit structure at the cost of a moderate increase in runtime. The main technical challenge lies in ensuring that the more involved portions of the recursive circuit decomposition and error analysis from [CC21] still hold in higher dimensions, which requires small modifications to the analysis in some places.


title: Constructing Generic Effective Field Theory for All Masses and Spins

abstract: We fully solve the long-standing problem of operator basis construction for fields with any masses and spins. Based on the on-shell method, we propose a novel method to systematically construct a complete set of lowest dimensional amplitude bases at any given dimension through semi-standard Young tableaus of Lorentz subgroup $SU(2)_r$ and global symmetry $U(N)$ ($N$ is the number of external legs), which can be directly mapped into physical operator bases. We first construct a complete set of monomial bases whose dimension is not the lowest and a redundant set of bases that always contains a complete set of amplitude bases with the lowest dimension. Then we decompose the bases of the redundant set into the monomial bases from low to high dimension and eliminate the linear correlation bases. Finally, the bases with the lowest dimension can be picked up. We also propose a matrix projection method to construct the massive amplitude bases involving identical particles. The operator bases of a generic massive effective field theory can be efficiently constructed by the computer programs. A complete set of four-vector operators at dimensions up to six is presented.


title: Flat tori with large Laplacian eigenvalues in dimensions up to eight

abstract: We consider the optimization problem of maximizing the $k$-th Laplacian eigenvalue, $\lambda_{k}$, over flat $d$-dimensional tori of fixed volume. For $k=1$, this problem is equivalent to the densest lattice sphere packing problem. For larger $k$, this is equivalent to the NP-hard problem of finding the $d$-dimensional (dual) lattice with longest $k$-th shortest lattice vector. As a result of extensive computations, for $d \leq 8$, we obtain a sequence of flat tori, $T_{k,d}$, each of volume one, such that the $k$-th Laplacian eigenvalue of $T_{k,d}$ is very large; for each (finite) $k$ the $k$-th eigenvalue exceeds the value in (the $k\to \infty$ asymptotic) Weyl's law by a factor between 1.54 and 2.01, depending on the dimension. Stationarity conditions are derived and numerically verified for $T_{k,d}$ and we describe the degeneration of the tori as $k \to \infty$.


title: Spatial decay of discretely self-similar solutions to the Navier-Stokes equations

abstract: Forward self-similar and discretely self-similar weak solutions of the Navier-Stokes equations are known to exist globally in time for large self-similar and discretely self-similar initial data and are known to be regular outside of a space-time paraboloid. In this paper, we establish spatial decay rates for such solutions which hold in the region of regularity provided the initial data has locally sub-critical regularity away from the origin. In particular, we (1) lower the H\"older regularity of the data required to obtain an optimal decay rate for the nonlinear part of the flow compared to the existing literature, (2) establish new decay rates without logarithmic corrections for some smooth data, (3) provide new decay rates for solutions with rough data, and, as an application of our decay rates, (4) provide new upper bounds on how rapidly potentially {non-unique}, scaling invariant local energy solutions can separate away from the origin.


title: Remarks on the use of objective probabilities in Bell-CHSH inequalities

abstract: The violation of Bell inequalities is often interpreted as showing that, if hidden variables exist, they must be contextual and non local. But they can also be explained questioning the probability space employed, or the validity of the Kolmogorov axioms. In this article we explore the additional constrains which can be deduced from two widely used objetive probability theories: frequentism and propensities. One of the strongest objections in the deduction of one version of Bell inequalities goes about the probability space, which assumes the existence of values for the output of the experiment in each run, while only two of the four values can be measured each time, making them counterfactual. It is shown that frequentism rejects the possibility of using counterfactual situations, while long-run propensities allow their use. In this case the introduction of locality and contextuality does not help to explain the violation, and an alternative explanation could point to a failure of the probability. Single case propensities were designed to associate probabilities to single events, but they need to be conditional to the whole universe, and do not have a clear link with the observed relative frequencies. It heavily limits their use.


title: Effect of minimal length uncertainty on neutrino oscillation

abstract: In this paper, we study the effect of the minimal length on neutrino oscillation in a static magnetic field. In the framework of the generalized uncertainty principle, we reformulate the Hamiltonian for a relativistic neutrino moving in a magnetic field oriented along the z-direction of Cartesian coordinates. Using the modified energy spectrum, we obtain the oscillation probability for different neutrino flavors. In addition, we obtain the energy differences for the neutrino-mass eigenstates. We find that the energy and energy difference depend on the minimal length parameter {\alpha}, and the energy difference becomes independent of {\alpha} when the magnetic field is not present. In addition, we find that the modified probability of oscillation differs from the usual probability of oscillation if a magnetic field is present. Using the current experimental result, we estimate the upper bound on the deformation parameter and the minimal length, and find that the upper bound on the minimal length scale is less than the electroweak scale. If the minimal length is at Planck scale, the minimal length formalism leads to the same result as a quantum theory of gravity with an $SU(2)_{L} \times U(1)$ effective invariant dimension-5 Lagrangian including neutrino and Higgs fields.


title: Skewness as a test of dark energy perturbations

abstract: We investigate the role played by dark energy perturbations in the skewness $S_3$ of large-scale matter distribution. We consider a two-fluid universe composed by matter and dark energy, with perturbations in both components, and we estimate numerically the skewness of the matter density field as a function of the dark energy parameters. We characterize today's $S_3$ value for quintessence and phantom dark energy cosmologies as well as its dependence on the matter density parameter $\Omega_{m0}$ and the dark energy sound speed $c^2_s$ with accurate numerical fitting. These fits can be used to test cosmology against future high quality data on large scale structure.


title: Negative result about the construction of genuinely entangled subspaces from unextendible product bases

abstract: Unextendible product bases (UPBs) provide a versatile tool with various applications across different areas of quantum information theory. Their comprehensive characterization is thus of great importance and has been a subject of vital interest for over two decades now. An open question asks about the existence of UPBs, which are genuinely unextendible, i.e., they are not extendible even with biproduct vectors. In other words, the problem is to verify whether there exist genuinely entangled subspaces (GESs), subspaces composed solely of genuinely multiparty entangled states, complementary to UPBs. We solve this problem in the negative for many sizes of UPBs in different multipartite scenarios. In particular, in the all-important case of equal local dimensions, we show that there are always forbidden cardinalities for such UPBs, including the minimal ones corresponding to GESs of the maximal dimensions.


title: On the Complexity of Some Variations of Sorting by Transpositions

abstract: One of the main challenges in Computational Biology is to find the evolutionary distance between two organisms. In the field of comparative genomics, one way to estimate such distance is to find a minimum cost sequence of rearrangements (large scale mutations) needed to transform one genome into another, which is called the rearrangement distance. In the past decades, these problems were studied considering many types of rearrangements (such as reversals, transpositions, transreversals, and revrevs) and considering the same weight for all rearrangements, or different weights depending on the types of rearrangements. The complexity of the problems involving reversals, transpositions, and both rearrangements is known, even though the hardness proof for the problem combining reversals and transpositions was recently given. In this paper, we enhance the knowledge for these problems by proving that models involving transpositions alongside reversals, transreversals, and revrevs are NP-hard, considering weights w1 for reversals and w2 for the other rearrangements such that w2/w1 is less than or equal to 1.5. In addition, we address a cost function related to the number of fragmentations caused by a rearrangement, proving that the problem of finding a minimum cost sorting sequence, considering the fragmentation cost function with some restrictions, is NP-hard for transpositions and the combination of reversals and transpositions.


title: Programmable Interface for Statistical & Simulation Models (PRISM): Towards Greater Accessibility of Clinical and Healthcare Decision Models

abstract: Background: Increasingly, decision-making in healthcare relies on computer models, be it clinical prediction models at point of care or decision-analytic models at the policymaking level. Given the important role models play in both contexts, their structure and implementation be rigorously scrutinized. The ability to interrogate input/output associations without facing barriers can improve quality assurance mechanisms while satisfying privacy/confidentiality concerns and facilitating the integration of models into decision-making. This paper reports on the development of Programmable Interface for Statistical & Simulation Models (PRISM), a cloud-based platform for model accessibility. Methods: PRISM emphasizes two main principles: 1) minimal specifications on the side of model developer to make the model fit for cloud hosting, and 2) making client access completely independent of the resource requirement and software dependencies of the model. The server architecture integrates a RESTful Application Programming Interface (API) infrastructure, JSON for data transfer, a routing layer for access management, container technology for management of computer resources and package dependencies, and the capacity for synchronous or asynchronous model calls. Results: We discuss the architecture, the minimal API standards that enable a universal language for access to such models, the underlying server infrastructure, and the standards used for data transfer. An instance of PRISM is available as a service via the Peer Models Network http://peermodelsnetwork.com. Through a series of case studies, we demonstrate how interrogating models becomes possible in standardized fashion, in a way that is irrespective of the specifics of any model. Conclusions: We have developed a publicly accessible platform and minimalist standards that facilitate model accessibility for both clinical and policy models.


title: Virtual Maps for Autonomous Exploration of Cluttered Underwater Environments

abstract: We consider the problem of autonomous mobile robot exploration in an unknown environment, taking into account a robot's coverage rate, map uncertainty, and state estimation uncertainty. This paper presents a novel exploration framework for underwater robots operating in cluttered environments, built upon simultaneous localization and mapping (SLAM) with imaging sonar. The proposed system comprises path generation, place recognition forecasting, belief propagation and utility evaluation using a virtual map, which estimates the uncertainty associated with map cells throughout a robot's workspace. We evaluate the performance of this framework in simulated experiments, showing that our algorithm maintains a high coverage rate during exploration while also maintaining low mapping and localization error. The real-world applicability of our framework is also demonstrated on an underwater remotely operated vehicle (ROV) exploring a harbor environment.


title: Vision Models Are More Robust And Fair When Pretrained On Uncurated Images Without Supervision

abstract: Discriminative self-supervised learning allows training models on any random group of internet images, and possibly recover salient information that helps differentiate between the images. Applied to ImageNet, this leads to object centric features that perform on par with supervised features on most object-centric downstream tasks. In this work, we question if using this ability, we can learn any salient and more representative information present in diverse unbounded set of images from across the globe. To do so, we train models on billions of random images without any data pre-processing or prior assumptions about what we want the model to learn. We scale our model size to dense 10 billion parameters to avoid underfitting on a large data size. We extensively study and validate our model performance on over 50 benchmarks including fairness, robustness to distribution shift, geographical diversity, fine grained recognition, image copy detection and many image classification datasets. The resulting model, not only captures well semantic information, it also captures information about artistic style and learns salient information such as geolocations and multilingual word embeddings based on visual content only. More importantly, we discover that such model is more robust, more fair, less harmful and less biased than supervised models or models trained on object centric datasets such as ImageNet.


title: Vectorization of a thread-parallel Jacobi singular value decomposition method

abstract: The eigenvalue decomposition (EVD) of (a batch of) Hermitian matrices of order two has a role in many numerical algorithms, of which the one-sided Jacobi method for the singular value decomposition (SVD) is the prime example. In this paper the batched EVD is vectorized, with a vector-friendly data layout and the AVX-512 SIMD instructions of Intel CPUs, alongside other key components of a real and a complex OpenMP-parallel Jacobi-type SVD method, inspired by the sequential xGESVJ routines from LAPACK. These vectorized building blocks should be portable to other platforms that support similar vector operations. Unconditional numerical reproducibility is guaranteed for the batched EVD, sequential or threaded, and for the column transformations, that are, like the scaled dot-products, presently sequential but can be threaded if nested parallelism is desired. No avoidable overflow of the results can occur with the proposed EVD or the whole SVD. The measured accuracy of the proposed EVD often surpasses that of the xLAEV2 routines from LAPACK. While the batched EVD outperforms the matching sequence of xLAEV2 calls, speedup of the parallel SVD is modest but can be improved and is already beneficial with enough threads. Regardless of their number, the proposed SVD method gives identical results, but of somewhat lower accuracy than xGESVJ.


title: Alliance Makes Difference? Maximizing Social Welfare in Cross-Silo Federated Learning

abstract: As one of the typical settings of Federated Learning (FL), cross-silo FL allows organizations to jointly train an optimal Machine Learning (ML) model. In this case, some organizations may try to obtain the global model without contributing their local training power, lowering the social welfare. In this paper, we model the interactions among organizations in cross-silo FL as a public goods game and theoretically prove that there exists a social dilemma where the maximum social welfare is not achieved in Nash equilibrium. To overcome this dilemma, we employ the Multi-player Multi-action Zero-Determinant (MMZD) strategy to maximize the social welfare. With the help of the MMZD, an individual organization can unilaterally control the social welfare without extra cost. Since the MMZD strategy can be adopted by all organizations, we further study the case of multiple organizations jointly adopting the MMZD strategy to form an MMZD Alliance (MMZDA). We prove that the MMZDA strategy can strengthen the control of the maximum social welfare. Experimental results validate that the MMZD strategy is effective in obtaining the maximum social welfare and the MMZDA can achieve a larger maximum value.


title: Learning-Enabled Robust Control with Noisy Measurements

abstract: We present a constructive approach to bounded $\ell_2$-gain adaptive control with noisy measurements for linear time-invariant scalar systems with uncertain parameters belonging to a finite set. The gain bound refers to the closed-loop system, including the learning procedure. The approach is based on forward dynamic programming to construct a finite-dimensional information state consisting of $\mathcal H_\infty$-observers paired with a recursively computed performance metric. We do not assume prior knowledge of a stabilizing controller.


title: Transition from localized to mean field behaviour of cascading failures in the fiber bundle model on complex networks

abstract: We study the failure process of fiber bundles on complex networks focusing on the effect of the degree of disorder of fibers' strength on the transition from localized to mean field behaviour. Starting from a regular square lattice we apply the Watts-Strogatz rewiring technique to introduce long range random connections in the load transmission network and analyze how the ultimate strength of the bundle and the statistics of the size of failure cascades change when the rewiring probability is gradually increased. Our calculations revealed that the degree of strength disorder of nodes of the network has a substantial effect on the localized to mean field transition. In particular, we show that the transition sets on at a finite value of the rewiring probability, which shifts to higher values as the degree of disorder is reduced. The transition is limited to a well defined range of disorder, so that there exists a threshold disorder of nodes' strength below which the randomization of the network structure does not provide any improvement neither of the overall load bearing capacity nor of the cascade tolerance of the system. At low strength disorder the fully random network is the most stable one, while at high disorder best cascade tolerance is obtained at a lower structural randomness. Based on the interplay of the network structure and strength disorder we construct an analytical argument which provides a reasonable description of the numerical findings.


title: Investigating the low-flux states in six Intermediate Polars

abstract: We present optical photometry of six intermediate polars that exhibit transitions to a low-flux state. For four of these systems, DW Cnc, V515 And, V1223 Sgr and RX J2133.7+5107, we are able to perform timing analysis in and out of the low states. We find that, for DW Cnc and V515 And, the dominant periodicities in the light curves change as the flux decreases, indicating a change in the sources' accretion properties as they transition to the low state. For V1223 Sgr we find that the variability is almost completely quenched at the lowest flux, but do not find evidence for a changing accretion geometry. For RX J2133.7+5107, the temporal properties do not change in the low state, but we do see a period of enhanced accretion that is coincident with increased variability on the beat frequency, which we do not associate with a change in the accretion mechanisms in the system.


title: Preference Learning in School Choice Problems

abstract: In school choice, students make decisions based on their expectations of particular schools' suitability, and the decision to gather information about schools is influenced by the acceptance odds determined by the mechanism in place. We study a school choice model where students can obtain information about their preferences by incurring a cost. We demonstrate greater homogeneity in rank-order reports and reduced information acquisition under the Deferred-Acceptance (DA) mechanism, resulting in an increased reliance on random tie-breaking and ultimately inefficient outcomes. Thus, it is critical for the DA mechanism to have easy access to school information in order to maintain its efficiency.


title: Plasmoid-fed prominence formation (PF$^2$) during flux rope eruption

abstract: We report a new, plasmoid-fed scenario for the formation of an eruptive prominence (PF$^2$), involving reconnection and condensation. We use grid-adaptive resistive two-and-a-half-dimensional magnetohydrodynamic (MHD) simulations in a chromosphere-to-corona setup to resolve this plasmoid-fed scenario. We study a pre-existing flux rope (FR) in the low corona that suddenly erupts due to catastrophe, which also drives a fast shock above the erupting FR. A current sheet (CS) forms underneath the erupting FR, with chromospheric matter squeezed into it. The plasmoid instability occurs and multiple magnetic islands appear in the CS once the Lundquist number reaches $\sim 3.5\times 10^{4}$. The remnant chromospheric matter in the CS is then transferred to the FR by these newly formed magnetic islands. The dense and cool mass transported by the islands accumulates in the bottom of the FR, thereby forming a prominence during the eruption phase. More coronal plasma continuously condenses into the prominence due to the thermal instability as the FR rises. Due to the fine structure brought in by the PF$^2$ process, the model naturally forms filament threads, aligned above the polarity inversion line. Synthetic views at our resolution of $15 \mathrm{km}$ show many details that may be verified in future high-resolution observations.


title: Posterior Predictive Propensity Scores and $p$-Values

abstract: \citet{Rosenbaum83ps} introduced the notion of the propensity score and discussed its central role in causal inference with observational studies. Their paper, however, caused a fundamental incoherence with an early paper by \citet{Rubin78}, which showed that the propensity score does not play any role in the Bayesian analysis of unconfounded observational studies if the priors on the propensity score and outcome models are independent. Despite the serious efforts made in the literature, it is generally difficult to reconcile these contradicting results. We offer a simple approach to incorporating the propensity score in Bayesian causal inference based on the posterior predictive $p$-value. To motivate a simple procedure, we focus on the model with the strong null hypothesis of no causal effects for any units whatsoever. Computationally, the proposed posterior predictive $p$-value equals the classic $p$-value based on the Fisher randomization test averaged over the posterior predictive distribution of the propensity score. Moreover, using the studentized doubly robust estimator as the test statistic, the proposed $p$-value inherits the doubly robust property and is also asymptotically valid for testing the weak null hypothesis of zero average causal effect. Perhaps surprisingly, this Bayesianly motivated $p$-value can have better frequentist's finite-sample performance than the frequentist's $p$-value based on the asymptotic approximation especially when the propensity scores can take extreme values.


title: Accurate Ground-based Astrometry of Naked-eye Stars: The United States Naval Observatory Bright-Star Astrometric Database

abstract: We present the United States Naval Observatory (USNO) Bright-Star Astrometric Database (UBAD), a current-epoch high-accuracy astrometric catalog. The catalog consists of 364 bright northern hemisphere stars, including all but five such stars with either $V < 3.5$ or with $I < 3.2$ and $V < 6$, as well as a large fraction of slightly fainter stars; 36 of the brightest catalog stars are not included in Gaia Early Data Release 3 (EDR3). Observations were conducted with the USNO, Flagstaff Station, Kaj Strand 61-inch Astrometric Reflector. Target stars were imaged through a small 12.5-magnitude neutral-density spot, while the remainder of the stars in the field of view were unattenuated. This allowed for unsaturated images of the bright target stars to be calibrated directly against much fainter reference stars from Gaia EDR3. The median position errors are 1.9 mas in both right ascension and declination at the catalog epoch of 2017.0, with 90% of catalog stars having errors less than 2.6 mas; systematic errors are 1 -- 3 mas. Combining UBAD observations with Hipparcos-2 positions yields proper motions with median errors of 0.045 and 0.049 mas year$^{-1}$ in right ascension and declination, respectively, with 90% of stars having errors less than 0.1 mas year$^{-1}$; systematic errors are about 0.1 mas year$^{-1}$. Single-frame accuracy for positions of the target stars is typically 5 -- 6 mas. Gaia EDR3 astrometry for these bright stars, which are heavily saturated in the Gaia observations, is validated over the magnitude range $2 \lesssim G \lesssim 6$.


title: CAREER: A Foundation Model for Labor Sequence Data

abstract: Labor economists regularly analyze employment data by fitting predictive models to small, carefully constructed longitudinal survey datasets. Although machine learning methods offer promise for such problems, these survey datasets are too small to take advantage of them. In recent years large datasets of online resumes have also become available, providing data about the career trajectories of millions of individuals. However, standard econometric models cannot take advantage of their scale or incorporate them into the analysis of survey data. To this end we develop CAREER, a foundation model for job sequences. CAREER is first fit to large, passively-collected resume data and then fine-tuned to smaller, better-curated datasets for economic inferences. We fit CAREER to a dataset of 24 million job sequences from resumes, and adjust it on small longitudinal survey datasets. We find that CAREER forms accurate predictions of job sequences, outperforming econometric baselines on three widely-used economics datasets. We further find that CAREER can be used to form good predictions of other downstream variables. For example, incorporating CAREER into a wage model provides better predictions than the econometric models currently in use.


title: The Quarks of Attention

abstract: Attention plays a fundamental role in both natural and artificial intelligence systems. In deep learning, attention-based neural architectures, such as transformer architectures, are widely used to tackle problems in natural language processing and beyond. Here we investigate the fundamental building blocks of attention and their computational properties. Within the standard model of deep learning, we classify all possible fundamental building blocks of attention in terms of their source, target, and computational mechanism. We identify and study three most important mechanisms: additive activation attention, multiplicative output attention (output gating), and multiplicative synaptic attention (synaptic gating). The gating mechanisms correspond to multiplicative extensions of the standard model and are used across all current attention-based deep learning architectures. We study their functional properties and estimate the capacity of several attentional building blocks in the case of linear and polynomial threshold gates. Surprisingly, additive activation attention plays a central role in the proofs of the lower bounds. Attention mechanisms reduce the depth of certain basic circuits and leverage the power of quadratic activations without incurring their full cost.


title: Fuzzy Pooling

abstract: Convolutional Neural Networks (CNNs) are artificial learning systems typically based on two operations: convolution, which implements feature extraction through filtering, and pooling, which implements dimensionality reduction. The impact of pooling in the classification performance of the CNNs has been highlighted in several previous works, and a variety of alternative pooling operators have been proposed. However, only a few of them tackle with the uncertainty that is naturally propagated from the input layer to the feature maps of the hidden layers through convolutions. In this paper we present a novel pooling operation based on (type-1) fuzzy sets to cope with the local imprecision of the feature maps, and we investigate its performance in the context of image classification. Fuzzy pooling is performed by fuzzification, aggregation and defuzzification of feature map neighborhoods. It is used for the construction of a fuzzy pooling layer that can be applied as a drop-in replacement of the current, crisp, pooling layers of CNN architectures. Several experiments using publicly available datasets show that the proposed approach can enhance the classification performance of a CNN. A comparative evaluation shows that it outperforms state-of-the-art pooling approaches.


title: Text-Based Action-Model Acquisition for Planning

abstract: Although there have been approaches that are capable of learning action models from plan traces, there is no work on learning action models from textual observations, which is pervasive and much easier to collect from real-world applications compared to plan traces. In this paper we propose a novel approach to learning action models from natural language texts by integrating Constraint Satisfaction and Natural Language Processing techniques. Specifically, we first build a novel language model to extract plan traces from texts, and then build a set of constraints to generate action models based on the extracted plan traces. After that, we iteratively improve the language model and constraints until we achieve the convergent language model and action models. We empirically exhibit that our approach is both effective and efficient.


title: Composite excitonic states in doped semiconductors

abstract: We present a theoretical model of composite excitonic states in doped semiconductors. Many-body interactions between a photoexcited electron-hole pair and the electron gas are integrated into a computationally tractable few-body problem, solved by the variational method. We focus on electron-doped ML-MoSe$_2$ and ML-WSe$_2$ due to the contrasting character of their conduction bands. In both cases, the core of the composite is a tightly-bound trion (two electrons and valence-band hole), surrounded by a region depleted of electrons. The composite in ML-WSe$_2$ further includes a satellite electron with different quantum numbers. The theory is general and can be applied to semiconductors with various energy-band properties, allowing one to calculate their excitonic states and to quantify the interaction with the Fermi sea.


title: Hexcitons and oxcitons in monolayer WSe$_2$

abstract: In the archetypal monolayer semiconductor WSe$_2$, the distinct ordering of spin-polarized valleys (low-energy pockets) in the conduction band allows for studies of not only simple neutral excitons and charged excitons (i.e., trions), but also more complex many-body states that are predicted at higher electron densities. We discuss magneto-optical measurements of electron-rich WSe$_2$ monolayers, and interpret the spectral lines that emerge at high electron doping as optical transitions of 6-body exciton states ("hexcitons") and 8-body exciton states ("oxcitons"). These many-body states emerge when a photoexcited electron-hole pair interacts simultaneously with multiple Fermi seas, each having distinguishable spin and valley quantum numbers. In addition, we identify the energies of primary and satellite optical transitions of hexcitons in the photoluminescence spectrum.


title: How to Fill the Optimum Set? Population Gradient Descent with Harmless Diversity

abstract: Although traditional optimization methods focus on finding a single optimal solution, most objective functions in modern machine learning problems, especially those in deep learning, often have multiple or infinite numbers of optima. Therefore, it is useful to consider the problem of finding a set of diverse points in the optimum set of an objective function. In this work, we frame this problem as a bi-level optimization problem of maximizing a diversity score inside the optimum set of the main loss function, and solve it with a simple population gradient descent framework that iteratively updates the points to maximize the diversity score in a fashion that does not hurt the optimization of the main loss. We demonstrate that our method can efficiently generate diverse solutions on a variety of applications, including text-to-image generation, text-to-mesh generation, molecular conformation generation and ensemble neural network training.


title: Generic nonadditivity of quantum capacity in simple channels

abstract: Determining capacities of quantum channels is a fundamental question in quantum information theory. Despite having rigorous coding theorems quantifying the flow of information across quantum channels, their capacities are poorly understood due to super-additivity effects. Studying these phenomena is important for deepening our understanding of quantum information, yet simple and clean examples of super-additive channels are scarce. Here we study a family of channels called platypus channels. Its simplest member, a qutrit channel, is shown to display super-additivity of coherent information when used jointly with a variety of qubit channels. Higher-dimensional family members display super-additivity of quantum capacity together with an erasure channel. Subject to the "spin-alignment conjecture" introduced in the companion paper [IEEE Trans. Inf. Theory 69(6), pp. 3825-3849, 2023; arXiv:2202.08380], our results on super-additivity of quantum capacity extend to lower-dimensional channels as well as larger parameter ranges. In particular, super-additivity occurs between two weakly additive channels each with large capacity on their own, in stark contrast to previous results. Remarkably, a single, novel transmission strategy achieves super-additivity in all examples. Our results show that super-additivity is much more prevalent than previously thought. It can occur across a wide variety of channels, even when both participating channels have large quantum capacity.


title: Turning many-body problems to few-body ones in photoexcited semiconductors using the stochastic variational method in momentum space, SVM-k

abstract: We develop an efficient computational technique to calculate composite excitonic states in photoexcited semiconductors through the stochastic variational method (SVM). Many-body interactions between an electron gas and the excitonic state are embodied in the problem through Fermi holes in the conduction band, introduced when electrons are pulled out of the Fermi sea to bind the photoexcited electron-hole pair. We consider the direct Coulomb interaction between distinguishable particles in the complex, the exchange-induced band-gap renormalization effect, and electron-hole exchange interaction between an electron and its conduction-band hole. We provide analytical expressions for potential matrix elements, using a technique that allows us to circumvent the difficulty imposed by the occupation of low-energy electron states in the conduction band. We discuss the computational steps one should implement in order to perform the calculation, and how to extract kinetic energies of individual particles in the complex, average inter-particle distances, and density distributions.


title: Tetrons, pexcitons, and hexcitons in monolayer transition-metal dichalcogenides

abstract: We present a comprehensive theoretical analysis of composite excitonic states in doped transition-metal dichalcogenide monolayers. Making use of the pair distribution function, we introduce a method to include the effect of screening in the interaction between the excitonic state and Fermi-sea electrons. Employing the screened potential, we study tetrons and pexcitons in which a trion is bound to one and two holes in the conduction band, respectively. A conduction-band hole denotes the lack of electrons with certain quantum numbers at the vicinity of the trion. We then analyze the hexciton complex wherein conduction-band holes facilitate binding of three electrons with one valence-band hole. We introduce a simple model from which one can readily calculate the binding energy of the third (satellite) electron in the hexciton complex. Finally, we compare the simulated results with experiment and point out the shortcomings and successes of the analysis.


title: The platypus of the quantum channel zoo

abstract: Understanding quantum channels and the strange behavior of their capacities is a key objective of quantum information theory. Here we study a remarkably simple, low-dimensional, single-parameter family of quantum channels with exotic quantum information-theoretic features. As the simplest example from this family, we focus on a qutrit-to-qutrit channel that is intuitively obtained by hybridizing together a simple degradable channel and a completely useless qubit channel. Such hybridizing makes this channel's capacities behave in a variety of interesting ways. For instance, the private and classical capacity of this channel coincide and can be explicitly calculated, even though the channel does not belong to any class for which the underlying information quantities are known to be additive. Moreover, the quantum capacity of the channel can be computed explicitly, given a clear and compelling conjecture is true. This "spin alignment conjecture," which may be of independent interest, is proved in certain special cases and additional numerical evidence for its validity is provided. Finally, we generalize the qutrit channel in two ways, and the resulting channels and their capacities display similarly rich behavior. In the companion paper [Phys. Rev. Lett. 130, 200801 (2023); arXiv:2202.08377], we further show that the qutrit channel demonstrates superadditivity when transmitting quantum information jointly with a variety of assisting channels, in a manner unknown before.


title: Improving Performance Bounds for Weighted Round-Robin Schedulers under Constrained Cross-Traffic

abstract: Weighted round robin (WRR) is an effective, yet particularly easy-to-implement packet scheduler. A slight modification in the implementation of WRR, interleaved weighted round robin, has been proposed as an enhancement of the initial version and has been recently investigated. Network calculus is a versatile framework to model and analyze such network schedulers. By means of this, one can derive theoretical upper bounds on network performance metrics, such as delay or backlog. In our previous work, we derive performance bounds by showing that both round-robin variants belong to a class called bandwidth-sharing policy; however, the proofs are incomplete and thus, we cannot conclude that the round-robin schedulers are bandwidth-sharing policies (under variable packet sizes).To that end, in the subsequent erratum, we introduce so-called resource-segregating policies and show the round-robin schedulers to be members of this class. We first present our original work, as published in [CNS22-1], and then the erratum correcting the previously mentioned shortcoming. In our erratum, we provide slightly worse delay bounds compared to [CNS22-1]; yet, across all our experiments, they significantly outperform the state of the art.


title: The relative class number one problem for function fields, I

abstract: We reduce the classification of finite extensions of function fields (of curves over finite fields) with the same class number to a finite computation; complete this computation in all cases except when both curves have base field $\mathbb{F}_2$ and genus $>1$; and give a conjectural answer in the remaining cases. The conjecture will be resolved in subsequent papers.


title: Morse Graphs: Topological Tools for Analyzing the Global Dynamics of Robot Controllers

abstract: Understanding the global dynamics of a robot controller, such as identifying attractors and their regions of attraction (RoA), is important for safe deployment and synthesizing more effective hybrid controllers. This paper proposes a topological framework to analyze the global dynamics of robot controllers, even data-driven ones, in an effective and explainable way. It builds a combinatorial representation representing the underlying system's state space and non-linear dynamics, which is summarized in a directed acyclic graph, the Morse graph. The approach only probes the dynamics locally by forward propagating short trajectories over a state-space discretization, which needs to be a Lipschitz-continuous function. The framework is evaluated given either numerical or data-driven controllers for classical robotic benchmarks. It is compared against established analytical and recent machine learning alternatives for estimating the RoAs of such controllers. It is shown to outperform them in accuracy and efficiency. It also provides deeper insights as it describes the global dynamics up to the discretization's resolution. This allows to use the Morse graph to identify how to synthesize controllers to form improved hybrid solutions or how to identify the physical limitations of a robotic system.


title: Limitations of Neural Collapse for Understanding Generalization in Deep Learning

abstract: The recent work of Papyan, Han, & Donoho (2020) presented an intriguing "Neural Collapse" phenomenon, showing a structural property of interpolating classifiers in the late stage of training. This opened a rich area of exploration studying this phenomenon. Our motivation is to study the upper limits of this research program: How far will understanding Neural Collapse take us in understanding deep learning? First, we investigate its role in generalization. We refine the Neural Collapse conjecture into two separate conjectures: collapse on the train set (an optimization property) and collapse on the test distribution (a generalization property). We find that while Neural Collapse often occurs on the train set, it does not occur on the test set. We thus conclude that Neural Collapse is primarily an optimization phenomenon, with as-yet-unclear connections to generalization. Second, we investigate the role of Neural Collapse in feature learning. We show simple, realistic experiments where training longer leads to worse last-layer features, as measured by transfer-performance on a downstream task. This suggests that neural collapse is not always desirable for representation learning, as previously claimed. Finally, we give preliminary evidence of a "cascading collapse" phenomenon, wherein some form of Neural Collapse occurs not only for the last layer, but in earlier layers as well. We hope our work encourages the community to continue the rich line of Neural Collapse research, while also considering its inherent limitations.


title: Monte Carlo studies of the Blume-Capel model on nonregular two- and three-dimensional lattices: Phase diagrams, tricriticality, and critical exponents

abstract: We perform Monte Carlo simulations, combining both the Wang-Landau and the Metropolis algorithms, to investigate the phase diagrams of the Blume-Capel model on different types of nonregular lattices (Lieb lattice (LL), decorated triangular lattice (DTL), and decorated simple cubic lattice (DSC)). The nonregular character of the lattices induces a double transition (reentrant behavior) in the region of the phase diagram at which the nature of the phase transition changes from first-order to second-order. A physical mechanism underlying this reentrance is proposed. The large-scale Monte Carlo simulations are performed with the finite-size scaling analysis to compute the critical exponents and the critical Binder cumulant for three different values of the anisotropy $\Delta/J \in \big\{ 0, 1, 1.34 \textrm{ (for LL)}, 1.51 \textrm{ (for DTL and DSC)} \big\}$, showing thus no deviation from the standard Ising universality class in two and three dimensions. We report also the location of the tricritical point to considerable precision: ($\Delta_t/J=1.3457(1)$; $k_B T_t/J=0.309(2)$), ($\Delta_t/J=1.5766(1)$; $k_B T_t/J=0.481(2)$), and ($\Delta_t/J=1.5933(1)$; $k_B T_t/J=0.569(4)$) for LL, DTL, and DSC, respectively.


title: Laplacian operator on statistical manifold

abstract: In this paper, we define a Laplacian operator on a statistical manifold, called the vector Laplacian. This vector Laplacian incorporates information from the Amari-Chentsov tensor. We derive a formula for the vector Laplacian. We also give two applications using the heat kernel associated with the vector Laplacian.


title: TROPHY: Trust Region Optimization Using a Precision Hierarchy

abstract: We present an algorithm to perform trust-region-based optimization for nonlinear unconstrained problems. The method selectively uses function and gradient evaluations at different floating-point precisions to reduce the overall energy consumption, storage, and communication costs; these capabilities are increasingly important in the era of exascale computing. In particular, we are motivated by a desire to improve computational efficiency for massive climate models. We employ our method on two examples: the CUTEst test set and a large-scale data assimilation problem to recover wind fields from radar returns. Although this paper is primarily a proof of concept, we show that if implemented on appropriate hardware, the use of mixed-precision can significantly reduce the computational load compared with fixed-precision solvers.


title: Generalizable Information Theoretic Causal Representation

abstract: It is evidence that representation learning can improve model's performance over multiple downstream tasks in many real-world scenarios, such as image classification and recommender systems. Existing learning approaches rely on establishing the correlation (or its proxy) between features and the downstream task (labels), which typically results in a representation containing cause, effect and spurious correlated variables of the label. Its generalizability may deteriorate because of the unstability of the non-causal parts. In this paper, we propose to learn causal representation from observational data by regularizing the learning procedure with mutual information measures according to our hypothetical causal graph. The optimization involves a counterfactual loss, based on which we deduce a theoretical guarantee that the causality-inspired learning is with reduced sample complexity and better generalization ability. Extensive experiments show that the models trained on causal representations learned by our approach is robust under adversarial attacks and distribution shift.


title: On solutions of codimension-one $A$-hypergeometric systems

abstract: By a codimension-one system we mean a system whose lattice of relations has rank one. We consider codimension-one $A$-hypergeometric systems and explicitly construct some of the logarithmic series solutions at the origin. When the parameter vector $\beta$ is nonresonant we obtain a full set of logarithmic series solutions at the origin by this procedure. We also determine when a codimension-one system with nonresonant parameter can have maximal unipotent monodromy at the origin.


title: An analogue of the Robin inequality of the second type for odd integers

abstract: In this paper we give a variant of the Robin inequality which states that $\frac{\sigma(n)}{n} \leq \frac{e^\gamma}{2} \log\log n + \frac{0.7398\cdots}{\log\log n}$ for any odd integer $n \geq 3$.


title: Graph Masked Autoencoders with Transformers

abstract: Recently, transformers have shown promising performance in learning graph representations. However, there are still some challenges when applying transformers to real-world scenarios due to the fact that deep transformers are hard to train from scratch and the quadratic memory consumption w.r.t. the number of nodes. In this paper, we propose Graph Masked Autoencoders (GMAEs), a self-supervised transformer-based model for learning graph representations. To address the above two challenges, we adopt the masking mechanism and the asymmetric encoder-decoder design. Specifically, GMAE takes partially masked graphs as input, and reconstructs the features of the masked nodes. The encoder and decoder are asymmetric, where the encoder is a deep transformer and the decoder is a shallow transformer. The masking mechanism and the asymmetric design make GMAE a memory-efficient model compared with conventional transformers. We show that, when serving as a conventional self-supervised graph representation model, GMAE achieves state-of-the-art performance on both the graph classification task and the node classification task under common downstream evaluation protocols. We also show that, compared with training in an end-to-end manner from scratch, we can achieve comparable performance after pre-training and fine-tuning using GMAE while simplifying the training process.


title: Charge-Swapping Q-balls in a Logarithmic Potential and Affleck-Dine condensate fragmentation

abstract: We study charge-swapping Q-balls, a kind of composite Q-ball where positive and negative charges co-exist and swap with time, in models with a logarithmic potential that arises naturally in supersymmetric extensions of the Standard Model. We show that charge-swapping Q-balls can be copiously generated in the Affleck-Dine fragmentation process in the early universe. We find that the charge-swapping Q-balls with the logarithmic potential are extremely stable. By performing long time, parallelized lattice simulations with absorbing boundary conditions, we find that the lifetimes of such objects with low multipoles are at least $4.6 \times 10^5/m$ in 3+1D and $2.5 \times 10^7/m$ in 2+1D, where $m$ is the mass scale of the scalar field. We also chart the attractor basin of the initial conditions to form these charge-swapping Q-balls.


title: Robust Reinforcement Learning via Genetic Curriculum

abstract: Achieving robust performance is crucial when applying deep reinforcement learning (RL) in safety critical systems. Some of the state of the art approaches try to address the problem with adversarial agents, but these agents often require expert supervision to fine tune and prevent the adversary from becoming too challenging to the trainee agent. While other approaches involve automatically adjusting environment setups during training, they have been limited to simple environments where low-dimensional encodings can be used. Inspired by these approaches, we propose genetic curriculum, an algorithm that automatically identifies scenarios in which the agent currently fails and generates an associated curriculum to help the agent learn to solve the scenarios and acquire more robust behaviors. As a non-parametric optimizer, our approach uses a raw, non-fixed encoding of scenarios, reducing the need for expert supervision and allowing our algorithm to adapt to the changing performance of the agent. Our empirical studies show improvement in robustness over the existing state of the art algorithms, providing training curricula that result in agents being 2 - 8x times less likely to fail without sacrificing cumulative reward. We include an ablation study and share insights on why our algorithm outperforms prior approaches.


title: Dynamical phases in growing populations: understanding recovery from bottlenecks

abstract: Since steep declines in a population's size also typically alter its composition, population bottlenecks are considered highly important for evolution. However, despite such significance, the mechanisms governing the impact of a given population bottleneck remain poorly understood. In this context, we show that long-term post-bottleneck outcomes can depend crucially on the rate at which a diminished population grows whilst recovering. That is, two otherwise identical populations, each having undergone the same bottleneck, can fixate on dramatically different demographics due to different rates of post-bottleneck growth. This behaviour is moreover shown to change non-trivially with different levels of mutation. The underlying mechanism can be traced to intrinsic fluctuations whose standard deviation decreases in time with the inverse square root of the growing population size. Crucially, such decreasing fluctuations couple to the underlying dynamics and result in distinct regimes of (non-equilibrium) demographic behaviour, delimited by abrupt transitions at critical population sizes. Seen through this lens, growth and mutation alter evolutionary outcomes by changing the duration and character of dynamical phases; a feature that we speculate may be of generic importance across many systems.


title: SWIM: Selective Write-Verify for Computing-in-Memory Neural Accelerators

abstract: Computing-in-Memory architectures based on non-volatile emerging memories have demonstrated great potential for deep neural network (DNN) acceleration thanks to their high energy efficiency. However, these emerging devices can suffer from significant variations during the mapping process i.e., programming weights to the devices), and if left undealt with, can cause significant accuracy degradation. The non-ideality of weight mapping can be compensated by iterative programming with a write-verify scheme, i.e., reading the conductance and rewriting if necessary. In all existing works, such a practice is applied to every single weight of a DNN as it is being mapped, which requires extensive programming time. In this work, we show that it is only necessary to select a small portion of the weights for write-verify to maintain the DNN accuracy, thus achieving significant speedup. We further introduce a second derivative based technique SWIM, which only requires a single pass of forward and backpropagation, to efficiently select the weights that need write-verify. Experimental results on various DNN architectures for different datasets show that SWIM can achieve up to 10x programming speedup compared with conventional full-blown write-verify while attaining a comparable accuracy.


title: Augment with Care: Contrastive Learning for Combinatorial Problems

abstract: Supervised learning can improve the design of state-of-the-art solvers for combinatorial problems, but labelling large numbers of combinatorial instances is often impractical due to exponential worst-case complexity. Inspired by the recent success of contrastive pre-training for images, we conduct a scientific study of the effect of augmentation design on contrastive pre-training for the Boolean satisfiability problem. While typical graph contrastive pre-training uses label-agnostic augmentations, our key insight is that many combinatorial problems have well-studied invariances, which allow for the design of label-preserving augmentations. We find that label-preserving augmentations are critical for the success of contrastive pre-training. We show that our representations are able to achieve comparable test accuracy to fully-supervised learning while using only 1% of the labels. We also demonstrate that our representations are more transferable to larger problems from unseen domains. Our code is available at https://github.com/h4duan/contrastive-sat.


title: Superrotation of Titan's stratosphere driven by the radiative heating of the haze layer

abstract: Titan's stratosphere has been observed in a superrotation state, where the atmosphere rotates many times faster than the surface does. Another characteristics of Titan's atmosphere is the presence of thick haze layer. In this paper, we performed numerical experiments using a General Circulation Model (GCM), to explore the effects of the haze layer on the stratospheric superrotation. We employed a semi-gray radiation model of Titan's atmosphere following McKay et al. (1999), which takes account of the sunlight absorption by haze particles. The phase change of methane or the seasonal changes were not taken into account. Our model with the radiation parameters tuned for Titan yielded the global eastward wind around the equator with larger velocities at higher altitudes except at around 70 km after $10^5$ Earth days. Although the atmosphere is not in an equilibrium state, the zonal wind profiles is approximately consistent with the observed one. Analysis on our experiments suggests that the quasi-stationary stratospheric superrotation is maintained by the balance between the meridional circulation decoupled from the surface, and the eddies that transport angular momentum equatorward. This is different from, but similar to, the so-called Gierasch mechanism, in which momentum is supplied from the surface. This structure may explain the no-wind region at about $80$ km in altitude.


title: Modelling stars with Gaussian Process Regression: Augmenting Stellar Model Grid

abstract: Grid-based modelling is widely used for estimating stellar parameters. However, stellar model grid is sparse because of the computational cost. This paper demonstrates an application of a machine-learning algorithm using the Gaussian Process (GP) Regression that turns a sparse model grid onto a continuous function. We train GP models to map five fundamental inputs (mass, equivalent evolutionary phase, initial metallicity, initial helium fraction, and the mixing-length parameter) to observable outputs (effective temperature, surface gravity, radius, surface metallicity, and stellar age). We test the GP predictions for the five outputs using off-grid stellar models and find no obvious systematic offsets, indicating good accuracy in predictions.As a further validation, we apply these GP models to characterise 1,000 fake stars. Inferred masses and ages determined with GP models well recover true values within one standard deviation. An important consequence of using GP-based interpolation is that stellar ages are more precise than those estimated with the original sparse grid because of the full sampling of fundamental inputs.


title: Shift-Memory Network for Temporal Scene Segmentation

abstract: Semantic segmentation has achieved great accuracy in understanding spatial layout. For real-time tasks based on dynamic scenes, we extend semantic segmentation in temporal domain to enhance the spatial accuracy with motion. We utilize a shift-mode network over streaming input to ensure zero-latency output. For the data overlap under shifting network, this paper identifies repeated computation in fixed periods across network layers. To avoid this redundancy, we derive a Shift-Memory Network (SMN) from encoding-decoding baseline to reuse the network values without accuracy loss. Trained in patch-mode, the SMN extracts the network parameters for SMN to perform inference promptly in compact memory. We segment dynamic scenes from 1D scanning input and 2D video. The experiments of SMN achieve equivalent accuracy as shift-mode but in faster inference speeds and much smaller memory. This will facilitate semantic segmentation in real-time application on edge devices.


title: Real Time Motion Planning Using Constrained Iterative Linear Quadratic Regulator for On-Road Self-Driving

abstract: Collision avoidance is one of the most challenging tasks people need to consider for developing the self-driving technology. In this paper we propose a new spatiotemporal motion planning algorithm that efficiently solves a constrained nonlinear optimal control problem using the iterative linear quadratic regulator (iLQR), which takes into account the uncertain driving behaviors of the traffic vehicles and minimizes the collision risks between the self-driving vehicle (referred to as the "ego" vehicle) and the traffic vehicles such that the ego vehicle is able to maintain sufficiently large distances to all the surrounding vehicles for achieving the desired collision avoidance maneuver in traffic. To this end, we introduce the concept of the "collision polygon" for computing the minimum distances between the ego vehicle and the traffic vehicles, and provide two different solutions for designing the constraints of the motion planning problem by properly modeling the behaviors of the traffic vehicles in order to evaluate the collision risk. Finally, the iLQR motion planning algorithm is validated in multiple real-time tasks for collision avoidance using both a simulator and a level-3 autonomous driving test platform.


title: Multi-Modal Fusion in Contact-Rich Precise Tasks via Hierarchical Policy Learning

abstract: Combined visual and force feedback play an essential role in contact-rich robotic manipulation tasks. Current methods focus on developing the feedback control around a single modality while underrating the synergy of the sensors. Fusing different sensor modalities is necessary but remains challenging. A key challenge is to achieve an effective multi-modal and generalized control scheme to novel objects with precision. This paper proposes a practical multi-modal sensor fusion mechanism using hierarchical policy learning. To begin with, we use a self-supervised encoder that extracts multi-view visual features and a hybrid motion/force controller that regulates force behaviors. Next, the multi-modality fusion is simplified by hierarchical integration of the vision, force, and proprioceptive data in the reinforcement learning (RL) algorithm. Moreover, with hierarchical policy learning, the control scheme can exploit the visual feedback limits and explore the contribution of individual modality in precise tasks. Experiments indicate that robots with the control scheme could assemble objects with 0.25mm clearance in simulation. The system could be generalized to widely varied initial configurations and new shapes. Experiments validate that the simulated system can be robustly transferred to reality without fine-tuning.


title: Federated Stochastic Gradient Descent Begets Self-Induced Momentum

abstract: Federated learning (FL) is an emerging machine learning method that can be applied in mobile edge systems, in which a server and a host of clients collaboratively train a statistical model utilizing the data and computation resources of the clients without directly exposing their privacy-sensitive data. We show that running stochastic gradient descent (SGD) in such a setting can be viewed as adding a momentum-like term to the global aggregation process. Based on this finding, we further analyze the convergence rate of a federated learning system by accounting for the effects of parameter staleness and communication resources. These results advance the understanding of the Federated SGD algorithm, and also forges a link between staleness analysis and federated computing systems, which can be useful for systems designers.


title: Moderate deviations for fully coupled multiscale weakly interacting particle systems

abstract: We consider a collection of fully coupled weakly interacting diffusion processes moving in a two-scale environment. We study the moderate deviations principle of the empirical distribution of the particles' positions in the combined limit as the number of particles grow to infinity and the time-scale separation parameter goes to zero simultaneously. We make use of weak convergence methods, which provide a convenient representation for the moderate deviations rate function in terms of an effective mean field control problem. We rigorously obtain equivalent representations for the moderate deviations rate function in an appropriate "negative Sobolev" form, which is reminiscent of the large deviations rate function for the empirical measure of weakly interacting diffusions obtained in the 1987 seminal paper by Dawson-G\"{a}rtner. In the course of the proof we obtain related ergodic theorems and we rigorously study the regularity of Poisson type of equations associated to McKean-Vlasov problems, both of which are topics of independent interest.


title: Global-in-time existence of weak solutions for Vlasov-Manev-Fokker-Planck system

abstract: We consider the Vlasov-Manev-Fokker-Planck (VMFP) system in three dimensions, which differs from the Vlasov-Poisson-Fokker-Planck in that it has the gravitational potential of the form $-1/r - 1/r^2$ instead of the Newtonian one. For the VMFP system, we establish the global-in-time existence of weak solutions. The proof extends to several related kinetic systems.


title: UV counterpart of an X-ray ultra-fast outflow in IRAS 17020+4544

abstract: We report on the discovery of a UV absorption counterpart of a low-ionization X-ray ultra-fast outflow (UFO) in the Narrow-Line Seyfert-1 galaxy IRAS 17020+4544. This UV signature of the UFO is seen as a narrow and blueshifted Lyman-alpha absorption feature in the far-UV spectrum, taken with the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST). The Lyman-alpha feature is found to be outflowing with a velocity of -23430 km/s (0.078 c). We carry out high-resolution UV spectroscopy and photoionization modeling to study the UFO that is seen in the HTS/COS spectrum. The results of our modeling show that the UV UFO corresponds to a low-ionization, low-velocity component of the X-ray UFO found previously with XMM-Newton's Reflection Grating Spectrometer (RGS). The other higher-velocity and higher-ionization components of the X-ray UFOs are not significantly detected in the HST/COS spectrum, consistent with predictions of our photoionization calculations. The multiple ionization and velocity components of the UFOs in IRAS 17020+4544 suggest a scenario where a powerful primary UFO entrains and shocks the ambient medium, resulting in formation of weaker secondary UFO components, such as the one found in the UV band.


title: The Unboxing Experience: Exploration and Design of Initial Interactions Between Children and Social Robots

abstract: Social robots are increasingly introduced into children's lives as educational and social companions, yet little is known about how these products might best be introduced to their environments. The emergence of the "unboxing" phenomenon in media suggests that introduction is key to technology adoption where initial impressions are made. To better understand this phenomenon toward designing a positive unboxing experience in the context of social robots for children, we conducted three field studies with families of children aged 8 to 13: (1) an exploratory free-play activity ($n=12$); (2) a co-design session ($n=11$) that informed the development of a prototype box and a curated unboxing experience; and (3) a user study ($n=9$) that evaluated children's experiences. Our findings suggest the unboxing experience of social robots can be improved through the design of a creative aesthetic experience that engages the child socially to guide initial interactions and foster a positive child-robot relationship.


title: AutoScore-Ordinal: An interpretable machine learning framework for generating scoring models for ordinal outcomes

abstract: Background: Risk prediction models are useful tools in clinical decision-making which help with risk stratification and resource allocations and may lead to a better health care for patients. AutoScore is a machine learning-based automatic clinical score generator for binary outcomes. This study aims to expand the AutoScore framework to provide a tool for interpretable risk prediction for ordinal outcomes. Methods: The AutoScore-Ordinal framework is generated using the same 6 modules of the original AutoScore algorithm including variable ranking, variable transformation, score derivation (from proportional odds models), model selection, score fine-tuning, and model evaluation. To illustrate the AutoScore-Ordinal performance, the method was conducted on electronic health records data from the emergency department at Singapore General Hospital over 2008 to 2017. The model was trained on 70% of the data, validated on 10% and tested on the remaining 20%. Results: This study included 445,989 inpatient cases, where the distribution of the ordinal outcome was 80.7% alive without 30-day readmission, 12.5% alive with 30-day readmission, and 6.8% died inpatient or by day 30 post discharge. Two point-based risk prediction models were developed using two sets of 8 predictor variables identified by the flexible variable selection procedure. The two models indicated reasonably good performance measured by mean area under the receiver operating characteristic curve (0.785 and 0.793) and generalized c-index (0.737 and 0.760), which were comparable to alternative models. Conclusion: AutoScore-Ordinal provides an automated and easy-to-use framework for development and validation of risk prediction models for ordinal outcomes, which can systematically identify potential predictors from high-dimensional data.


title: Multivariate Time Series Forecasting with Dynamic Graph Neural ODEs

abstract: Multivariate time series forecasting has long received significant attention in real-world applications, such as energy consumption and traffic prediction. While recent methods demonstrate good forecasting abilities, they have three fundamental limitations. (i) Discrete neural architectures: Interlacing individually parameterized spatial and temporal blocks to encode rich underlying patterns leads to discontinuous latent state trajectories and higher forecasting numerical errors. (ii) High complexity: Discrete approaches complicate models with dedicated designs and redundant parameters, leading to higher computational and memory overheads. (iii) Reliance on graph priors: Relying on predefined static graph structures limits their effectiveness and practicability in real-world applications. In this paper, we address all the above limitations by proposing a continuous model to forecast $\textbf{M}$ultivariate $\textbf{T}$ime series with dynamic $\textbf{G}$raph neural $\textbf{O}$rdinary $\textbf{D}$ifferential $\textbf{E}$quations ($\texttt{MTGODE}$). Specifically, we first abstract multivariate time series into dynamic graphs with time-evolving node features and unknown graph structures. Then, we design and solve a neural ODE to complement missing graph topologies and unify both spatial and temporal message passing, allowing deeper graph propagation and fine-grained temporal information aggregation to characterize stable and precise latent spatial-temporal dynamics. Our experiments demonstrate the superiorities of $\texttt{MTGODE}$ from various perspectives on five time series benchmark datasets.


title: Million.js: A Fast Compiler-Augmented Virtual DOM for the Web

abstract: Interactive web applications created with declarative JavaScript User Interface (UI) libraries have increasingly dominated the modern internet. However, existing libraries are primarily made for run-time execution, and rely on the user to load and render web applications. This led us to create Million.js, a fast compiler-augmented virtual Document Object Model (DOM) for the web. Million.js reduces load time and time-to-interactive by creating a compiler to compute interactive regions of a web application before the user visits the page. The virtual DOM run-time optimizes interactive content through compiler flags, compute batching, scheduling, and reactive data primitives to achieve optimal performance. When benchmarked against the most popular virtual DOM libraries, Million.js resulted in 133% to 300% faster rendering and 2347\% faster load. In a real-world web application with both comparative benchmarks and an informal user study, Million.js loaded 35.11% faster after migrating from React. The findings show that web applications have the potential to be orders of magnitude faster through JavaScript UI libraries that use Million.js.


title: Non-linear Transport by Bethe Bound States

abstract: We consider non-linear ballistic spin transport in the XXZ spin chain and derive an analytical result for the non-linear Drude weight $D^{(3)}$ at infinite temperatures. In contrast to the linear Drude weight $D^{(1)}$, we find that the result not only depends on anisotropy but also on the string length of the quasiparticles transporting the spin current. Our result provides further insights into transport by quasiparticles and raises questions about Luttinger liquid universality.


title: On the local rank of fusion systems

abstract: In this paper we define the notion of local rank for fusion systems so as to reformulate the Alperin's weight conjecture in the framework of block fusion systems following the work by Kn\"orr and Robinson.


title: The Scale Invariant Vacuum Paradigm: main results and current progress

abstract: We present a summary of the main results within the Scale Invariant Vacuum (SIV) paradigm as related to the Weyl Integrable Geometry. After a brief review of the mathematical framework, we will highlight the main results related to inflation within the SIV [9], the growth of the density fluctuations [8], and the application of the SIV to scale-invariant dynamics of Galaxies, MOND, Dark Matter, and the Dwarf Spheroidals [7]. The connection of the weak-field SIV results to the un-proper time parametrization within the re-parametrization paradigm is also discussed [14].


title: Entropic Associative Memory for Manuscript Symbols

abstract: Manuscript symbols can be stored, recognized and retrieved from an entropic digital memory that is associative and distributed but yet declarative; memory retrieval is a constructive operation, memory cues to objects not contained in the memory are rejected directly without search, and memory operations can be performed through parallel computations. Manuscript symbols, both letters and numerals, are represented in Associative Memory Registers that have an associated entropy. The memory recognition operation obeys an entropy trade-off between precision and recall, and the entropy level impacts on the quality of the objects recovered through the memory retrieval operation. The present proposal is contrasted in several dimensions with neural networks models of associative memory. We discuss the operational characteristics of the entropic associative memory for retrieving objects with both complete and incomplete information, such as severe occlusions. The experiments reported in this paper add evidence on the potential of this framework for developing practical applications and computational models of natural memory.


title: FPIC: A Novel Semantic Dataset for Optical PCB Assurance

abstract: Outsourced printed circuit board (PCB) fabrication necessitates increased hardware assurance capabilities. Several assurance techniques based on automated optical inspection (AOI) have been proposed that leverage PCB images acquired using digital cameras. We review state-of-the-art AOI techniques and observe a strong, rapid trend toward machine learning (ML) solutions. These require significant amounts of labeled ground truth data, which is lacking in the publicly available PCB data space. We contribute the FICS PCB Image Collection (FPIC) dataset to address this need. Additionally, we outline new hardware security methodologies enabled by our data set.


title: Continuity Postulates and Solvability Axioms in Economic Theory and in Mathematical Psychology: A Consolidation of the Theory of Individual Choice

abstract: This paper presents four theorems that connect continuity postulates in mathematical economics to solvability axioms in mathematical psychology, and ranks them under alternative supplementary assumptions. Theorem 1 connects notions of continuity (full, separate, Wold, weak Wold, Archimedean, mixture) with those of solvability (restricted, unrestricted) under the completeness and transitivity of a binary relation. Theorem 2 uses the primitive notion of a separately-continuous function to answer the question when an analogous property on a relation is fully continuous. Theorem 3 provides a portmanteau theorem on the equivalence between restricted solvability and various notions of continuity under weak monotonicity. Finally, Theorem 4 presents a variant of Theorem 3 that follows Theorem 1 in dispensing with the dimensionality requirement and in providing partial equivalences between solvability and continuity notions. These theorems are motivated for their potential use in representation theorems.


title: Massive Galaxy Mergers Have Distinctive Global HI Profiles

abstract: The global 21 cm HI emission-line profile of a galaxy encodes valuable information on the spatial distribution and kinematics of the neutral atomic gas. Galaxy interactions significantly influence the HI disk and imprint observable features on the integrated HI line profile. In this work, we study the neutral atomic gas properties of galaxy mergers selected from the Great Observatories All-sky LIRG Survey. The HI spectra come from new observations with the Five-hundred-meter Aperture Spherical Telescope and from a collection of archival data. We quantify the HI profile of the mergers with a newly developed method that uses the curve-of-growth of the line profile. Using a control sample of non-merger galaxies carefully selected to match the stellar mass of the merger sample, we show that mergers have a larger proportion of single-peaked HI profiles, as well as a greater tendency for the HI central velocity to deviate from the systemic optical velocity of the galaxy. By contrast, the HI profiles of mergers are not significantly more asymmetric than those of non-mergers.


title: Retrieval-Augmented Reinforcement Learning

abstract: Most deep reinforcement learning (RL) algorithms distill experience into parametric behavior policies or value functions via gradient updates. While effective, this approach has several disadvantages: (1) it is computationally expensive, (2) it can take many updates to integrate experiences into the parametric model, (3) experiences that are not fully integrated do not appropriately influence the agent's behavior, and (4) behavior is limited by the capacity of the model. In this paper we explore an alternative paradigm in which we train a network to map a dataset of past experiences to optimal behavior. Specifically, we augment an RL agent with a retrieval process (parameterized as a neural network) that has direct access to a dataset of experiences. This dataset can come from the agent's past experiences, expert demonstrations, or any other relevant source. The retrieval process is trained to retrieve information from the dataset that may be useful in the current context, to help the agent achieve its goal faster and more efficiently. he proposed method facilitates learning agents that at test-time can condition their behavior on the entire dataset and not only the current state, or current trajectory. We integrate our method into two different RL agents: an offline DQN agent and an online R2D2 agent. In offline multi-task problems, we show that the retrieval-augmented DQN agent avoids task interference and learns faster than the baseline DQN agent. On Atari, we show that retrieval-augmented R2D2 learns significantly faster than the baseline R2D2 agent and achieves higher scores. We run extensive ablations to measure the contributions of the components of our proposed method.


title: Neural Marionette: Unsupervised Learning of Motion Skeleton and Latent Dynamics from Volumetric Video

abstract: We present Neural Marionette, an unsupervised approach that discovers the skeletal structure from a dynamic sequence and learns to generate diverse motions that are consistent with the observed motion dynamics. Given a video stream of point cloud observation of an articulated body under arbitrary motion, our approach discovers the unknown low-dimensional skeletal relationship that can effectively represent the movement. Then the discovered structure is utilized to encode the motion priors of dynamic sequences in a latent structure, which can be decoded to the relative joint rotations to represent the full skeletal motion. Our approach works without any prior knowledge of the underlying motion or skeletal structure, and we demonstrate that the discovered structure is even comparable to the hand-labeled ground truth skeleton in representing a 4D sequence of motion. The skeletal structure embeds the general semantics of possible motion space that can generate motions for diverse scenarios. We verify that the learned motion prior is generalizable to the multi-modal sequence generation, interpolation of two poses, and motion retargeting to a different skeletal structure.


title: High-Dimensional Time-Varying Coefficient Estimation

abstract: In this paper, we develop a novel high-dimensional time-varying coefficient estimation method, based on high-dimensional Ito diffusion processes. To account for high-dimensional time-varying coefficients, we first estimate local (or instantaneous) coefficients using a time-localized Dantzig selection scheme under a sparsity condition, which results in biased local coefficient estimators due to the regularization. To handle the bias, we propose a debiasing scheme, which provides well-performing unbiased local coefficient estimators. With the unbiased local coefficient estimators, we estimate the integrated coefficient, and to further account for the sparsity of the coefficient process, we apply thresholding schemes. We call this Thresholding dEbiased Dantzig (TED). We establish asymptotic properties of the proposed TED estimator. In the empirical analysis, we apply the TED procedure to analyzing high-dimensional factor models using high-frequency data.


title: Time-Correlated Sparsification for Efficient Over-the-Air Model Aggregation in Wireless Federated Learning

abstract: Federated edge learning (FEEL) is a promising distributed machine learning (ML) framework to drive edge intelligence applications. However, due to the dynamic wireless environments and the resource limitations of edge devices, communication becomes a major bottleneck. In this work, we propose time-correlated sparsification with hybrid aggregation (TCS-H) for communication-efficient FEEL, which exploits jointly the power of model compression and over-the-air computation. By exploiting the temporal correlations among model parameters, we construct a global sparsification mask, which is identical across devices, and thus enables efficient model aggregation over-the-air. Each device further constructs a local sparse vector to explore its own important parameters, which are aggregated via digital communication with orthogonal multiple access. We further design device scheduling and power allocation algorithms for TCS-H. Experiment results show that, under limited communication resources, TCS-H can achieve significantly higher accuracy compared to the conventional top-K sparsification with orthogonal model aggregation, with both i.i.d. and non-i.i.d. data distributions.


title: Some identities involving degenerate r-Stirling numbers

abstract: Recently, authors studied the unsigned degenerate r-Stirling number of the first kind and the degenerate r-Stirling number of the second kind, respectively of which are the degenerate versions of the unsigned r-Stirling numbers of the first kind and those of the r-Stirling numbers of the second kind. The aim of this paper is to derive some identities involving such special numbers from the inverse relations for them.


title: Euler's scheme of Mckean-Vlasov SDEs with non-Lipschitz coefficients

abstract: In this paper, we show the strong well-posedness of Mckean-Vlasov SDEs with non-Lipschitz coefficients. Moreover, propagation of chaos and the convergence rate for Euler's scheme of Mckean-Vlasov SDEs are also obtained.


title: Chord-Conditioned Melody Harmonization with Controllable Harmonicity

abstract: Melody harmonization has long been closely associated with chorales composed by Johann Sebastian Bach. Previous works rarely emphasised chorale generation conditioned on chord progressions, and there has been a lack of focus on assistive compositional tools. In this paper, we first designed a music representation that encoded chord symbols for chord conditioning, and then proposed DeepChoir, a melody harmonization system that can generate a four-part chorale for a given melody conditioned on a chord progression. With controllable harmonicity, users can control the extent of harmonicity for generated chorales. Experimental results reveal the effectiveness of the music representation and the controllability of DeepChoir.


title: Higher-order Topological Phases of Magnons in van der Waals Honeycomb Ferromagnets

abstract: We theoretically propose a second-order topological magnon insulator by stacking the van der Waals honeycomb ferromagnets with antiferromagnetic interlayer coupling. The system exhibits Z$_{2}$ topological phase, protected by pseudo-time-reversal symmetry (PTRS). An easy-plane anisotropy term breaks PTRS and destroys the topological phase. Nevertheless, it respects a magnetic two-fold rotational symmetry which protects a second-order topological phase with corner modes in bilayer and hinge modes along stacking direction. Moreover, an introduced staggered interlayer coupling establishes a Z$_{2}$$\times$Z topology, giving rise to gapped topological surface modes carrying non-zero Chern numbers. Consequently, chiral hinge modes propagate along the horizontal hinges in a cuboid geometry and are robust against disorders. Our work bridges the higher-order topology and magnons in van der Waals platforms, and could be used for constructing topological magnonic devices.


title: The log canonical threshold and rational singularities

abstract: We show that if $f$ is a nonzero, noninvertible function on a smooth complex variety $X$ and $J_f$ is the Jacobian ideal of $f$, then ${\rm lct}(f,J_f^2)>1$ if and only if the hypersurface defined by $f$ has rational singularities. Moreover, if it does not have rational singularities, then ${\rm lct}(f,J_f^2)={\rm lct}(f)$. We give two proofs, one relying on arc spaces and one that goes through the inequality $\widetilde{\alpha}(f)\geq{\rm lct}(f,J_f^2)$, where $\widetilde{\alpha}(f)$ is the minimal exponent of $f$. In the case of a polynomial over $\overline{\mathbf{Q}}$, we also prove an analogue of this latter inequality, with $\widetilde{\alpha}(f)$ replaced by the motivic oscillation index ${\rm moi}(f)$. We also show a part of Igusa's strong monodromy conjecture, for poles larger than $-{\rm lct}(f,J_f^2)$. We end with a discussion of lct-maximal ideals: these are ideals $I$ with the property that ${\rm lct}(I)<{\rm lct}(J)$ for every $J$ with $I\subsetneq J$.


title: Synthetic Control As Online Linear Regression

abstract: This paper notes a simple connection between synthetic control and online learning. Specifically, we recognize synthetic control as an instance of Follow-The-Leader (FTL). Standard results in online convex optimization then imply that, even when outcomes are chosen by an adversary, synthetic control predictions of counterfactual outcomes for the treated unit perform almost as well as an oracle weighted average of control units' outcomes. Synthetic control on differenced data performs almost as well as oracle weighted difference-in-differences, potentially making it an attractive choice in practice. We argue that this observation further supports the use of synthetic control estimators in comparative case studies.


title: Weak-odd chromatic index of special digraph classes

abstract: Give a digraph $D=(V(D),A(D))$, let $\partial^+_D(v)=\{vw|w\in N^+_D(v)\}$ and $\partial^-_D(v)=\{uv|u\in N^-_D(v)\}$ be semi-cuts of $v$. A mapping $\varphi:A(D)\rightarrow [k]$ is called a weak-odd $k$-edge coloring of $D$ if it satisfies the condition: for each $v\in V(D)$, there is at least one color with an odd number of occurrences on each non-empty semi-cut of $v$. We call the minimum integer $k$ the weak-odd chromatic index of $D$. When limit to 2 colors, use $def(D)$ to denote the defect of $D$, the minimum number of vertices in $D$ at which the above condition is not satisfied. In this paper, we give a descriptive characterization about the weak-odd chromatic index and the defect of semicomplete digraphs and extended tournaments, which generalize results of tournaments to broader classes. And we initiated the study of weak-odd edge covering on digraphs.


title: Hidden local symmetry breaking in a kagome-lattice magnetic Weyl semimetal

abstract: Exploring the relationship between intriguing physical properties and structural complexity is a central topic in studying modern functional materials. Co$_{3}$Sn$_{2}$S$_{2}$, a new discovered kagome-lattice magnetic Weyl semimetal, has triggered intense interest owing to the intimate coupling between topological semimetallic states and peculiar magnetic properties. However, the origins of the magnetic phase separation and spin glass state below $T_{C}$ in this ordered compound are two unresolved yet important puzzles in understanding its magnetism. Here, we report the discovery of local symmetry breaking surprisingly co-emerges with the onset of ferromagnetic order in Co$_{3}$Sn$_{2}$S$_{2}$, by a combined use of neutron total scattering and half polarized neutron diffraction. The mismatch of local and average symmetries occurs below $T_{C}$, indicating that Co$_{3}$Sn$_{2}$S$_{2}$ evolves to an intrinsically lattice disordered system when the ferromagnetic order is established. The local symmetry breaking with intrinsic lattice disorder provides new understandings to the puzzling magnetic properties. Our density function theory calculation indicates that the local symmetry breaking is expected to reorient local ferromagnetic moments, unveiling the existence of the ferromagnetic instability associated with the lattice instability. Furthermore, DFT calculation unveils that the local symmetry breaking could affect the Weyl property by breaking mirror plane. Our findings highlight the fundamentally important role that the local symmetry breaking plays in advancing our understanding on the magnetic and topological properties in Co$_{3}$Sn$_{2}$S$_{2}$, which may draw the attention to explore the overlooked local symmetry breaking in Co$_{3}$Sn$_{2}$S$_{2}$, its derivatives, and more broadly in other topological Dirac/Weyl semimetals and kagome-lattice magnets.


title: CHEX: Multiversion Replay with Ordered Checkpoints

abstract: In scientific computing and data science disciplines, it is often necessary to share application workflows and repeat results. Current tools containerize application workflows, and share the resulting container for repeating results. These tools, due to containerization, do improve sharing of results. However, they do not improve the efficiency of replay. In this paper, we present the multiversion replay problem which arises when multiple versions of an application are containerized, and each version must be replayed to repeat results. To avoid executing each version separately, we develop CHEX, which checkpoints program state and determines when it is permissible to reuse program state across versions. It does so using system call-based execution lineage. Our capability to identify common computations across versions enables us to consider optimizing replay using an in-memory cache, based on a checkpoint-restore-switch system. We show the multiversion replay problem is NP-hard, and propose efficient heuristics for it. CHEX reduces overall replay time by sharing common computations but avoids storing a large number of checkpoints. We demonstrate that CHEX maintains lightweight package sharing, and improves the total time of multiversion replay by 50% on average.