title: Simulating the Evolution of Signaling Signatures during CART-Cell -- Tumor Cell Interactions

abstract: Immunotherapies have been proven to have significant therapeutic efficacy in the treatment of cancer. The last decade has seen adoptive cell therapies, such as chimeric antigen receptor T-cell (CART-cell) therapy, gain FDA approval against specific cancers. Additionally, there are numerous clinical trials ongoing investigating additional designs and targets. Nevertheless, despite the excitement and promising potential of CART-cell therapy, response rates to therapy vary greatly between studies, patients, and cancers. There remains an unmet need to develop computational frameworks that more accurately predict CART-cell function and clinical efficacy. Here we present a coarse-grained model simulated with logical rules that demonstrates the evolution of signaling signatures following the inter-action between CART-cells and tumor cells and allows for in silico based prediction of CART-cell functionality prior to experimentation.


title: On the valence of logharmonic polynomials

abstract: Investigating a problem posed by W. Hengartner (2000), we study the maximal valence (number of preimages of a prescribed point in the complex plane) of logharmonic polynomials, i.e., complex functions that take the form $f(z) = p(z) \overline{q(z)}$ of a product of an analytic polynomial $p(z)$ of degree $n$ and the complex conjugate of another analytic polynomial $q(z)$ of degree $m$. In the case $m=1$, we adapt an indirect technique utilizing anti-holomorphic dynamics to show that the valence is at most $3n-1$. This confirms a conjecture of Bshouty and Hengartner (2000). Using a purely algebraic method based on Sylvester resultants, we also prove a general upper bound for the valence showing that for each $n,m \geq 1$ the valence is at most $n^2+m^2$. This improves, for every choice of $n,m \geq 1$, the previously established upper bound $(n+m)^2$ based on Bezout's theorem. We also consider the more general setting of polyanalytic polynomials where we show that this latter result can be extended under a nondegeneracy assumption.


title: When does humoral memory enhance infection?

abstract: Antibodies and humoral memory are key components of the adaptive immune system. We consider and computationally model mechanisms by which humoral memory present at baseline might instead increase infection load; we refer to this effect as EI-HM (enhancement of infection by humoral memory). We first consider antibody dependent enhancement (ADE) in which antibody enhances the growth of the pathogen, typically a virus, and typically at intermediate "Goldilocks" levels of antibody. Our ADE model reproduces ADE in vitro and enhancement of infection in vivo from passive antibody transfer. But notably the simplest implementation of our ADE model never results in EI-HM. Adding complexity, by making the cross-reactive antibody much less neutralizing than the de novo generated antibody or by including a sufficiently strong non-antibody immune response, allows for ADE-mediated EI-HM. We next consider the possibility that cross-reactive memory causes EI-HM by crowding out a possibly superior de novo immune response. We show that, even without ADE, EI-HM can occur when the cross-reactive response is both less potent and "directly" (i.e. independently of infection load) suppressive with regard to the de novo response. In this case adding a non-antibody immune response to our computational model greatly reduces or completely eliminates EI-HM, which suggests that "crowding out" is unlikely to cause substantial EI-HM. Hence, our results provide examples in which simple models give qualitatively opposite results compared to models with plausible complexity. Our results may be helpful in interpreting and reconciling disparate experimental findings, especially from dengue, and for vaccination.


title: Neonatal Face and Facial Landmark Detection from Video Recordings

abstract: This paper explores automated face and facial landmark detection of neonates, which is an important first step in many video-based neonatal health applications, such as vital sign estimation, pain assessment, sleep-wake classification, and jaundice detection. Utilising three publicly available datasets of neonates in the clinical environment, 366 images (258 subjects) and 89 (66 subjects) were annotated for training and testing, respectively. Transfer learning was applied to two YOLO-based models, with input training images augmented with random horizontal flipping, photo-metric colour distortion, translation and scaling during each training epoch. Additionally, the re-orientation of input images and fusion of trained deep learning models was explored. Our proposed model based on YOLOv7Face outperformed existing methods with a mean average precision of 84.8% for face detection, and a normalised mean error of 0.072 for facial landmark detection. Overall, this will assist in the development of fully automated neonatal health assessment algorithms.


title: Graph and backward asymptotics of the tent map

abstract: The tent map family is arguably the simplest 1-parametric family of maps with non-trivial dynamics and it is still an active subject of research. In recent works the second author, jointly with J. Yorke, studied the graph and backward limits of S-unimodal maps. In this article we generalize those results to tent-like unimodal maps. By tent-like here we mean maps that share fundamental properties that characterize tent maps, namely unimodal maps without wandering intervals nor attracting cycles and whose graph has a finite number of nodes.


title: CRL+: A Novel Semi-Supervised Deep Active Contrastive Representation Learning-Based Text Classification Model for Insurance Data

abstract: Financial sector and especially the insurance industry collect vast volumes of text on a daily basis and through multiple channels (their agents, customer care centers, emails, social networks, and web in general). The information collected includes policies, expert and health reports, claims and complaints, results of surveys, and relevant social media posts. It is difficult to effectively extract label, classify, and interpret the essential information from such varied and unstructured material. Therefore, the Insurance Industry is among the ones that can benefit from applying technologies for the intelligent analysis of free text through Natural Language Processing (NLP). In this paper, CRL+, a novel text classification model combining Contrastive Representation Learning (CRL) and Active Learning is proposed to handle the challenge of using semi-supervised learning for text classification. In this method, supervised (CRL) is used to train a RoBERTa transformer model to encode the textual data into a contrastive representation space and then classify using a classification layer. This (CRL)-based transformer model is used as the base model in the proposed Active Learning mechanism to classify all the data in an iterative manner. The proposed model is evaluated using unstructured obituary data with objective to determine the cause of the death from the data. This model is compared with the CRL model and an Active Learning model with the RoBERTa base model. The experiment shows that the proposed method can outperform both methods for this specific task.


title: Learning Dynamical Systems by Leveraging Data from Similar Systems

abstract: We consider the problem of learning the dynamics of a linear system when one has access to data generated by an auxiliary system that shares similar (but not identical) dynamics, in addition to data from the true system. We use a weighted least squares approach, and provide finite sample error bounds of the learned model as a function of the number of samples and various system parameters from the two systems as well as the weight assigned to the auxiliary data. We show that the auxiliary data can help to reduce the intrinsic system identification error due to noise, at the price of adding a portion of error that is due to the differences between the two system models. We further provide a data-dependent bound that is computable when some prior knowledge about the systems, such as upper bounds on noise levels and model difference, is available. This bound can also be used to determine the weight that should be assigned to the auxiliary data during the model training stage.


title: Inefficiency of CFMs: hedging perspective and agent-based simulations

abstract: We investigate whether the fee income from trades on the CFM is sufficient for the liquidity providers to hedge away the exposure to market risk. We first analyse this problem through the lens of continuous-time financial mathematics and derive an upper bound for not-arbitrage fee income that would make CFM efficient and liquidity provision fair. We then evaluate our findings by performing multi-agent simulations by varying CFM fees, market volatility, and rate of arrival of liquidity takers. We observe that, on average, fee income generated from liquidity provision is insufficient to compensate for market risk.


title: On the stabilization of the topological complexity of graph braid groups

abstract: We establish the first nontrivial lower bound on the (higher) topological complexity of the unordered configuration spaces of a general graph. As an application, we show that, for most graphs, the topological complexity eventually stabilizes at its maximal possible value, a direct analogue of a stability phenomenon in the ordered setting first conjectured by Farber. We estimate the stable range in terms of the number of trivalent vertices.


title: Inequalities for sections and projections of convex bodies

abstract: This article belongs to the area of geometric tomography, which is the study of geometric properties of solids based on data about their sections and projections. We describe a new direction in geometric tomography where different volumetric results are considered in a more general setting, with volume replaced by an arbitrary measure. Surprisingly, such a general approach works for a number of volumetric results. In particular, we discuss the Busemann-Petty problem on sections of convex bodies for arbitrary measures and the slicing problem for arbitrary measures. We present generalizations of these questions to the case of functions. A number of generalizations of questions related to projections, such as the problem of Shephard, are also discussed as well as some questions in discrete tomography.


title: Addressing Systematics in the Traceback Age of the $\beta$ Pictoris Moving Group

abstract: We characterize the impact of several sources of systematic errors on the computation of the traceback age of the $\beta$ Pictoris Moving Group ($\beta$PMG). We find that uncorrected gravitational redshift and convective blueshift bias absolute radial velocity measurements by $\sim$ 0.6 kms${}^{-1}$, which leads to erroneously younger traceback ages by $\sim$ 2 Myr. Random errors on parallax, proper motion, and radial velocity measurements lead to an additional bias of $\sim$ 1.5 Myr on traceback ages. Contamination of astrometric and kinematic data by kinematic outliers and unresolved multiple systems in the full input sample of 76 members and candidates of $\beta$PMG also erroneously lowers traceback ages by ${\sim}$ 3 Myr. We apply our new numerical traceback analysis tool to a core sample of 25 carefully vetted members of $\beta$PMG using Gaia Data Release 3 (DR3) data products and other kinematic surveys. Our method yields a corrected age of 20.4 $\pm$ 2.5 Myr, bridging the gap between kinematic ages (11$-$19 Myr) and other age-dating methods, such as isochrones and lithium depletion boundary (20$-$26 Myr). We explore several association size metrics that can track the spatial extent of $\beta$PMG over time, and we determine that minimizing the variance along the heliocentric curvilinear coordinate $\xi^{\prime}$ (i.e., toward the Galactic Center) offers the least random and systematic errors, due to the wider UVW space velocity dispersion of members of $\beta$PMG along the U-axis, which tends to maximize the spatial growth of the association along the $\xi^{\prime}$-axis over time.


title: Symbolic Quantum Simulation with Quasimodo

