title: Compositional Separation of Control Flow and Data Flow

abstract: Every Model of High-Level Computation (MHC) has an underlying composition mechanism for combining simple computation devices into more complex ones. Composition can be done by (explicitly or implicitly) defining control flow, data flow or any combination thereof. Control flow specifies the order in which individual computations are activated, whereas data flow defines how data is exchanged among them. Unfortunately, traditional MHCs either mix data and control or only consider one dimension explicitly, which makes it difficult to reason about data flow and control flow separately. Reasoning about these dimensions orthogonally is a crucial desideratum for optimisation, maintainability and verification purposes. In this paper, we introduce a novel MHC that explicitly treats data flow and control flow as separate dimensions, while providing modularity. As the model is rooted in category theory, it provides category-theoretic operations for compositionally constructing sequential, parallel, branchial or iterative composites. Compositionality entails that a composite exhibits the same properties as its respective constituents, including separation of concerns and modularity. We conclude the paper by demonstrating how our proposal can be used to model high-level computations in two different application domains: software engineering and artificial intelligence.


title: Bifurcation and periodic solutions to neuroscience models with a small parameter

abstract: The existence of periodic solutions is proven for some neuroscience models with a small parameter. Moreover, the stability of such solutions is investigated, as well. The results are based on a theoretical research dealing with the functional differential equation with parameters $$ \dot{x}(t)=L(\tau) x_t + \varepsilon f(t, x_t), $$ where $L: \mathbb{R}_+\rightarrow \mathcal{L}(C; \mathbb{R})$ and $f: \mathbb{R} \times C \rightarrow \mathbb{R}$ are, respectively, linear and nonlinear operators, and $\varepsilon>0$ is a small enough parameter. The theoretical results are applied to a Parkinson's disease model, where the obtained conclusions are illustrated by numerical simulations.


title: Orbital-selective effect of spin reorientation on the Dirac fermions in a three-dimensional kagome ferromagnet Fe$_3$Ge

abstract: Kagome magnets provide a fascinating platform for the realization of correlated topological quantum phases under various magnetic ground states. However, the intricate effect of the magnetic spin configurations on the characteristic electronic structure directly from the kagome lattice layer remains still elusive. Here, utilizing angle-resolved photoemission spectroscopy and density functional theory calculations, we report the spectroscopic evidence for the spin-reorientation effect of a kagome ferromagnet Fe$_3$Ge, which is composed only of the kagome planes. There are two kinds of kagome-derived Dirac fermions due to the structural three-dimensionality -- one is less dispersive ($k_z$ $\sim$ 0) and the other disperses linearly ($k_z$ $\sim$ $\pi$). As the Fe moments cant from the $c$ axis into the $ab$ plane upon cooling, the Dirac fermion in $k_z$ $\sim$ 0 plane with a mixture of the Fe-$3d_{xy}$ and Fe-$3d_{x^2-y^2}$ components evolves from gapped into nearly gapless, while the Dirac cone in $k_z$ $\sim$ $\pi$ plane mainly of the Fe-$3d_{x^2-y^2}$ orbital character remains intact, suggesting that the effect of spin reorientation on the Dirac fermions has an orbital selectivity. Our unambiguous observations provide a feasible route to design and manipulate the mass of Dirac fermions for realizing the novel quantum phases. We also perform comparative studies between the non-charge-ordered Fe$_3$Ge and its sibling compound FeGe, a newly established charge-density-wave kagome magnet, the results suggest that the orbital-selective van Hove singularities near the Fermi level play an indispensable part in driving the charge order on a magnetic kagome lattice.


title: Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

abstract: Understanding the impact of extractant functionalization in solvent extraction is essential to guide the development of better separations processes. Traditionally, computational extractant design uses electronic structure calculations to determine the metal binding energy of the lowest energy state. Although highly accurate, this approach does not account for all the relevant physics encountered under experimental conditions, such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical MD simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. We generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the head group and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the head group has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal binding energetics based on experimentally reported distribution ratios. This study demonstrates the feasibility of using molecular dynamics simulations with advance sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.


title: Diagonal operators, $q$-Whittaker functions and rook theory

abstract: We discuss the problem posed by Bender, Coley, Robbins and Rumsey of enumerating the number of subspaces which have a given profile with respect to a linear operator over the finite field $\mathbb{F}_q$. We solve this problem in the case where the operator is diagonalizable. The solution leads us to a new class of polynomials $b_{\mu\nu}(q)$ indexed by pairs of integer partitions. These polynomials have several interesting specializations and can be expressed as positive sums over semistandard tableaux. We present a new correspondence between set partitions and semistandard tableaux. A close analysis of this correspondence reveals the existence of several new set partition statistics which generate the polynomials $b_{\mu\nu}(q)$; each such statistic arises from a Mahonian statistic on multiset permutations. The polynomials $b_{\mu\nu}(q)$ are also given a description in terms of coefficients in the monomial expansion of $q$-Whittaker symmetric functions which are specializations of Macdonald polynomials. We express the Touchard--Riordan generating polynomial for chord diagrams by number of crossings in terms of $q$-Whittaker functions. We also introduce a class of $q$-Stirling numbers defined in terms of the polynomials $b_{\mu\nu}(q)$ and present connections with $q$-rook theory in the spirit of Garsia and Remmel.


title: Expressive dynamics models with nonlinear injective readouts enable reliable recovery of latent features from neural activity

abstract: The advent of large-scale neural recordings has enabled new methods to discover the computational mechanisms of neural circuits by understanding the rules that govern how their state evolves over time. While these \textit{neural dynamics} cannot be directly measured, they can typically be approximated by low-dimensional models in a latent space. How these models represent the mapping from latent space to neural space can affect the interpretability of the latent representation. We show that typical choices for this mapping (e.g., linear or MLP) often lack the property of injectivity, meaning that changes in latent state are not obligated to affect activity in the neural space. During training, non-injective readouts incentivize the invention of dynamics that misrepresent the underlying system and the computation it performs. Combining our injective Flow readout with prior work on interpretable latent dynamics models, we created the Ordinary Differential equations autoencoder with Injective Nonlinear readout (ODIN), which captures latent dynamical systems that are nonlinearly embedded into observed neural activity via an approximately injective nonlinear mapping. We show that ODIN can recover nonlinearly embedded systems from simulated neural activity, even when the nature of the system and embedding are unknown. Additionally, ODIN enables the unsupervised recovery of underlying dynamical features (e.g., fixed points) and embedding geometry. When applied to biological neural recordings, ODIN can reconstruct neural activity with comparable accuracy to previous state-of-the-art methods while using substantially fewer latent dimensions. Overall, ODIN's accuracy in recovering ground-truth latent features and ability to accurately reconstruct neural activity with low dimensionality make it a promising method for distilling interpretable dynamics that can help explain neural computation.


title: On the status of DELL systems

abstract: A detailed review of the $p,q$-duality for Calogero system and its generalizations is given. For the first time, we present some of elliptic-trigonometric Hamiltonians dual to the elliptic Ruijsenaars Hamiltonians (i.e. trigonometric-elliptic ones), and explain their relations to the bi-elliptic Koroteev-Shakirov (KS) model. The most interesting self-dual double-elliptic (DELL) system remains a mystery, but we provide a clearer formulation of the problem and describe the steps that are still to be done.


title: Tunable Circular Photogalvanic and Photovoltaic Effect in 2D Tellurium with Different Chirality

abstract: Chirality arises from the asymmetry of matters, where two counterparts are the mirror image of each other. The interaction between circular-polarization light and quantum materials is enhanced in chiral space groups due to the structural chirality. Tellurium (Te) possesses the simplest chiral crystal structure, with Te atoms covalently bonded into a spiral atomic chain (left- or right-handed) with a periodicity of three. Here, we investigate the tunable circular photo-electric responses in 2D Te field-effect transistor with different chirality, including the longitudinal circular photogalvanic effect induced by the radial spin texture (electron-spin polarization parallel to the electron momentum direction) and the circular photovoltaic induced by the chiral crystal structure (helical Te atomic chains). Our work demonstrates the controllable manipulation of the chirality degree of freedom in materials.


title: 4$f$ electron temperature driven ultrafast electron localization

abstract: Valence transitions in strongly correlated electron systems are caused by orbital hybridization and Coulomb interactions between localized and delocalized electrons. The transition can be triggered by changes in the electronic structure and is sensitive to temperature variations, applications of magnetic fields, and physical or chemical pressure. Launching the transition by photoelectric fields can directly excite the electronic states and thus provides an ideal platform to study the correlation among electrons on ultrafast timescales. The EuNi$_2$(Si$_{0.21}$Ge$_{0.79}$)$_2$ mixed-valence metal is an ideal material to investigate the valence transition of the Eu ions via the amplified orbital hybridization by the photoelectric field on sub-picosecond timescales. A direct view on the 4$f$ electron occupancy of the Eu ions is required to understand the microscopic origin of the transition. Here we probe the 4$f$ electron states of EuNi$_2$(Si$_{0.21}$Ge$_{0.79}$)$_2$ at the sub-ps timescale after photoexcitation by X-ray absorption spectroscopy across the Eu $M_5$-absorption edge. The observed spectral changes due to the excitation indicate a population change of total angular momentum multiplet states $J$ = 0, 1, 2, and 3 of Eu$^{3+}$, and the Eu$^{2+}$ $J$ = 7/2 multiplet state caused by an increase in 4$f$ electron temperature that results in a 4$f$ localization process. This electronic temperature increase combined with fluence-dependent screening accounts for the strongly non-linear effective valence change. The data allow us to extract a time-dependent determination of an effective temperature of the 4$f$ shell, which is also of great relevance in the understanding of metallic systems' properties, such as the ultrafast demagnetization of ferromagnetic rare-earth intermetallics and their all-optical magnetization switching.


title: In operando cryo-STEM of pulse-induced charge density wave switching in TaS$_2$

