title: Statistical Analysis of the Radial Evolution of the Solar Winds between 0.1 and 1 au, and their Semi-empirical Iso-poly Fluid Modeling

abstract: Statistical classification of the Helios solar wind observations into several populations sorted by bulk speed has revealed an outward acceleration of the wind. The faster the wind is, the smaller is this acceleration in the 0.3 - 1 au radial range (Maksimovic et al. 2020). In this article we show that recent measurements from the Parker Solar Probe (PSP) are compatible with an extension closer to the Sun of the latter Helios classification. For instance the well established bulk speed/proton temperature (u,Tp) correlation and bulk speed/electron temperature (u,Te) anti-correlation, together with the acceleration of the slowest winds, are verified in PSP data. We also model the combined PSP and Helios data, using empirical Parker-like models for which the solar wind undergoes an "iso-poly" expansion: isothermal in the corona, then polytropic at distances larger than the sonic point radius. The polytropic indices are derived from the observed temperature and density gradients. Our modelling reveals that the electron thermal pressure has a major contribution in the acceleration process of slow and intermediate winds (in the range of 300-500 km/s at 1 au), over a broad range of distances and that the global (electron and protons) thermal energy, alone, is able to explain the acceleration profiles. Moreover, we show that the very slow solar wind requires in addition to the observed pressure gradients, another source of acceleration.


title: Big Learning

abstract: Recent advances in big/foundation models reveal a promising path for deep learning, where the roadmap steadily moves from big data to big models to (the newly-introduced) big learning. Specifically, the big learning exhaustively exploits the information inherent in its large-scale complete/incomplete training data, by simultaneously modeling many/all joint/conditional/marginal data distributions across potentially diverse domains, with one universal foundation model. We reveal that big learning ($i$) underlies most existing foundation models, ($ii$) is equipped with extraordinary flexibilities for complete/incomplete training data and trustworthy data tasks, ($iii$) is capable of delivering all joint/conditional/marginal data capabilities with one universal model, and ($iv$) unifies conventional machine learning paradigms and enables their flexible cooperations, manifested as a universal learning paradigm. Diverse experiments are carried out to validate the effectiveness of the presented big learning.


title: An investigation of Hertzian contact in soft materials using photoelastic tomography

abstract: Hertzian contact of a rigid sphere and a highly deformable soft solid is investigated using integrated photoelasticity. The experiments are performed by pressing a styrene sphere of 15 mm diameter against a 44 x 44 x 47 mm$^3$ cuboid made of 5% wt. gelatin, inside a circular polariscope, and with a range of forces. The emerging light rays are processed by considering that the retardation of each ray carries the cumulative effect of traversing the contact-induced axisymmetric stress field. Then, assuming Hertzian theory is valid, the retardation is analytically calculated for each ray and compared to the experimental one. Furthermore, a finite element model of the process introduces the effect of finite displacements and strains. Beyond the qualitative comparison of the retardation fields, the experimental, theoretical, and numerical results are quantitatively compared in terms of the maximum equivalent stress, surface displacement, and contact radius dimensions. A favorable agreement is found at lower force levels, where the assumptions of Hertz theory hold, whereas deviations are observed at higher force levels. A major discovery of this work is that at the maximum equivalent stress location, all three components of principal stress can be determined experimentally, and show satisfactory agreement with theoretical and numerical ones in our measurement range. This provides valuable insight into Hertzian contact problems since the maximum equivalent stress controls the initiation of plastic deformation or failure. The measured displacement and contact radii also reasonably agree with the theoretical and numerical ones. Finally, the limitations that arise due to the linearization of this problem are explored.


title: Reproducing sensory induced hallucinations via neural fields

abstract: Understanding sensory-induced cortical patterns in the primary visual cortex V1 is an important challenge both for physiological motivations and for improving our understanding of human perception and visual organisation. In this work, we focus on pattern formation in the visual cortex when the cortical activity is driven by a geometric visual hallucination-like stimulus. In particular, we present a theoretical framework for sensory-induced hallucinations which allows one to reproduce novel psychophysical results such as the MacKay effect (Nature, 1957) and the Billock and Tsou experiences (PNAS, 2007).


title: Interaction Pattern Disentangling for Multi-Agent Reinforcement Learning

abstract: Deep cooperative multi-agent reinforcement learning has demonstrated its remarkable success over a wide spectrum of complex control tasks. However, recent advances in multi-agent learning mainly focus on value decomposition while leaving entity interactions still intertwined, which easily leads to over-fitting on noisy interactions between entities. In this work, we introduce a novel interactiOn Pattern disenTangling (OPT) method, to disentangle the entity interactions into interaction prototypes, each of which represents an underlying interaction pattern within a subgroup of the entities. OPT facilitates filtering the noisy interactions between irrelevant entities and thus significantly improves generalizability as well as interpretability. Specifically, OPT introduces a sparse disagreement mechanism to encourage sparsity and diversity among discovered interaction prototypes. Then the model selectively restructures these prototypes into a compact interaction pattern by an aggregator with learnable weights. To alleviate the training instability issue caused by partial observability, we propose to maximize the mutual information between the aggregation weights and the history behaviors of each agent. Experiments on single-task, multi-task and zero-shot benchmarks demonstrate that the proposed method yields results superior to the state-of-the-art counterparts. Our code is available at https://github.com/liushunyu/OPT.


title: Discovery of five cyano derivatives of propene with the QUIJOTE line survey

abstract: We report the discovery of five cyano derivatives of propene towards TMC-1 with the QUIJOTE line survey: $trans$ and $cis$-crotononitrile ($t$-CH$_3$CHCHCN, $c$-CH$_3$CHCHCN), methacrylonitrile (CH$_2$C(CH$_3$)CN), and $gauche$ and $cis$-allyl cyanide ($g$-CH$_2$CHCH$_2$CN and $c$-CH$_2$CHCH$_2$CN). The observed transitions allowed us to derive a common rotational temperature of 7$\pm$1 K for all them. The derived column densities are N($t$-CH$_3$CHCHCN)=(5$\pm$0.5)$\times$10$^{10}$ cm$^{-2}$, N($c$-CH$_3$CHCHCN)=(1.3$\pm$0.2)$\times$10$^{11}$ cm$^{-2}$, N(CH$_2$C(CH$_3$)CN)=(1.0$\pm$0.1)$\times$10$^{11}$ cm$^{-2}$, N($g$-CH$_2$CHCH$_2$CN)=(8.0$\pm$0.8)$\times$10$^{10}$ cm$^{-2}$, and N($c$-CH$_2$CHCH$_2$CN)=(7.0$\pm$0.7)$\times$10$^{10}$ cm$^{-2}$, respectively. The abundance of cyano-propene relative to that of propene is thus $\sim$10$^{-2}$, which is considerably lower than those of other cyano derivatives of abundant hydrocarbons. Upper limits are obtained for two ethynyl derivatives of propene ($E$ and $Z$-CH$_3$CHCHCCH).


title: Privacy Preservation by Local Design in Cooperative Networked Control Systems

abstract: In this paper, we study the privacy preservation problem in a cooperative networked control system working for the task of LQG control. The system consists of a user and a server: the user owns the state process, the server provides computation capability, and the user employs the server to compute control inputs for it. To enable the server's computation, the user needs to provide the measurements of the process states to the server. However, the user regards the states as privacy and wants the server to have deviated knowledge of the state estimates rather than the true value. Regarding that, we propose a novel local design of privacy preservation, which creates a deviation in the server's knowledge of the state estimates from the true value. Meanwhile, we propose an associated privacy metric to measure the deviation and analyze the privacy performance accordingly. Furthermore, we analyze the service loss in the LQG control performance caused by the privacy scheme. Finally, we study the performance trade-off between the privacy preservation and the LQG control, where the proposed optimization problems are solved by numerical methods efficiently.


title: Towards tangled finite element analysis over partially inverted hexahedral elements

abstract: If a finite element mesh contains concave elements, it is said to tangled. Tangled meshes can occur during mesh generation, mesh optimization, and large deformation simulations, and will lead to erroneous results during finite element analysis. Recently, the authors introduced the tangled finite element method (TFEM) to accurately handle tangled 2D concave quadrilateral elements. In this paper, TFEM is extended to 3D hexahedral elements. We demonstrate that TFEM leads to accurate results over tangled hexahedral elements, requiring minimal changes to existing FEM framework.


title: Predicting binding energies of astrochemically relevant molecules via machine learning

abstract: The behaviour of molecules in space is to a large extent governed by where they freeze out or sublimate. The molecular binding energy is thus an important parameter for many astrochemical studies. This parameter is usually determined with time-consuming experiments, computationally expensive quantum chemical calculations, or the inexpensive, but inaccurate, linear addition method. In this work we propose a new method based on machine learning for predicting binding energies that is accurate, yet computationally inexpensive. A machine learning model based on Gaussian Process Regression is created and trained on a database of binding energies of molecules collected from laboratory experiments presented in the literature. The molecules in the database are categorized by their features, such as mono- or multilayer coverage, binding surface, functional groups, valence electrons, and H-bond acceptors and donors. The performance of the model is assessed with five-fold and leave-one-molecule-out cross validation. Predictions are generally accurate, with differences between predicted and literature binding energies values of less than $\pm$20\%. The validated model is used to predict the binding energies of twenty one molecules that have recently been detected in the interstellar medium, but for which binding energy values are not known. A simplified model is used to visualize where the snowlines of these molecules would be located in a protoplanetary disk. This work demonstrates that machine learning can be employed to accurately and rapidly predict binding energies of molecules. Machine learning complements current laboratory experiments and quantum chemical computational studies. The predicted binding energies will find use in the modelling of astrochemical and planet-forming environments.


title: Anomalous Fano factor as a signature of Bogoliubov Fermi surfaces