abstract: The simulation of quantum circuits on classical computers is an important problem in quantum computing. Such simulation requires representations of distributions over very large sets of basis vectors, and recent work has used symbolic data-structures such as Binary Decision Diagrams (BDDs) for this purpose. In this tool paper, we present Quasimodo, an extensible, open-source Python library for symbolic simulation of quantum circuits. Quasimodo is specifically designed for easy extensibility to other backends. Quasimodo allows simulations of quantum circuits, checking properties of the outputs of quantum circuits, and debugging quantum circuits. It also allows the user to choose from among several symbolic data-structures -- both unweighted and weighted BDDs, and a recent structure called Context-Free-Language Ordered Binary Decision Diagrams (CFLOBDDs) -- and can be easily extended to support other symbolic data-structures.


title: One-parameter families of conformal mappings of the half-plane onto polygonal domains with several slits

abstract: Among various methods of finding accessory parameters in the Schwarz-Christoffel integrals, Kufarev's method, based on the Loewner differential equation, plays an important role. It is used for describing one-parameter families of functions that conformally map a canonical domain onto a polygon with a slit the endpoint of which moves along a polygonal line starting from a boundary point. We present a modification of Kufarev's method for the case of several slits, the lengths of which have depend of each other in a certain way. We justify the method and find a system of ODEs describing the dynamics of accessory parameters. We also present the results of numerical calculations which confirm the efficiency of our method.


title: Fuzzing Automatic Differentiation in Deep-Learning Libraries

abstract: Deep learning (DL) has attracted wide attention and has been widely deployed in recent years. As a result, more and more research efforts have been dedicated to testing DL libraries and frameworks. However, existing work largely overlooked one crucial component of any DL system, automatic differentiation (AD), which is the basis for the recent development of DL. To this end, we propose $\nabla$Fuzz, the first general and practical approach specifically targeting the critical AD component in DL libraries. Our key insight is that each DL library API can be abstracted into a function processing tensors/vectors, which can be differentially tested under various execution scenarios (for computing outputs/gradients with different implementations). We have implemented $\nabla$Fuzz as a fully automated API-level fuzzer targeting AD in DL libraries, which utilizes differential testing on different execution scenarios to test both first-order and high-order gradients, and also includes automated filtering strategies to remove false positives caused by numerical instability. We have performed an extensive study on four of the most popular and actively-maintained DL libraries, PyTorch, TensorFlow, JAX, and OneFlow. The result shows that $\nabla$Fuzz substantially outperforms state-of-the-art fuzzers in terms of both code coverage and bug detection. To date, $\nabla$Fuzz has detected 173 bugs for the studied DL libraries, with 144 already confirmed by developers (117 of which are previously unknown bugs and 107 are related to AD). Remarkably, $\nabla$Fuzz contributed 58.3% (7/12) of all high-priority AD bugs for PyTorch and JAX during a two-month period. None of the confirmed AD bugs were detected by existing fuzzers.


title: The X-ray synchrotron rims in Cassiopeia A narrow with energy

abstract: Some young supernova remnants exhibit thin filaments of X-ray synchrotron radiation coinciding with the forward shock due to accelerated electrons interacting with the local magnetic field. The two main models accounting for the radial brightness evolution of these filaments differ in their prediction of the narrowing (or not) of the filaments with increasing photon energy. In this paper, we report our observation of such a narrowing of the synchrotron filaments in Cassiopeia A at X-ray energies, and how this finding could help in understanding the mechanisms at stake in their formation. We used a new blind source separation method on the 1 Ms Chandra observation of Cassiopeia A, in order to obtain detailed and unpolluted images of the synchrotron emission in three energy bands. We then extracted the profiles of several filaments at the forward shock and the reverse shock to estimate and compare their widths. We find that there is indeed a narrowing with energy of the synchrotron filaments both at the forward and at the reverse shocks in Cassiopeia A. The energy dependency of this narrowing seems stronger at high energy, which is indicative of a damping effect, confirmed by radio observations.


title: Solving the Transmission Problem for Open Wave-Guides, I Fundamental Solutions and Integral Equations

abstract: We introduce a layer potential representation for the solution of the transmission problem defined by two dielectric channels, or open wave-guides, meeting along the straight-line interface, $\{x_1=0\}.$ The main observation is that the outgoing fundamental solution for the operator $\Delta +k_1^2+q(x_2),$ acting on functions defined in ${\mathbb R}^2,$ is easily constructed using the Fourier transform in the $x_1$-variable and the elementary theory of ordinary differential equations. These fundamental solutions can then be used to represent the solution to the transmission problem in half planes. The transmission boundary conditions lead to integral equations along the intersection of the half planes, which, in our normalization, is the $x_2$-axis. We show that, in appropriate Banach spaces, these integral equations are Fredholm equations of second kind, which are therefore generically solvable. We analyze the representation of the guided modes in our formulation.


title: Consider or Choose? The Role and Power of Consideration Sets

abstract: Consideration sets play a crucial role in discrete choice modeling, where customers are commonly assumed to go through a two-stage decision making process. Specifically, customers are assumed to form consideration sets in the first stage and then use a second-stage choice mechanism to pick the product with the highest utility from the consideration sets. Recent studies mostly aim to propose more powerful choice mechanisms based on advanced non-parametric models to improve prediction accuracy. In contrast, this paper takes a step back from exploring more complex second-stage choice mechanisms and instead focus on how effectively we can model customer choice relying only on the first-stage consideration set formation. To this end, we study a class of nonparametric choice models that is only specified by a distribution over consideration sets and has a bounded rationality interpretation. We denote it as the consideration set model. Intriguingly, we show that this class of choice models can be characterized by the axiom of symmetric demand cannibalization, which enables complete statistical identification. We further consider the model's downstream assortment planning as an application. We first present an exact description of the optimal assortment, proving that it is revenue-ordered based on the blocks defined by the consideration sets. Despite this compelling structure, we establish that the assortment optimization problem under this model is NP-hard even to approximate. This result shows that accounting for consideration sets in the model inevitably results in inapproximability in assortment planning, even though the consideration set model uses the simplest possible uniform second-stage choice mechanism. Finally, using a real-world dataset, we show the tremendous power of the first-stage consideration sets when modeling customers' decision-making processes.


title: MedDiff: Generating Electronic Health Records using Accelerated Denoising Diffusion Model

abstract: Due to patient privacy protection concerns, machine learning research in healthcare has been undeniably slower and limited than in other application domains. High-quality, realistic, synthetic electronic health records (EHRs) can be leveraged to accelerate methodological developments for research purposes while mitigating privacy concerns associated with data sharing. The current state-of-the-art model for synthetic EHR generation is generative adversarial networks, which are notoriously difficult to train and can suffer from mode collapse. Denoising Diffusion Probabilistic Models, a class of generative models inspired by statistical thermodynamics, have recently been shown to generate high-quality synthetic samples in certain domains. It is unknown whether these can generalize to generation of large-scale, high-dimensional EHRs. In this paper, we present a novel generative model based on diffusion models that is the first successful application on electronic health records. Our model proposes a mechanism to perform class-conditional sampling to preserve label information. We also introduce a new sampling strategy to accelerate the inference speed. We empirically show that our model outperforms existing state-of-the-art synthetic EHR generation methods.


title: BLO spaces associated with Laguerre polynomials expansions

abstract: In this paper we introduce spaces of $\textup{BLO}$-type related to Laguerre polynomial expansions. We consider the probability measure on $(0,\infty)$ defined by $d\gamma_\alpha(x)=\frac{2}{\Gamma(\alpha+1)}e^{-x^2}x^{2\alpha+1}dx$ with $\alpha>-\frac12$. For every $a>0$, the space $\textup{BLO}_a((0,\infty),\gamma_\alpha)$ consists of all those measurable functions defined on $(0,\infty)$ having bounded lower oscillation with respect to $\gamma_\alpha$ over an admissible family $\mathcal{B}_a$ of intervals in $(0,\infty)$. The space $\textup{BLO}_a((0,\infty),\gamma_\alpha)$ is a subspace of the space $\textup{BMO}_a((0,\infty),\gamma_\alpha)$ of bounded mean oscillation functions with respect to $\gamma_\alpha$ and $\mathcal{B}_a$. The natural $a$-local centered maximal function defined by $\gamma_\alpha$ is bounded from $\textup{BMO}_a((0,\infty),\gamma_\alpha)$ into $\textup{BLO}_a((0,\infty),\gamma_\alpha)$. We prove that the maximal operator, the $\rho$-variation and the oscillation operators associated with local truncations of the Riesz transforms in the Laguerre setting are bounded from $L^\infty((0,\infty),\gamma_\alpha)$ into $\textup{BLO}_a((0,\infty),\gamma_\alpha)$. Also, we obtain a similar result for the maximal operator of local truncations for spectral Laplace transform type multipliers.


title: On Computable Online Learning

abstract: We initiate a study of computable online (c-online) learning, which we analyze under varying requirements for "optimality" in terms of the mistake bound. Our main contribution is to give a necessary and sufficient condition for optimal c-online learning and show that the Littlestone dimension no longer characterizes the optimal mistake bound of c-online learning. Furthermore, we introduce anytime optimal (a-optimal) online learning, a more natural conceptualization of "optimality" and a generalization of Littlestone's Standard Optimal Algorithm. We show the existence of a computational separation between a-optimal and optimal online learning, proving that a-optimal online learning is computationally more difficult. Finally, we consider online learning with no requirements for optimality, and show, under a weaker notion of computability, that the finiteness of the Littlestone dimension no longer characterizes whether a class is c-online learnable with finite mistake bound. A potential avenue for strengthening this result is suggested by exploring the relationship between c-online and CPAC learning, where we show that c-online learning is as difficult as improper CPAC learning.


title: Mitigating Bias in Visual Transformers via Targeted Alignment