abstract: The charge density wave (CDW) material 1T-TaS$_2$ exhibits a pulse-induced insulator-to-metal transition, which shows promise for next-generation electronics such as memristive memory and neuromorphic hardware. However, the rational design of TaS$_2$ devices is hindered by a poor understanding of the switching mechanism, the pulse-induced phase, and the influence of material defects. Here, we operate a 2-terminal TaS$_2$ device within a scanning transmission electron microscope (STEM) at cryogenic temperature, and directly visualize the changing CDW structure with nanoscale spatial resolution and down to 300 {\mu}s temporal resolution. We show that the pulse-induced transition is driven by Joule heating, and that the pulse-induced state corresponds to nearly commensurate and incommensurate CDW phases, depending on the applied voltage amplitude. With our in operando cryo-STEM experiments, we directly correlate the CDW structure with the device resistance, and show that dislocations significantly impact device performance. This work resolves fundamental questions of resistive switching in TaS$_2$ devices critical for engineering reliable and scalable TaS$_2$ electronics.


title: Self-Similar Outflows at the Source of the Fast Solar Wind: A Smoking Gun of Multiscale Impulsive Reconnection?

abstract: We present results of a quantitative analysis of structured plasma outflows above a polar coronal hole observed by the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory spacecraft. In a 6-hour interval of continuous high-cadence SDO/AIA images, we identified more than 2300 episodes of small-scale plasma flows in the polar corona. The mean upward flow speed measured by the surfing transform technique (Uritsky et al., 2013) is estimated to be 122 $\pm$ 34 \kms, which is comparable to the local sound speed. The typical recurrence period of the flow episodes is 10 to 30 minutes, and the mean duration and transverse size of each episode are about 3-5 min and 3-4 Mm, respectively. The largest identifiable episodes last for tens of minutes and reach widths up to $40$ Mm. For the first time, we demonstrate that the polar coronal-hole outflows obey a family of power-law probability distributions characteristic of impulsive interchange magnetic reconnection. Turbulent photospheric driving may play a crucial role in releasing magnetically confined plasma onto open field. The estimated occurrence rate of the detected self-similar coronal outflows is sufficient for them to make a dominant contribution to the fast-wind mass and energy fluxes and to account for the wind's small-scale structure.


title: Surface Casimir densities on branes orthogonal to the boundary of anti-de Sitter spacetime

abstract: We investigate the vacuum expectation value of the surface energy-momentum tensor (SEMT) for a scalar field with general curvature coupling in the geometry of two branes orthogonal to the boundary of anti-de Sitter (AdS) spacetime. For Robin boundary conditions on the branes, the SEMT is decomposed into the contributions corresponding to the self-energies of the branes and the parts induced by the presence of the second brane. The renormalization is required for the first parts only and for the corresponding regularization the generalized zeta function method is employed. The induced SEMT is finite and is free from renormalization umbiguities. For an observer living on the brane, the corresponding equation of state is of the cosmological constant type. Depending on the boundary conditions and on the separation between the branes, the surface energy densities can be either positive or negative. The energy density induced on the brane vanishes in special cases of Dirichlet and Neumann boundary conditions on that brane. The effect of gravity on the induced SEMT is essential at separations between the branes of the order or larger than the curvature radius for AdS spacetime. In the large separation limit the decay of the SEMT, as a function of the proper separation, follows a power law for both massless and massive fields. For parallel plates in Minkowski bulk and for massive fields the fall-off of the corresponding expectation value is exponential.


title: Design and Implementation of DC-to-5~MHz Wide-Bandwidth High-Power High-Fidelity Converter

abstract: Advances in power electronics have made it possible to achieve high power levels, e.g., reaching GW in grids, or alternatively high output bandwidths, e.g., beyond MHz in communication. Achieving both simultaneously, however, remains challenging. Various applications, ranging from efficient multichannel wireless power transfer to cutting-edge medical and neuroscience applications, are demanding both high power and wide bandwidth. Conventional inverters can achieve high power and high quality at grid or specific frequency ranges but lose their fidelity when reaching higher output frequencies. Resonant circuits can promise a high output frequency but only a narrow bandwidth. We overcome the hardware challenges by combining gallium-nitride (GaN) transistors with modular cascaded double-H bridge circuits and control that can manage typical timing and balancing issues. We developed a lightweight embedded control solution that includes an improved look-up-table digital synthesizer and a novel adaptive-bias-elimination nearest-level modulation. This solution effectively solves the conflict between a high power level and high output bandwidth and can--in contrast to previous approaches--in principle be scaled in both dimensions. Our prototype exhibits a frequency range from DC to 5 MHz with <18% total voltage distortion across the entire frequency spectrum, while achieving a power level of >5 kW. We conducted tests by sweeping the output frequency and two channel-mixing trials, which included a practical magnetogenetics-oriented stimulation pulse and an entertaining trial to reproduce the famous Arecibo message with the current spectrum.


title: Solving the Pulsar Equation using Physics-Informed Neural Networks

abstract: In this study, Physics-Informed Neural Networks (PINNs) are skilfully applied to explore a diverse range of pulsar magneto-spheric models, specifically focusing on axisymmetric cases. The study successfully reproduced various axisymmetric models found in the literature, including those with non-dipolar configurations, while effectively characterizing current sheet features. Energy losses in all studied models were found to exhibit reasonable similarity, differing by no more than a factor of three from the classical dipole case. This research lays the groundwork for a reliable elliptic Partial Differential Equation solver tailored for astrophysical problems. Based on these findings, we foresee that the utilization of PINNs will become the most efficient approach in modelling three-dimensional magnetospheres. This methodology shows significant potential and facilitates an effortless generalization, contributing to the advancement of our understanding of pulsar magnetospheres.


title: Bose-Einstein condensation and muon production in ultra-high energy cosmic ray particle collisions

abstract: Collisions of cosmic ray particles with ultra-high initial energies with nuclei in the atmosphere open a wide room for appearing of the novel dynamical features for multiparticle production processes. In particular, the pion-lasing behavior driven by Bose-Einstein condensation would result in the shift to larger multiplicities and, as consequence, could provide, in general, the enhanced yield of cosmic muons. In the present work the critical value of the space charged particle density for onset of Bose-Einstein condensation of the boson (pion) wave-packets into the same wave-packet state is estimated within the model with complete multiparticle symmetrization for the energy domain corresponded to the ultra-high energy cosmic rays (UHECR). Energy dependence of mean density of charged pions is evaluated for the cases of absent of the Bose-Einstein effects and for presence of laser-like behavior of pions. The possible influence of the Bose-Einstein condensation is discussed for the muon production in UHECR particle collisions with the atmosphere.


title: Quantum transport in a multi-path graphene Aharonov-Bohm inteferometer

abstract: We investigate the quantum transport dynamics of electrons in a multi-path Aharonov-Bohm interferometer comprising several parallel graphene nanoribbons. At low magnetic field strengths, the conductance displays a complex oscillatory behavior stemming from the interference of electron wave functions from different paths, reminiscent of the diffraction grating in optics. With increasing magnetic field strength, certain nanoribbons experience transport blockade, leading to conventional Aharonov-Bohm oscillations arising from two-path interference. We also discuss the impact of edge effects and the influence of finite temperature. Our findings offer valuable insights for experimental investigations of quantum transport in multi-path devices and their potential application for interferometry and quantum sensing.


title: On Computationally Efficient Learning of Exponential Family Distributions

abstract: We consider the classical problem of learning, with arbitrary accuracy, the natural parameters of a $k$-parameter truncated \textit{minimal} exponential family from i.i.d. samples in a computationally and statistically efficient manner. We focus on the setting where the support as well as the natural parameters are appropriately bounded. While the traditional maximum likelihood estimator for this class of exponential family is consistent, asymptotically normal, and asymptotically efficient, evaluating it is computationally hard. In this work, we propose a novel loss function and a computationally efficient estimator that is consistent as well as asymptotically normal under mild conditions. We show that, at the population level, our method can be viewed as the maximum likelihood estimation of a re-parameterized distribution belonging to the same class of exponential family. Further, we show that our estimator can be interpreted as a solution to minimizing a particular Bregman score as well as an instance of minimizing the \textit{surrogate} likelihood. We also provide finite sample guarantees to achieve an error (in $\ell_2$-norm) of $\alpha$ in the parameter estimation with sample complexity $O({\sf poly}(k)/\alpha^2)$. Our method achives the order-optimal sample complexity of $O({\sf log}(k)/\alpha^2)$ when tailored for node-wise-sparse Markov random fields. Finally, we demonstrate the performance of our estimator via numerical experiments.


title: Just-in-Time autotuning

abstract: Performance portability is a major concern on current architectures. One way to achieve it is by using autotuning. In this paper, we are presenting how we exten ded a just-in-time compilation infrastructure to introduce autotuning capabiliti es triggered at run-time. When a function is executed, the first iterations optimize it, and once the best solution has been found, it is used for subsequent calls to the function. This just-in-time autotuning infrastructure is relevant for optimizing computation kernels that will be called numerous times with similar parameters through the execution, re-optimizes kernels when they are called with other parameters, and the programmer can obtain the optimal parameters to use them for other kernels. We present an experimental performance evaluation of our approach. Compiling the code introduces an overhead on the first iterations, and this overhead is compensated for during subsequent iterations. We also determined that the optimum found seems stable and accurate.


title: Down the Toxicity Rabbit Hole: A Novel Framework to Bias Audit Large Language Models

abstract: This paper makes three contributions. First, it presents a generalizable, novel framework dubbed \textit{toxicity rabbit hole} that iteratively elicits toxic content from a wide suite of large language models. Spanning a set of 1,266 identity groups, we first conduct a bias audit of \texttt{PaLM 2} guardrails presenting key insights. Next, we report generalizability across several other models. Through the elicited toxic content, we present a broad analysis with a key emphasis on racism, antisemitism, misogyny, Islamophobia, homophobia, and transphobia. Finally, driven by concrete examples, we discuss potential ramifications.