abstract: Noise spectroscopy is a key technique to investigate the nature and dynamics of charge carriers in superconductors. The recently discovered superconducting hybrids with Bogoliubov Fermi surfaces exhibit a particularly intriguing and rich charge dynamics, as their charge carriers consist of both Cooper pairs and an extensive number of Bogoliubov quasiparticles. Motivated by this, we compute the noise spectra of Bogoliubov Fermi surfaces and identify their key signatures in the differential conductance and the Fano factor. Specifically, we consider a semiconductor/superconductor hybrid device with an in-plane magnetic field, which exhibits several Bogoliubov Fermi surfaces. The number and orientation of the Bogoliubov Fermi surfaces in this device can be readily controlled by the applied magnetic field, which in turn alters the noise signal. In particular, we find that the Fano factor exhibits a reduced value, substantially lower than two, whenever the charge dynamics is governed by a large number of Bogoliubov quasiparticles. Using experimentally relevant parameters, we make a number of specific predictions for the noise spectra, that can be used as direct evidence of Bogoliubov Fermi surfaces. In particular, we find that the Fano factor as a function of magnetic field and spin-orbit coupling exhibits characteristic discontinuities at the transition lines that separate phases with different number of Bogoliubov Fermi surfaces.


title: Continuous Nakayama Representations

abstract: We introduce continuous analogues of Nakayama algebras. In particular, we introduce the notion of (pre-)Kupisch functions, which play a role as Kupisch series of Nakayama algebras, and view continuous Nakayama representations as a special type of representation of $\mathbb{R}$ or $\mathbb{S}^1$. We investigate equivalences and connectedness of the categories of Nakayama representations. Specifically, we prove that orientation-preserving homeomorphisms on $\mathbb{R}$ and on $\mathbb{S}^1$ induce equivalences between these categories. Connectedness is characterized by a special type of points called separation points determined by (pre-)Kupisch functions. We also construct an exact embedding from the category of finite-dimensional representations for any finite-dimensional Nakayama algebra, to a category of continuous Nakayama representaitons.


title: Ensemble random forest filter: An alternative to the ensemble Kalman filter for inverse modeling

abstract: The ensemble random forest filter (ERFF) is presented as an alternative to the ensemble Kalman filter (EnKF) for the purpose of inverse modeling. The EnKF is a data assimilation approach that forecasts and updates parameter estimates sequentially in time as observations are being collected. The updating step is based on the experimental covariances computed from an ensemble of realizations and the updates are given as linear combinations of the differences between observations and forecasted system state values. The ERFF replaces the linear combination in the update step with a non-linear function represented by a random forest. In this way, the non-linear relationships between the parameters to be updated and the observations can be captured and a better update produced. The ERFF is demonstrated for the purpose of log-conductivity identification from piezometric head observations in a number of scenarios with varying degrees of heterogeneity (log-conductivity variances going from 1 up to 6.25 (ln m/d)2), number of realizations in the ensemble (50 or 100), and number of piezometric head observations (18 or 36). In all scenarios, the ERFF works well, being able to reconstruct the log-conductivity spatial heterogeneity while matching the observed piezometric heads at selected control points. For benchmarking purposes the ERFF is compared to the restart EnKF to find that the ERFF is superior to the EnKF for the number of ensemble realizations used (small in typical EnKF applications). Only when the number of realizations grows to 500, the restart EnKF is able to match the performance of the ERFF, albeit at triple the computational cost.


title: Temporal patterns of reciprocity in communication networks

abstract: Human communication, the essence of collective social phenomena ranging from small-scale organizations to worldwide online platforms, features intense reciprocal interactions between members in order to achieve stability, cohesion, and cooperation in social networks. While high levels of reciprocity are well known in aggregated communication data, temporal patterns of reciprocal information exchange have received far less attention. Here we propose measures of reciprocity based on the time ordering of interactions and explore them in data from multiple communication channels, including calls, messaging and social media. By separating each channel into reciprocal and non-reciprocal temporal networks, we find persistent trends that point to the distinct roles of one-on-one exchange versus information broadcast. We implement several null models of communication activity, which identify memory, a higher tendency to repeat interactions with past contacts, as a key source of reciprocity. When adding memory to a model of activity-driven, time-varying networks, we reproduce the levels of reciprocity seen in empirical data. Our work adds to the theoretical understanding of the emergence of reciprocity in human communication systems, hinting at the mechanisms behind the formation of norms in social exchange and large-scale cooperation.


title: Two long-period transiting exoplanets on eccentric orbits: NGTS-20 b (TOI-5152 b) and TOI-5153 b

abstract: Long-period transiting planets provide the opportunity to better understand the formation and evolution of planetary systems. Their atmospheric properties remain largely unaltered by tidal or radiative effects of the host star, and their orbital arrangement reflects a different, and less extreme, migrational history compared to close-in objects. The sample of long-period exoplanets with well determined masses and radii is still limited, but a growing number of long-period objects reveal themselves in the TESS data. Our goal is to vet and confirm single transit planet candidates detected in the TESS space-based photometric data through spectroscopic and photometric follow up observations with ground-based instruments. We use the Next Generation Transit Survey (NGTS) to photometrically monitor the candidates in order to observe additional transits. We report the discovery of two massive, warm Jupiter-size planets, one orbiting the F8-type star TOI-5153 and the other orbiting the G1-type star NGTS-20 (=TOI-5152). From our spectroscopic analysis, both stars are metal-rich with a metallicity of 0.12 and 0.15, respectively. Follow-up radial velocity observations were carried out with CORALIE, CHIRON, FEROS, and HARPS. TOI-5153 hosts a 20.33 day period planet with a planetary mass of 3.26 (+-0.18) Mj, a radius of 1.06 (+-0.04) Rj , and an orbital eccentricity of 0.091 (+-0.026). NGTS-20 b is a 2.98 (+-0.16) Mj planet with a radius of 1.07 (+-0.04) Rj on an eccentric (0.432 +- 0.023) orbit with an orbital period of 54.19 days. Both planets are metal-enriched and their heavy element content is in line with the previously reported mass-metallicity relation for gas giants. Both warm Jupiters orbit moderately bright host stars making these objects valuable targets for follow-up studies of the planetary atmosphere and measurement of the spin-orbit angle of the system.


title: A Mask Attention Interaction and Scale Enhancement Network for SAR Ship Instance Segmentation

abstract: Most of existing synthetic aperture radar (SAR) ship in-stance segmentation models do not achieve mask interac-tion or offer limited interaction performance. Besides, their multi-scale ship instance segmentation performance is moderate especially for small ships. To solve these problems, we propose a mask attention interaction and scale enhancement network (MAI-SE-Net) for SAR ship instance segmentation. MAI uses an atrous spatial pyra-mid pooling (ASPP) to gain multi-resolution feature re-sponses, a non-local block (NLB) to model long-range spa-tial dependencies, and a concatenation shuffle attention block (CSAB) to improve interaction benefits. SE uses a content-aware reassembly of features block (CARAFEB) to generate an extra pyramid bottom-level to boost small ship performance, a feature balance operation (FBO) to improve scale feature description, and a global context block (GCB) to refine features. Experimental results on two public SSDD and HRSID datasets reveal that MAI-SE-Net outperforms the other nine competitive models, better than the suboptimal model by 4.7% detec-tion AP and 3.4% segmentation AP on SSDD and by 3.0% detection AP and 2.4% segmentation AP on HRSID.


title: NOx emissions trends in hydrogen lean premixed flamelets at high strain rate

abstract: NO$_{\rm x}$ formation in lean premixed and highly-strained pure hydrogen-air flamelets is investigated numerically. Lean conditions are established at an equivalence ratio of 0.7. Detailed-chemistry, one-dimensional simulations are performed on a reactants-to-products counter-flow configuration with an applied strain rate ranging from $a=100 \, {\rm s}^{-1}$ to $a=10000 \, {\rm s}^{-1}$ and the \texttt{GRI3.0} mechanism. Following a similar setup, two-dimensional direct numerical simulations are also conducted for representative strain rates of $2000 \, {\rm s}^{-1}$ and $5000 \, {\rm s}^{-1}$. Both solutions show a decreasing NO$_{\rm x}$ trend as the applied strain rate is increased. This decreasing emission outcome is highlighted for the first time in this study for lean pure-hydrogen flamelets. A deep analysis of the 2D solution underlines that there is no production of NO$_{\rm x}$ in the second dimension, thus proving that the emission trend is not a result of a setup preconditioning, but is instead a direct physical effect of stretch on the flame. Furthermore, a detailed analysis of the NO$_{\rm x}$ formation pathways at $a=2000 \, {\rm s}^{-1}$ and $a=5000 \, {\rm s}^{-1}$ is performed. Thermal NO$_{\rm x}$ and NNH pathways are shown to both contribute significantly to the total NO$_{\rm x}$ production. While the NNH route contribution is roughly constant at different strain rates, a significant decrease is observed along the thermal NO$_{\rm x}$ route. Overall, results show that lean and highly-strained hydrogen flames experience a significant decrease of NO$_{\rm x}$. This property is discussed and analysed in the paper.


title: The Atacama Large Aperture Submillimeter Telescope: Key science drivers

abstract: The Atacama Large Aperture Submillimeter Telescope (AtLAST) is a concept for a 50m class single-dish telescope that will provide high sensitivity, fast mapping of the (sub-)millimeter sky. Expected to be powered by renewable energy sources, and to be constructed in the Atacama desert in the 2030s, AtLAST's suite of up to six state-of-the-art instruments will take advantage of its large field of view and high throughput to deliver efficient continuum and spectroscopic observations of the faint, large-scale emission that eludes current facilities. Here we present the key science drivers for the telescope characteristics, and discuss constraints that the transformational science goals place on future instrumentation.


title: Using Quantile Forecasts for Dynamic Equivalents of Active Distribution Grids under Uncertainty