abstract: As transformer architectures become increasingly prevalent in computer vision, it is critical to understand their fairness implications. We perform the first study of the fairness of transformers applied to computer vision and benchmark several bias mitigation approaches from prior work. We visualize the feature space of the transformer self-attention modules and discover that a significant portion of the bias is encoded in the query matrix. With this knowledge, we propose TADeT, a targeted alignment strategy for debiasing transformers that aims to discover and remove bias primarily from query matrix features. We measure performance using Balanced Accuracy and Standard Accuracy, and fairness using Equalized Odds and Balanced Accuracy Difference. TADeT consistently leads to improved fairness over prior work on multiple attribute prediction tasks on the CelebA dataset, without compromising performance.


title: Sentiment analysis and opinion mining on educational data: A survey

abstract: Sentiment analysis AKA opinion mining is one of the most widely used NLP applications to identify human intentions from their reviews. In the education sector, opinion mining is used to listen to student opinions and enhance their learning-teaching practices pedagogically. With advancements in sentiment annotation techniques and AI methodologies, student comments can be labelled with their sentiment orientation without much human intervention. In this review article, (1) we consider the role of emotional analysis in education from four levels: document level, sentence level, entity level, and aspect level, (2) sentiment annotation techniques including lexicon-based and corpus-based approaches for unsupervised annotations are explored, (3) the role of AI in sentiment analysis with methodologies like machine learning, deep learning, and transformers are discussed, (4) the impact of sentiment analysis on educational procedures to enhance pedagogy, decision-making, and evaluation are presented. Educational institutions have been widely invested to build sentiment analysis tools and process their student feedback to draw their opinions and insights. Applications built on sentiment analysis of student feedback are reviewed in this study. Challenges in sentiment analysis like multi-polarity, polysemous, negation words, and opinion spam detection are explored and their trends in the research space are discussed. The future directions of sentiment analysis in education are discussed.


title: Kinodynamic Rapidly-exploring Random Forest for Rearrangement-Based Nonprehensile Manipulation

abstract: Rearrangement-based nonprehensile manipulation still remains as a challenging problem due to the high-dimensional problem space and the complex physical uncertainties it entails. We formulate this class of problems as a coupled problem of local rearrangement and global action optimization by incorporating free-space transit motions between constrained rearranging actions. We propose a forest-based kinodynamic planning framework to concurrently search in multiple problem regions, so as to enable global exploration of the most task-relevant subspaces, while facilitating effective switches between local rearranging actions. By interleaving dynamic horizon planning and action execution, our framework can adaptively handle real-world uncertainties. With extensive experiments, we show that our framework significantly improves the planning efficiency and manipulation effectiveness while being robust against various uncertainties.


title: Robust trajectory optimisation for transitions of tiltwing VTOL aircraft

abstract: We propose a method to generate robust and optimal trajectories for the transition of a tiltwing Vertical Take-Off and Landing (VTOL) aircraft leveraging concepts from convex optimisation, tube-based nonlinear Model Predictive Control (MPC) and Difference of Convex (DC) functions decomposition. The approach relies on computing DC decompositions of dynamic models in order to exploit convexity properties and develop a tractable robust optimisation that solves a sequence of convex programs converging to a local optimum of the trajectory generation problem. The algorithm developed is applied to an Urban Air Mobility case study. The resulting solutions are robust to approximation errors in dynamic models and provide safe trajectories for aggressive transition manoeuvres at constant altitude.


title: Disentangling Learning Representations with Density Estimation

abstract: Disentangled learning representations have promising utility in many applications, but they currently suffer from serious reliability issues. We present Gaussian Channel Autoencoder (GCAE), a method which achieves reliable disentanglement via flexible density estimation of the latent space. GCAE avoids the curse of dimensionality of density estimation by disentangling subsets of its latent space with the Dual Total Correlation (DTC) metric, thereby representing its high-dimensional latent joint distribution as a collection of many low-dimensional conditional distributions. In our experiments, GCAE achieves highly competitive and reliable disentanglement scores compared with state-of-the-art baselines.


title: Towards Model-Agnostic Federated Learning over Networks

abstract: We present a model-agnostic federated learning method for networks of heterogeneous data and models. The network structure reflects similarities between the (statistics of) local datasets and, in turn, their associated local("personal") models. Our method is an instance of empirical risk minimization, with the regularization term derived from the network structure of data. In particular, we require well-connected local models, forming clusters, to yield similar predictions on a common test set. The proposed method allows for a wide range of local models. The only restriction on these local models is that they allow for efficient implementation of regularized empirical risk minimization (training). For a wide range of models, such implementations are available in high-level programming libraries including scikit-learn, Keras or PyTorch.


title: Exotic Baryons in Hot Neutron Stars

abstract: We study the nuclear isentropic equation of state for a stellar matter composed of nucleons, hyperons, and $\Delta$-resonances. We investigate different snapshots of the evolution of a neutron star, from its birth as a lepton-rich protoneutron star in the aftermath of a supernova explosion to a lepton-poor regime when the star starts cooling to a catalyzed configuration. We use a relativistic model within the mean-field approximation to describe the hot stellar matter and adopt density-dependent couplings adjusted by the DDME2 parameterization. We use baryon-meson couplings for the spin-$1/2$ baryonic octet and spin-$3/2$ decuplet determined in a unified manner relying on $\text{SU}(6)$ and $\text{SU}(3)$ symmetry arguments. We observe that $\Lambda$ is the dominant exotic particle in the star at different entropies for both neutrino-free and neutrino-trapped stellar matter. For a fixed entropy, the inclusion of new particles (hyperons and/or delta resonances) in the stellar matter decreases the temperature. Also, an increase in entropy per baryon ($1\;\text{to}\; 2$) with decreasing lepton number density ($0.4\;\text{to}\; 0.2$) leads to an increase in stellar radii and a decrease in its mass due to neutrino diffusion. In the neutrino transparent matter, the radii decrease from entropy per baryon $2$ to $T\,=\,0$ without a significant change in stellar mass.


title: Spectrographs and Spectroscopists for the Sandage Test

abstract: The redshift drift is a small, dynamic change in the redshift of objects following the Hubble flow. Its measurement provides a direct, real-time, model-independent mapping of the expansion rate of the Universe. It is fundamentally different from other cosmological probes: instead of mapping our (present-day) past light-cone, it directly compares different past light-cones. Being independent of any assumptions on gravity, geometry or clustering, it directly tests the pillars of the Lambda-CDM paradigm. Recent theoretical studies have uncovered unique synergies with other cosmological probes, including the characterization of the physical properties of dark energy. At the time of the original proposal by Sandage (1962) the expected change in the redshift of objects at cosmological distances appeared to be exceedingly small for reasonable observing times and beyond technological capabilities. In the last decades progress in the spectrographs (e.g. ESPRESSO), in the collecting area of telescopes and in the samples of cosmic beacons, enabled by new datasets and new machine-learning-based selections, have drastically changed the situation, bringing the Redshift Drift Grail within reach. As a consequence, this measurement is a flagship objective of the Extremely Large Telescope (ELT), specifically of its high-resolution spectrograph, ANDES.


title: Machine learning detects multiplicity of the first stars in stellar archaeology data

abstract: In unveiling the nature of the first stars, the main astronomical clue is the elemental compositions of the second generation of stars, observed as extremely metal-poor (EMP) stars, in our Milky Way Galaxy. However, no observational constraint was available on their multiplicity, which is crucial for understanding early phases of galaxy formation. We develop a new data-driven method to classify observed EMP stars into mono- or multi-enriched stars with Support Vector Machines. We also use our own nucleosynthesis yields of core-collapse supernovae with mixing-fallback that can explain many of observed EMP stars. Our method predicts, for the first time, that $31.8\% \pm 2.3\%$ of 462 analyzed EMP stars are classified as mono-enriched. This means that the majority of EMP stars are likely multi-enriched, suggesting that the first stars were born in small clusters. Lower metallicity stars are more likely to be enriched by a single supernova, most of which have high carbon enhancement. We also find that Fe, Mg. Ca, and C are the most informative elements for this classification. In addition, oxygen is very informative despite its low observability. Our data-driven method sheds a new light on solving the mystery of the first stars from the complex data set of Galactic archaeology surveys.


title: Thermodynamics of blackbody radiation in nonlinear electrodynamics

abstract: We study the blackbody properties and the thermodynamic equilibrium quantities of a photon gas in the framework of nonlinear electrodynamics. In this vein, we take into account the photon propagation in a uniform external magnetic field in the weak field approximation, where an angular anisotropic energy density distribution appears in the frequency spectrum. The particular case when the photon propagates perpendicular to the background magnetic field is also discussed, which allows us to probe the strong field regime. We then derive a modified blackbody spectral distribution and the Stefan-Boltzmann law in this situation. Considerations about Wien's displacement law and the Rayleigh-Jeans formula are contemplated as well. Deviations from the thermodynamic quantities at thermal equilibrium such as energy, pressure, entropy, and heat capacity densities are obtained from the Helmholtz free energy. As an application, we study three nonlinear electrodynamics, namely, the Euler-Heisenberg, the generalized Born-Infeld, and the logarithmic electrodynamics. Possible implications on stellar systems with strong magnetic fields such as magnetars are discussed.

categories: cond-mat.stat-mech astro-ph.HE gr-qc hep-th


title: Channelformer: Attention based Neural Solution for Wireless Channel Estimation and Effective Online Training