title: Kubo formula for dc conductivity: generalization to systems with spin-orbit coupling

abstract: We revise the Kubo formula for the electric dc conductivity in the presence of spin-orbit coupling (SOC). We discover that each velocity operator that enters this formula differs from $\partial H/\partial \boldsymbol p$, where $H$ is the Hamiltonian and $\boldsymbol p$ is the canonical momentum. Moreover, we find an additional contribution to the Hall dc conductivity from noncommuting coordinates that is missing in the conventional Kubo-Streda formula. This contribution originates from the "electron-positron" matrix elements of the velocity and position operators. We argue that the widely used Rashba model does in fact provide a finite anomalous Hall dc conductivity in the metallic regime (in the noncrossing approximation) if SOC-corrections to the velocity and position operators are properly taken into account. While we focus on the response of the charge current to the electric field, linear response theories of other SOC-related effects should be modified similarly.


title: The trigger system for the CSR external-target experiment

abstract: A trigger system has been designed and implemented for the HIRFL-CSR external target experiment (CEE), the spectrometer for studying nuclear matter properties with heavy ion collisions in the GeV energy region. The system adopts master-slave structure and serial data transmission mode using optical fiber to deal with different types of detectors and long-distance signal transmission. The trigger logic can be accessed based on command register and controlled by a remote computer. The overall field programmable gate array (FPGA) logic can be flexibly reconfigured online to match the physical requirements of the experiment. The trigger system has been tested in beam experiment. It is demonstrated that the trigger system functions correctly and meets the physical requirements of CEE.


title: C4CAM: A Compiler for CAM-based In-memory Accelerators

abstract: Machine learning and data analytics applications increasingly suffer from the high latency and energy consumption of conventional von Neumann architectures. Recently, several in-memory and near-memory systems have been proposed to remove this von Neumann bottleneck. Platforms based on content-addressable memories (CAMs) are particularly interesting due to their efficient support for the search-based operations that form the foundation for many applications, including K-nearest neighbors (KNN), high-dimensional computing (HDC), recommender systems, and one-shot learning among others. Today, these platforms are designed by hand and can only be programmed with low-level code, accessible only to hardware experts. In this paper, we introduce C4CAM, the first compiler framework to quickly explore CAM configurations and to seamlessly generate code from high-level TorchScript code. C4CAM employs a hierarchy of abstractions that progressively lowers programs, allowing code transformations at the most suitable abstraction level. Depending on the type and technology, CAM arrays exhibit varying latencies and power profiles. Our framework allows analyzing the impact of such differences in terms of system-level performance and energy consumption, and thus supports designers in selecting appropriate designs for a given application.


title: Radiology-Llama2: Best-in-Class Large Language Model for Radiology

abstract: This paper introduces Radiology-Llama2, a large language model specialized for radiology through a process known as instruction tuning. Radiology-Llama2 is based on the Llama2 architecture and further trained on a large dataset of radiology reports to generate coherent and clinically useful impressions from radiological findings. Quantitative evaluations using ROUGE metrics on the MIMIC-CXR and OpenI datasets demonstrate that Radiology-Llama2 achieves state-of-the-art performance compared to other generative language models, with a Rouge-1 score of 0.4834 on MIMIC-CXR and 0.4185 on OpenI. Additional assessments by radiology experts highlight the model's strengths in understandability, coherence, relevance, conciseness, and clinical utility. The work illustrates the potential of localized language models designed and tuned for specialized domains like radiology. When properly evaluated and deployed, such models can transform fields like radiology by automating rote tasks and enhancing human expertise.


title: Verifiable Reinforcement Learning Systems via Compositionality

abstract: We propose a framework for verifiable and compositional reinforcement learning (RL) in which a collection of RL subsystems, each of which learns to accomplish a separate subtask, are composed to achieve an overall task. The framework consists of a high-level model, represented as a parametric Markov decision process, which is used to plan and analyze compositions of subsystems, and of the collection of low-level subsystems themselves. The subsystems are implemented as deep RL agents operating under partial observability. By defining interfaces between the subsystems, the framework enables automatic decompositions of task specifications, e.g., reach a target set of states with a probability of at least 0.95, into individual subtask specifications, i.e. achieve the subsystem's exit conditions with at least some minimum probability, given that its entry conditions are met. This in turn allows for the independent training and testing of the subsystems. We present theoretical results guaranteeing that if each subsystem learns a policy satisfying its subtask specification, then their composition is guaranteed to satisfy the overall task specification. Conversely, if the subtask specifications cannot all be satisfied by the learned policies, we present a method, formulated as the problem of finding an optimal set of parameters in the high-level model, to automatically update the subtask specifications to account for the observed shortcomings. The result is an iterative procedure for defining subtask specifications, and for training the subsystems to meet them. Experimental results demonstrate the presented framework's novel capabilities in environments with both full and partial observability, discrete and continuous state and action spaces, as well as deterministic and stochastic dynamics.


title: AGMDT: Virtual Staining of Renal Histology Images with Adjacency-Guided Multi-Domain Transfer

abstract: Renal pathology, as the gold standard of kidney disease diagnosis, requires doctors to analyze a series of tissue slices stained by H&E staining and special staining like Masson, PASM, and PAS, respectively. These special staining methods are costly, time-consuming, and hard to standardize for wide use especially in primary hospitals. Advances of supervised learning methods have enabled the virtually conversion of H&E images into special staining images, but achieving pixel-to-pixel alignment for training remains challenging. In contrast, unsupervised learning methods regarding different stains as different style transfer domains can utilize unpaired data, but they ignore the spatial inter-domain correlations and thus decrease the trustworthiness of structural details for diagnosis. In this paper, we propose a novel virtual staining framework AGMDT to translate images into other domains by avoiding pixel-level alignment and meanwhile utilizing the correlations among adjacent tissue slices. We first build a high-quality multi-domain renal histological dataset where each specimen case comprises a series of slices stained in various ways. Based on it, the proposed framework AGMDT discovers patch-level aligned pairs across the serial slices of multi-domains through glomerulus detection and bipartite graph matching, and utilizes such correlations to supervise the end-to-end model for multi-domain staining transformation. Experimental results show that the proposed AGMDT achieves a good balance between the precise pixel-level alignment and unpaired domain transfer by exploiting correlations across multi-domain serial pathological slices, and outperforms the state-of-the-art methods in both quantitative measure and morphological details.


title: Wurtzite vs rock-salt MnSe epitaxy: electronic and altermagnetic properties

abstract: Newly discovered altermagnets are magnetic materials exhibiting both compensated magnetic order, similar to antiferromagnets, and simultaneous non-relativistic spin-splitting of the bands, akin to ferromagnets. This characteristic arises from the specific symmetry operations that connect the spin sublattices. In this report, we show with ab initio calculations that the semiconductive MnSe exhibits altermagnetic spin-splitting in the wurtzite phase as well as a critical temperature well above room temperature. It is the first material from such space group identified to possess altermagnetic properties. Furthermore, we demonstrate experimentally through structural characterization techniques that it is possible to obtain thin films of both the intriguing wurtzite phase of MnSe and the more common rock-salt MnSe using molecular beam epitaxy on GaAs substrates. The choice of buffer layers plays a crucial role in determining the resulting phase and consequently extends the array of materials available for the physics of altermagnetism.


title: Accelerating Defect Predictions in Semiconductors Using Graph Neural Networks

abstract: Here, we develop a framework for the prediction and screening of native defects and functional impurities in a chemical space of Group IV, III-V, and II-VI zinc blende (ZB) semiconductors, powered by crystal Graph-based Neural Networks (GNNs) trained on high-throughput density functional theory (DFT) data. Using an innovative approach of sampling partially optimized defect configurations from DFT calculations, we generate one of the largest computational defect datasets to date, containing many types of vacancies, self-interstitials, anti-site substitutions, impurity interstitials and substitutions, as well as some defect complexes. We applied three types of established GNN techniques, namely Crystal Graph Convolutional Neural Network (CGCNN), Materials Graph Network (MEGNET), and Atomistic Line Graph Neural Network (ALIGNN), to rigorously train models for predicting defect formation energy (DFE) in multiple charge states and chemical potential conditions. We find that ALIGNN yields the best DFE predictions with root mean square errors around 0.3 eV, which represents a prediction accuracy of 98 % given the range of values within the dataset, improving significantly on the state-of-the-art. Models are tested for different defect types as well as for defect charge transition levels. We further show that GNN-based defective structure optimization can take us close to DFT-optimized geometries at a fraction of the cost of full DFT. DFT-GNN models enable prediction and screening across thousands of hypothetical defects based on both unoptimized and partially-optimized defective structures, helping identify electronically active defects in technologically-important semiconductors.


title: Unveiling the potential of large language models in generating semantic and cross-language clones

abstract: Semantic and Cross-language code clone generation may be useful for code reuse, code comprehension, refactoring and benchmarking. OpenAI's GPT model has potential in such clone generation as GPT is used for text generation. When developers copy/paste codes from Stack Overflow (SO) or within a system, there might be inconsistent changes leading to unexpected behaviours. Similarly, if someone possesses a code snippet in a particular programming language but seeks equivalent functionality in a different language, a semantic cross-language code clone generation approach could provide valuable assistance. In this study, using SemanticCloneBench as a vehicle, we evaluated how well the GPT-3 model could help generate semantic and cross-language clone variants for a given fragment.We have comprised a diverse set of code fragments and assessed GPT-3s performance in generating code variants.Through extensive experimentation and analysis, where 9 judges spent 158 hours to validate, we investigate the model's ability to produce accurate and semantically correct variants. Our findings shed light on GPT-3's strengths in code generation, offering insights into the potential applications and challenges of using advanced language models in software development. Our quantitative analysis yields compelling results. In the realm of semantic clones, GPT-3 attains an impressive accuracy of 62.14% and 0.55 BLEU score, achieved through few-shot prompt engineering. Furthermore, the model shines in transcending linguistic confines, boasting an exceptional 91.25% accuracy in generating cross-language clones