abstract: While distribution networks (DNs) turn from consumers to active and responsive intelligent DNs, the question of how to represent them in large-scale transmission network (TN) studies is still under investigation. The standard approach that uses aggregated models for the inverter-interfaced generation and conventional load models introduces significant errors to the dynamic modeling that can lead to instabilities. This paper presents a new approach based on quantile forecasting to represent the uncertainty originating in DNs at the TN level. First, we aquire a required rich dataset employing Monte Carlo simulations of a DN. Then, we use machine learning (ML) algorithms to not only predict the most probable response but also intervals of potential responses with predefined confidence. These quantile methods represent the variance in DN responses at the TN level. The results indicate excellent performance for most ML techniques. The tuned quantile equivalents predict accurate bands for the current at the TN/DN-interface, and tests with unseen TN conditions indicate robustness. A final assessment that compares the MC trajectories against the predicted intervals highlights the potential of the proposed method.


title: Estimating States and Model Uncertainties Jointly by a Sparsity Promoting UKF

abstract: State estimation when only a partial model of a considered system is available remains a major challenge in many engineering fields. This work proposes a joint, square-root unscented Kalman filter to estimate states and model uncertainties simultaneously by linear combinations of physics-motivated library functions. Using a sparsity promoting approach, a selection of those linear combinations is chosen and thus an interpretable model can be extracted. Results indicate a small estimation error compared to a traditional square-root unscented Kalman filter and exhibit the enhancement of physically meaningful models.


title: RePFormer: Refinement Pyramid Transformer for Robust Facial Landmark Detection

abstract: This paper presents a Refinement Pyramid Transformer (RePFormer) for robust facial landmark detection. Most facial landmark detectors focus on learning representative image features. However, these CNN-based feature representations are not robust enough to handle complex real-world scenarios due to ignoring the internal structure of landmarks, as well as the relations between landmarks and context. In this work, we formulate the facial landmark detection task as refining landmark queries along pyramid memories. Specifically, a pyramid transformer head (PTH) is introduced to build both homologous relations among landmarks and heterologous relations between landmarks and cross-scale contexts. Besides, a dynamic landmark refinement (DLR) module is designed to decompose the landmark regression into an end-to-end refinement procedure, where the dynamically aggregated queries are transformed to residual coordinates predictions. Extensive experimental results on four facial landmark detection benchmarks and their various subsets demonstrate the superior performance and high robustness of our framework.


title: Synthesis and Superconductivity in yttrium superhydrides under high pressure

abstract: The flourishing rare earth superhydrides are a class of recently discovered materials that possess near-room-temperature superconductivity at high pressures, opening a new era of superconductivity research at high pressures. Among these superhydrides, yttrium superhydrides attracted great interest owing to their abundance of stoichiometries and excellent superconductivities. Here, we carried out a comprehensive study of yttrium superhydrides in a wide pressure range of 145-300 GPa. We successfully synthesized a series of superhydrides with the compositions of YH4, YH6, YH7, and YH9, and reported their superconducting transition temperatures of 82 K at 167 GPa, 218 K at 165 GPa, 29 K at 162 GPa, and 230 K at 300 GPa, respectively, which were evidenced by a sharp drop of resistivity. The structure and superconductivity of YH4, which was taken as a representative example, were also examined by X-ray diffraction measurements and the suppression of the superconductivity under external magnetic fields, respectively. Clathrate YH10 as a candidate of room-temperature superconductor was not synthesized within the studied pressure and temperature ranges of up to 300 GPa and 2000 K, respectively. The current work created a detailed platform for further searching room-temperature superconductors in polynary yttrium-based superhydrides.


title: Accurate Kohn-Sham auxiliary system from the ground state density of solids

abstract: The Kohn-Sham (KS) system is an auxiliary system whose effective potential is unknown in most cases. It is in principle determined by the ground state density, and it has been found numerically for some low-dimensional systems by inverting the KS equations starting from a given accurate density. For solids, only approximate results are available. In this work, we determine accurate exchange correlation (xc) potentials for Si and NaCl using the ground state densities obtained from Auxiliary Field Quantum Monte Carlo calculations. We show that these xc potentials can be rationalized as an ensemble of environment-adapted functions of the local density. The KS band structure can be obtained with high accuracy. The true KS band gap turns out to be larger than the prediction of the local density approximation, but significantly smaller than the measurable photoemission gap, which confirms previous estimates. Finally, our findings show that the conjecture that very different xc potentials can lead to very similar densities and other KS observables is true also in solids, which questions the meaning of details of the potentials and, at the same time, confirms the stability of the KS system.


title: Towards Semantic Communication Protocols: A Probabilistic Logic Perspective

abstract: Classical medium access control (MAC) protocols are interpretable, yet their task-agnostic control signaling messages (CMs) are ill-suited for emerging mission-critical applications. By contrast, neural network (NN) based protocol models (NPMs) learn to generate task-specific CMs, but their rationale and impact lack interpretability. To fill this void, in this article we propose, for the first time, a semantic protocol model (SPM) constructed by transforming an NPM into an interpretable symbolic graph written in the probabilistic logic programming language (ProbLog). This transformation is viable by extracting and merging common CMs and their connections while treating the NPM as a CM generator. By extensive simulations, we corroborate that the SPM tightly approximates its original NPM while occupying only 0.02% memory. By leveraging its interpretability and memory-efficiency, we demonstrate several SPM-enabled applications such as SPM reconfiguration for collision-avoidance, as well as comparing different SPMs via semantic entropy calculation and storing multiple SPMs to cope with non-stationary environments.


title: nlfem: A flexible 2d Fem Code for Nonlocal Convection-Diffusion and Mechanics

abstract: In this work we present the mathematical foundation of an assembly code for finite element approximations of nonlocal models with compactly supported, weakly singular kernels. We demonstrate the code on a nonlocal diffusion model in various configurations and on a two-dimensional bond-based peridynamics model. The code nlfem is published under the MIT License and can be freely downloaded.


title: Transport of currents and geometric Rademacher-type theorems

abstract: The transport of many kinds of singular structures in a medium, such as vortex points/lines/sheets in fluids, dislocation loops in crystalline plastic solids, or topological singularities in magnetism, can be expressed in terms of the geometric (Lie) transport equation \[ \frac{\mathrm{d}}{\mathrm{d} t} T_t + \mathcal{L}_{b_t} T_t = 0 \] for a time-indexed family of integral or normal $k$-currents $t \mapsto T_t$ in $\mathbb{R}^d$. Here, $b_t$ is the driving vector field and $\mathcal{L}_{b_t} T_t$ is the Lie derivative of $T_t$ with respect to $b_t$. Written in coordinates for different values of $k$, this PDE encompasses the classical transport equation ($k = d$), the continuity equation ($k = 0$), as well as the equations for the transport of dislocation lines in crystals ($k = 1$) and membranes in liquids ($k =d-1$). The top-dimensional and bottom-dimensional cases have received a great deal of attention in connection with the DiPerna--Lions and Ambrosio theories of Regular Lagrangian Flows. On the other hand, very little is rigorously known at present in the intermediate-dimensional cases. This work develops the theory of the geometric transport equation for arbitrary $k$ and in the case of boundaryless currents $T_t$, covering in particular existence and uniqueness of solutions, structure theorems, rectifiability, and a number of Rademacher-type differentiability results. The latter yield, given an absolutely continuous (in time) path $t \mapsto T_t$, the existence almost everywhere of a ''geometric derivative'', namely a driving vector field $b_t$. This subtle question turns out to be intimately related to the critical set of the evolution, a new notion introduced in this work, which is closely related to Sard's theorem and concerns singularities that are ''smeared out in time''. Our differentiability results are sharp, which we demonstrate through an explicit example.


title: Bernstein-Kouchnirenko-Khovanskii with a symmetry

abstract: A generic polynomial f(x,y,z) with a prescribed Newton polytope defines a symmetric spatial curve f(x,y,z)=f(y,x,z)=0. We study its geometry: the number, degree and genus of its irreducible components, the number and type of singularities, etc. and discuss to what extent these results generalize to higher dimension and more complicated symmetries. As an application, we characterize generic one-parameter families of complex univariate polynomials, whose Galois group is a complete symmetric group.


title: Isospectral graphs via spectral bracketing

abstract: In this article, we develop a perturbative technique to construct families of non-isomorphic discrete graphs which are isospectral for the standard (also called normalised) Laplacian and its signless version. We use vertex contractions as a graph perturbation and spectral bracketing with auxiliary graphs which have certain eigenvalues with high multiplicity. There is no need to know explicitly the eigenvalues or eigenfunctions of the corresponding graphs. We illustrate the method by presenting several families of examples of isospectral graphs including fuzzy complete bipartite graphs and subdivision graphs obtained from the previous examples. All the examples constructed turn out to be also isospectral for the standard (Kirchhoff) Laplacian on the associated equilateral metric graph.


title: Search for resonant $WZ \rightarrow \ell\nu \ell^{\prime}\ell^{\prime}$ production in proton$-$proton collisions at $\mathbf{\sqrt{s} = 13}$ TeV with the ATLAS detector

abstract: A search for a $WZ$ resonance, in the fully leptonic final state (electrons and muons), is performed using 139 fb$^{-1}$ of data collected at a centre-of-mass energy of 13 TeV by the ATLAS detector at the Large Hadron Collider. The results are interpreted in terms of a singly charged Higgs boson of the Georgi$-$Machacek model, produced by $WZ$ fusion, and of a Heavy Vector Triplet, with the resonance produced by $WZ$ fusion or the Drell$-$Yan process. No significant excess over the Standard Model predictions is observed and limits are set on the production cross-section times branching ratio as a function of the resonance mass for these processes.


title: On the Universality of Random Persistence Diagrams

abstract: One of the most elusive challenges within the area of topological data analysis is understanding the distribution of persistence diagrams. Despite much effort, this is still largely an open problem. In this paper, we present a series of novel conjectures regarding the behavior of persistence diagrams arising from random point-clouds. We claim that these diagrams obey a universal probability law, and include an explicit expression as a candidate for what this law is. We back these conjectures with an exhaustive set of experiments, including both simulated and real data. We demonstrate the power of these conjectures by proposing a new hypothesis testing framework for individual features within persistence diagrams.


title: BAST: Binaural Audio Spectrogram Transformer for Binaural Sound Localization