abstract: In this paper, we propose an encoder-decoder neural architecture (called Channelformer) to achieve improved channel estimation for orthogonal frequency-division multiplexing (OFDM) waveforms in downlink scenarios. The self-attention mechanism is employed to achieve input precoding for the input features before processing them in the decoder. In particular, we implement multi-head attention in the encoder and a residual convolutional neural architecture as the decoder, respectively. We also employ a customized weight-level pruning to slim the trained neural network with a fine-tuning process, which reduces the computational complexity significantly to realize a low complexity and low latency solution. This enables reductions of up to 70\% in the parameters, while maintaining an almost identical performance compared with the complete Channelformer. We also propose an effective online training method based on the fifth generation (5G) new radio (NR) configuration for the modern communication systems, which only needs the available information at the receiver for online training. Using industrial standard channel models, the simulations of attention-based solutions show superior estimation performance compared with other candidate neural network methods for channel estimation.


title: Unsupervised Learning of Initialization in Deep Neural Networks via Maximum Mean Discrepancy

abstract: Despite the recent success of stochastic gradient descent in deep learning, it is often difficult to train a deep neural network with an inappropriate choice of its initial parameters. Even if training is successful, it has been known that the initial parameter configuration may negatively impact generalization. In this paper, we propose an unsupervised algorithm to find good initialization for input data, given that a downstream task is d-way classification. We first notice that each parameter configuration in the parameter space corresponds to one particular downstream task of d-way classification. We then conjecture that the success of learning is directly related to how diverse downstream tasks are in the vicinity of the initial parameters. We thus design an algorithm that encourages small perturbation to the initial parameter configuration leads to a diverse set of d-way classification tasks. In other words, the proposed algorithm ensures a solution to any downstream task to be near the initial parameter configuration. We empirically evaluate the proposed algorithm on various tasks derived from MNIST with a fully connected network. In these experiments, we observe that our algorithm improves average test accuracy across most of these tasks, and that such improvement is greater when the number of labelled examples is small.


title: Adaptive State-Dependent Diffusion for Derivative-Free Optimization

abstract: This paper develops and analyzes a stochastic derivative-free optimization strategy. A key feature is the state-dependent adaptive variance. We prove global convergence in probability with algebraic rate and give the quantitative results in numerical examples. A striking fact is that convergence is achieved without explicit information of the gradient and even without comparing different objective function values as in established methods such as the simplex method and simulated annealing. It can otherwise be compared to annealing with state-dependent temperature.


title: Thermodynamic properties of the noncommutative quantum Hall effect with anomalous magnetic moment

abstract: In this paper, we study the thermodynamic properties of the noncommutative quantum Hall effect (NCQHE) with anomalous magnetic moment (AMM) for both relativistic and nonrelativistic cases in the high temperatures regime. Thus, we use the canonical ensemble for a set of $N$-particles in contact with a thermal bath. Next, we explicitly determine the thermodynamic properties of our interest, namely: the Helmholtz free energy, the entropy, the mean energy, and the heat capacity. In order to perform the calculations, we work with the Euler-MacLaurin formula to construct the partition function of the system. In that way, we plotted the graphs of thermodynamic properties as a function of temperature for six different values of the magnetic field and of the NC parameters. As a result, we note that the Helmholtz free energy decreases with the temperature, increases with the NC parameters, and can decrease or increase for certain values of the magnetic field, white that the entropy increases with the temperature, decreases with the NC parameters, and can decrease or increase for certain values of the magnetic field. Besides, the mean energy increases linearly with the temperature and its values for the relativistic case are twice of the nonrelativistic case, consequently, the heat capacity for the relativistic case is twice of the nonrelativistic case, where both are constants, and therefore, satisfying the so-called Dulong-Petit law. Finally, we also verify that there is no influence of the AMM on the thermodynamic properties of the system.


title: Bandwidth Control and Symmetry Breaking in a Mott-Hubbard Correlated Metal

abstract: In Mott materials strong electron correlation yields a spectrum of complex electronic structures. Recent synthesis advancements open realistic opportunities for harnessing Mott physics to design transformative devices. However, a major bottleneck in realizing such devices remains the lack of control over the electron correlation strength. This stems from the complexity of the electronic structure, which often veils the basic mechanisms underlying the correlation strength. Here, we present control of the correlation strength by tuning the degree of orbital overlap using picometer-scale lattice engineering. We illustrate how bandwidth control and concurrent symmetry breaking can govern the electronic structure of a correlated $SrVO_3$ model system. We show how tensile and compressive biaxial strain oppositely affect the $SrVO_3$ in-plane and out-of-plane orbital occupancy, resulting in the partial alleviation of the orbital degeneracy. We derive and explain the spectral weight redistribution under strain and illustrate how high tensile strain drives the system towards a Mott insulating state. Implementation of such concepts will drive correlated electron phenomena closer towards new solid state devices and circuits. These findings therefore pave the way for understanding and controlling electron correlation in a broad range of functional materials, driving this powerful resource for novel electronics closer towards practical realization.


title: Measuring the Privacy Leakage via Graph Reconstruction Attacks on Simplicial Neural Networks (Student Abstract)

abstract: In this paper, we measure the privacy leakage via studying whether graph representations can be inverted to recover the graph used to generate them via graph reconstruction attack (GRA). We propose a GRA that recovers a graph's adjacency matrix from the representations via a graph decoder that minimizes the reconstruction loss between the partial graph and the reconstructed graph. We study three types of representations that are trained on the graph, i.e., representations output from graph convolutional network (GCN), graph attention network (GAT), and our proposed simplicial neural network (SNN) via a higher-order combinatorial Laplacian. Unlike the first two types of representations that only encode pairwise relationships, the third type of representation, i.e., SNN outputs, encodes higher-order interactions (e.g., homological features) between nodes. We find that the SNN outputs reveal the lowest privacy-preserving ability to defend the GRA, followed by those of GATs and GCNs, which indicates the importance of building more private representations with higher-order node information that could defend the potential threats, such as GRAs.


title: Near-Optimal Adversarial Reinforcement Learning with Switching Costs

abstract: Switching costs, which capture the costs for changing policies, are regarded as a critical metric in reinforcement learning (RL), in addition to the standard metric of losses (or rewards). However, existing studies on switching costs (with a coefficient $\beta$ that is strictly positive and is independent of $T$) have mainly focused on static RL, where the loss distribution is assumed to be fixed during the learning process, and thus practical scenarios where the loss distribution could be non-stationary or even adversarial are not considered. While adversarial RL better models this type of practical scenarios, an open problem remains: how to develop a provably efficient algorithm for adversarial RL with switching costs? This paper makes the first effort towards solving this problem. First, we provide a regret lower-bound that shows that the regret of any algorithm must be larger than $\tilde{\Omega}( ( H S A )^{1/3} T^{2/3} )$, where $T$, $S$, $A$ and $H$ are the number of episodes, states, actions and layers in each episode, respectively. Our lower bound indicates that, due to the fundamental challenge of switching costs in adversarial RL, the best achieved regret (whose dependency on $T$ is $\tilde{O}(\sqrt{T})$) in static RL with switching costs (as well as adversarial RL without switching costs) is no longer achievable. Moreover, we propose two novel switching-reduced algorithms with regrets that match our lower bound when the transition function is known, and match our lower bound within a small factor of $\tilde{O}( H^{1/3} )$ when the transition function is unknown. Our regret analysis demonstrates the near-optimal performance of them.


title: A Near-Optimal Algorithm for Safe Reinforcement Learning Under Instantaneous Hard Constraints

abstract: In many applications of Reinforcement Learning (RL), it is critically important that the algorithm performs safely, such that instantaneous hard constraints are satisfied at each step, and unsafe states and actions are avoided. However, existing algorithms for ''safe'' RL are often designed under constraints that either require expected cumulative costs to be bounded or assume all states are safe. Thus, such algorithms could violate instantaneous hard constraints and traverse unsafe states (and actions) in practice. Therefore, in this paper, we develop the first near-optimal safe RL algorithm for episodic Markov Decision Processes with unsafe states and actions under instantaneous hard constraints and the linear mixture model. It not only achieves a regret $\tilde{O}(\frac{d H^3 \sqrt{dK}}{\Delta_c})$ that tightly matches the state-of-the-art regret in the setting with only unsafe actions and nearly matches that in the unconstrained setting, but is also safe at each step, where $d$ is the feature-mapping dimension, $K$ is the number of episodes, $H$ is the number of steps in each episode, and $\Delta_c$ is a safety-related parameter. We also provide a lower bound $\tilde{\Omega}(\max\{dH \sqrt{K}, \frac{H}{\Delta_c^2}\})$, which indicates that the dependency on $\Delta_c$ is necessary. Further, both our algorithm design and regret analysis involve several novel ideas, which may be of independent interest.


title: Efficient Planning in Combinatorial Action Spaces with Applications to Cooperative Multi-Agent Reinforcement Learning

abstract: A practical challenge in reinforcement learning are combinatorial action spaces that make planning computationally demanding. For example, in cooperative multi-agent reinforcement learning, a potentially large number of agents jointly optimize a global reward function, which leads to a combinatorial blow-up in the action space by the number of agents. As a minimal requirement, we assume access to an argmax oracle that allows to efficiently compute the greedy policy for any Q-function in the model class. Building on recent work in planning with local access to a simulator and linear function approximation, we propose efficient algorithms for this setting that lead to polynomial compute and query complexity in all relevant problem parameters. For the special case where the feature decomposition is additive, we further improve the bounds and extend the results to the kernelized setting with an efficient algorithm.


title: ER network heterogeneity guides diffusive transport and kinetics

abstract: The endoplasmic reticulum (ER) is a dynamic network of interconnected sheets and tubules that orchestrates the distribution of lipids, ions, and proteins throughout the cell. The impact of its complex, dynamic morphology on its function as an intracellular transport hub remains poorly understood. To elucidate the functional consequences of ER network structure and dynamics, we quantify how the heterogeneity of the peripheral ER in COS7 cells affects diffusive protein transport. In vivo imaging of photoactivated ER membrane proteins demonstrates their non-uniform spreading to adjacent regions, in a manner consistent with simulations of diffusing particles on extracted network structures. Using a minimal network model to represent tubule rearrangements, we demonstrate that ER network dynamics are sufficiently slow to have little effect on diffusive protein transport. Furthermore, stochastic simulations reveal a novel consequence of ER network heterogeneity: the existence of 'hot spots' where sparse diffusive reactants are more likely to find one another. Intriguingly, ER exit sites are disproportionately found in these highly accessible regions. Combining in vivo experiments with analytic calculations, quantitative image analysis, and computational modeling, we demonstrate how structure guides diffusive protein transport and reactions in the ER.