title: Random matrix statistics and safety rest areas on interstates in the United States

abstract: We analyze physical spacings between locations of safety rest areas on interstates in the United States. We show normalized safety rest area spacings on major interstates exhibit Wigner surmise statistics, which align with the eigenvalue spacings for the Gaussian Unitary Ensemble from random matrix theory as well as the one-dimensional gas interactions via the Coulomb potential. We identify economic and geographic regional traits at the state level that exhibit Poissonian statistics, which become more pronounced with increased geographical obstacles in interstate travel. Other regional filters (e.g., historical or political) produced results that did not diverge substantially from the overall Wigner surmise model.

categories: physics.soc-ph math-ph math.MP math.PR stat.AP


title: Suppression of lift-up effect in the 3D Boussinesq equations around a stably stratified Couette flow

abstract: In this paper, we establish linear enhanced dissipation results for the three-dimensional Boussinesq equations around a stably stratified Couette flow, in the viscous and thermally diffusive setting. The dissipation rates are faster compared to those observed in the homogeneous Navier-Stokes equations, in light of the interplay between velocity and temperature, driven by buoyant forces. Our approach involves introducing a change of variables grounded in a Fourier space symmetrization framework. This change elucidates the energy structure inherent in the system. Specifically, we handle non-streaks modes through an augmented energy functional, while streaks modes are amenable to explicit solutions. This explicit treatment reveals the oscillatory nature of shear modes, providing the elimination of the well-known three-dimensional instability mechanism known as the ``lift-up effect''.


title: Symmetric Stair Preconditioning of Linear Systems for Parallel Trajectory Optimization

abstract: There has been a growing interest in parallel strategies for solving trajectory optimization problems. One key step in many algorithmic approaches to trajectory optimization is the solution of moderately-large and sparse linear systems. Iterative methods are particularly well-suited for parallel solves of such systems. However, fast and stable convergence of iterative methods is reliant on the application of a high-quality preconditioner that reduces the spread and increase the clustering of the eigenvalues of the target matrix. To improve the performance of these approaches, we present a new parallel-friendly symmetric stair preconditioner. We prove that our preconditioner has advantageous theoretical properties when used in conjunction with iterative methods for trajectory optimization such as a more clustered eigenvalue spectrum. Numerical experiments with typical trajectory optimization problems reveal that as compared to the best alternative parallel preconditioner from the literature, our symmetric stair preconditioner provides up to a 34% reduction in condition number and up to a 25% reduction in the number of resulting linear system solver iterations.


title: Tail Gini Functional under Asymptotic Independence

abstract: Tail Gini functional is a measure of tail risk variability for systemic risks, and has many applications in banking, finance and insurance. Meanwhile, there is growing attention on aymptotic independent pairs in quantitative risk management. This paper addresses the estimation of the tail Gini functional under asymptotic independence. We first estimate the tail Gini functional at an intermediate level and then extrapolate it to the extreme tails. The asymptotic normalities of both the intermediate and extreme estimators are established. The simulation study shows that our estimator performs comparatively well in view of both bias and variance. The application to measure the tail variability of weekly loss of individual stocks given the occurence of extreme events in the market index in Hong Kong Stock Exchange provides meaningful results, and leads to new insights in risk management.


title: Efficient Inference on High-Dimensional Linear Models with Missing Outcomes

abstract: This paper is concerned with inference on the regression function of a high-dimensional linear model when outcomes are missing at random. We propose an estimator which combines a Lasso pilot estimate of the regression function with a bias correction term based on the weighted residuals of the Lasso regression. The weights depend on estimates of the missingness probabilities (propensity scores) and solve a convex optimization program that trades off bias and variance optimally. Provided that the propensity scores can be pointwise consistently estimated at in-sample data points, our proposed estimator for the regression function is asymptotically normal and semi-parametrically efficient among all asymptotically linear estimators. Furthermore, the proposed estimator keeps its asymptotic properties even if the propensity scores are estimated by modern machine learning techniques. We validate the finite-sample performance of the proposed estimator through comparative simulation studies and the real-world problem of inferring the stellar masses of galaxies in the Sloan Digital Sky Survey.


title: Into the Mystic: ALMA ACA observations of the Mystic Mountains in Carina

abstract: We present new observations of the Mystic Mountains cloud complex in the Carina Nebula using the ALMA Atacama Compact Array (ACA) to quantify the impact of strong UV radiation on the structure and kinematics of the gas. Our Band~6 observations target CO, $^{13}$CO, and C$^{18}$O; we also detect DCN J=3-2 and $^{13}$CS J=5-4. A dendrogram analysis reveals that the Mystic Mountains are a coherent structure, with continuous emission over $-$10.5 km s$^{-1}$ $<$ v < $-$2 km s$^{-1}$. We perform multiple analyses to isolate non-thermal motions in the Mystic Mountains including computing the turbulent driving parameter, $b$, which indicates whether compressive or solenoidal modes dominate. Each analysis yields values similar to other pillars in Carina that have been observed in a similar way but are subject to an order of magnitude less intense ionizing radiation. We find no clear correlation between the velocity or turbulent structure of the gas and the incident radiation, in contrast to other studies targeting different regions of Carina. This may reflect differences in the initial densities of regions that go on to collapse into pillars and those that still look like clouds or walls in the present day. Pre-existing over-densities that enable pillar formation may also explain why star formation in the pillars appears more evolved (from the presence of jets) than in other heavily-irradiated but non-pillar-like regions. High resolution observations of regions subject to an array of incident radiation are required to test this hypothesis.


title: Limit theorems for critical faces above the vanishing threshold

abstract: We investigate convergence of point processes associated with critical faces for a \v{C}ech filtration built over a homogeneous Poisson point process in the $d$-dimensional flat torus. The convergence of our point process is established in terms of the $\mathcal M_0$-topology, when the connecting radius of a \v{C}ech complex decays to $0$, so slowly that critical faces are even less likely to occur than those in the regime of threshold for homological connectivity. We also obtain a series of limit theorems for positive and negative critical faces, all of which are considerably analogous to those for critical faces.


title: Can the Parker Solar Probe Detect a CME-flare Current Sheet?

abstract: A current sheet (CS) is the central structure in the disrupting magnetic configuration during solar eruptions. More than 90\% of the free magnetic energy (the difference between the energy in the non-potential magnetic field and that in the potential one) stored in the coronal magnetic field beforehand is converted into heating and kinetic energy of the plasma, as well as accelerating charged particles, by magnetic reconnection occurring in the CS. However, the detailed physical properties and fine structures of the CS are still unknown since there is no relevant information obtained via in situ detections. The Parker Solar Probe (PSP) may provide us such information should it traverse a CS in the eruption. The perihelion of PSP's final orbit is located at about 10 solar radii from the center of the Sun, so it can observe the CS at a very close distance, or even traverses the CS, which provides us a unique opportunity to look into fine properties and structures of the CS, helping reveal the detailed physics of large-scale reconnection that was impossible before. We evaluate the probability that PSP can traverse a CS, and examine the orbit of a PSP-like spacecraft that has the highest probability to traverse a CS.


title: Concentration properties of theta lifts on orthogonal groups

abstract: Let $n>m\geqslant 1$ be integers with $n+m\geqslant 4$ even. We prove the existence of Maass forms with large sup norms on anisotropic ${\rm O}(n,m)$, by combining a counting argument with a new period relation showing that a certain orthogonal period on ${\rm O}(n,m)$ distinguishes theta lifts from ${\rm Sp}_{2m}$. This generalizes a method of Rudnick and Sarnak in the rank one case, when $m = 1$. Our lower bound is naturally expressed as a ratio of the Plancherel measures for the groups ${\rm O}(n,m)$ and ${\rm Sp}_{2m}(\mathbb{R})$, up to logarithmic factors, and strengthens the lower bounds of our previous paper for such groups. In the case of odd-dimensional hyperbolic spaces, the growth exponent we obtain improves on a result of Donnelly, and is optimal under the purity conjecture of Sarnak.


title: Negative spectrum of Schr\"odinger Operators with Rapidly Oscillating Potentials

abstract: For Schr\"odinger operators with potentials that are asymptotically homogeneous of degree $-2$, the size of the coupling determines whether it has finite or infinitely many negative eigenvalues. In the latter case the asymptotic accumulation of these eigenvalues at zero has been determined by Kirsch and Simon. A similar regime occurs for potentials which are not asymptotically monotone, but oscillatory. In this case when the ratio between the amplitude and frequency of oscillation is asymptotically homogeneous of degree $-1$ the coupling determines the finiteness of the negative spectrum. We present a new proof of this fact by making use of a ground state representation. As a consequence of this approach we derive an asymptotic formula analogous to that of Kirsch and Simon.


title: Early time dynamics far from equilibrium via holography

abstract: We investigate the early time dynamics of heavy ion collisions studying the time evolution of the energy-momentum tensor as well as energy-momentum correlations within a uniformly thermalizing holographic QGP. From these quantities, we suggest a far-from equilibrium definition of shear viscosity, which is a crucial property of QCD matter as it significantly determines the generation of elliptic flow already at early times. During an exemplary initial heating phase of the holographic QGP the shear viscosity of entropy density ratio decreases down to 60%, followed by an overshoot to 110% of the near-equilibrium value, $\eta/s=1/(4\pi)$. Implications for the QCD QGP are discussed. Subsequently, we consider a holographic QGP which is Bjorken-expanding. Its energy-momentum tensor components have a known hydrodynamic attractor to which all time evolutions collapse independent of the initial conditions. Based on this, we propose a definition for a far from equilibrium speed of sound, and analytically compute its hydrodynamic attractor. Subjecting this Bjorken-expanding plasma to an external magnetic field and an axial chemical potential, we study the chiral magnetic effect far from equilibrium.