abstract: Accurate sound localization in a reverberation environment is essential for human auditory perception. Recently, Convolutional Neural Networks (CNNs) have been utilized to model the binaural human auditory pathway. However, CNN shows barriers in capturing the global acoustic features. To address this issue, we propose a novel end-to-end Binaural Audio Spectrogram Transformer (BAST) model to predict the sound azimuth in both anechoic and reverberation environments. Two modes of implementation, i.e. BAST-SP and BAST-NSP corresponding to BAST model with shared and non-shared parameters respectively, are explored. Our model with subtraction interaural integration and hybrid loss achieves an angular distance of 1.29 degrees and a Mean Square Error of 1e-3 at all azimuths, significantly surpassing CNN based model. The exploratory analysis of the BAST's performance on the left-right hemifields and anechoic and reverberation environments shows its generalization ability as well as the feasibility of binaural Transformers in sound localization. Furthermore, the analysis of the attention maps is provided to give additional insights on the interpretation of the localization process in a natural reverberant environment.


title: Accelerating Material Design with the Generative Toolkit for Scientific Discovery

abstract: With the growing availability of data within various scientific domains, generative models hold enormous potential to accelerate scientific discovery. They harness powerful representations learned from datasets to speed up the formulation of novel hypotheses with the potential to impact material discovery broadly. We present the Generative Toolkit for Scientific Discovery (GT4SD). This extensible open-source library enables scientists, developers, and researchers to train and use state-of-the-art generative models to accelerate scientific discovery focused on material design.


title: String Diagrams for Layered Explanations

abstract: We propose a categorical framework to reason about scientific explanations: descriptions of a phenomenon meant to translate it into simpler terms, or into a context that has been already understood. Our motivating examples come from systems biology, electrical circuit theory, and concurrency. We demonstrate how three explanatory models in these seemingly diverse areas can be all understood uniformly via a graphical calculus of layered props. Layered props allow for a compact visual presentation of the same phenomenon at different levels of precision, as well as the translation between these levels. Notably, our approach allows for partial explanations, that is, for translating just one part of a diagram while keeping the rest of the diagram untouched. Furthermore, our approach paves the way for formal reasoning about counterfactual models in systems biology.


title: The possibility of $X(4014)$ as a $D^\ast \bar{D}^\ast$ molecular state

abstract: Within the framework of the local hidden gauge approach, we have studied the near-threshold interaction of the $D^* \bar{D}^*$ channel with quantum numbers $I(J^{PC}) = 0(0^{++})$, $0(2^{++})$, $1(0^{++})$, and $1(2^{++})$, respectively. The contact interaction is taken into account alone, since it is expected to give the dominant contribution near the threshold. One pole, which is found in the case of the quantum numbers $I(J^{PC}) = 0(0^{++})$, could be associated to the $X(4014)$ recently observed by the Belle Collaboration. Thus, we suggest that the $X(4014)$ may be an $D^* \bar{D}^*$ molecular state with $I(J^{PC}) = 0(0^{++})$, and more precise information about $X(4014)$, such as resonance parameters and decay properties, could be useful to shed light on its structure.


title: Random Colorings in Manifolds

abstract: We develop a general method for constructing random manifolds and submanifolds in arbitrary dimensions. The method is based on associating colors to the vertices of a triangulated manifold, as in recent work for curves in 3-dimensional space by Sheffield and Yadin (2014). We determine conditions on which submanifolds can arise, in terms of Stiefel-Whitney classes and other properties. We then consider the random submanifolds that arise from randomly coloring the vertices. Since this model generates submanifolds, it allows for studying properties and using tools that are not available in processes that produce general random subcomplexes. The case of 3 colors in a triangulated 3-ball gives rise to random knots and links. In this setting, we answer a question raised by de Crouy-Chanel and Simon (2019), showing that the probability of generating an unknot decays exponentially. In the general case of k colors in d-dimensional manifolds, we investigate the random submanifolds of different codimensions, as the number of vertices in the triangulation grows. We compute the expected Euler characteristic, and discuss relations to homological percolation and other topological properties. Finally, we explore a method to search for solutions to topological problems by generating random submanifolds. We describe computer experiments that search for a low-genus surface in the 4-dimensional ball whose boundary is a given knot in the 3-dimensional sphere.


title: Memory-free Online Change-point Detection: A Novel Neural Network Approach