categories: physics.bio-ph cond-mat.soft q-bio.QM q-bio.SC


title: Parallel Derandomization for Coloring

abstract: Graph coloring problems are among the most fundamental problems in parallel and distributed computing, and have been studied extensively in both settings. In this context, designing efficient deterministic algorithms for these problems has been found particularly challenging. In this work we consider this challenge, and design a novel framework for derandomizing algorithms for coloring-type problems in the Massively Parallel Computation (MPC) model with sublinear space. We give an application of this framework by showing that a recent $(degree+1)$-list coloring algorithm by Halldorsson et al. (STOC'22) in the LOCAL model of distributed computation can be translated to the MPC model and efficiently derandomized. Our algorithm runs in $O(\log \log \log n)$ rounds, which matches the complexity of the state of the art algorithm for the $(\Delta + 1)$-coloring problem.


title: Et Tu Certifications: Robustness Certificates Yield Better Adversarial Examples

abstract: In guaranteeing the absence of adversarial examples in an instance's neighbourhood, certification mechanisms play an important role in demonstrating neural net robustness. In this paper, we ask if these certifications can compromise the very models they help to protect? Our new \emph{Certification Aware Attack} exploits certifications to produce computationally efficient norm-minimising adversarial examples $74 \%$ more often than comparable attacks, while reducing the median perturbation norm by more than $10\%$. While these attacks can be used to assess the tightness of certification bounds, they also highlight that releasing certifications can paradoxically reduce security.


title: Covariate Adjustment in Experiments with Matched Pairs

abstract: This paper studies inference on the average treatment effect in experiments in which treatment status is determined according to "matched pairs" and it is additionally desired to adjust for observed, baseline covariates to gain further precision. By a "matched pairs" design, we mean that units are sampled i.i.d. from the population of interest, paired according to observed, baseline covariates and finally, within each pair, one unit is selected at random for treatment. Importantly, we presume that not all observed, baseline covariates are used in determining treatment assignment. We study a broad class of estimators based on a "doubly robust" moment condition that permits us to study estimators with both finite-dimensional and high-dimensional forms of covariate adjustment. We find that estimators with finite-dimensional, linear adjustments need not lead to improvements in precision relative to the unadjusted difference-in-means estimator. This phenomenon persists even if the adjustments are interacted with treatment; in fact, doing so leads to no changes in precision. However, gains in precision can be ensured by including fixed effects for each of the pairs. Indeed, we show that this adjustment is the "optimal" finite-dimensional, linear adjustment. We additionally study two estimators with high-dimensional forms of covariate adjustment based on the LASSO. For each such estimator, we show that it leads to improvements in precision relative to the unadjusted difference-in-means estimator and also provide conditions under which it leads to the "optimal" nonparametric, covariate adjustment. A simulation study confirms the practical relevance of our theoretical analysis, and the methods are employed to reanalyze data from an experiment using a "matched pairs" design to study the effect of macroinsurance on microenterprise.


title: The failure of semiclassical approach in the dissipative fully-connected Ising model

abstract: We solve the fully-connected Ising model in the presence of dissipation and time-periodic field, with the corresponding Lindblad equation having a time-periodic Liouvillian. The dynamics of the magnetizations is studied by using both the semiclassical approach and the numerical simulation with the help of permutation symmetry. The semiclassical approach shows a transition from the periodic response for small field amplitude to the chaotic dynamics for large amplitude. The trajectory of the magnetizations and the Lyapunov exponents are calculated, which support the existence of a chaotic phase. But in the exact numerical simulation, the response is periodic for both small and large amplitude. The scaling analysis of Floquet Liouvillian spectrum confirms the periodic response in the thermodynamic limit. The semiclassical approximation is found to fail as the field amplitude is large.


title: On the relative isoperimetric problem for the cube

abstract: In this article, we solve the relative isoperimetric problem in $[0,1]^3$ for subsets whose boundaries are in the union of a finite number of coordinate-parallel hyperplanes. Up to isometries of the cube or sets of measure $0$, the minimizers are of the form $[0,\epsilon]^3$, $[0,\epsilon]^2 \times [0,1]$, or $[0,\epsilon] \times [0,1]^2$ for some $\epsilon > 0$. This should be compared to the conjectured minimizers for the unconstrained relative isoperimetric problem in $[0,1]^3$, which are (up to isometries and sets of measure $0$) of the form $\left( B^3(\epsilon) \right) \cap [0,1]^3$, $\left( B^2(\epsilon) \times [0,1] \right) \cap [0,1]^3$, or $[0,\epsilon] \times [0,1]^2$ for some $\epsilon > 0$. Here, $B^k(\epsilon)$ is the closed ball in $\mathbb{R}^k$ of radius $\epsilon$ centered at the origin.


title: Privacy-Preserving Representation Learning for Text-Attributed Networks with Simplicial Complexes

abstract: Although recent network representation learning (NRL) works in text-attributed networks demonstrated superior performance for various graph inference tasks, learning network representations could always raise privacy concerns when nodes represent people or human-related variables. Moreover, standard NRLs that leverage structural information from a graph proceed by first encoding pairwise relationships into learned representations and then analysing its properties. This approach is fundamentally misaligned with problems where the relationships involve multiple points, and topological structure must be encoded beyond pairwise interactions. Fortunately, the machinery of topological data analysis (TDA) and, in particular, simplicial neural networks (SNNs) offer a mathematically rigorous framework to learn higher-order interactions between nodes. It is critical to investigate if the representation outputs from SNNs are more vulnerable compared to regular representation outputs from graph neural networks (GNNs) via pairwise interactions. In my dissertation, I will first study learning the representations with text attributes for simplicial complexes (RT4SC) via SNNs. Then, I will conduct research on two potential attacks on the representation outputs from SNNs: (1) membership inference attack, which infers whether a certain node of a graph is inside the training data of the GNN model; and (2) graph reconstruction attacks, which infer the confidential edges of a text-attributed network. Finally, I will study a privacy-preserving deterministic differentially private alternating direction method of multiplier to learn secure representation outputs from SNNs that capture multi-scale relationships and facilitate the passage from local structure to global invariant features on text-attributed networks.


title: SF-SGL: Solver-Free Spectral Graph Learning from Linear Measurements

abstract: This work introduces a highly-scalable spectral graph densification framework (SGL) for learning resistor networks with linear measurements, such as node voltages and currents. We show that the proposed graph learning approach is equivalent to solving the classical graphical Lasso problems with Laplacian-like precision matrices. We prove that given $O(\log N)$ pairs of voltage and current measurements, it is possible to recover sparse $N$-node resistor networks that can well preserve the effective resistance distances on the original graph. In addition, the learned graphs also preserve the structural (spectral) properties of the original graph, which can potentially be leveraged in many circuit design and optimization tasks. To achieve more scalable performance, we also introduce a solver-free method (SF-SGL) that exploits multilevel spectral approximation of the graphs and allows for a scalable and flexible decomposition of the entire graph spectrum (to be learned) into multiple different eigenvalue clusters (frequency bands). Such a solver-free approach allows us to more efficiently identify the most spectrally-critical edges for reducing various ranges of spectral embedding distortions. Through extensive experiments for a variety of real-world test cases, we show that the proposed approach is highly scalable for learning sparse resistor networks without sacrificing solution quality. We also introduce a data-driven EDA algorithm for vectorless power/thermal integrity verifications to allow estimating worst-case voltage/temperature (gradient) distributions across the entire chip by leveraging a few voltage/temperature measurements.


title: Structural changes induced by electric currents in a single crystal of Pr$_2$CuO$_4$

abstract: We demonstrate a novel approach to the structural and electronic property modification of perovskites, focusing on Pr$_2$CuO$_4$, an undoped parent compound of a class of electron-doped copper-oxide superconductors. Currents were passed parallel or perpendicular to the copper-oxygen layers with the voltage ramped up until a rapid drop in the resistivity was achieved, a process referred to as "flash". The current was then further increased tenfold in current-control mode. This state was quenched by immersion into liquid nitrogen. Flash can drive many compounds into different atomic structures with new properties, whereas the quench freezes them into a long-lived state. Single-crystal neutron diffraction of as-grown and modified Pr$_2$CuO$_4$ revealed a $\sqrt{10}$x$\sqrt{10}$ superlattice due to oxygen-vacancy order. The diffraction peak intensities of the superlattice of the modified sample were significantly enhanced relative to the pristine sample. Raman-active phonons in the modified sample were considerably sharper. Measurements of electrical resistivity, magnetization and two-magnon Raman scattering indicate that the modification affected only the Pr-O layers, but not the Cu-O planes. These results point to enhanced oxygen-vacancy order in the modified samples well beyond what can be achieved without passing electrical current. Our work opens a new avenue toward electric field/quench control of structure and properties of layered perovskite oxides.


title: Machine Learning Capability: A standardized metric using case difficulty with applications to individualized deployment of supervised machine learning

abstract: Model evaluation is a critical component in supervised machine learning classification analyses. Traditional metrics do not currently incorporate case difficulty. This renders the classification results unbenchmarked for generalization. Item Response Theory (IRT) and Computer Adaptive Testing (CAT) with machine learning can benchmark datasets independent of the end-classification results. This provides high levels of case-level information regarding evaluation utility. To showcase, two datasets were used: 1) health-related and 2) physical science. For the health dataset a two-parameter IRT model, and for the physical science dataset a polytonomous IRT model, was used to analyze predictive features and place each case on a difficulty continuum. A CAT approach was used to ascertain the algorithms' performance and applicability to new data. This method provides an efficient way to benchmark data, using only a fraction of the dataset (less than 1%) and 22-60x more computationally efficient than traditional metrics. This novel metric, termed Machine Learning Capability (MLC) has additional benefits as it is unbiased to outcome classification and a standardized way to make model comparisons within and across datasets. MLC provides a metric on the limitation of supervised machine learning algorithms. In situations where the algorithm falls short, other input(s) are required for decision-making.