title: Holographic Tensor Networks with Bulk Gauge Symmetries

abstract: Tensor networks are useful toy models for understanding the structure of entanglement in holographic states and reconstruction of bulk operators within the entanglement wedge. They are, however, constrained to only prepare so-called "fixed-area states" with flat entanglement spectra, limiting their utility in understanding general features of holographic entanglement. Here, we overcome this limitation by constructing a variant of random tensor networks that enjoys bulk gauge symmetries. Our model includes a gauge theory on a general graph, whose gauge-invariant states are fed into a random tensor network. We show that the model satisfies the quantum-corrected Ryu-Takayanagi formula with a nontrivial area operator living in the center of a gauge-invariant algebra. We also demonstrate nontrivial, n-dependent contributions to the R\'enyi entropy and R\'enyi mutual information from this area operator, a feature shared by general holographic states.


title: Non-constant ground configurations in the disordered ferromagnet

abstract: The disordered ferromagnet is a disordered version of the ferromagnetic Ising model in which the coupling constants are non-negative quenched random. A ground configuration is an infinite-volume configuration whose energy cannot be reduced by finite modifications. It is a long-standing challenge to ascertain whether the disordered ferromagnet on the $\mathbb{Z}^D$ lattice admits non-constant ground configurations. We answer this affirmatively in dimensions $D\ge 4$, when the coupling constants are sampled independently from a sufficiently concentrated distribution. The obtained ground configurations are further shown to be translation-covariant with respect to $\mathbb{Z}^{D-1}$ translations of the disorder. Our result is proved by showing that the finite-volume interface formed by Dobrushin boundary conditions is localized, and converges to an infinite-volume interface. This may be expressed in purely combinatorial terms, as a result on the fluctuations of certain minimal cutsets in the lattice $\mathbb{Z}^D$ endowed with independent edge capacities.


title: Exploring Non-additive Randomness on ViT against Query-Based Black-Box Attacks

abstract: Deep Neural Networks can be easily fooled by small and imperceptible perturbations. The query-based black-box attack (QBBA) is able to create the perturbations using model output probabilities of image queries requiring no access to the underlying models. QBBA poses realistic threats to real-world applications. Recently, various types of robustness have been explored to defend against QBBA. In this work, we first taxonomize the stochastic defense strategies against QBBA. Following our taxonomy, we propose to explore non-additive randomness in models to defend against QBBA. Specifically, we focus on underexplored Vision Transformers based on their flexible architectures. Extensive experiments show that the proposed defense approach achieves effective defense, without much sacrifice in performance.


title: Attention De-sparsification Matters: Inducing Diversity in Digital Pathology Representation Learning

abstract: We propose DiRL, a Diversity-inducing Representation Learning technique for histopathology imaging. Self-supervised learning techniques, such as contrastive and non-contrastive approaches, have been shown to learn rich and effective representations of digitized tissue samples with limited pathologist supervision. Our analysis of vanilla SSL-pretrained models' attention distribution reveals an insightful observation: sparsity in attention, i.e, models tends to localize most of their attention to some prominent patterns in the image. Although attention sparsity can be beneficial in natural images due to these prominent patterns being the object of interest itself, this can be sub-optimal in digital pathology; this is because, unlike natural images, digital pathology scans are not object-centric, but rather a complex phenotype of various spatially intermixed biological components. Inadequate diversification of attention in these complex images could result in crucial information loss. To address this, we leverage cell segmentation to densely extract multiple histopathology-specific representations, and then propose a prior-guided dense pretext task for SSL, designed to match the multiple corresponding representations between the views. Through this, the model learns to attend to various components more closely and evenly, thus inducing adequate diversification in attention for capturing context rich representations. Through quantitative and qualitative analysis on multiple tasks across cancer types, we demonstrate the efficacy of our method and observe that the attention is more globally distributed.


title: LEAP Hand: Low-Cost, Efficient, and Anthropomorphic Hand for Robot Learning

abstract: Dexterous manipulation has been a long-standing challenge in robotics. While machine learning techniques have shown some promise, results have largely been currently limited to simulation. This can be mostly attributed to the lack of suitable hardware. In this paper, we present LEAP Hand, a low-cost dexterous and anthropomorphic hand for machine learning research. In contrast to previous hands, LEAP Hand has a novel kinematic structure that allows maximal dexterity regardless of finger pose. LEAP Hand is low-cost and can be assembled in 4 hours at a cost of 2000 USD from readily available parts. It is capable of consistently exerting large torques over long durations of time. We show that LEAP Hand can be used to perform several manipulation tasks in the real world -- from visual teleoperation to learning from passive video data and sim2real. LEAP Hand significantly outperforms its closest competitor Allegro Hand in all our experiments while being 1/8th of the cost. We release detailed assembly instructions, the Sim2Real pipeline and a development platform with useful APIs on our website at https://leap-hand.github.io/


title: Learning Disentangled Avatars with Hybrid 3D Representations

abstract: Tremendous efforts have been made to learn animatable and photorealistic human avatars. Towards this end, both explicit and implicit 3D representations are heavily studied for a holistic modeling and capture of the whole human (e.g., body, clothing, face and hair), but neither representation is an optimal choice in terms of representation efficacy since different parts of the human avatar have different modeling desiderata. For example, meshes are generally not suitable for modeling clothing and hair. Motivated by this, we present Disentangled Avatars~(DELTA), which models humans with hybrid explicit-implicit 3D representations. DELTA takes a monocular RGB video as input, and produces a human avatar with separate body and clothing/hair layers. Specifically, we demonstrate two important applications for DELTA. For the first one, we consider the disentanglement of the human body and clothing and in the second, we disentangle the face and hair. To do so, DELTA represents the body or face with an explicit mesh-based parametric 3D model and the clothing or hair with an implicit neural radiance field. To make this possible, we design an end-to-end differentiable renderer that integrates meshes into volumetric rendering, enabling DELTA to learn directly from monocular videos without any 3D supervision. Finally, we show that how these two applications can be easily combined to model full-body avatars, such that the hair, face, body and clothing can be fully disentangled yet jointly rendered. Such a disentanglement enables hair and clothing transfer to arbitrary body shapes. We empirically validate the effectiveness of DELTA's disentanglement by demonstrating its promising performance on disentangled reconstruction, virtual clothing try-on and hairstyle transfer. To facilitate future research, we also release an open-sourced pipeline for the study of hybrid human avatar modeling.


title: Geometric semimetals and their simulation in synthetic matter

abstract: Topological semimetals, such as the Weyl and Dirac semimetals, represent one of the most active research fields in modern condensed matter physics. The peculiar physical properties of these systems mainly originate from their underlying symmetries, emergent relativistic dispersion, and band topology. In this Letter, we present a different class of gapless systems in three dimensions, dubbed $\textit{geometric semimetals}$. These semimetals are protected by the generalized chiral and rotation symmetries, but are topologically trivial. Nevertheless, we show that their band geometry is nontrivial, as evidenced by the nonzero quantum metric trace with possible quantization. The possible realization in synthetic-matter experiments is also discussed.


title: DESI luminous red galaxy samples for cross-correlations

abstract: We present two galaxy samples, selected from DESI Legacy Imaging Surveys (LS) DR9, with approximately 20,000 square degrees of coverage and spectroscopic redshift distributions designed for cross-correlations such as with CMB lensing, galaxy lensing, and the Sunyaev-Zel'dovich effect. The first sample is identical to the DESI Luminous Red Galaxy (LRG) sample, and the second sample is an extended LRG sample with 2-3 times the DESI LRG density. We present the improved photometric redshifts, tomographic binning and their spectroscopic redshift distributions and imaging systematics weights, and magnification bias coefficients. The catalogs and related data products will be made publicly available. The cosmological constraints using this sample and Planck lensing maps are presented in a companion paper. We also make public the new set of general-purpose photometric redshifts trained using DESI spectroscopic redshifts, which are used in this work, for all galaxies in LS DR9.


title: Connecting Everyday Objects with the Metaverse: A Unified Recognition Framework

abstract: The recent Facebook rebranding to Meta has drawn renewed attention to the metaverse. Technology giants, amongst others, are increasingly embracing the vision and opportunities of a hybrid social experience that mixes physical and virtual interactions. As the metaverse gains in traction, it is expected that everyday objects may soon connect more closely with virtual elements. However, discovering this "hidden" virtual world will be a crucial first step to interacting with it in this new augmented world. In this paper, we address the problem of connecting physical objects with their virtual counterparts, especially through connections built upon visual markers. We propose a unified recognition framework that guides approaches to the metaverse access points. We illustrate the use of the framework through experimental studies under different conditions, in which an interactive and visually attractive decoration pattern, an Artcode, is used as the approach to enable the connection. This paper will be of interest to, amongst others, researchers working in Interaction Design or Augmented Reality who are seeking techniques or guidelines for augmenting physical objects in an unobtrusive, complementary manner.


title: Correlations of Baryon and Charge Stopping in Heavy Ion Collisions

abstract: Baryon numbers are carried by valence quarks in the standard QCD picture of the baryon structure, while some theory proposed an alternative baryon number carrier, a non-perturbative Y-shaped configuration of the gluon field, called the baryon junction in the 1970s. However, neither of the theories has been verified experimentally. It was recently suggested to search for the baryon junction by investigating the correlation of net-charge and net-baryon yields at midrapidity in heavy-ion collisions. This paper presents studies of such correlations in collisions of various heavy ions from Oxygen to Uranium with the UrQMD Monte Carlo model. The UrQMD model implements valence quark transport as the primary means of charge and baryon stopping at midrapidity. Detailed study are also carried out for isobaric $_{40}^{96}\mathrm{Zr}$ + $_{40}^{96}\mathrm{Zr}$ and $_{44}^{96}\mathrm{Ru}$ + $_{44}^{96}\mathrm{Ru}$ collisions. We found a universal trend of the charge stopping with respect to the baryon stopping, and that the charge stopping is always more than the baryon stopping. This study provides a model baseline in valence quark transport for what is expected in net-charge and net-baryon yields at midrapidity of relativistic heavy-ion collisions.