abstract: Change-point detection (CPD), which detects abrupt changes in the data distribution, is recognized as one of the most significant tasks in time series analysis. Despite the extensive literature on offline CPD, unsupervised online CPD still suffers from major challenges, including scalability, hyperparameter tuning, and learning constraints. To mitigate some of these challenges, in this paper, we propose a novel deep learning approach for unsupervised online CPD from multi-dimensional time series, named Adaptive LSTM-Autoencoder Change-Point Detection (ALACPD). ALACPD exploits an LSTM-autoencoder-based neural network to perform unsupervised online CPD. It continuously adapts to the incoming samples without keeping the previously received input, thus being memory-free. We perform an extensive evaluation on several real-world time series CPD benchmarks. We show that ALACPD, on average, ranks first among state-of-the-art CPD algorithms in terms of quality of the time series segmentation, and it is on par with the best performer in terms of the accuracy of the estimated change-points. The implementation of ALACPD is available online on Github\footnote{\url{https://github.com/zahraatashgahi/ALACPD}}.


title: A law of adversarial risk, interpolation, and label noise

abstract: In supervised learning, it has been shown that label noise in the data can be interpolated without penalties on test accuracy. We show that interpolating label noise induces adversarial vulnerability, and prove the first theorem showing the relationship between label noise and adversarial risk for any data distribution. Our results are almost tight if we do not make any assumptions on the inductive bias of the learning algorithm. We then investigate how different components of this problem affect this result, including properties of the distribution. We also discuss non-uniform label noise distributions; and prove a new theorem showing uniform label noise induces nearly as large an adversarial risk as the worst poisoning with the same noise rate. Then, we provide theoretical and empirical evidence that uniform label noise is more harmful than typical real-world label noise. Finally, we show how inductive biases amplify the effect of label noise and argue the need for future work in this direction.


title: Active Learning-based Isolation Forest (ALIF): Enhancing Anomaly Detection in Decision Support Systems

abstract: The detection of anomalous behaviours is an emerging need in many applications, particularly in contexts where security and reliability are critical aspects. While the definition of anomaly strictly depends on the domain framework, it is often impractical or too time consuming to obtain a fully labelled dataset. The use of unsupervised models to overcome the lack of labels often fails to catch domain specific anomalies as they rely on general definitions of outlier. This paper suggests a new active learning based approach, ALIF, to solve this problem by reducing the number of required labels and tuning the detector towards the definition of anomaly provided by the user. The proposed approach is particularly appealing in the presence of a Decision Support System (DSS), a case that is increasingly popular in real-world scenarios. While it is common that DSS embedded with anomaly detection capabilities rely on unsupervised models, they don't have a way to improve their performance: ALIF is able to enhance the capabilities of DSS by exploiting the user feedback during common operations. ALIF is a lightweight modification of the popular Isolation Forest that proved superior performances with respect to other state-of-art algorithms in a multitude of real anomaly detection datasets.


title: ControlBurn: Nonlinear Feature Selection with Sparse Tree Ensembles

abstract: ControlBurn is a Python package to construct feature-sparse tree ensembles that support nonlinear feature selection and interpretable machine learning. The algorithms in this package first build large tree ensembles that prioritize basis functions with few features and then select a feature-sparse subset of these basis functions using a weighted lasso optimization criterion. The package includes visualizations to analyze the features selected by the ensemble and their impact on predictions. Hence ControlBurn offers the accuracy and flexibility of tree-ensemble models and the interpretability of sparse generalized additive models. ControlBurn is scalable and flexible: for example, it can use warm-start continuation to compute the regularization path (prediction error for any number of selected features) for a dataset with tens of thousands of samples and hundreds of features in seconds. For larger datasets, the runtime scales linearly in the number of samples and features (up to a log factor), and the package support acceleration using sketching. Moreover, the ControlBurn framework accommodates feature costs, feature groupings, and $\ell_0$-based regularizers. The package is user-friendly and open-source: its documentation and source code appear on https://pypi.org/project/ControlBurn/ and https://github.com/udellgroup/controlburn/.


title: Search by triplet: An efficient local track reconstruction algorithm for parallel architectures

abstract: Millions of particles are collided every second at the LHCb detector placed inside the Large Hadron Collider at CERN. The particles produced as a result of these collisions pass through various detecting devices which will produce a combined raw data rate of up to 40 Tbps by 2021. These data will be fed through a data acquisition system which reconstructs individual particles and filters the collision events in real time. This process will occur in a heterogeneous farm employing exclusively off-the-shelf CPU and GPU hardware, in a two stage process known as High Level Trigger. The reconstruction of charged particle trajectories in physics detectors, also referred to as track reconstruction or tracking, determines the position, charge and momentum of particles as they pass through detectors. The Vertex Locator subdetector (VELO) is the closest such detector to the beamline, placed outside of the region where the LHCb magnet produces a sizable magnetic field. It is used to reconstruct straight particle trajectories which serve as seeds for reconstruction of other subdetectors and to locate collision vertices. The VELO subdetector will detect up to 1000 million particles every second, which need to be reconstructed in real time in the High Level Trigger. We present Search by triplet, an efficient track reconstruction algorithm. Our algorithm is designed to run efficiently across parallel architectures. We extend on previous work and explain the algorithm evolution since its inception. We show the scaling of our algorithm under various situations, and analyze its amortized time in terms of complexity for each of its constituent parts and profile its performance. Our algorithm is the current state-of-the-art in VELO track reconstruction on SIMT architectures, and we qualify its improvements over previous results.


title: Neutrino oscillations in gravitational fields and astrophysical applications

abstract: The neutrino propagation and oscillations in various gravitational fields are studied. First, we consider the neutrino scattering off a black hole accounting for the neutrino spin precession. Then, we study the evolution of flavor neutrinos in stochastic gravitational waves. The astrophysical applications of the obtained results are considered.


title: An Approach to Ensure Fairness in News Articles

abstract: Recommender systems, information retrieval, and other information access systems present unique challenges for examining and applying concepts of fairness and bias mitigation in unstructured text. This paper introduces Dbias, which is a Python package to ensure fairness in news articles. Dbias is a trained Machine Learning (ML) pipeline that can take a text (e.g., a paragraph or news story) and detects if the text is biased or not. Then, it detects the biased words in the text, masks them, and recommends a set of sentences with new words that are bias-free or at least less biased. We incorporate the elements of data science best practices to ensure that this pipeline is reproducible and usable. We show in experiments that this pipeline can be effective for mitigating biases and outperforms the common neural network architectures in ensuring fairness in the news articles.


title: Sensemaking and Scientific Modeling: Intertwined processes analyzed in the context of physics problem solving

abstract: Researchers in physics education have advocated both for including modeling in science classrooms as well as promoting student engagement with sensemaking. These two processes facilitate the generation of new knowledge by connecting to one's existing ideas. Despite being two distinct processes, modeling is often described as sensemaking of the physical world. In the current work, we provide an explicit, framework-based analysis of the intertwining between modeling and sensemaking by analyzing think-aloud interviews of two students solving a physics problem. While one student completes the task, the other abandons their approach. The case studies reveal that particular aspects of modeling and sensemaking processes co-occur. For instance, the priming on the `given' information from the problem statement constituted the students' engagement with their mental models, and their attempts to resolve inconsistencies in understanding involved the use of external representations. We find that barriers experienced in modeling can inhibit students' sustained sensemaking. These results suggest ways for future research to support students' sensemaking in physics by promoting modeling practices.


title: Bistochastic privacy

abstract: We introduce a new privacy model relying on bistochastic matrices, that is, matrices whose components are nonnegative and sum to 1 both row-wise and column-wise. This class of matrices is used to both define privacy guarantees and a tool to apply protection on a data set. The bistochasticity assumption happens to connect several fields of the privacy literature, including the two most popular models, k-anonymity and differential privacy. Moreover, it establishes a bridge with information theory, which simplifies the thorny issue of evaluating the utility of a protected data set. Bistochastic privacy also clarifies the trade-off between protection and utility by using bits, which can be viewed as a natural currency to comprehend and operationalize this trade-off, in the same way than bits are used in information theory to capture uncertainty. A discussion on the suitable parameterization of bistochastic matrices to achieve the privacy guarantees of this new model is also provided.


title: Local Irregularity Conjecture vs. cacti

abstract: A graph is locally irregular if the degrees of the end-vertices of every edge are distinct. An edge coloring of a graph G is locally irregular if every color induces a locally irregular subgraph of G. A colorable graph G is any graph which admits a locally irregular edge coloring. The locally irregular chromatic index X'irr(G) of a colorable graph G is the smallest number of colors required by a locally irregular edge coloring of G. The Local Irregularity Conjecture claims that all colorable graphs require at most 3 colors for locally irregular edge coloring. Recently, it has been observed that the conjecture does not hold for the bow-tie graph B [7]. Cacti are important class of graphs for this conjecture since B and all non-colorable graphs are cacti. In this paper we show that for every colorable cactus graph G != B it holds that X'irr(G) <= 3. This makes us to believe that B is the only colorable graph with X'irr(B) > 3, and consequently that B is the only counterexample to the Local Irregularity Conjecture.


title: Microscopic origin of anomalous interlayer exciton transport in van der Waals heterostructures

abstract: Van der Waals heterostructures constitute a platform for investigating intriguing many-body quantum phenomena. In particular, transition-metal dichalcogenide (TMD) hetero-bilayers host long-lived interlayer excitons which exhibit permanent out-of-plane dipole moments. Here, we develop a microscopic theory for interlayer exciton-exciton interactions including both the dipolar nature of interlayer excitons as well as their fermionic substructure, which gives rise to an attractive fermionic exchange. We find that these interactions contribute to a drift force resulting in highly non-linear exciton propagation at elevated densities in the MoSe$_2$-WSe$_2$ heterostructure. We show that the propagation can be tuned by changing the number of hBN spacers between the TMD layers or by adjusting the dielectric environment. In particular, although counter-intuitive, we reveal that interlayer excitons in free-standing samples propagate slower than excitons in hBN-encapsulated TMDs - due to an enhancement of the net Coulomb-drift with stronger environmental screening. Overall, our work contributes to a better microscopic understanding of the interlayer exciton transport in technologically promising atomically thin semiconductors.


title: A Geometric Condition for Uniqueness of Fr\'echet Means of Persistence Diagrams

abstract: The Fr\'echet mean is an important statistical summary and measure of centrality of data; it has been defined and studied for persistent homology captured by persistence diagrams. However, the complicated geometry of the space of persistence diagrams implies that the Fr\'echet mean for a given set of persistence diagrams is not necessarily unique, which prohibits theoretical guarantees for empirical means with respect to population means. In this paper, we derive a variance expression for a set of persistence diagrams exhibiting a multi-matching between the persistence points known as a grouping. Moreover, we propose a condition for groupings, which we refer to as flatness: sets of persistence diagrams that exhibit flat groupings give rise to unique Fr\'echet means. We derive a finite sample convergence result for general groupings, which results in convergence for Fr\'echet means if the groupings are flat. Finally, we interpret flat groupings in a recently-proposed general framework of Fr\'echet means in Alexandrov geometry. Together with recent results from Alexandrov geometry, this allows for the first derivation of a finite sample convergence rate for sets of persistence diagrams and lays the ground for viability of the Fr\'echet mean as a practical statistical summary of persistent homology.


title: Shrouded black holes in Einstein-Gauss-Bonnet gravity

abstract: We study black holes in a modified gravity scenario involving a scalar field quadratically coupled to the Gauss-Bonnet invariant. The scalar is assumed to be in a spontaneously broken phase at spatial infinity due to a bare Higgs-like potential. For a proper choice of sign, the non-minimal coupling to gravity leads to symmetry restoration near the black hole horizon, prompting the development of the scalar wall in its vicinity. The wall thickness depends on the bare mass of the scalar and can be much smaller than the Schwarzschild radius. In a weakly coupled regime, the quadratic coupling to the Gauss-Bonnet invariant effectively becomes linear, and no walls are formed. We find approximate analytical solutions for the scalar field in the test field regime, and obtain numerically static black hole solutions within this setup. We discuss cosmological implications of the model and show that it is fully consistent with the existence of an inflationary stage, unlike the spontaneous scalarization scenario assuming the opposite sign of the non-minimal coupling to gravity. Our model predicts the speed of gravitational waves to be extremely close to unity, - in a comfortable agreement with the observation of the GW170817 event and its electromagnetic counterpart.


title: High Performance Simulation for Scalable Multi-Agent Reinforcement Learning

abstract: Multi-agent reinforcement learning experiments and open-source training environments are typically limited in scale, supporting tens or sometimes up to hundreds of interacting agents. In this paper we demonstrate the use of Vogue, a high performance agent based model (ABM) framework. Vogue serves as a multi-agent training environment, supporting thousands to tens of thousands of interacting agents while maintaining high training throughput by running both the environment and reinforcement learning (RL) agents on the GPU. High performance multi-agent environments at this scale have the potential to enable the learning of robust and flexible policies for use in ABMs and simulations of complex systems. We demonstrate training performance with two newly developed, large scale multi-agent training environments. Moreover, we show that these environments can train shared RL policies on time-scales of minutes and hours.


title: Complementarity relations of a delayed-choice quantum eraser in a quantum circuit

abstract: We propose a quantum circuit that emulates a delayed-choice quantum eraser via bipartite entanglement with the extension that the degree of entanglement between the two paired quantons is adjustable. This provides a broader setting to test complementarity relations between interference visibility and which-way distinguishability in the scenario that the which-way information is obtained through entanglement without direct contact with the quantum state for interference. The visibility-distinguishability relations are investigated from three perspectives that differ in how the which-way information is taken into consideration. These complementarity relations can be understood in terms of entropic uncertainty relations in the information-theoretic framework and the triality relation that incorporates single-particle and bipartite properties. We then perform experiments on the quantum computers provided by the IBM Quantum platform to verify the theoretical predictions. We also apply the delay gate to delay the measurement of the which-way information to affirm that the measurement can be made truly in the "delayed-choice" manner.


title: Classical Limit of Higher-Spin String Amplitudes

abstract: It has been shown that a special set of three-point amplitudes between two massive spinning states and a graviton reproduces the linearised stress-energy tensor for a Kerr black hole in the classical limit. In this work we revisit this result and compare it to the analysis of the amplitudes describing the interaction of leading Regge states of the open and closed superstring. We find an all-spin result for the classical limit of two massive spinning states interacting with a photon or graviton. This result differs from Kerr and instead matches the current four-vector and the stress-energy tensor generated by a classical string coupled to electromagnetism and gravity respectively. For the superstring amplitudes, contrary to the black-hole case, we find that the spin to infinity limit is necessary to reproduce the classical spin multipoles.


title: Self-Adjusting Linear Networks with Ladder Demand Graph

abstract: Self-adjusting networks (SANs) have the ability to adapt to communication demand by dynamically adjusting the workload (or demand) embedding, i.e., the mapping of communication requests into the network topology. SANs can thus reduce routing costs for frequently communicating node pairs by paying a cost for adjusting the embedding. This is particularly beneficial when the demand has structure, which the network can adapt to. Demand can be represented in the form of a demand graph, which is defined by the set of network nodes (vertices) and the set of pairwise communication requests (edges). Thus, adapting to the demand can be interpreted by embedding the demand graph to the network topology. This can be challenging both when the demand graph is known in advance (offline) and when it revealed edge-by-edge (online). The difficulty also depends on whether we aim at constructing a static topology or a dynamic (self-adjusting) one that improves the embedding as more parts of the demand graph are revealed. Yet very little is known about these self-adjusting embeddings. In this paper, the network topology is restricted to a line and the demand graph to a ladder graph, i.e., a $2^n$ grid, including all possible subgraphs of the ladder. We present an online self-adjusting network that matches the known lower bound asymptotically and is $12$-competitive in terms of request cost. As a warm up result, we present an asymptotically optimal algorithm for the cycle demand graph. We also present an oracle-based algorithm for an arbitrary demand graph that has a constant overhead.


title: Computing Electronic Correlation Energies using Linear Depth Quantum Circuits

abstract: Efficient computation of molecular energies is an exciting application of quantum computing for quantum chemistry, but current noisy intermediate-scale quantum (NISQ) devices can only execute shallow circuits, limiting existing variational quantum algorithms, which require deep entangling quantum circuit ansatzes to capture correlations, to small molecules. Here we demonstrate a variational NISQ-friendly algorithm that generates a set of mean-field Hartree-Fock (HF) ansatzes using multiple shallow circuits with depth linear in the number of qubits to estimate electronic correlation energies via perturbation theory up to the second order. We tested the algorithm on several small molecules, both with classical simulations including noise models and on cloud quantum processors, showing that it not only reproduces the equilibrium molecular energies but it also captures the perturbative electronic correlation effects at longer bond distances. As fidelities of quantum processors continue to improve our algorithm will enable the study of larger molecules compared to other approaches requiring higher-order polynomial circuit depth.


title: Fast and Furious dynamo action in the anisotropic dynamo

abstract: In the limit of large magnetic Reynolds numbers, it is shown that a smooth differential rotation can lead to fast dynamo action, provided that the electrical conductivity or magnetic permeability is anisotropic. If the shear is infinite, for example between two rotating solid bodies, the anisotropic dynamo becomes furious, meaning that the magnetic growth rate increases toward infinity with an increasing magnetic Reynolds number.


title: Anomalous transport in periodic photonic chains with designed loss

abstract: Here we show that a coherent random walk in a perfectly periodic chain of bosonic modes with designed loss can exhibit a variety of different anomalous transfer regimes in dependence on the initial state of the chain. In particular, for any given finite initial time-interval there is a set of initial states leading to a hyperballistic transport regime. Also, there are initial states allowing one to achieve a subdiffusive regime or even localization for a given time-interval, or change an asymptotic long-time diffusion rate. We show how these anomalous transport regimes can be practically realized in a laser-written network of single-mode waveguides in balk glass or how a planar system of coupled single-mode waveguides can be realized with an integrated photonic platform.


title: Learning with Muscles: Benefits for Data-Efficiency and Robustness in Anthropomorphic Tasks

abstract: Humans are able to outperform robots in terms of robustness, versatility, and learning of new tasks in a wide variety of movements. We hypothesize that highly nonlinear muscle dynamics play a large role in providing inherent stability, which is favorable to learning. While recent advances have been made in applying modern learning techniques to muscle-actuated systems both in simulation as well as in robotics, so far, no detailed analysis has been performed to show the benefits of muscles when learning from scratch. Our study closes this gap and showcases the potential of muscle actuators for core robotics challenges in terms of data-efficiency, hyperparameter sensitivity, and robustness.


title: Nonlinear three-state quantum walks

abstract: The dynamics of a three-state quantum walk with amplitude-dependent phase shifts is investigated. We consider two representative inputs whose linear evolution is known to display either full dispersion of the wave packet or intrinsic localization on the initial position. The nonlinear counterpart presents much more involved dynamics featuring self-trapping, solitonic pulses, radiation, and chaotic-like behavior. We show that nonlinearity leads to a metastable self-trapped wavepacket component that radiates in the long-time regime with the survival probability $\propto t^{-1/2}$. A sudden dynamical transition from such metastable state to the point when the radiation process is triggered is found for a set of parameters.


title: Black and Gray Box Learning of Amplitude Equations: Application to Phase Field Systems

abstract: We present a data-driven approach to learning surrogate models for amplitude equations, and illustrate its application to interfacial dynamics of phase field systems. In particular, we demonstrate learning effective partial differential equations describing the evolution of phase field interfaces from full phase field data. We illustrate this on a model phase field system, where analytical approximate equations for the dynamics of the phase field interface (a higher order eikonal equation and its approximation, the Kardar-Parisi-Zhang (KPZ) equation) are known. For this system, we discuss data-driven approaches for the identification of equations that accurately describe the front interface dynamics. When the analytical approximate models mentioned above become inaccurate, as we move beyond the region of validity of the underlying assumptions, the data-driven equations outperform them. In these regimes, going beyond black-box identification, we explore different approaches to learn data-driven corrections to the analytically approximate models, leading to effective gray box partial differential equations.

categories: stat.ML cs.LG math.DS nlin.AO physics.comp-ph


title: Outliers, Dynamics, and the Independence Postulate

abstract: We show that outliers occur almost surely in computable dynamics over infinite sequences. Ever greater outliers can be found as the number of visited states increases. We show the Independence Postulate explains how outliers are found in the physical world. We generalize the outliers theorem to uncomputable sampling methods.


title: BPS Correlators for $\text{AdS}_3/\text{CFT}_2$

abstract: The BPS correlators of the symmetric product orbifold $\text{Sym}_N(\mathbb{T}^4)$ are reproduced from the dual worldsheet theory describing strings on $\text{AdS}_3\times {\rm S}^3\times \mathbb{T}^4$ with minimal ($k=1$) NS-NS flux. More specifically, we show that the worldsheet duals of the symmetric orbifold BPS states can be identified with their lift to the covering surface, thereby making the matching of the correlators essentially manifest. We also argue that the argument can be generalised to arbitrary descendants, using suitable DDF operators on the worldsheet.


title: Communication Acceleration of Local Gradient Methods via an Accelerated Primal-Dual Algorithm with Inexact Prox

abstract: Inspired by a recent breakthrough of Mishchenko et al (2022), who for the first time showed that local gradient steps can lead to provable communication acceleration, we propose an alternative algorithm which obtains the same communication acceleration as their method (ProxSkip). Our approach is very different, however: it is based on the celebrated method of Chambolle and Pock (2011), with several nontrivial modifications: i) we allow for an inexact computation of the prox operator of a certain smooth strongly convex function via a suitable gradient-based method (e.g., GD, Fast GD or FSFOM), ii) we perform a careful modification of the dual update step in order to retain linear convergence. Our general results offer the new state-of-the-art rates for the class of strongly convex-concave saddle-point problems with bilinear coupling characterized by the absence of smoothness in the dual function. When applied to federated learning, we obtain a theoretically better alternative to ProxSkip: our method requires fewer local steps ($O(\kappa^{1/3})$ or $O(\kappa^{1/4})$, compared to $O(\kappa^{1/2})$ of ProxSkip), and performs a deterministic number of local steps instead. Like ProxSkip, our method can be applied to optimization over a connected network, and we obtain theoretical improvements here as well.