title: Catch Planner: Catching High-Speed Targets in the Flight

abstract: Catching high-speed targets in the flight is a complex and typical highly dynamic task. In this paper, we propose Catch Planner, a planning-with-decision scheme for catching. For sequential decision making, we propose a policy search method based on deep reinforcement learning. In order to make catching adaptive and flexible, we propose a trajectory optimization method to jointly optimize the highly coupled catching time and terminal state while considering the dynamic feasibility and safety. We also propose a flexible constraint transcription method to catch targets at any reasonable attitude and terminal position bias. The proposed Catch Planner provides a new paradigm for the combination of learning and planning and is integrated on the quadrotor designed by ourselves, which runs at 100hz on the onboard computer. Extensive experiments are carried out in real and simulated scenes to verify the robustness of the proposed method and its expansibility when facing a variety of high-speed flying targets.


title: The Radio to GeV Afterglow of GRB 221009A

abstract: GRB 221009A ($z=0.151$) is one of the closest known long $\gamma$-ray bursts (GRBs). Its extreme brightness across all electromagnetic wavelengths provides an unprecedented opportunity to study a member of this still-mysterious class of transients in exquisite detail. We present multi-wavelength observations of this extraordinary event, spanning 15 orders of magnitude in photon energy from radio to $\gamma$-rays. We find that the data can be partially explained by a forward shock (FS) from a highly-collimated relativistic jet interacting with a low-density wind-like medium. Under this model, the jet's beaming-corrected kinetic energy ($E_K \sim 4\times10^{50}$ erg) is typical for the GRB population. The radio and mm data provide strong limiting constraints on the FS model, but require the presence of an additional emission component. From equipartition arguments, we find that the radio emission is likely produced by a small amount of mass ($\lesssim6\times10^{-7} M_\odot$) moving relativistically ($\Gamma\gtrsim9$) with a large kinetic energy ($\gtrsim10^{49}$ erg). However, the temporal evolution of this component does not follow prescriptions for synchrotron radiation from a single power-law distribution of electrons (e.g. in a reverse shock or two-component jet), or a thermal electron population, perhaps suggesting that one of the standard assumptions of afterglow theory is violated. GRB 221009A will likely remain detectable with radio telescopes for years to come, providing a valuable opportunity to track the full lifecycle of a powerful relativistic jet.


title: Verification of Distributed Artificial Intelligence Systems in Bioinformatics

abstract: Software is a great enabler for a number of projects that otherwise would be impossible to perform. Such projects include Space Exploration, Weather Modeling, Genome Projects, and many others. It is critical that software aiding these projects does what it is expected to do. In the terminology of software engineering, software that corresponds to requirements, that is does what it is expected to do is called correct. Checking the correctness of software has been the focus of a great deal of research in the area of software engineering. Practitioners in the field in which software is applied quite often do not assign much value to checking this correctness. Yet, as software systems become larger, potentially combined with distributed subsystems written by different authors, such verification becomes even more important. Concurrent, distributed systems are prone to dangerous errors due to different speeds of execution of their components such as deadlocks, race conditions, or violation of project-specific properties. This project describes an application of a static analysis method called model checking to verification of a distributed system for the Bioinformatics process. In it, we evaluate the efficiency of the model checking approach to the verification of combined processes with an increasing number of concurrently executed steps. We show that our experimental results correspond to analytically derived expectations. We also highlight the importance of static analysis to combined processes in the Bioinformatics field.


title: Cosmological gravitational particle production of massive spin-2 particles

abstract: The phenomenon of cosmological gravitational particle production (CGPP) is expected to occur during the period of inflation and the transition into a hot big bang cosmology. Particles may be produced even if they only couple directly to gravity, and so CGPP provides a natural explanation for the origin of dark matter. In this work we study the gravitational production of massive spin-2 particles assuming two different couplings to matter. We evaluate the full system of mode equations, including the helicity-0 modes, and by solving them numerically we calculate the spectrum and abundance of massive spin-2 particles that results from inflation on a hilltop potential. We conclude that CGPP might provide a viable mechanism for the generation of massive spin-2 particle dark matter during inflation, and we identify the favorable region of parameter space in terms of the spin-2 particle's mass and the reheating temperature. As a secondary product of our work, we identify the conditions under which such theories admit ghost or gradient instabilities, and we thereby derive a generalization of the Higuchi bound to Friedmann-Robertson-Walker (FRW) spacetimes.


title: The Re-Label Method For Data-Centric Machine Learning

abstract: In industry deep learning application, our manually labeled data has a certain number of noisy data. To solve this problem and achieve more than 90 score in dev dataset, we present a simple method to find the noisy data and re-label the noisy data by human, given the model predictions as references in human labeling. In this paper, we illustrate our idea for a broad set of deep learning tasks, includes classification, sequence tagging, object detection, sequence generation, click-through rate prediction. The dev dataset evaluation results and human evaluation results verify our idea.


title: Second order fractional mean-field SDEs with singular kernels and measure initial data

abstract: In this paper we establish the local and global well-posedness of weak and strong solutions to second order fractional mean-field SDEs with singular/distribution interaction kernels and measure initial value, where the kernel can be Newton or Coulomb potential, Riesz potential, Biot-Savart law, etc. Moreover, we also show the stability, smoothness and the short time singularity and large time decay estimates of the distribution density. Our results reveal a phenomenon that for {\it nonlinear} mean-field equations, the regularity of the initial distribution could balance the singularity of the kernel. The precise relationship between the singularity of kernels and the regularity of initial values are calculated, which belongs to the subcritical regime in the scaling sense. In particular, our results provide a microscopic probabilistic explanation and establish a unified treatment for many physical models such as the fractional Vlasov-Poisson-Fokker-Planck system, the vorticity formulation of 2D-fractal Navier-Stokes equations, surface quasi-geostrophic models, fractional porous medium equation with viscosity, etc.


title: The SWAP Filter: A Simple Azimuthally Varying Radial Filter for Wide-Field EUV Solar Images

abstract: We present the SWAP Filter: an azimuthally varying, radial normalizing filter specifically developed for EUV images of the solar corona, named for the Sun Watcher with Active Pixels and Image Processing (SWAP) instrument on the Project for On-Board Autonomy 2 spacecraft. We discuss the origins of our technique, its implementation and key user-configurable parameters, and highlight its effects on data via a series of examples. We discuss the filter's strengths in a data environment in which wide field-of-view observations that specifically target the low signal-to-noise middle corona are newly available and expected to grow in the coming years.


title: Understanding the role of entropy in high entropy oxides

abstract: The field of high entropy oxides (HEOs) flips traditional materials science paradigms on their head by seeking to understand what properties arise in the presence of profound configurational disorder. This disorder, which originates from multiple elements sharing a single lattice site, can take on a kaleidoscopic character due to the vast numbers of possible elemental combinations. High configurational disorder appears to imbue some HEOs with functional properties that far surpass their non-disordered analogs. While experimental discoveries abound, efforts to characterize the true magnitude of the configurational entropy and understand its role in stabilizing new phases and generating superior functional properties have lagged behind. Understanding the role of configurational disorder in existing HEOs is the crucial link to unlocking the rational design of new HEOs with targeted properties. In this Perspective, we attempt to establish a framework for articulating and beginning to address these questions in pursuit of a deeper understanding of the true role of entropy in HEOs.


title: Optimized Hybrid Focal Margin Loss for Crack Segmentation

abstract: Many loss functions have been derived from cross-entropy loss functions such as large-margin softmax loss and focal loss. The large-margin softmax loss makes the classification more rigorous and prevents overfitting. The focal loss alleviates class imbalance in object detection by down-weighting the loss of well-classified examples. Recent research has shown that these two loss functions derived from cross entropy have valuable applications in the field of image segmentation. However, to the best of our knowledge, there is no unified formulation that combines these two loss functions so that they can not only be transformed mutually, but can also be used to simultaneously address class imbalance and overfitting. To this end, we subdivide the entropy-based loss into the regularizer-based entropy loss and the focal-based entropy loss, and propose a novel optimized hybrid focal loss to handle extreme class imbalance and prevent overfitting for crack segmentation. We have evaluated our proposal in comparison with three crack segmentation datasets (DeepCrack-DB, CRACK500 and our private PanelCrack dataset). Our experiments demonstrate that the focal margin component can significantly increase the IoU of cracks by 0.43 on DeepCrack-DB and 0.44 on our PanelCrack dataset, respectively.


title: Detection of Gravitational memory effect in LISA using triggers from ground-based detectors