title: Local Spectral Optimisation for Robin Problems with Negative Boundary Parameter on Quadrilaterals

abstract: We investigate the Robin eigenvalue problem for the Laplacian with negative boundary parameter on quadrilateral domains of fixed area. In this paper, we prove that the square is a local maximiser of the first eigenvalue with respect to the Hausdorff metric. We also provide asymptotic results relating to the optimality of the square for extreme values of the Robin parameter.


title: Matrix-analytic methods for the evolution of species trees, gene trees, and their reconciliation

abstract: We consider the reconciliation problem, in which the task is to find a mapping of a gene tree into a species tree, so as to maximize the likelihood of such fitting, given the available data. We describe a model for the evolution of the species tree, a subfunctionalisation model for the evolution of the gene tree, and provide an algorithm to compute the likelihood of the reconciliation. We derive our results using the theory of matrix-analytic methods and describe efficient algorithms for the computation of a range of useful metrics. We illustrate the theory with examples and provide the physical interpretations of the discussed quantities, with a focus on the practical applications of the theory to incomplete data.


title: Noisy Demkov-Kunike model

abstract: The Demkov-Kunike (DK) model, characterized by a time-dependent Rabi coupling $J~\text{sech}(t/T)$ and on-site detuning $\Delta_0+\Delta_1\tanh(t/T)$, has one of the most general forms of an exactly solvable two-state quantum system, and, therefore, it provides a paradigm for coherent manipulations of a qubit's quantum state. Despite its extensive applications in the noise-free cases, the exploration of the noisy DK model remains limited. Here, we extend the coherent DK model to take into account of a noisy coupling term $J\rightarrow J_{\text{noisy}}(t)$. We consider colored Markovian noise sources represented by the telegraph noise and Gaussian noise. We present exact solutions for the survival probability $Q^{\text{noisy}}_{\text{DK}}$ of the noisy DK model, namely the probability of the system to remain in its initial state. For the slow telegraph noise, we identify parameter regimes where the survival probability $Q^{\text{noisy}}_{\text{DK}}$ is suppressed rather than enhanced by noise. In contrast, for slow Gaussian noise, the noise always enhances the survival probability $Q^{\text{noisy}}_{\text{DK}}$, due to the absorption of noise quanta across the energy gap. This study not only complements the existing research on the noisy Landau-Zener model, but also provides valuable insights for the control of two-level quantum systems.


title: Quantized Non-Volatile Nanomagnetic Synapse based Autoencoder for Efficient Unsupervised Network Anomaly Detection

abstract: In the autoencoder based anomaly detection paradigm, implementing the autoencoder in edge devices capable of learning in real-time is exceedingly challenging due to limited hardware, energy, and computational resources. We show that these limitations can be addressed by designing an autoencoder with low-resolution non-volatile memory-based synapses and employing an effective quantized neural network learning algorithm. We propose a ferromagnetic racetrack with engineered notches hosting a magnetic domain wall (DW) as the autoencoder synapses, where limited state (5-state) synaptic weights are manipulated by spin orbit torque (SOT) current pulses. The performance of anomaly detection of the proposed autoencoder model is evaluated on the NSL-KDD dataset. Limited resolution and DW device stochasticity aware training of the autoencoder is performed, which yields comparable anomaly detection performance to the autoencoder having floating-point precision weights. While the limited number of quantized states and the inherent stochastic nature of DW synaptic weights in nanoscale devices are known to negatively impact the performance, our hardware-aware training algorithm is shown to leverage these imperfect device characteristics to generate an improvement in anomaly detection accuracy (90.98%) compared to accuracy obtained with floating-point trained weights. Furthermore, our DW-based approach demonstrates a remarkable reduction of at least three orders of magnitude in weight updates during training compared to the floating-point approach, implying substantial energy savings for our method. This work could stimulate the development of extremely energy efficient non-volatile multi-state synapse-based processors that can perform real-time training and inference on the edge with unsupervised data.


title: Lambert series in analytic number theory

abstract: Annotated bibliography of 18th, 19th, and early 20th century works involving Lambert series. A tour of 19th and early 20th century analytic number theory.


title: Statistical Applications in Pain Management: A Review

abstract: Chronic pain is a widespread and debilitating condition that affects millions of individuals worldwide. Recent research has unveiled a connection between chronic pain and epigenetic processes, shedding light on the complex interplay between genetics, environment, and pain perception. This review synthesizes findings from several studies exploring the relationship between epigenetic modifications, pain intensity, disability, and various other factors in individuals with chronic pain. The studies encompass a range of chronic pain conditions, including chronic low back pain, knee osteoarthritis pain, and pain in the context of underlying health conditions. The review highlights key findings and implications from each study and discusses their collective contributions to our understanding of the epigenetic mechanisms underpinning chronic pain.


title: Solar Atmospheric Heating Due to Small-scale Events in an Emerging Flux Region

abstract: We investigate the thermal, kinematic and magnetic structure of small-scale heating events in an emerging flux region (EFR). We use high-resolution multi-line observations (including Ca II 8542~\AA, Ca II K, and Fe I 6301~\AA line pair) of an EFR located close to the disk center from the CRISP and CHROMIS instruments at the Swedish 1-m Solar Telescope. We perform non-LTE inversions of multiple spectral lines to infer the temperature, velocity, and magnetic field structure of the heating events. Additionally, we use the data-driven Coronal Global Evolutionary Model to simulate the evolution of the 3D magnetic field configuration above the events and understand their dynamics. Furthermore, we analyze the differential emission measure to gain insights into the heating of the coronal plasma in the EFR. Our analysis reveals the presence of numerous small-scale heating events in the EFR, primarily located at polarity inversion lines of bipolar structures. These events not only heat the lower atmosphere but also significantly heat the corona. The data-driven simulations, along with the observed enhancement of currents and Poynting flux, suggest that magnetic reconnection in the lower atmosphere is likely responsible for the observed heating at these sites.


title: Narrowing the Gap between Supervised and Unsupervised Sentence Representation Learning with Large Language Model

abstract: Sentence Representation Learning (SRL) is a fundamental task in Natural Language Processing (NLP), with the Contrastive Learning of Sentence Embeddings (CSE) being the mainstream technique due to its superior performance. An intriguing phenomenon in CSE is the significant performance gap between supervised and unsupervised methods, with their only difference lying in the training data. Previous works attribute this performance gap to differences in two representation properties (alignment and uniformity). However, since alignment and uniformity only measure the results, they fail to answer "What aspects of the training data contribute to the performance gap?" and "How can the performance gap be narrowed?", In this paper, we conduct empirical experiments to answer these "What" and "How" questions. We first answer the "What" question by thoroughly comparing the behavior of supervised and unsupervised CSE during their respective training processes. From the comparison, we identify the similarity pattern as a key factor to the performance gap, and introduce a metric, called Relative Fitting Difficulty (RFD), to measure the complexity of the similarity pattern. Then, based on the insights gained from the "What" question, we tackle the "How" question by increasing the pattern complexity of the training data. We achieve this by leveraging the In-Context Learning (ICL) capability of the Large Language Model (LLM) to generate data that simulates complex patterns. By utilizing the hierarchical patterns in the LLM-generated data, we effectively narrow the gap between supervised and unsupervised CSE. We release our codes and appendix at https://github.com/BDBC-KG-NLP/NGCSE.


title: Photonsphere, shadow, quasinormal modes, and greybody bounds of non-rotating Simpson-Visser black hole

abstract: In this manuscript, we study photonsphere, shadow, quasinormal modes, Hawking temperature, and greybody bounds of a non-rotating Simpson-Visser black hole which is a regular black hole. We observe that though the radius of the photonsphere does depend on the Simpson-Visser parameter $\alpha$, the shadow radius is independent of it. The shadow radius is found to be equal to that for Schwarzschild black hole. We, then, study quasinormal frequencies of the Simpson-Visser black hole for scalar and electromagnetic perturbations with the help of $6$th order WKB method. We tabulate values of quasinormal frequencies for various values of $\alpha$, angular momentum $\ell$, and overtone number $n$. We also graphically show the dependence of real and imaginary parts of quasinormal frequency on $\alpha$ and $\ell$. Additionally, We study the convergence of the WKB method for various values of pair $(n,\ell)$. Finally, we shed light on the dependence of the Hawking temperature on the parameter $\alpha$ and the dependence of greybody bounds on $\alpha$ and $\ell$.


title: Multimodal Outcomes in N-of-1 Trials: Combining Unsupervised Learning and Statistical Inference

abstract: N-of-1 trials are randomized multi-crossover trials in single participants with the purpose of investigating the possible effects of one or more treatments. Research in the field of N-of-1 trials has primarily focused on scalar outcomes. However, with the increasing use of digital technologies, we propose to adapt this design to multimodal outcomes, such as audio, video, or image data or also sensor measurements, that can easily be collected by the trial participants on their personal mobile devices. We present here a fully automated approach for analyzing multimodal N-of-1 trials by combining unsupervised deep learning models with statistical inference. First, we train an autoencoder on all images across all patients to create a lower-dimensional embedding. In the second step, the embeddings are reduced to a single dimension by projecting on the first principal component, again using all images. Finally, we test on an individual level whether treatment and non-treatment periods differ with respect to the component. We apply our proposed approach to a published series of multimodal N-of-1 trials of 5 participants who tested the effect of creams on acne captured through images over 16 days. We compare several parametric and non-parametric statistical tests, and we also compare the results to an expert analysis that rates the pictures directly with respect to their acne severity and applies a t-test on these scores. The results indicate a treatment effect for one individual in the expert analysis. This effect was replicated with the proposed unsupervised pipeline. In summary, our proposed approach enables the use of novel data types in N-of-1 trials while avoiding the need for manual labels. We anticipate that this can be the basis for further explorations of valid and interpretable approaches and their application in clinical multimodal N-of-1 trials.