title: Ab-Initio Simulation of Field Evaporation

abstract: A new simulation approach of field evaporation is presented. The model combines classical electrostatics with molecular dynamics (MD) simulations. Unlike previous atomic-level simulation approaches, our method does not rely on an evaporation criterion based on thermal activation theory, instead, electric-field-induced forces on atoms are explicitly calculated and added to the interatomic forces. Atoms then simply move according to the laws of classical molecular dynamics and are "evaporated" when the external force overcomes interatomic bonding. This approach thus makes no ad-hoc assumptions concerning evaporation fields and criteria, which makes the simulation fully physics-based and "ab-initio" apart from the interatomic potential. As proof of principle, we perform simulations to determine material dependent critical voltages which allow assessing the evaporation fields and the corresponding steady-state tip shapes in different metals. We also extract critical evaporation fields in elemental metals and sublimation energies in a high entropy alloy to have a more direct comparison with tabulated values. In contrast to previous approaches, we show that our method is able to successfully reproduce the enhanced zone lines observed in experimental field desorption patterns. We also demonstrate the need for careful selection of the interatomic potential by a comparative study for the example of Cu-Ni alloys.


title: Reactive Neural Path Planning with Dynamic Obstacle Avoidance in a Condensed Configuration Space

abstract: We present a biologically inspired approach for path planning with dynamic obstacle avoidance. Path planning is performed in a condensed configuration space of a robot generated by self-organizing neural networks (SONN). The robot itself and static as well as dynamic obstacles are mapped from the Cartesian task space into the configuration space by precomputed kinematics. The condensed space represents a cognitive map of the environment, which is inspired by place cells and the concept of cognitive maps in mammalian brains. The generation of training data as well as the evaluation were performed on a real industrial robot accompanied by simulations. To evaluate the reactive collision-free online planning within a changing environment, a demonstrator was realized. Then, a comparative study regarding sample-based planners was carried out. So we could show that the robot is able to operate in dynamically changing environments and re-plan its motion trajectories within impressing 0.02 seconds, which proofs the real-time capability of our concept.


title: Detection of Furigana Text in Images