abstract: The LIGO-Virgo-Kagra (LVK) collaboration has detected gravitational waves from 90 Compact Binary Coalescences. In addition to fortifying the linearized theory of General Relativity (GR), the statistical ensemble of detections also provides prospects of detecting nonlinear effects predicted by GR, one such prediction being the nonlinear gravitational memory effect. For detected stellar and intermediate mass compact binaries, the induced strain from the memory effect is one or two orders below the detector noise background. Additionally, since most of the energy is radiated at merger the strain induced by the memory effect resembles a step function at the merger time. These facts motivate the idea of coherently stacking up data streams from recorded GW events at these merger times so that the cumulative memory strain is detected with a sufficient SNR. GW detectors essentially record the integrated strain response at time scales of the round trip light travel time, making future space-based long arm interferometers like LISA ideal for detecting the memory effect at low frequencies. In this paper, we propose a method that uses the event catalog of ground-based detectors and searches for corresponding memory strains in the LISA data stream. Given LVK's O3 science run catalog, we use scaling arguments and assumptions on the source population models to estimate the run time required for LISA to accumulate a memory SNR of 5, using triggers from current and future ground-based detectors. Finally, we extend these calculations for using beyond LISA missions like ALIA, AMIGO, and Folkner to detect the gravitational memory effect. The results for LISA indicate a possible detection of the memory effect within the 10 year LISA mission lifetime and the corresponding results for beyond LISA missions are even more promising.


title: Designing Workflows for Materials Characterization

abstract: Experimental science is enabled by the combination of synthesis, imaging, and functional characterization. Synthesis of a new material is typically followed by a set of characterization methods aiming to provide feedback for optimization or discover fundamental mechanisms. However, the sequence of synthesis and characterization methods and their interpretation, or research workflow, has traditionally been driven by human intuition and is highly domain specific. Here we explore concepts of scientific workflows that emerge at the interface between theory, characterization, and imaging. We discuss the criteria by which these workflows can be constructed for special cases of multi-resolution structural imaging and structural and functional characterization. Some considerations for theory-experiment workflows are provided. We further pose that the emergence of user facilities and cloud labs disrupt the classical progression from ideation, orchestration, and execution stages of workflow development and necessitate development of universal frameworks for workflow design, including universal hyper-languages describing laboratory operation, reward functions and their integration between domains, and policy development for workflow optimization. These tools will enable knowledge-based workflow optimization, enable lateral instrumental networks, sequential and parallel orchestration of characterization between dissimilar facilities, and empower distributed research.


title: FDD Massive MIMO Without CSI Feedback

abstract: Transmitter channel state information (CSIT) is indispensable for the spectral efficiency gains offered by massive multiple-input multiple-output (MIMO) systems. In a frequency-division-duplexing (FDD) massive MIMO system, CSIT is typically acquired through downlink channel estimation and user feedback, but as the number of antennas increases, the overhead for CSI training and feedback per user grows, leading to a decrease in spectral efficiency. In this paper, we show that, using uplink pilots in FDD, the downlink sum spectral efficiency gain with perfect downlink CSIT is achievable when the number of antennas at a base station is infinite under some mild channel conditions. The key idea showing our result is the mean squared error-optimal downlink channel reconstruction method using uplink pilots, which exploits the geometry reciprocity of uplink and downlink channels. We also present a robust downlink precoding method harnessing the reconstructed channel with the error covariance matrix. Our system-level simulations show that our proposed precoding method can attain comparable sum spectral efficiency to zero-forcing precoding with perfect downlink CSIT, without CSI training and feedback.


title: Balanced Allocations in Batches: The Tower of Two Choices

abstract: In balanced allocations, the goal is to place $m$ balls into $n$ bins, so as to minimize the gap (difference of max to average load). The One-Choice process places each ball to a bin sampled independently and uniformly at random. The Two-Choice process places balls in the least loaded of two sampled bins. Finally, the $(1+\beta)$-process mixes these processes, meaning each ball is allocated using Two-Choice with probability $\beta\in(0,1)$, and using One-Choice otherwise. Despite Two-Choice being optimal in the sequential setting, it has been observed in practice that it does not perform well in a parallel environment, where load information may be outdated. Following [BCEFN12], we study such a parallel setting where balls are allocated in batches of size $b$, and balls within the same batch are allocated with the same strategy and based on the same load information. For small batch sizes $b\in[n,n\log n]$, it was shown in [LS22a] that Two-Choice achieves an asymptotically optimal gap among all processes with a constant number of samples. In this work, we focus on larger batch sizes $b\in[n\log n,n^3]$. It was proved in [LS22c] that Two-Choice leads to a gap of $\Theta(b/n)$. As our main result, we prove that the gap reduces to $O(\sqrt{(b/n)\cdot\log n})$, if one runs the $(1+\beta)$-process with an appropriately chosen $\beta$ (in fact this result holds for a larger class of processes). This not only proves the phenomenon that Two-Choice is not the best (leading to the formation of "towers" over previously light bins), but also that mixing two processes (One-Choice and Two-Choice) leads to a process which achieves a gap that is asymptotically smaller than both. We also derive a matching lower bound of $\Omega(\sqrt{(b/n)\cdot\log n})$ for any allocation process, which demonstrates that the above $(1+\beta)$-process is asymptotically optimal.


title: Discovering interpretable Lagrangian of dynamical systems from data

abstract: A complete understanding of physical systems requires models that are accurate and obeys natural conservation laws. Recent trends in representation learning involve learning Lagrangian from data rather than the direct discovery of governing equations of motion. The generalization of equation discovery techniques has huge potential; however, existing Lagrangian discovery frameworks are black-box in nature. This raises a concern about the reusability of the discovered Lagrangian. In this article, we propose a novel data-driven machine-learning algorithm to automate the discovery of interpretable Lagrangian from data. The Lagrangian are derived in interpretable forms, which also allows the automated discovery of conservation laws and governing equations of motion. The architecture of the proposed framework is designed in such a way that it allows learning the Lagrangian from a subset of the underlying domain and then generalizing for an infinite-dimensional system. The fidelity of the proposed framework is exemplified using examples described by systems of ordinary differential equations and partial differential equations where the Lagrangian and conserved quantities are known.


title: Minimal entropy production in the presence of anisotropic fluctuations

abstract: Anisotropy in temperature, chemical potential, or ion concentration, provides the fuel that feeds dynamical processes that sustain life. At the same time, anisotropy is a root cause of incurred losses manifested as entropy production. In this work we consider a rudimentary model of an overdamped stochastic thermodynamic system in an anisotropic temperature heat bath, and study minimum entropy production when driving the system between thermodynamic states in finite time. While entropy production in isotropic temperature environments can be expressed in terms of the length (in the Wasserstein-2 metric) traversed by the thermodynamic state of the system, anisotropy complicates substantially the mechanism of entropy production since, besides dissipation, seepage of energy between ambient anisotropic heat sources by way of the system dynamics is often a major contributing factor. A key result of the paper is to show that in the presence of anisotropy, minimization of entropy production can once again be expressed via a modified Optimal Mass Transport (OMT) problem. However, in contrast to the isotropic situation that leads to a classical OMT problem and a Wasserstein length, entropy production may not be identically zero when the thermodynamic state remains unchanged (unless one has control over non-conservative forces); this is due to the fact that maintaining a Non-Equilibrium Steady-State (NESS) incurs an intrinsic entropic cost that can be traced back to a seepage of heat between heat baths. As alluded to, NESSs represent hallmarks of life, since living matter by necessity operates far from equilibrium. Therefore, the question studied herein, to characterize minimal entropy production in anisotropic environments, appears of central importance in biological processes and on how such processes may have evolved to optimize for available usage of resources.


title: Cartesian closed varieties I: the classification theorem

abstract: In 1990, Johnstone gave a syntactic characterisation of the equational theories whose associated varieties are cartesian closed. Among such theories are all unary theories -- whose models are sets equipped with an action by a monoid M -- and all hyperaffine theories -- whose models are sets with an action by a Boolean algebra B. We improve on Johnstone's result by showing that an equational theory is cartesian closed just when its operations have a unique hyperaffine-unary decomposition. It follows that any non-degenerate cartesian closed variety is a variety of sets equipped with compatible actions by a monoid M and a Boolean algebra B; this is the classification theorem of the title.


title: Cartesian closed varieties II: links to operator algebra

abstract: This paper is the second in a series investigating cartesian closed varieties. In first of these, we showed that every non-degenerate finitary cartesian variety is a variety of sets equipped with an action by a Boolean algebra B and a monoid M which interact to form what we call a matched pair [B|M]. In this paper, we show that such pairs [B|M] are equivalent to Boolean restriction monoids and also to ample source-\'etale topological categories; these are generalisations of the Boolean inverse monoids and ample \'etale topological groupoids used by operator algebraists to encode structures such as Cuntz and Cuntz--Krieger C*-algebras, Leavitt path algebras and the C*-algebras associated to self-similar group actions. We explain and illustrate these links, and begin the programme of understanding how topological and algebraic properties of such groupoids can be understood from the logical perspective of the associated varieties.


title: Bargain hunting in a Coxeter group