title: Searching for sbottom LSP at the LHC

abstract: Assuming that the sbottom is the lightest supersymmetric particle (LSP), we carry out an analysis of the relevant signals expected at the LHC. The discussion is established in the framework of the $\mu\nu$SSM, where the presence of $R$-parity violating couplings involving right-handed neutrinos solves simultaneously the $\mu$-problem and the accommodation of neutrino masses and mixing angles. The sbottoms are pair produced at the LHC, decaying to a lepton and a top quark or a neutrino and a bottom quark. The decays can be prompt or displaced, depending on the regions of the parameter space of the model. We focus the analysis on the right sbottom LSP, since the left sbottom is typically heavier than the left stop because of the D-term contribution. We compare the predictions of this scenario with ATLAS and CMS searches for prompt and long-lived particles. To analyze the parameter space we sample the $\mu\nu$SSM for a right sbottom LSP, paying special attention to reproduce the current experimental data on neutrino and Higgs physics, as well as flavor observables. For displaced (prompt) decays, our results translate into lower limits on the mass of the right sbottom LSP of about $1041$ GeV ($1070$ GeV). The largest possible value found for the decay length is about $3.5$ mm.


title: Opportunistic Reflection in Reconfigurable Intelligent Surface-Assisted Wireless Networks

abstract: This paper focuses on multiple-access protocol design in a wireless network assisted by multiple reconfigurable intelligent surfaces (RISs). By extending the existing approaches in single-user or single-RIS cases, we present two benchmark schemes for this multi-user multi-RIS scenario. Inspecting their shortcomings, a simple but efficient method coined opportunistic multi-user reflection (OMUR) is proposed. The key idea is to opportunistically select the best user as the anchor for optimizing the RISs, and non-orthogonally transmitting all users' signals simultaneously. A simplified version of OMUR exploiting random phase shifts is also proposed to avoid the complexity of RIS channel estimation.


title: Quantum multi-secret sharing scheme with access structures and cheat identification

abstract: This work proposes a $d$-dimensional quantum multi-secret sharing scheme with a cheat detection mechanism. The dealer creates multiple secrets and distributes the shares of these secrets using multi-access structures and a monotone span program. The dealer detects the cheating of each participant using the Black box's cheat detection mechanism. To detect the participants' deceit, the dealer distributes secret shares' shadows derived from a randomly invertible matrix $X$ to the participants, stored in the black box. The Black box identifies the participant's deceitful behavior during the secret recovery phase. Only honest participants authenticated by the Black box acquire their secret shares to recover the multiple secrets. After the Black box cheating verification, the participants reconstruct the secrets by utilizing the unitary operations and quantum Fourier transform. The proposed protocol is reliable in preventing attacks from eavesdroppers and participants. The scheme's efficiency is demonstrated in different noise environments: dit-flip noise, $d$-phase-flip noise, and amplitude-damping noise, indicating its robustness in practical scenarios. The proposed protocol provides greater versatility, security, and practicality.


title: Sensitivity Analysis for Quantiles of Hidden Biases in Matched Observational Studies

abstract: In matched observational studies, the inferred causal conclusions pretending that matching has taken into account all confounding can be sensitive to unmeasured confounding. In such cases, a sensitivity analysis is often conducted, which investigates whether the observed association between treatment and outcome is due to effects caused by the treatment or it is due to hidden confounding. In general, a sensitivity analysis tries to infer the minimum amount of hidden biases needed in order to explain away the observed association between treatment and outcome, assuming that the treatment has no effect. If the needed bias is large, then the treatment is likely to have significant effects. The Rosenbaum sensitivity analysis is a modern approach for conducting sensitivity analysis for matched observational studies. It investigates what magnitude the maximum of the hidden biases from all matched sets needs to be in order to explain away the observed association, assuming that the treatment has no effect. However, such a sensitivity analysis can be overly conservative and pessimistic, especially when the investigators believe that some matched sets may have exceptionally large hidden biases. In this paper, we generalize Rosenbaum's framework to conduct sensitivity analysis on quantiles of hidden biases from all matched sets, which are more robust than the maximum. Moreover, we demonstrate that the proposed sensitivity analysis on all quantiles of hidden biases is simultaneously valid and is thus a free lunch added to the conventional sensitivity analysis. The proposed approach works for general outcomes, general matched studies and general test statistics. Finally, we demonstrate that the proposed sensitivity analysis also works for bounded null hypotheses as long as the test statistic satisfies certain properties. An R package implementing the proposed method is also available online.


title: Widely Interpretable Semantic Representation: Frameless Meaning Representation for Broader Applicability

abstract: This paper presents a novel semantic representation, WISeR, that overcomes challenges for Abstract Meaning Representation (AMR). Despite its strengths, AMR is not easily applied to languages or domains without predefined semantic frames, and its use of numbered arguments results in semantic role labels, which are not directly interpretable and are semantically overloaded for parsers. We examine the numbered arguments of predicates in AMR and convert them to thematic roles that do not require reference to semantic frames. We create a new corpus of 1K English dialogue sentences annotated in both WISeR and AMR. WISeR shows stronger inter-annotator agreement for beginner and experienced annotators, with beginners becoming proficient in WISeR annotation more quickly. Finally, we train a state-of-the-art parser on the AMR 3.0 corpus and a WISeR corpus converted from AMR 3.0. The parser is evaluated on these corpora and our dialogue corpus. The WISeR model exhibits higher accuracy than its AMR counterpart across the board, demonstrating that WISeR is easier for parsers to learn.


title: Moments of $L$-functions via the relative trace formula

abstract: We prove an asymptotic formula for the second moment of the $\mathrm{GL}(n)\times\mathrm{GL}(n+1)$ Rankin--Selberg central $L$-values $L(1/2,\Pi\otimes\pi)$, where $\pi$ is a fixed cuspidal representation of $\mathrm{GL}(n)$ that is tempered and unramified at every place, while $\Pi$ varies over a family of automorphic representations of $\mathrm{PGL}(n+1)$ ordered by (archimedean or non-archimedean) conductor. As another application of our method, we prove the existence of infinitely many cuspidal representations $\Pi$ of $\mathrm{PGL}(n+1)$ such that $L(1/2,\Pi\otimes\pi_1)$ and $L(1/2,\Pi\otimes\pi_2)$ do not vanish simultaneously where $\pi_1$ and $\pi_2$ are cuspidal representations of $\mathrm{GL}(n)$ that are unramified and tempered at every place and have trivial central characters.


title: Action Segmentation Using 2D Skeleton Heatmaps and Multi-Modality Fusion

abstract: This paper presents a 2D skeleton-based action segmentation method with applications in fine-grained human activity recognition. In contrast with state-of-the-art methods which directly take sequences of 3D skeleton coordinates as inputs and apply Graph Convolutional Networks (GCNs) for spatiotemporal feature learning, our main idea is to use sequences of 2D skeleton heatmaps as inputs and employ Temporal Convolutional Networks (TCNs) to extract spatiotemporal features. Despite lacking 3D information, our approach yields comparable/superior performances and better robustness against missing keypoints than previous methods on action segmentation datasets. Moreover, we improve the performances further by using both 2D skeleton heatmaps and RGB videos as inputs. To our best knowledge, this is the first work to utilize 2D skeleton heatmap inputs and the first work to explore 2D skeleton+RGB fusion for action segmentation.


title: Neutrino mass models at $\mu$TRISTAN

abstract: We study the prospects of probing neutrino mass models at the newly proposed antimuon collider $\mu$TRISTAN, involving $\mu^+e^-$ scattering at $\sqrt{s}= 346$ GeV and $\mu^+\mu^+$ scattering at $\sqrt{s}= 2$ TeV. We show that $\mu$TRISTAN is uniquely sensitive to leptophilic neutral and doubly-charged scalars naturally occurring in various neutrino mass models, such as Zee, Zee-Babu, cocktail, and type-II seesaw models, over a wide range of mass and coupling values, well beyond the current experimental constraints. It also allows for the possibility to correlate the collider signals with neutrino mixing parameters and charged lepton flavor violating observables.


title: Bootstrapping the Stochastic Resonance

abstract: Stochastic resonance is a phenomenon where a noise of appropriate intensity enhances the input signal strength. In this work, by employing the recently developed convex optimization methods in the context of dynamical systems and stochastic processes, we derive rigorous two-sided bounds on the expected power at the input signal frequency for the prototypical example of stochastic resonance, the double-well potential with periodic forcing and Gaussian white noise.


title: New bounds on light millicharged particles from the tip of the red-giant branch

abstract: Stellar energy loss is a sensitive probe of light, weakly coupled dark sectors, including ones containing millicharged particles (MCPs). The emission of MCPs can affect stellar evolution, and therefore can alter the observed properties of stellar populations. In this work, we improve upon the accuracy of existing stellar limits on MCPs by self-consistently modelling (1) the MCP emission rate, accounting for all relevant in-medium effects and production channels, and (2) the evolution of stellar interiors (including backreactions from MCP emission) using the MESA stellar evolution code. We find MCP emission leads to significant brightening of the tip of the red-giant branch. Based on photometric observations of 15 globular clusters whose bolometric magnitudes are inferred using parallaxes from Gaia astrometry, we obtain robust bounds on the existence of MCPs with masses below 100 keV.


title: The integrated galaxy-wide stellar initial mass function over the radial acceleration range of early-type galaxies