abstract: Furigana are pronunciation notes used in Japanese writing. Being able to detect these can help improve optical character recognition (OCR) performance or make more accurate digital copies of Japanese written media by correctly displaying furigana. This project focuses on detecting furigana in Japanese books and comics. While there has been research into the detection of Japanese text in general, there are currently no proposed methods for detecting furigana. We construct a new dataset containing Japanese written media and annotations of furigana. We propose an evaluation metric for such data which is similar to the evaluation protocols used in object detection except that it allows groups of objects to be labeled by one annotation. We propose a method for detection of furigana that is based on mathematical morphology and connected component analysis. We evaluate the detections of the dataset and compare different methods for text extraction. We also evaluate different types of images such as books and comics individually and discuss the challenges of each type of image. The proposed method reaches an F1-score of 76\% on the dataset. The method performs well on regular books, but less so on comics, and books of irregular format. Finally, we show that the proposed method can improve the performance of OCR by 5\% on the manga109 dataset. Source code is available via \texttt{\url{https://github.com/nikolajkb/FuriganaDetection}}


title: CoSIm: Commonsense Reasoning for Counterfactual Scene Imagination

abstract: As humans, we can modify our assumptions about a scene by imagining alternative objects or concepts in our minds. For example, we can easily anticipate the implications of the sun being overcast by rain clouds (e.g., the street will get wet) and accordingly prepare for that. In this paper, we introduce a new task/dataset called Commonsense Reasoning for Counterfactual Scene Imagination (CoSIm) which is designed to evaluate the ability of AI systems to reason about scene change imagination. In this task/dataset, models are given an image and an initial question-response pair about the image. Next, a counterfactual imagined scene change (in textual form) is applied, and the model has to predict the new response to the initial question based on this scene change. We collect 3.5K high-quality and challenging data instances, with each instance consisting of an image, a commonsense question with a response, a description of a counterfactual change, a new response to the question, and three distractor responses. Our dataset contains various complex scene change types (such as object addition/removal/state change, event description, environment change, etc.) that require models to imagine many different scenarios and reason about the changed scenes. We present a baseline model based on a vision-language Transformer (i.e., LXMERT) and ablation studies. Through human evaluation, we demonstrate a large human-model performance gap, suggesting room for promising future work on this challenging counterfactual, scene imagination task. Our code and dataset are publicly available at: https://github.com/hyounghk/CoSIm


title: MeerKAT radio observations of the neutron star low-mass X-ray binary Cen X-4 at low accretion rates

abstract: Centaurus X-4 (Cen X-4) is a relatively nearby neutron star low-mass X-ray binary that showed outbursts in 1969 and 1979, but has not shown a full outburst since. Due to its proximity and sustained period of quiescence, it is a prime target to study the coupling between accretion and jet ejection in quiescent neutron star low-mass X-ray binaries. Here, we present four MeerKAT radio observations at 1.3 GHz of Cen X-4, combined with NICER and Swift X-ray monitoring. During the first and most sensitive observation, Cen X-4 was in a fully quiescent X-ray state. The three later and shorter observations targeted a brief period of faint X-ray activity in January 2021, which has been referred to as a 'mis-fired' outburst. Cen X-4 is not detected in any of the four MeerKAT observations. We place these radio non-detections on the X-ray -- radio luminosity diagram, improving the constraints on the correlation between the two luminosities from earlier quiescent radio studies. We confirm that Cen X-4 is radio fainter than the transitional milli-second pulsar PSR J1023+0038 at the same X-ray luminosity. We discuss the radio behaviour of accreting neutron stars at low X-ray luminosity more generally and finally comment on future observing campaigns.


title: The diagonal behaviour of the one-particle Coulombic density matrix

abstract: We obtain bounds for all derivatives of the non-relativistic Coulombic one-particle density matrix $\gamma(x, y)$ near the diagonal $x = y$.


title: Interacting topological Dirac magnons

abstract: In this work, we study the magnon-magnon interaction effect in typical honeycomb ferromagnets consisting of van der Waals-bonded stacks of honeycomb layers, e.g., chromium trihalides CrX3 (X = F, Cl, Br, and I), that display two spin-wave modes (Dirac magnon). Using Green's function formalism with the presence of the Dzyaloshinskii-Moriya interaction, we obtain a spinor Dyson equation up to the second-order approximation by the cluster expansion method. Numerical calculations show prominent renormalizations of the single-particle spectrum. Furthermore, we propose a tunable renormalization effect using a parametric magnon amplification scheme. By amplifying the magnon population at different k points, the enabled renormalization effect not only reshapes the band structure but also modifies the Berry curvature distribution. Our work demonstrates the interplay between band geometry, interactions, and the external light field in the bosonic system and can potentially lead to new insights into the properties of magnon-based spintronic devices.


title: Power Response and Modelling Aspects of Power Electronic Loads in Case of Voltage Drops

abstract: In this paper, the power response of power electronic loads in case of voltage drops are measured and their dynamics are analysed. Based on this, dynamic simulation models are derived which can be used for voltage stability investigations. For this, four loads with different power factor techniques are considered. In addition, the impact of the grid impedance and input filter on the power response are measured in the laboratory. Based on the measurements, the simulation models are described. It is also outlined under which aspects the components of the loads are normally dimensioned if no detailed information is available. A comparison with the measurements demonstrates that the simulation models capture the main dynamics. At the end of the paper, the load models are compared to a constant power load in a short-term voltage stability use case. The results indicate that the power electronic loads have a more positive influence on short-term voltage stability in case of voltage drops. Overall, the contributions of the paper are the identification of the basic power dynamics of power electronic loads for different voltage drops and a subsequent derivation of suitable simulation load models for voltage stability investigations.


title: Thin-shell wormholes in N-dimensional F(R) gravity

abstract: In this work, spherically symmetric thin-shell wormholes with a conformally invariant Maxwell field for $N$-dimensional $F(R)$ gravity and constant scalar curvature $R$ are built. Two cases are considered: symmetric wormholes and asymmetric ones in the scalar curvature. Their stability under radial perturbations is analyzed, finding stable solutions made of exotic matter for a given range of the parameters.


title: Geometric blow-up of a dynamic Turing instability in the Swift-Hohenberg equation

abstract: We present a rigorous analysis of the slow passage through a Turing bifurcation in the Swift-Hohenberg equation using a novel approach based on geometric blow-up. We show that the formally derived multiple scales ansatz which is known from classical modulation theory can be adapted for use in the fast-slow setting, by reformulating it as a blow-up transformation. This leads to dynamically simpler modulation equations posed in the blown-up space, via a formal procedure which directly extends the established approach to the time-dependent setting. The modulation equations take the form of non-autonomous Ginzburg-Landau equations, which can be analysed within the blow-up. The asymptotics of solutions in weighted Sobelev spaces are given in two different cases: (i) A symmetric case featuring a delayed loss of stability, and (ii) A second case in which the symmetry is broken by a source term. In order to characterise the dynamics of the Swift-Hohenberg equation itself we derive rigorous estimates on the error of the dynamic modulation approximation. These estimates are obtained by bounding weak solutions to an evolution equation for the error which is also posed in the blown-up space. Using the error estimates obtained, we are able to infer the asymptotics of a large class of solutions to the dynamic Swift-Hohenberg equation. We provide rigorous asymptotics for solutions in both cases (i) and (ii). We also prove the existence of the delayed loss of stability in the symmetric case (i), and provide a lower bound for the delay time.


title: Do minor interactions trigger star formation in galaxy pairs?

abstract: We analyze the galaxy pairs in a set of volume limited samples from the SDSS to study the effects of minor interactions on the star formation rate (SFR) and colour of galaxies. We carefully design control samples of the isolated galaxies by matching the stellar mass and redshift of the minor pairs. The SFR distributions and colour distributions in the minor pairs differ from their controls at $>99\%$ significance level. We also simultaneously match the control galaxies in stellar mass, redshift and local density to assess the role of the environment. The null hypothesis can be rejected at $>99\%$ confidence level even after matching the environment. Our analysis shows a quenching in the minor pairs where the degree of quenching decreases with the increasing pair separation and plateaus beyond 50 kpc. We also prepare a sample of minor pairs with $H_{\alpha}$ line information. We calculate the SFR of these galaxies using the $H_{\alpha}$ line and repeat our analysis. We observe a quenching in the $H_{\alpha}$ sample too. We find that the majority of the minor pairs are quiescent systems that could be quenched due to minor interactions. Combining data from the Galaxy Zoo and Galaxy Zoo2, we find that only $\sim 1\%$ galaxies have a dominant bulge, $4\%-7\%$ galaxies host a bar, and $5\%-10\%$ galaxies show the AGN activity in minor pairs. This indicates that the presence of bulge, bar and AGN activity plays an insignificant role in quenching the galaxies in minor pairs. The more massive companion satisfies the criteria for mass quenching in most of the minor pairs. We propose that the stripping and starvation likely caused the quenching in the less massive companion at a later stage of evolution.


title: A club guessing toolbox I

abstract: Club guessing principles were introduced by Shelah as a weakening of Jensen's diamond. Most spectacularly, they were used to prove Shelah's ZFC bound on the power of the first singular cardinal. These principles have found many other applications: in cardinal arithmetic and PCF theory; in the construction of combinatorial objects on uncountable cardinals such as Jonsson algebras, strong colourings, Souslin trees, and pathological graphs; to the non-existence of universals in model theory; to the non-existence of forcing axioms at higher uncountable cardinals; and many more. In this paper, the first part of a series, we survey various forms of club-guessing that have appeared in the literature, and then systematically study the various ways in which a club-guessing sequences can be improved, especially in the way the frequency of guessing is calibrated. We include an expository section intended for those unfamiliar with club-guessing and which can be read independently of the rest of the article.


title: Boundary and domain wall theories of 2d generalized quantum double model

abstract: The generalized quantum double lattice realization of 2d topological orders based on Hopf algebras is discussed in this work. Both left-module and right-module constructions are investigated. The ribbon operators and the classification of topological excitations based on the representations of the quantum double of Hopf algebras are discussed. To generalize the model to a 2d surface with boundaries and surface defects, we present a systematic construction of the boundary Hamiltonian and domain wall Hamiltonian. The algebraic data behind the gapped boundary and domain wall are comodule algebras and bicomodule algebras. The topological excitations in the boundary and domain wall are classified by bimodules over these algebras. The ribbon operator realization of boundary-bulk duality is also discussed. Finally, via the Hopf tensor network representation of the quantum many-body states, we solve the ground state of the model in the presence of the boundary and domain wall.

categories: quant-ph cond-mat.str-el math-ph math.MP


title: Fast high-fidelity gates for galvanically-coupled fluxonium qubits using strong flux modulation

abstract: Long coherence times, large anharmonicity and robust charge-noise insensitivity render fluxonium qubits an interesting alternative to transmons. Recent experiments have demonstrated record coherence times for low-frequency fluxonia. Here, we propose a galvanic-coupling scheme with flux-tunable $\textit{XX}$ coupling. To implement a high-fidelity entangling $\sqrt{i\mathrm{SWAP}}$ gate, we modulate the strength of this coupling and devise variable-time identity gates to synchronize required single-qubit operations. Both types of gates are implemented using strong ac flux drives, lasting for only a few drive periods. We employ a theoretical framework capable of capturing qubit dynamics beyond the rotating-wave approximation (RWA) as required for such strong drives. We predict an open-system fidelity of $F>0.999$ for the $\sqrt{i\mathrm{SWAP}}$ gate under realistic conditions.


title: Weakly bounded cohomology classes and a counterexample to a conjecture of Gromov

abstract: We exhibit a finitely presented group whose second cohomology contains a weakly bounded, but not bounded, class. As an application, we disprove a long-standing conjecture of Gromov about bounded primitives of differential forms on universal covers of closed manifolds.


title: A new compact active turbulence generator for premixed combustion: Non-reacting flow characteristics

abstract: A new compact active turbulence generator is developed, tested, and characterized, which extends the capabilities of such generators used in turbulent premixed combustion research. The generator is composed of two blades that resemble the shape of two bow-ties. Hot-wire anemometry and high-speed imaging are performed to characterize the non-reacting flow produced by the generator and the blades dynamics, respectively. Two mean bulk flow velocities of 5.0 and 7.0 m/s are examined. For comparison purposes, in addition to the developed generator, tests are also performed for a free jet as well as one and two perforated plates. The results show that the centerline root mean square velocity fluctuations can become about 1.8~m/s, which is several folds larger than that in the literature for reacting flows generated by active grids. For the newly developed device, the power-law decay of the one-dimensional kinetic energy is -1.0 and -1.3 for the mean bulk flow velocities of 5.0 and 7.0 m/s, respectively. The normalized energy dissipation rate is smallest for the newly developed device, while the energy dissipation rate is the largest. The spectral analysis of the velocity data does not show dominant frequencies equal to the blades rotation frequencies; and, the one dimensional kinetic energy and dissipation spectra follow -5/3 and 1/3 power-law relations. It is shown that, the small eddies produced by the newly developed device dissipate the turbulent kinetic energy faster than those corresponding to the rest of the tested turbulence generators.


title: The induced saturation problem for posets

abstract: For a fixed poset $P$, a family $\mathcal F$ of subsets of $[n]$ is induced $P$-saturated if $\mathcal F$ does not contain an induced copy of $P$, but for every subset $S$ of $[n]$ such that $ S\not \in \mathcal F$, $P$ is an induced subposet of $\mathcal F \cup \{S\}$. The size of the smallest such family $\mathcal F$ is denoted by $\text{sat}^* (n,P)$. Keszegh, Lemons, Martin, P\'alv\"olgyi and Patk\'os [Journal of Combinatorial Theory Series A, 2021] proved that there is a dichotomy of behaviour for this parameter: given any poset $P$, either $\text{sat}^* (n,P)=O(1)$ or $\text{sat}^* (n,P)\geq \log _2 n$. In this paper we improve this general result showing that either $\text{sat}^* (n,P)=O(1)$ or $\text{sat}^* (n,P) \geq \min\{ 2 \sqrt{n}, n/2+1\}$. Our proof makes use of a Tur\'an-type result for digraphs. Curiously, it remains open as to whether our result is essentially best possible or not. On the one hand, a conjecture of Ivan states that for the so-called diamond poset $\Diamond$ we have $\text{sat}^* (n,\Diamond)=\Theta (\sqrt{n})$; so if true this conjecture implies our result is tight up to a multiplicative constant. On the other hand, a conjecture of Keszegh, Lemons, Martin, P\'alv\"olgyi and Patk\'os states that given any poset $P$, either $\text{sat}^* (n,P)=O(1)$ or $\text{sat}^* (n,P)\geq n+1$. We prove that this latter conjecture is true for a certain class of posets $P$.


title: Unconventional Collective Resonance as Nonlinear Mechanism of Ectopic Activity in Excitable Media

abstract: Many physical, chemical and biological processes rely on intrinsic oscillations to employ resonance responses to external stimuli of certain frequency. Such resonance phenomena in biological systems are typically explained by one of two mechanisms: either a classical linear resonance of harmonic oscillator, or entrainment and phase locking of nonlinear limit cycle oscillators subjected to periodic forcing. Here, we discover a nonlinear mechanism, which does not require intrinsic oscillations. Instead, the resonant frequency dependence arises from coupling between an excitable and a monostable region of the medium. This composite system is endowed with emergent bistability between a stable steady state and stable spatiotemporal oscillations. The resonant transition from stable state to oscillatory state is induced by waves of particular frequency travelling through the medium. This transition to the spatiotemporal oscillatory state requires accumulation of multiple waves, resulting in the exclusion of lower frequencies. The cutting off of high frequencies is realized by damping of wave amplitude in the monostable zone and then by activating amplitude sensitive dynamics in the monostable units. We demonstrate this new resonance mechanism in a simplistic reaction-diffusion model. Also, we reveal this collective resonance mechanism in in-vitro experiments and detailed biophysical simulations representing a major type of arrhythmia. We further demonstrate, both experimentally and theoretically, that the ongoing spatiotemporal oscillations, such as ectopic activity in cardiac tissue, can be stopped by travelling waves of high frequency. Overall, we claim the universality of this resonance mechanism in a broad class of nonlinear biophysical systems. Specifically, we hypothesize that such phenomena could be found in neuronal systems as an alternative to traditional resonant processes.


title: Generalization-Memorization Machines

abstract: Classifying the training data correctly without over-fitting is one of the goals in machine learning. In this paper, we propose a generalization-memorization mechanism, including a generalization-memorization decision and a memory modeling principle. Under this mechanism, error-based learning machines improve their memorization abilities of training data without over-fitting. Specifically, the generalization-memorization machines (GMM) are proposed by applying this mechanism. The optimization problems in GMM are quadratic programming problems and could be solved efficiently. It should be noted that the recently proposed generalization-memorization kernel and the corresponding support vector machines are the special cases of our GMM. Experimental results show the effectiveness of the proposed GMM both on memorization and generalization.


title: Dreamento: an open-source dream engineering toolbox for sleep EEG wearables

abstract: We introduce Dreamento (Dream engineering toolbox), an open-source Python package for dream engineering using sleep electroencephalography (EEG) wearables. Dreamento main functions are (1) real-time recording, monitoring, analysis, and sensory stimulation, and (2) offline post-processing of the resulting data, both in a graphical user interface (GUI). In real-time, Dreamento is capable of (1) data recording, visualization, and navigation, (2) power-spectrum analysis, (3) automatic sleep scoring, (4) sensory stimulation (visual, auditory, tactile), (5) establishing text-to-speech communication, and (6) managing annotations of automatic and manual events. The offline functions aid in post-processing the acquired data with features to reformat the wearable data and integrate it with non-wearable recorded modalities such as electromyography (EMG). While Dreamento was primarily developed for (lucid) dreaming studies, its applications can be extended to other areas of sleep research such as closed-loop auditory stimulation and targeted memory reactivation.


title: 2n-Stream Conservative Scattering

abstract: We show how to use matrix methods of quantum mechanics to efficiently and accurately calculate axially symmetric radiation transfer in clouds, with conservative scattering of arbitrary anisotropy. Analyses of conservative scattering, where the single scattering albedo is $\tilde\omega =1$ and no energy is exchanged between the radiation and scatterers, began with work by Schwarzschild, Milne, Eddington and others on radiative transfer in stars. There the scattering is isotropic or nearly so. It has been difficult to extend traditional methods to highly anisotropic scattering, like that of sunlight in Earth's clouds. The $2n$-stream method described here is a practical way to handle highly anisotropic, conservative scattering. The basic ideas of the $2n$-stream method are an extension of Wick's seminal work on transport of thermal neutrons by isotropic scattering to scattering with arbitrary anisotropy. How to do this for finite absorption and $\tilde\omega<1$ was described in our previous paper (arXiv:2205.09713v2). But those methods fail for conservative scattering, when $\tilde \omega = 1$. Here we show that minor modifications to the fundamental $2n$-scattering theory for $\tilde\omega <1$ make it suitable for $\tilde\omega = 1$.


title: Analytic Nullstellens\"atze and the model theory of valued fields

abstract: We present a uniform framework for establishing Nullstellens\"atze for power series rings using quantifier elimination results for valued fields. As an application we obtain Nullstellens\"atze for $p$-adic power series (both formal and convergent) analogous to R\"uckert's complex and Risler's real Nullstellensatz, as well as a $p$-adic analytic version of Hilbert's 17th Problem. Analogous statements for restricted power series, both real and $p$-adic, are also considered.


title: Simple master equations for describing driven systems subject to classical non-Markovian noise

abstract: Driven quantum systems subject to non-Markovian noise are typically difficult to model even if the noise is classical. We present a systematic method based on generalized cumulant expansions for deriving a time-local master equation for such systems. This master equation has an intuitive form that directly parallels a standard Lindblad equation, but contains several surprising features: the combination of driving and non-Markovianity results in effective time-dependent dephasing rates that can be negative, and the noise can generate Hamiltonian renormalizations even though it is classical. We analyze in detail the highly relevant case of a Rabi-driven qubit subject to various kinds of non-Markovian noise including $1/f$ fluctuations, finding an excellent agreement between our master equation and numerically-exact simulations over relevant timescales. The approach outlined here is more accurate than commonly employed phenomenological master equations which ignore the interplay between driving and noise.


title: Perturbation of Systems with a First Integral: Motion on the Reeb Graph

abstract: We consider the long-time behavior of systems close to a system with a smooth first integral. Under certain assumptions, the limiting behavior, to some extent, turns out to be universal: it is determined by the first integral, the deterministic perturbation, and the initial point. Furthermore, it is the same for a broad class of noises. In particular, the long-time behavior of a deterministic system can, in a sense, be stochastic and stochastic systems can have a reduced stochasticity. The limiting distribution is calculated explicitly.