abstract: Petersen and Tenner defined the depth statistic for Coxeter group elements which, in the symmetric group, can be described in terms of a cost function on transpositions. We generalize that cost function to the other classical (finite and affine) Weyl groups, letting the cost of an individual reflection $t$ be the distance between the integers transposed by $t$ in the combinatorial representation of the group (\`a la Eriksson and Eriksson). Arbitrary group elements then have a well-defined cost, obtained by minimizing the sum of the transposition costs among all factorizations of the element. We show that the cost of arbitrary elements can be computed directly from the elements themselves using a simple, intrinsic formula.


title: Extended Decoupled Anisotropic Solutions in $f(\mathcal{R},\mathcal{T},\mathcal{R}_{\gamma\chi}\mathcal{T}^{\gamma\chi})$ Gravity

abstract: In this paper, we consider static spherical structure to develop some anisotropic solutions by employing the extended gravitational decoupling scheme in the background of $f(\mathcal{R},\mathcal{T},\mathcal{R}_{\gamma\chi}\mathcal{T}^{\gamma\chi})$ gravity, where $\mathcal{R}$ and $\mathcal{T}$ indicate the Ricci scalar and trace of the energy-momentum tensor, respectively. We transform both radial as well as temporal metric functions and apply them on the field equations that produce two different sets corresponding to the decoupling parameter $\xi$. The first set is associated with isotropic distribution, i.e., modified Krori-Barua solution. The second set is influenced from anisotropic factor and contains unknowns which are determined by taking some constraints. The impact of decoupling parameter is then analyzed on the obtained physical variables and anisotropy. We also investigate energy conditions and some other parameters such as mass, compactness and redshift graphically. It is found that our solution corresponding to pressure-like constraint shows stable behavior throughout in this gravity for the considered range of $\xi$.


title: Neural Architecture Search via Two Constant Shared Weights Initialisations

abstract: In recent years, zero-cost metrics are gaining ground in neural architecture search (NAS). There metrics allow finding the optimal neural network for a given task faster and with a lesser computational load than conventional NAS methods. Equally important is that they also shed some light on the internal workings of neural architectures. This paper presents a zero-cost metric that highly correlated with the train set accuracy across the NAS-Bench-101, NAS-Bench-201 and NAS-Bench-NLP benchmark datasets. We evaluate a neural achitecture's potential based on the outputs' statistics after two constant shared weights initialisations. For this, we only use an unlabelled mini-batch of data. We observe that the dispersion of the outputs between two initialisations positively correlates with trained accuracy. The correlation further improves when we normalise dispersion by average output magnitude. The resulting metric, epsilon, does not require gradients computation and unbinds the NAS procedure from training hyperparameters, loss metrics and human-labelled data. Our method is easy to integrate within existing NAS algorithms and takes a fraction of a second to evaluate a single network. The code supporting this study can be found on GitHub at https://github.com/egracheva/epsinas.


title: Bayesian Non-parametric Hidden Markov Model for Agile Radar Pulse Sequences Streaming Analysis

abstract: Multi-function radars (MFRs) are sophisticated types of sensors with the capabilities of complex agile inter-pulse modulation implementation and dynamic work mode scheduling. The developments in MFRs pose great challenges to modern electronic reconnaissance systems or radar warning receivers for recognition and inference of MFR work modes. To address this issue, this paper proposes an online processing framework for parameter estimation and change point detection of MFR work modes. At first, this paper designed a fully-conjugate Bayesian non-parametric hidden Markov model with a designed prior distribution (agile BNP-HMM) to represent the MFR pulse agility characteristics. The proposed model allows fully-variational Bayesian inference. Then, the proposed framework is constructed by two main parts. The first part is the agile BNP-HMM model for automatically inferring the number of HMM hidden states and emission distribution of the corresponding hidden states. An estimation error lower bound on performance is derived and the proposed algorithm is shown to be close to the bound. The second part utilizes the streaming Bayesian updating to facilitate computation, and designed an online work mode change detection framework based upon a weighted sequential probability ratio test. We demonstrate that the proposed framework is consistently highly effective and robust to baseline methods on diverse simulated data-sets.


title: Updated Constraints on Primordial Black Hole Evaporation

abstract: The Hawking evaporation process, leading to the production of detectable particle species, constrains the abundance of light black holes, presumably of primordial origin. Here, we reconsider and correct constraints from soft gamma-ray observations, including of the gamma-ray line, at 511 keV, produced by electron-positron pair-annihilation, where positrons originate from black hole evaporation. First, we point out that the INTEGRAL detection of the Large Magellanic Cloud provides one of the strongest bounds attainable with present observations; and that future MeV gamma-ray telescopes, such as GECCO, will greatly enhance such constraints. Second, we discuss issues with previous limits from the isotropic flux at 511 keV and we provide updated, robust constraints from recent measurements of the diffuse Galactic soft gamma-ray emission and from the isotropic soft gamma-ray background.


title: Intrinsic auxeticity and mechanical anisotropy of Si9C15 siligraphene

abstract: The graphene-like two-dimensional (2D) silicon carbide or siligraphene has attracted remarkable attentions, owing to its fascinating physical properties. Nevertheless, the first high-quality siligraphene, i.e., monolayer Si9C15 was synthesised very recently, which exhibits an excellent semiconducting behaviour. In this work, we investigate the mechanical properties of Si9C15 siligraphene by using atomistic simulations including density functional theory (DFT) calculations and molecular dynamics (MD) simulations. Both methods confirm the existence of intrinsic negative Poisson's ratios in Si9C15 siligraphene, which, as illustrated by MD simulations, is attributed to the tension-induced de-wrinkling behaviours of its intrinsic rippled configuration. Different de-wrinkling behaviours are observed in different directions of Si9C15 siligraphene, which result in the anisotropy of its auxetic property. The fracture property of Si9C15 siligraphene is similarly anisotropic, but relatively large fracture strains are observed in different orientations, indicating the stretchability of Si9C15 siligraphene. The stretchability together with the strain-sensitive bandgap of Si9C15 siligraphene observed in DFT calculations indicates the effectiveness of strain engineering in modulating its electronic property. The combination of unique auxetic property, excellent mechanical property and tunable electronic property may render Si9C15 siligraphene a novel 2D material with multifunctional applications.


title: Simulation-to-reality UAV Fault Diagnosis with Deep Learning

abstract: Accurate diagnosis of propeller faults is crucial for ensuring the safe and efficient operation of quadrotors. Training a fault classifier using simulated data and deploying it on a real quadrotor is a cost-effective and safe approach. However, the simulation-to-reality gap often leads to poor performance of the classifier when applied in real flight. In this work, we propose a deep learning model that addresses this issue by utilizing newly identified features (NIF) as input and utilizing domain adaptation techniques to reduce the simulation-to-reality gap. In addition, we introduce an adjusted simulation model that generates training data that more accurately reflects the behavior of real quadrotors. The experimental results demonstrate that our proposed approach achieves an accuracy of 96\% in detecting propeller faults. To the best of our knowledge, this is the first reliable and efficient method for simulation-to-reality fault diagnosis of quadrotor propellers.


title: Geometry of Score Based Generative Models

abstract: In this work, we look at Score-based generative models (also called diffusion generative models) from a geometric perspective. From a new view point, we prove that both the forward and backward process of adding noise and generating from noise are Wasserstein gradient flow in the space of probability measures. We are the first to prove this connection. Our understanding of Score-based (and Diffusion) generative models have matured and become more complete by drawing ideas from different fields like Bayesian inference, control theory, stochastic differential equation and Schrodinger bridge. However, many open questions and challenges remain. One problem, for example, is how to decrease the sampling time? We demonstrate that looking from geometric perspective enables us to answer many of these questions and provide new interpretations to some known results. Furthermore, geometric perspective enables us to devise an intuitive geometric solution to the problem of faster sampling. By augmenting traditional score-based generative models with a projection step, we show that we can generate high quality images with significantly fewer sampling-steps.


title: Spatiotemporal factor models for functional data with application to population map forecast

abstract: The proliferation of mobile devices has led to the collection of large amounts of population data. This situation has prompted the need to utilize this rich, multidimensional data in practical applications. In response to this trend, we have integrated functional data analysis (FDA) and factor analysis to address the challenge of predicting hourly population changes across various districts in Tokyo. Specifically, by assuming a Gaussian process, we avoided the large covariance matrix parameters of the multivariate normal distribution. In addition, the data were both time and spatially dependent between districts. To capture these characteristics, a Bayesian factor model was introduced, which modeled the time series of a small number of common factors and expressed the spatial structure through factor loading matrices. Furthermore, the factor loading matrices were made identifiable and sparse to ensure the interpretability of the model. We also proposed a Bayesian shrinkage method as a systematic approach for factor selection. Through numerical experiments and data analysis, we investigated the predictive accuracy and interpretability of our proposed method. We concluded that the flexibility of the method allows for the incorporation of additional time series features, thereby improving its accuracy.


title: Mathematical Modeling of Cyber Resilience

abstract: We identify quantitative characteristics of responses to cyber compromises that can be learned from repeatable, systematic experiments. We model a vehicle equipped with an autonomous cyber-defense system and which also has some inherent physical resilience features. When attacked by malware, this ensemble of cyber-physical features (i.e., "bonware") strives to resist and recover from the performance degradation caused by the malware's attack. We propose parsimonious continuous models, and develop stochastic models to aid in quantifying systems' resilience to cyber attacks.


title: Cosmology with the Galaxy Bispectrum Multipoles: Optimal Estimation and Application to BOSS Data

abstract: We present a framework for self-consistent cosmological analyses of the full-shape anisotropic bispectrum, including the quadrupole $(\ell=2)$ and hexadecapole $(\ell=4)$ moments. This features a novel window-free algorithm for extracting the latter quantities from data, derived using a maximum-likelihood prescription. Furthermore, we introduce a theoretical model for the bispectrum multipoles (which does not introduce new free parameters), and test both aspects of the pipeline on several high-fidelity mocks, including the PT Challenge suite of gigantic cumulative volume. This establishes that the systematic error is significantly below the statistical threshold, both for the measurement and modeling. As a realistic example, we extract the large-scale bispectrum multipoles from BOSS DR12 and analyze them in combination with the power spectrum data. Assuming a minimal $\Lambda$CDM model, with a BBN prior on the baryon density and a \textit{Planck} prior on $n_s$, we can extract the remaining cosmological parameters directly from the clustering data. The inclusion of the unwindowed higher-order $(\ell>0)$ large-scale bispectrum multipoles is found to moderately improve one-dimensional cosmological parameter posteriors (at the $5\%-10\%$ level), though these multipoles are detected only in three out of four BOSS data segments at $\approx 5\sigma$. Combining information from the power spectrum and bispectrum multipoles, the real space power spectrum, and the post-reconstructed BAO data, we find $H_0 = 68.2\pm 0.8~\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, $\Omega_m =0.33\pm 0.01$ and $\sigma_8 = 0.736\pm 0.033$ (the tightest yet found in perturbative full-shape analyses). Our estimate of the growth parameter $S_8=0.77\pm 0.04$ agrees with both weak lensing and CMB results.