abstract: The observed radial accelerations of 462 Early-type galaxies (ETGs) at their half-mass radii are discussed. They are compared to the baryonic masses of the same galaxies, which are derived from theoretical expectations for their stellar populations and cover a range from $\approx 10^4 \, {\rm M}_{\odot}$ to $\approx 10^{11} \, {\rm M}_{\odot}$. Both quantities are plotted against each other, and it is tested whether they lie (within errors) along theoretical radial acceleration relations (RARs). We choose the Newtonian RAR and two Milgromian, or MONDian RARs. At low radial accelerations (corresponding to low masses), the Newtonian RAR fails without non-baryonic dark matter, but the two MONDian ones may work, provided moderate out-of-equilibrium dynamics in some of the low-mass ETGs. However all three RARs fail at high accelerations (corresponding to high masses) if all ETGs have formed their stellar populations with the canonical stellar initial mass function (IMF). A much better agreement with the observations can however be accomplished, if the theory of the integrated galaxy-wide stellar initial mass functions (IGIMFs) is used instead. This is because the IGIMF-theory predicts the formation of an overabundance of stellar remnants during the lifetime of the massive ETGs. Thus their baryonic masses today are higher than they would be if the ETGs had formed with a canonical IMF. Also the masses of the stellar-mass black holes should be rather high, which would mean that most of them probably formed when the massive ETGs were not as metal-enriched as they are today. The IGIMF-approach confirms downsizing.


title: Planetesimals at DZ stars I: chondritic compositions and a massive accretion event

abstract: There is a wealth of evidence to suggest that planetary systems can survive beyond the main sequence. Most commonly, white dwarfs are found to be accreting material from tidally disrupted asteroids, whose bulk compositions are reflected by the metals polluting the stellar photospheres. While many examples are known, most lack the deep, high-resolution data required to detect multiple elements, and thus characterise the planetesimals that orbit them. Here, spectra of seven DZ white dwarfs observed with Keck HIRES are analysed, where up to nine metals are measured per star. Their compositions are compared against those of solar system objects, working in a Bayesian framework to infer or marginalise over the accretion history. All of the stars have been accreting primitive material, similar to chondrites, with hints of a Mercury-like composition at one star. The most polluted star is observed several Myr after its last major accretion episode, in which a Moon-sized object met its demise.


title: Anatomy of the eigenstates distribution: a quest for a genuine multifractality

abstract: Motivated by a series of recent works, an interest in multifractal phases has risen as they are believed to be present in the Many-Body Localized (MBL) phase and are of high demand in quantum annealing and machine learning. Inspired by the success of the RosenzweigPorter (RP) model with Gaussian-distributed hopping elements, several RP-like ensembles with the fat-tailed distributed hopping terms have been proposed, with claims that they host the desired multifractal phase. In the present work, we develop a general (graphical) approach allowing a self-consistent analytical calculation of fractal dimensions for a generic RP model and investigate what features of the RP Hamiltonians can be responsible for the multifractal phase emergence. We conclude that the only feature contributing to a genuine multifractality is the on-site energies' distribution, meaning that no random matrix model with a statistically homogeneous distribution of diagonal disorder and uncorrelated off-diagonal terms can host a multifractal phase.

categories: cond-mat.dis-nn cond-mat.stat-mech physics.data-an quant-ph


title: Holographic Thermal Observables and M2-branes

abstract: We use holography in conjunction with recent results from supersymmetric localization to compute certain thermal observables for 3d $\mathcal{N}=2$ holographic SCFTs arising on the worldvolume of $N$ M2-branes. We obtain results for the thermal free energy density on $S^1 \times \mathbb{R}^2$, the Casimir energy on $T^{2} \times \mathbb{R}$, and the three leading coefficients in the large temperature limit of the free energy on $S^1\times S^2$ valid to subleading order in the large $N$ limit. As a byproduct of our holographic analysis we also present a conjecture for the structure of the large temperature expansion of the thermal free energy of general 3d CFTs on $S^1\times S^2$.


title: Metal enrichment and evolution in four z > 6.5 quasar sightlines observed with JWST/NIRSpec

abstract: We present JWST/NIRSpec R~2700 spectra of four high-redshift quasars: VDES J0020-3653 (z = 6.860), DELS J0411-0907 (z = 6.825), UHS J0439+1634 (z = 6.519) and ULAS J1342+0928 (z = 7.535). The exquisite data quality, signal-to-noise ratio of 50-200, and large $0.86\!~\mu{\rm m}\le \lambda \le 5.5\!~\mu{\rm m}$ spectral coverage allows us to identify between 13 and 17 intervening and proximate metal absorption line systems in each quasar spectrum, with a total number of 61 absorption-line systems detected at 2.42<z<7.48 including the highest redshift intervening OI 1302 and MgII systems at z=7.37 and z=7.44. We investigate the evolution of the metal enrichment in the epoch of reionization at z>6 and find: i) A continued increase of the low-ionization OI, CII, and SiII incidence, ii) Decreasing high-ionization CIV and SiIV incidence with a transition from predominantly high- to low-ionization at $z\approx6.0$, and iii) a constant MgII incidence across all redshifts. The observations support a change in the ionization state of the intergalactic medium in the EoR rather than a change in metallicity. The abundance ratio of [Si/O] in five z>6 absorption systems show enrichment signatures produced by low-mass Pop III pair instability supernovae, and possibly Pop III hypernovae. In the Gunn-Peterson troughs we detect transmission spikes where Ly$\alpha$ photons can escape. From 22 absorption systems at z>5.7, only a single low-ionization system out of 13 lies within 2000 km/s from a spike, while four high-ionization systems out of nine lie within ~2000 km/s from a spike. This confirms that galaxies responsible for the heavy elements that are transported into the circumgalactic medium lie in predominantly in high-density, neutral environments, while lower density environments are ionized without being polluted by metals at $z\approx$ 6-7. [abridged]


title: Transient Corotating Clumps Around Adolescent Low-Mass Stars From Four Years of TESS

abstract: Complex periodic variables (CPVs) are stars that exhibit highly structured and periodic optical light curves. Previous studies have indicated that these stars are typically disk-free pre-main-sequence M dwarfs with rotation periods ranging from 0.2 to 2 days. To advance our understanding of these enigmatic objects, we conducted a blind search using TESS 2-minute data of 65,760 K and M dwarfs with $T$<16 mag and $d$<150 pc. We found 50 high-quality CPVs, and subsequently determined that most are members of stellar associations. Among the new discoveries are the brightest ($T$$\approx$9.5 mag), closest ($d$$\approx$20 pc), and oldest ($\approx$200 Myr) CPVs known. One exceptional object, LP 12-502, exhibited up to eight flux dips per cycle. Some of these dips coexisted with slightly different periods, and the shortest-duration dips precisely matched the expected timescale for transiting small bodies at the corotation radius. Broadly, our search confirms that CPVs are mostly young ($\lesssim$150 Myr) and low-mass ($\lesssim$0.4 $M_\odot$). The flux dips characteristic of the class have lifetimes of $\approx$100 cycles, although stellar flares seem to induce sudden dip collapse once every few months. The most plausible explanation for these phenomena remains corotating concentrations of gas or dust. The gas or dust is probably entrained by the star's magnetic field, and the sharp features could result from a multipolar field topology, a hypothesis supported by correspondences between the light curves of CPVs and of rapidly rotating B stars known to have multipolar magnetic fields.


title: Flows for Flows: Morphing one Dataset into another with Maximum Likelihood Estimation

abstract: Many components of data analysis in high energy physics and beyond require morphing one dataset into another. This is commonly solved via reweighting, but there are many advantages of preserving weights and shifting the data points instead. Normalizing flows are machine learning models with impressive precision on a variety of particle physics tasks. Naively, normalizing flows cannot be used for morphing because they require knowledge of the probability density of the starting dataset. In most cases in particle physics, we can generate more examples, but we do not know densities explicitly. We propose a protocol called flows for flows for training normalizing flows to morph one dataset into another even if the underlying probability density of neither dataset is known explicitly. This enables a morphing strategy trained with maximum likelihood estimation, a setup that has been shown to be highly effective in related tasks. We study variations on this protocol to explore how far the data points are moved to statistically match the two datasets. Furthermore, we show how to condition the learned flows on particular features in order to create a morphing function for every value of the conditioning feature. For illustration, we demonstrate flows for flows for toy examples as well as a collider physics example involving dijet events


title: Axion Paradigm with Color-Mediated Neutrino Masses

abstract: We propose a generalized Kim-Shifman-Vainshtein-Zakharov-type axion framework in which colored fermions and scalars act as two-loop Majorana neutrino-mass mediators. The global Peccei-Quinn symmetry under which exotic fermions are charged solves the strong CP problem. Within our general proposal, various setups can be distinguished by probing the axion-to-photon coupling at helioscopes and haloscopes. We also comment on axion dark-matter production in the early Universe.


title: Spectroscopic analysis of hot, massive stars in large spectroscopic surveys with de-idealised models

abstract: Upcoming large-scale spectroscopic surveys with e.g. WEAVE and 4MOST will provide thousands of spectra of massive stars, which need to be analysed in an efficient and homogeneous way. Usually, studies of massive stars are limited to samples of a few hundred objects which pushes current spectroscopic analysis tools to their limits because visual inspection is necessary to verify the spectroscopic fit. Often uncertainties are only estimated rather than derived and prior information cannot be incorporated without a Bayesian approach. In addition, uncertainties of stellar atmospheres and radiative transfer codes are not considered as a result of simplified, inaccurate or incomplete/missing physics or, in short, idealised physical models. Here, we address the question of "How to compare an idealised model of complex objects to real data?" with an empirical Bayesian approach and maximum a {\it posterior} approximations. We focus on application to large scale optical spectroscopic studies of complex astrophysical objects like stars. More specifically, we test and verify our methodology on samples of OB stars in 30 Doradus region of the Large Magellanic Clouds using a grid of FASTWIND model atmospheres. Our spectroscopic model de-idealisation analysis pipeline takes advantage of the statistics that large samples provide by determining the model error to account for the idealised stellar atmosphere models, which are included into the error budget. The pipeline performs well over a wide parameter space and derives robust stellar parameters with representative uncertainties.