title: Optimal Clearing Payments in a Financial Contagion Model

abstract: Financial networks are characterized by complex structures of mutual obligations. These obligations are fulfilled entirely or in part (when defaults occur) via a mechanism called clearing, which determines a set of payments that settle the claims by respecting rules such as limited liability, absolute priority, and proportionality (pro-rated payments). In the presence of shocks on the financial system, however, the clearing mechanism may lead to cascaded defaults and eventually to financial disaster. In this paper, we first study the clearing model under pro-rated payments of Eisenberg and Noe, and we derive novel necessary and sufficient conditions for the uniqueness of the clearing payments, valid for an arbitrary topology of the financial network. Then, we argue that the proportionality rule is one of the factors responsible for cascaded defaults, and that the overall system loss can be reduced if this rule is lifted. The proposed approach thus shifts the focus from the individual interest to the overall system's interest to control and contain adverse effects of cascaded failures, and we show that clearing payments in this setting can be computed by solving suitable convex optimization problems.


title: Fully Onboard AI-powered Human-Drone Pose Estimation on Ultra-low Power Autonomous Flying Nano-UAVs

abstract: Artificial intelligence-powered pocket-sized air robots have the potential to revolutionize the Internet-of-Things ecosystem, acting as autonomous, unobtrusive, and ubiquitous smart sensors. With a few cm$^{2}$ form-factor, nano-sized unmanned aerial vehicles (UAVs) are the natural befit for indoor human-drone interaction missions, as the pose estimation task we address in this work. However, this scenario is challenged by the nano-UAVs' limited payload and computational power that severely relegates the onboard brain to the sub-100 mW microcontroller unit-class. Our work stands at the intersection of the novel parallel ultra-low-power (PULP) architectural paradigm and our general development methodology for deep neural network (DNN) visual pipelines, i.e., covering from perception to control. Addressing the DNN model design, from training and dataset augmentation to 8-bit quantization and deployment, we demonstrate how a PULP-based processor, aboard a nano-UAV, is sufficient for the real-time execution (up to 135 frame/s) of our novel DNN, called PULP-Frontnet. We showcase how, scaling our model's memory and computational requirement, we can significantly improve the onboard inference (top energy efficiency of 0.43 mJ/frame) with no compromise in the quality-of-result vs. a resource-unconstrained baseline (i.e., full-precision DNN). Field experiments demonstrate a closed-loop top-notch autonomous navigation capability, with a heavily resource-constrained 27-gram Crazyflie 2.1 nano-quadrotor. Compared against the control performance achieved using an ideal sensing setup, onboard relative pose inference yields excellent drone behavior in terms of median absolute errors, such as positional (onboard: 41 cm, ideal: 26 cm) and angular (onboard: 3.7$^{\circ}$, ideal: 4.1$^{\circ}$).


title: A new method for constructing linear codes with small hulls

abstract: The hull of a linear code over finite fields is the intersection of the code and its dual, which was introduced by Assmus and Key. In this paper, we develop a method to construct linear codes with trivial hull ( LCD codes) and one-dimensional hull by employing the positive characteristic analogues of Gauss sums. These codes are quasi-abelian, and sometimes doubly circulant. Some sufficient conditions for a linear code to be an LCD code (resp. a linear code with one-dimensional hull) are presented. It is worth mentioning that we present a lower bound on the minimum distances of the constructed linear codes. As an application, using these conditions, we obtain some optimal or almost optimal LCD codes (resp. linear codes with one-dimensional hull) with respect to the online Database of Grassl.


title: Scalable Bayesian computation for crossed and nested hierarchical models

abstract: We develop sampling algorithms to fit Bayesian hierarchical models, the computational complexity of which scales linearly with the number of observations and the number of parameters in the model. We focus on crossed random effect and nested multilevel models, which are used ubiquitously in applied sciences. The posterior dependence in both classes is sparse: in crossed random effects models it resembles a random graph, whereas in nested multilevel models it is tree-structured. For each class we identify a framework for scalable computation, building on previous work. Methods for crossed models are based on extensions of appropriately designed collapsed Gibbs samplers, where we introduce the idea of local centering; while methods for nested models are based on sparse linear algebra and data augmentation. We provide a theoretical analysis of the proposed algorithms in some simplified settings, including a comparison with previously proposed methodologies and an average-case analysis based on random graph theory. Numerical experiments, including two challenging real data analyses on predicting electoral results and real estate prices, compare with off-the-shelf Hamiltonian Monte Carlo, displaying drastic improvement in performance.


title: Fluctuation control of non-thermal orbital order

abstract: Orbitally ordered states exhibit unique features which make them a promising platform for exploring the ultrafast dynamics of long-range order in solids: Their free energy typically has multiple discrete minima, and electric laser fields or selectively excited phonons can exert effective forces that may be used to steer the order parameter through these free energy landscapes. Moreover, their free energy strongly depends on fluctuations, and in some cases restoring forces close to a minimum are exclusively of entropic origin (order-by-disorder mechanisms). This can open pathways to control the dynamics of the order parameter via non-thermal fluctuations. In this work, we study the laser-induced non-equilibrium dynamics in a $120^\circ$ compass model, using time-dependent Ginzburg-Landau theory. We analyze protocols to switch the order parameter between equivalent configurations, with a focus on the interplay between the external force due to the driving field, and the non-thermal entropic forces. In particular, we find that remanent non-thermal fluctuations after some excitation can stabilize the high-symmetry phase even when the homogeneous potential has retrieved its low-temperature form, which facilitates laser-induced switching.


title: Supervisory Control of Multi-Agent Discrete-Event Systems with Partial Observation

abstract: In this paper we investigate multi-agent discrete-event systems with partial observation. The agents can be divided into several groups in each of which the agents have similar (isomorphic) state transition structures, and thus can be relabeled into the same template. Based on the template a scalable supervisor whose state size and computational cost are independent of the number of agents is designed for the case of partial observation. The scalable supervisor under partial observation does not need to be recomputed regardless of how many agents are added to or removed from the system. We generalize our earlier results to partial observation by proposing sufficient conditions for safety and maximal permissiveness of the scalable least restrictive supervisor on the template level. An example is provided to illustrate the proposed scalable supervisory synthesis.


title: Millisecond Magnetars

abstract: Two classes of X-ray/$\gamma$-ray sources, the Soft Gamma Repeaters and the Anomalous X-ray Pulsars have been identified with isolated, slowly spinning magnetars, neutron stars whose emission draws energy from their extremely strong magnetic field ($\sim 10^{15}-10^{16}$ G). Magnetars are believed to form with millisecond spin period and to represent an important fraction of the whole population of young neutron stars. Newborn magnetars can convert very quickly their rotational energy into electromagnetic and/or gravitational waves, by virtue of their strong magnetic fields and fast spins. This chapter provides a brief summary of astrophysical problems and scenarios in which millisecond magnetars are believed to play a key role: these include Gamma Ray Bursts, Supernovae, Gravitational Wave events and Fast Radio Bursts.


title: Asymptotic syzygies of secant varieties of curves

abstract: We prove that the minimal free resolution of the secant variety of a curve is asymptotically pure. As a corollary, we show that the Betti numbers of converge to a normal distribution.


title: Uniformization of compact foliated spaces by surfaces of hyperbolic type

abstract: We give a new proof of the uniformization theorem of the leaves of a lamination by surfaces of hyperbolic conformal type. We use a laminated version of the Ricci flow to prove the existence of a laminated Riemannian metric (smooth on the leaves, transversaly continuous) with leaves of constant Gaussian curvature equal to -1, which is conformally equivalent to the original metric.


title: Building Specifications in the Event-B Institution

abstract: This paper describes a formal semantics for the Event-B specification language using the theory of institutions. We define an institution for Event-B, EVT, and prove that it meets the validity requirements for satisfaction preservation and model amalgamation. We also present a series of functions that show how the constructs of the Event-B specification language can be mapped into our institution. Our semantics sheds new light on the structure of the Event-B language, allowing us to clearly delineate three constituent sub-languages: the superstructure, infrastructure and mathematical languages. One of the principal goals of our semantics is to provide access to the generic modularisation constructs available in institutions, including specification-building operators for parameterisation and refinement. We demonstrate how these features subsume and enhance the corresponding features already present in Event-B through a detailed study of their use in a worked example. We have implemented our approach via a parser and translator for Event-B specifications, EBtoEVT, which also provides a gateway to the Hets toolkit for heterogeneous specification.


title: Coupling fully resolved light particles with the Lattice Boltzmann method on adaptively refined grids

abstract: The simulation of geometrically resolved rigid particles in a fluid relies on coupling algorithms to transfer momentum both ways between the particles and the fluid. In this article, the fluid flow is modeled with a parallel Lattice Boltzmann method using adaptive grid refinement to improve numerical efficiency. The coupling with the particles is realized with the momentum exchange method. When implemented in plain form, instabilities may arise in the coupling when the particles are lighter than the fluid. The algorithm can be stabilized with a virtual mass correction specifically developed for the Lattice Boltzmann method. The method is analyzed for a wide set of physically relevant regimes, varying independently the body-to-fluid density ratio and the relative magnitude of inertial and viscous effects. These studies of a single rising particle exhibit periodic regimes of particle motion as well as chaotic behavior, as previously reported in the literature. The new method is carefully compared with available experimental and numerical results. This serves to validate the presented new coupled Lattice Boltzmann method and additionally it leads to new physical insight for some of the parameter settings.


title: Local and Global Existence of A Nonlocal Equation with A Singular Integral Drift Term

abstract: We study an initial value problem with fractional Laplacian and a singular drift term, and obtain local and global existence theorems similar to the results in Jourdain et al.(2005).


title: Localized Reduced Basis Additive Schwarz Methods

abstract: Reduced basis methods build low-rank approximation spaces for the solution sets of parameterized PDEs by computing solutions of the given PDE for appropriately selected snapshot parameters. Localized reduced basis methods reduce the offline cost of computing these snapshot solutions by instead constructing a global space from spatially localized less expensive problems. In the case of online enrichment, these local problems are iteratively solved in regions of high residual and correspond to subdomain solves in domain decomposition methods. We show in this note that indeed there is a close relationship between online-enriched localized reduced basis and domain decomposition methods by introducing a Localized Reduced Basis Additive Schwarz method (LRBAS), which can be interpreted as a locally adaptive multi-preconditioning scheme for the CG method.


title: Predicting Covid-19 EMS Incidents from Daily Hospitalization Trends

abstract: Introduction: The aim of our retrospective study was to quantify the impact of Covid-19 on the temporal distribution of Emergency Medical Services (EMS) demand in Travis County, Austin, Texas and propose a robust model to forecast Covid-19 EMS incidents. Methods: We analyzed the temporal distribution of EMS calls in the Austin-Travis County area between January 1st, 2019, and December 31st, 2020. Change point detection was performed to identify critical dates marking changes in EMS call distributions, and time series regression was applied for forecasting Covid-19 EMS incidents. Results: Two critical dates marked the impact of Covid-19 on the distribution of EMS calls: March 17th, when the daily number of non-pandemic EMS incidents dropped significantly, and May 13th, by which the daily number of EMS calls climbed back to 75% of the number in pre-Covid-19 time. The new daily count of the hospitalization of Covid-19 patients alone proves a powerful predictor of the number of pandemic EMS calls, with an r2 value equal to 0.85. In particular, for every 2.5 cases where EMS takes a Covid-19 patient to a hospital, one person is admitted. Conclusion: The mean daily number of non-pandemic EMS demand was significantly less than the period before Covid-19 pandemic. The number of EMS calls for Covid-19 symptoms can be predicted from the daily new hospitalization of Covid-19 patients. These findings may be of interest to EMS departments as they plan for future pandemics, including the ability to predict pandemic-related calls in an effort to adjust a targeted response.


title: Rotating deformed halo nuclei and shape decoupling effects

abstract: We explore the rotational feature of deformed halos by performing the angular momentum projection (AMP) on the ground state wave functions obtained from the deformed relativistic Hartree-Bogoliubov theory (DRHBc) in continuum. The DRHBc+AMP approach self-consistently describes the coupling between single particle bound states and the continuum not only in the ground state but also in rotational states. The rotational modes of deformed halos in $^{42}$Mg and $^{44}$Mg are investigated by studying properties of rotational states such as excitation energy, configuration, and density distribution. Our study demonstrates that the deformed halo structure persists from the ground state in the intrinsic frame to collective states. Especially, the typical behavior of shape decoupling effects in rotating deformed halo nuclei is revealed.


title: Gouy-phase-mediated propagation variations and revivals of transverse structure in vectorially structured light

abstract: Exploring the physics and potential applications of vectorially structured light with propagation-invariant transverse structures has benefited many areas of modern optics and photonics. In this paper, we investigate the non-eigen vector modes of paraxial light fields, focusing on the propagation variations and revivals of their transverse structures, including both spatial and polarization structures. We show that the physical mechanism behind the variations and revivals of their transverse structure is linked to the evolution of the intramodal phases between the constituting spatial modes. Such evolution originates from fractional Gouy phases, or rather, Geometric-phase difference between spatial modes with different orders under a same unitary transformation. This underlying principle, provides a general guideline for shaping vectorially structured light with custom propagation-evolution properties, and may also inspire a wide variety of new applications based on structured light.


title: Feedback from Pixels: Output Regulation via Learning-Based Scene View Synthesis

abstract: We propose a novel controller synthesis involving feedback from pixels, whereby the measurement is a high dimensional signal representing a pixelated image with Red-Green-Blue (RGB) values. The approach neither requires feature extraction, nor object detection, nor visual correspondence. The control policy does not involve the estimation of states or similar latent representations. Instead, tracking is achieved directly in image space, with a model of the reference signal embedded as required by the internal model principle. The reference signal is generated by a neural network with learning-based scene view synthesis capabilities. Our approach does not require an end-to-end learning of a pixel-to-action control policy. The approach is applied to a motion control problem, namely the longitudinal dynamics of a car-following problem. We show how this approach lend itself to a tractable stability analysis with associated bounds critical to establishing trustworthiness and interpretability of the closed-loop dynamics.


title: $p$-adic Directions of Primitive Vectors

abstract: Linnik type problems concern the distribution of projections of integral points on the unit sphere as their norm increases, and different generalizations of this phenomenon. Our work addresses a question of this type: we prove the uniform distribution of the projections of primitive $\mathbb{Z}^{2}$ points in the $p$-adic unit sphere, as their (real) norm tends to infinity. The proof is via counting lattice points in semi-simple $S$-arithmetic groups.


title: Simulation techniques for modified gravity

abstract: The standard paradigm of cosmology assumes General Relativity (GR) is a valid theory for gravity at scales in which it has not been properly tested. Developing novel tests of GR and its alternatives is crucial if we want to give strength to the model or find departures from GR in the data. Since alternatives to GR are usually defined through nonlinear equations, designing new tests for these theories implies a jump in complexity and thus, a need for refining the simulation techniques. We summarize existing techniques for dealing with modified gravity (MG) in the context of cosmological simulations. $N$-body codes for MG are usually based on standard gravity codes. We describe the required extensions, classifying the models not according to their original motivation, but by the numerical challenges that must be faced by numericists. MG models usually give rise to elliptic equations, for which multigrid techniques are well suited. Thus, we devote a large fraction of this review to describing this particular technique. Contrary to other reviews on multigrid methods, we focus on the specific techniques that are required to solve MG equations and describe useful tricks. Finally, we describe extensions for going beyond the static approximation and dealing with baryons.

categories: astro-ph.CO astro-ph.GA gr-qc physics.comp-ph


title: Accelerating SLIDE Deep Learning on Modern CPUs: Vectorization, Quantizations, Memory Optimizations, and More

abstract: Deep learning implementations on CPUs (Central Processing Units) are gaining more traction. Enhanced AI capabilities on commodity x86 architectures are commercially appealing due to the reuse of existing hardware and virtualization ease. A notable work in this direction is the SLIDE system. SLIDE is a C++ implementation of a sparse hash table based back-propagation, which was shown to be significantly faster than GPUs in training hundreds of million parameter neural models. In this paper, we argue that SLIDE's current implementation is sub-optimal and does not exploit several opportunities available in modern CPUs. In particular, we show how SLIDE's computations allow for a unique possibility of vectorization via AVX (Advanced Vector Extensions)-512. Furthermore, we highlight opportunities for different kinds of memory optimization and quantizations. Combining all of them, we obtain up to 7x speedup in the computations on the same hardware. Our experiments are focused on large (hundreds of millions of parameters) recommendation and NLP models. Our work highlights several novel perspectives and opportunities for implementing randomized algorithms for deep learning on modern CPUs. We provide the code and benchmark scripts at https://github.com/RUSH-LAB/SLIDE


title: Deep Label Fusion: A 3D End-to-End Hybrid Multi-Atlas Segmentation and Deep Learning Pipeline

abstract: Deep learning (DL) is the state-of-the-art methodology in various medical image segmentation tasks. However, it requires relatively large amounts of manually labeled training data, which may be infeasible to generate in some applications. In addition, DL methods have relatively poor generalizability to out-of-sample data. Multi-atlas segmentation (MAS), on the other hand, has promising performance using limited amounts of training data and good generalizability. A hybrid method that integrates the high accuracy of DL and good generalizability of MAS is highly desired and could play an important role in segmentation problems where manually labeled data is hard to generate. Most of the prior work focuses on improving single components of MAS using DL rather than directly optimizing the final segmentation accuracy via an end-to-end pipeline. Only one study explored this idea in binary segmentation of 2D images, but it remains unknown whether it generalizes well to multi-class 3D segmentation problems. In this study, we propose a 3D end-to-end hybrid pipeline, named deep label fusion (DLF), that takes advantage of the strengths of MAS and DL. Experimental results demonstrate that DLF yields significant improvements over conventional label fusion methods and U-Net, a direct DL approach, in the context of segmenting medial temporal lobe subregions using 3T T1-weighted and T2-weighted MRI. Further, when applied to an unseen similar dataset acquired in 7T, DLF maintains its superior performance, which demonstrates its good generalizability.


title: An operadic approach to substitution in Lie-Butcher series

abstract: The paper follows an operadic approach to provide a bialgebraic description of substitution for Lie-Butcher series. We first show how the well-known bialgebraic description for substitution in Butcher's $B$-series can be obtained from the pre-Lie operad. We then apply the same construction to the post-Lie operad to arrive at a bialgebra $\mathcal{Q}$. By considering a module over the post-Lie operad, we get a cointeraction between $\mathcal{Q}$ and the Hopf algebra $\mathcal{H}_N$ that describes composition for Lie-Butcher series. We use this coaction to describe substitution for Lie-Butcher series.


title: The influence of metallicity on the Leavitt Law from geometrical distances of Milky Way and Magellanic Clouds Cepheids

abstract: The Cepheid Period-Luminosity (PL) relation is the key tool for measuring astronomical distances and for establishing the extragalactic distance scale. In particular, the local value of the Hubble constant ($H_0$) strongly depends on Cepheid distance measurements. The recent Gaia Data Releases and other parallax measurements from the Hubble Space Telescope (HST) already enabled to improve the accuracy of the slope ($\alpha$) and intercept ($\beta$) of the PL relation. However, the dependence of this law on metallicity is still largely debated. In this paper, we combine three samples of Cepheids in the Milky Way (MW), the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC) in order to derive the metallicity term (hereafter $\gamma$) of the PL relation. The recent publication of extremely precise LMC and SMC distances based on late-type detached eclipsing binary systems (DEBs) provides a solid anchor for the Magellanic Clouds. In the MW, we adopt Cepheid parallaxes from the early third Gaia Data Release. We derive the metallicity effect in $V$, $I$, $J$, $H$, $K_S$, $W_{VI}$ and $W_{JK}$. In the $K_S$ band we report a metallicity effect of $-0.221 \pm 0.051$ mag/dex, the negative sign meaning that more metal-rich Cepheids are intrinsically brighter than their more metal-poor counterparts of the same pulsation period.


title: Sewer-ML: A Multi-Label Sewer Defect Classification Dataset and Benchmark

abstract: Perhaps surprisingly sewerage infrastructure is one of the most costly infrastructures in modern society. Sewer pipes are manually inspected to determine whether the pipes are defective. However, this process is limited by the number of qualified inspectors and the time it takes to inspect a pipe. Automatization of this process is therefore of high interest. So far, the success of computer vision approaches for sewer defect classification has been limited when compared to the success in other fields mainly due to the lack of public datasets. To this end, in this work we present a large novel and publicly available multi-label classification dataset for image-based sewer defect classification called Sewer-ML. The Sewer-ML dataset consists of 1.3 million images annotated by professional sewer inspectors from three different utility companies across nine years. Together with the dataset, we also present a benchmark algorithm and a novel metric for assessing performance. The benchmark algorithm is a result of evaluating 12 state-of-the-art algorithms, six from the sewer defect classification domain and six from the multi-label classification domain, and combining the best performing algorithms. The novel metric is a class-importance weighted F2 score, $\text{F}2_{\text{CIW}}$, reflecting the economic impact of each class, used together with the normal pipe F1 score, $\text{F}1_{\text{Normal}}$. The benchmark algorithm achieves an $\text{F}2_{\text{CIW}}$ score of 55.11% and $\text{F}1_{\text{Normal}}$ score of 90.94%, leaving ample room for improvement on the Sewer-ML dataset. The code, models, and dataset are available at the project page https://vap.aau.dk/sewer-ml/


title: Delta-kick Squeezing

abstract: We explore the possibility to overcome the standard quantum limit (SQL) in a free-fall atom interferometer using a Bose-Einstein condensate (BEC) in either of the two relevant cases of Bragg or Raman scattering light pulses. The generation of entanglement in the BEC is dramatically enhanced by amplifying the atom-atom interactions via the rapid action of an external trap focusing the matter-waves to significantly increase the atomic densities during a preparation stage -- a technique we refer to as delta-kick squeezing (DKS). The action of a second DKS operation at the end of the interferometry sequence allows to implement a non-linear readout scheme making the sub-SQL sensitivity highly robust against imperfect atom counting detection. We predict more than 30 dB of sensitivity gain beyond the SQL for the variance, assuming realistic parameters and $10^6$ atoms.


title: Bilinear Classes: A Structural Framework for Provable Generalization in RL

abstract: This work introduces Bilinear Classes, a new structural framework, which permit generalization in reinforcement learning in a wide variety of settings through the use of function approximation. The framework incorporates nearly all existing models in which a polynomial sample complexity is achievable, and, notably, also includes new models, such as the Linear $Q^*/V^*$ model in which both the optimal $Q$-function and the optimal $V$-function are linear in some known feature space. Our main result provides an RL algorithm which has polynomial sample complexity for Bilinear Classes; notably, this sample complexity is stated in terms of a reduction to the generalization error of an underlying supervised learning sub-problem. These bounds nearly match the best known sample complexity bounds for existing models. Furthermore, this framework also extends to the infinite dimensional (RKHS) setting: for the the Linear $Q^*/V^*$ model, linear MDPs, and linear mixture MDPs, we provide sample complexities that have no explicit dependence on the explicit feature dimension (which could be infinite), but instead depends only on information theoretic quantities.


title: Nonequilibrium control of thermal and mechanical changes in a levitated system

abstract: Fluctuation theorems are fundamental extensions of the second law of thermodynamics for small nonequilibrium systems. While work and heat are equally important forms of energy exchange, fluctuation relations have not been experimentally assessed for the generic situation of simultaneous mechanical and thermal changes. Thermal driving is indeed generally slow and more difficult to realize than mechanical driving. Here, we use feedback cooling techniques to implement fast and controlled temperature variations of an underdamped levitated microparticle that are one order of magnitude faster than the equilibration time. Combining mechanical and thermal control, we verify the validity of a fluctuation theorem that accounts for both contributions, well beyond the range of linear response theory. Our results allow the investigation of general far-from-equilibrium processes in microscopic systems that involve fast mechanical and thermal changes at the same time.


title: Experimental Realisation of a \pi/2 Vortex Mode Converter for Electrons Using a Spherical Aberration Corrector

abstract: In light optics, beams with orbital angular momentum (OAM) can be produced by employing a properly-tuned two-cylinder-lens arrangement, also called $\pi$/2 mode converter. It is not possible to convey this concept directly to the beam in an electron microscope due to the non-existence of cylinder lenses in commercial transmission electron microscope (TEM). A viable work-around are readily-available electron optical elements in the form of quadrupole lenses. In a proof-of-principle experiment in 2012, it has been shown that a single quadrupole in combination with a Hilbert phase plate produces a spatially-confined, transient vortex mode. Here, an analogue to an optical $\pi$/2 mode converter is realized by repurposing a CEOS DCOR probe corrector in an aberration corrected TEM in a way that it resembles a dual cylinder lens using two quadrupoles. In order to verify the presence of OAM in the output beam, a fork dislocation grating is used as an OAM analyser. The possibility to use magnetic quadrupole fields instead of, e.g., prefabricated fork dislocation gratings to produce electron beams carrying OAM enhances the beam brightness by almost an order of magnitude and delivers switchable high-mode purity vortex beams without unwanted side-bands.


title: Continental generalization of an AI system for clinical seizure recognition

abstract: Electroencephalogram (EEG) monitoring and objective seizure identification is an essential clinical investigation for some patients with epilepsy. Accurate annotation is done through a time-consuming process by EEG specialists. Computer-assisted systems for seizure detection currently lack extensive clinical utility due to retrospective, patient-specific, and/or irreproducible studies that result in low sensitivity or high false positives in clinical tests. We aim to significantly reduce the time and resources on data annotation by demonstrating a continental generalization of seizure detection that balances sensitivity and specificity. This is a prospective inference test of artificial intelligence on nearly 14,590 hours of adult EEG data from patients with epilepsy between 2011 and 2019 in a hospital in Sydney, Australia. The inference set includes patients with different types and frequencies of seizures across a wide range of ages and EEG recording hours. We validated our inference model in an AI-assisted mode with a human expert arbiter and a result review panel of expert neurologists and EEG specialists on 66 sessions to demonstrate achievement of the same performance with over an order-of-magnitude reduction in time. Our inference on 1,006 EEG recording sessions on the Australian dataset achieved 76.68% with nearly 56 [0, 115] false alarms per 24 hours on average, against legacy ground-truth annotations by human experts, conducted independently over nine years. Our pilot test of 66 sessions with a human arbiter, and reviewed ground truth by a panel of experts, confirmed an identical human performance of 92.19% with an AI-assisted system, while the time requirements reduce significantly from 90 to 7.62 minutes on average.


title: From Static to Dynamic Prediction: Wildfire Risk Assessment Based on Multiple Environmental Factors

abstract: Wildfire is one of the biggest disasters that frequently occurs on the west coast of the United States. Many efforts have been made to understand the causes of the increases in wildfire intensity and frequency in recent years. In this work, we propose static and dynamic prediction models to analyze and assess the areas with high wildfire risks in California by utilizing a multitude of environmental data including population density, Normalized Difference Vegetation Index (NDVI), Palmer Drought Severity Index (PDSI), tree mortality area, tree mortality number, and altitude. Moreover, we focus on a better understanding of the impacts of different factors so as to inform preventive actions. To validate our models and findings, we divide the land of California into 4,242 grids of 0.1 degrees $\times$ 0.1 degrees in latitude and longitude, and compute the risk of each grid based on spatial and temporal conditions. To verify the generalizability of our models, we further expand the scope of wildfire risk assessment from California to Washington without any fine tuning. By performing counterfactual analysis, we uncover the effects of several possible methods on reducing the number of high risk wildfires. Taken together, our study has the potential to estimate, monitor, and reduce the risks of wildfires across diverse areas provided that such environment data is available.


title: On the electron energy distribution function in the high power impulse magnetron sputtering discharge

abstract: We apply the Ionization Region Model (IRM) and the Orsay Boltzmann equation for ELectrons coupled with Ionization and eXcited states kinetics (OBELIX) model to study the electron kinetics of a high power impulse magnetron sputtering (HiPIMS) discharge. In the IRM the bulk (cold) electrons are assumed to exhibit a Maxwellian energy distribution and the secondary (hot) electrons, emitted from the target surface upon ion bombardment, are treated as a high energy tail, while in the OBELIX the electron energy distribution is calculated self-consistently using an isotropic Boltzmann equation. The two models are merged in the sense that the output from the IRM is used as an input for OBELIX. The temporal evolutions of the particle densities are found to agree very well between the two models. Furthermore, a very good agreement is demonstrated between the bi-Maxwellian electron energy distribution assumed by the IRM and the electron energy distribution calculated by the OBELIX model. It can therefore be concluded that assuming a bi-Maxwellian electron energy distribution, constituting a cold bulk electron group and a hot secondary electron group, is a good approximation for modeling the HiPIMS discharge.


title: A Unified Framework for IRS Enabled Wireless Powered Sensor Networks

abstract: This paper unveils the importance of intelligent reflecting surface (IRS) in a wireless powered sensor network (WPSN). Specifically, a multi-antenna power station (PS) employs energy beamforming to provide wireless charging for multiple Internet of Things (IoT) devices, which utilize the harvested energy to deliver their own messages to an access point (AP). Meanwhile, an IRS is deployed to enhance the performances of wireless energy transfer (WET) and wireless information transfer (WIT) by intelligently adjusting the phase shift of each reflecting element. To evaluate the performance of this IRS assisted WPSN, we are interested in maximizing its system sum throughput to jointly optimize the energy beamforming of the PS, the transmission time allocation, as well as the phase shifts of the WET and WIT phases. The formulated problem is not jointly convex due to the multiple coupled variables. To deal with its non-convexity, we first independently find the phase shifts of the WIT phase in closed-form. We further propose an alternating optimization (AO) algorithm to iteratively solve the sum throughput maximization problem. To be specific, a semidefinite programming (SDP) relaxation approach is adopted to design the energy beamforming and the time allocation for given phase shifts of WET phase, which is then optimized for given energy beamforming and time allocation. Moreover, we propose an AO low-complexity scheme to significantly reduce the computational complexity incurred by the SDP relaxation, where the optimal closed-form energy beamforming, time allocation, and phase shifts of the WET phase are derived. Finally, numerical results are demonstrated to validate the effectiveness of the proposed algorithm, and highlight the beneficial role of the IRS in comparison to the benchmark schemes.


title: Frobenius Numbers and Automatic Sequences

abstract: The Frobenius number $g(S)$ of a set $S$ of non-negative integers with $\gcd 1$ is the largest integer not expressible as a linear combination of elements of $S$. Given a sequence ${\bf s} = (s_i)_{i \geq 0}$, we can define the associated sequence $G_{\bf s} (i) = g(\{ s_i,s_{i+1},\ldots \})$. In this paper we compute $G_{\bf s} (i)$ for some classical automatic sequences: the evil numbers, the odious numbers, and the lower and upper Wythoff sequences. In contrast with the usual methods, our proofs are based largely on automata theory and logic.


title: Joint Parameter Discovery and Generative Modeling of Dynamic Systems

abstract: Given an unknown dynamic system such as a coupled harmonic oscillator with $n$ springs and point masses. We are often interested in gaining insights into its physical parameters, i.e. stiffnesses and masses, by observing trajectories of motion. How do we achieve this from video frames or time-series data and without the knowledge of the dynamics model? We present a neural framework for estimating physical parameters in a manner consistent with the underlying physics. The neural framework uses a deep latent variable model to disentangle the system physical parameters from canonical coordinate observations. It then returns a Hamiltonian parameterization that generalizes well with respect to the discovered physical parameters. We tested our framework with simple harmonic oscillators, $n=1$, and noisy observations and show that it discovers the underlying system parameters and generalizes well with respect to these discovered parameters. Our model also extrapolates the dynamics of the system beyond the training interval and outperforms a non-physically constrained baseline model. Our source code and datasets can be found at this URL: https://github.com/gbarber94/ConSciNet.


title: Diffuse optics for glaciology

abstract: Optical probing of glaciers has the potential for tremendous impact on environmental science. However, glacier ice is turbid, which prohibits the use of most established optical measurements for determining a glacier's interior structure. Here, we propose a method for determining the depth, scattering and absorption length based upon diffuse propagation of short optical pulses. Our model allows us to extract several characteristics of the glacier. Performing Monte Carlo simulations implementing Mie scattering and mixed boundary conditions, we show that the proposed approach should be feasible with current technology. The results suggest that optical properties and geometry of the glacier can be extracted from realistic measurements, which could be implemented with low cost and small footprint.


title: On $p$-adic $L$-functions for $GL_{2n}$ in finite slope Shalika families

abstract: In this paper, we prove new results on the geometry of the cuspidal eigenvariety for $\mathrm{GL}_{2n}$ over a totally real number field $F$ at classical points admitting Shalika models. We also construct $p$-adic $L$-functions over the eigenvariety around these points. Our proofs proceed in the opposite direction to established methods: rather than using the geometry of eigenvarieties to deduce results about $p$-adic $L$-functions, we instead show that -- via evaluation maps on parahoric overconvergent cohomology groups -- non-vanishing of a $p$-adic $L$-function implies smoothness of the eigenvariety at such points. More precisely, we attach a $p$-adic $L$-function to a non-critical refinement $\tilde\pi$ of a regular algebraic cuspidal automorphic representation $\pi$ of $\mathrm{GL}_{2n}/F$ which is spherical at $p$ and admits a Shalika model. This gives the first construction of $p$-adic $L$-functions in this generality beyond the $p$-ordinary setting. Further, when $\pi$ has regular weight and the corresponding $p$-adic Galois representation is irreducible, we show that the parabolic eigenvariety for $\mathrm{GL}_{2n}/F$ is \'etale at $\tilde\pi$ over an $([F:\mathbb{Q}]+1)$-dimensional weight space and contains a dense set of classical points admitting Shalika models. Finally, under a hypothesis on the local Shalika models at bad places which is empty for $\pi$ of level 1, we construct a $p$-adic $L$-function for the family.


title: Evaluating Global Blockage engineering parametrizations with LES

abstract: Under the term global blockage, the cumulative induction of wind turbines in a wind farm has been recently suspected to be responsible for observed overestimations of the energy yield in large-size wind farms. In this paper, the practice of modeling this effect after linear superposition of single turbine inductions, calculated with three of the most recent analytic models, is compared to Large-Eddy-Simulations of wind farms. We compare the models across two different farms, composed of 9 and 49 turbines, with two different heights of the atmospheric boundary layer, 300 and 500 m. The results show hat the differences between the analytical models are negligible while they substantially differ from the LES results. The linear superposition of induction consistently underestimates the velocity deficit in front of the farm with an error that increases as the wind farm size grows and the ABL height decreases. Also, when calculating the power output at the turbines of the farm, all the analytical models considered do not agree with the LES. These comparisons reveal that the farm interactions with the atmospheric boundary layer may highly outclass the turbine induction in determining the extent of the global blockage effect. Therefore, we present a first dimensional approach to the problem based on LES, aimed at simplifying its characterization.


title: IA Planner: Motion Planning Using Instantaneous Analysis for Autonomous Vehicle in the Dense Dynamic Scenarios on Highways

abstract: In dense and dynamic scenarios, planning a safe and comfortable trajectory is full of challenges when traffic participants are driving at high speed. The classic graph search and sampling methods first perform path planning and then configure the corresponding speed, which lacks a strategy to deal with the high-speed obstacles. Decoupling optimization methods perform motion planning in the S-L and S-T domains respectively. These methods require a large free configuration space to plan the lane change trajectory. In dense dynamic scenes, it is easy to cause the failure of trajectory planning and be cut in by others, causing slow driving speed and bring safety hazards. We analyze the collision relationship in the spatio-temporal domain, and propose an instantaneous analysis model which only analyzes the collision relationship at the same time. In the model, the collision-free constraints in 3D spatio-temporal domain is projected to the 2D space domain to remove redundant constraints and reduce computational complexity. Experimental results show that our method can plan a safe and comfortable lane-changing trajectory in dense dynamic scenarios. At the same time, it improves traffic efficiency and increases ride comfort.


title: Local fluid pressure perturbations inside faults matter in unexpectedly large injection-triggered earthquakes

abstract: Anticipating the maximum magnitude of injection-triggered earthquakes is highly valuable for the safe and efficient exploitation of geoenergies. The recent work by Li et al.(2021) reached the conclusion that unexpectedly large injection-triggered earthquakes are primarily caused by large pre-existing critical shear stresses on seismogenic faults. Also of great interest is the proposal of the ratio of fault slip to dilation as an index to anticipate the fault rupture energetics. However, their fluid injection experiments were conducted under fully drained conditions, where the fluid pressure distribution on the fault plane is always uniform. As already pointed out by the authors, local fluid pressure perturbations inside faults under locally undrained conditions also have the potential to cause unexpectedly large seismic events. Here we add more datasets and possible explanations for this viewpoint, which has not been explored in such detail by Li et al.(2021).


title: Performance Analysis of Deep Learning Workloads on a Composable System

abstract: A composable infrastructure is defined as resources, such as compute, storage, accelerators and networking, that are shared in a pool and that can be grouped in various configurations to meet application requirements. This freedom to 'mix and match' resources dynamically allows for experimentation early in the design cycle, prior to the final architectural design or hardware implementation of a system. This design provides flexibility to serve a variety of workloads and provides a dynamic co-design platform that allows experiments and measurements in a controlled manner. For instance, key performance bottlenecks can be revealed early on in the experimentation phase thus avoiding costly and time consuming mistakes. Additionally, various system-level topologies can be evaluated when experimenting with new System on Chip (SoCs) and new accelerator types. This paper details the design of an enterprise composable infrastructure that we have implemented and made available to our partners in the IBM Research AI Hardware Center (AIHC). Our experimental evaluations on the composable system give insights into how the system works and evaluates the impact of various resource aggregations and reconfigurations on representative deep learning benchmarks.


title: Copula Averaging for Tail Dependence in Insurance Claims Data

abstract: Analysing dependent risks is an important task for insurance companies. A dependency is reflected in the fact that information about one random variable provides information about the likely distribution of values of another random variable. Insurance companies in particular must investigate such dependencies between different lines of business and the effects that an extreme loss event, such as an earthquake or hurricane, has across multiple lines of business simultaneously. Copulas provide a popular model-based approach to analysing the dependency between risks, and the coefficient of tail dependence is a measure of dependence for extreme losses. Besides commonly used empirical estimators for estimating the tail dependence coefficient, copula fitting can lead to estimation of such coefficients directly or can verify their existence. Generally, a range of copula models is available to fit a data set well, leading to multiple different tail dependence results; a method based on Bayesian model averaging is designed to obtain a unified estimate of tail dependence. In this article, this model-based coefficient estimation method is illustrated through a variety of copula fitting approaches and results are presented for several simulated data sets and also a real general insurance loss data set.


title: Composite search of active particles in three-dimensional space based on non-directional cues

abstract: We theoretically address minimal search strategies of active, self-propelled particles towards hidden targets in three-dimensional space. The particles can sense if a target is close, e.g., by detecting signaling molecules released by a target, but they cannot deduce any directional cues. We focus on composite search strategies, where particles switch between extensive outer search and intensive inner search; inner search is started when the proximity of a target is detected and ends again when a certain inner search time has elapsed. In the simplest strategy, active particles move ballistically during outer search, and transiently reduce their directional persistence during inner search. In a second, adaptive strategy, particles exploit a dynamic scattering effect by reducing directional persistence only outside a well-defined target zone. These two search strategies require only minimal information processing capabilities and a single binary or tertiary internal state, respectively, yet increases the rate of target encounter substantially. The optimal inner search time scales as a power-law with exponent -2/3 with target density, reflecting a trade-off between exploration and exploitation.


title: Dynamical renormalization group for mode-coupling field theories with solenoidal constraint

abstract: The recent inflow of empirical data about the collective behaviour of strongly correlated biological systems has brought field theory and the renormalization group into the biophysical arena. Experiments on bird flocks and insect swarms show that social forces act on the particles' velocity through the generator of its rotations, namely the spin, indicating that mode-coupling field theories are necessary to reproduce the correct dynamical behaviour. Unfortunately, a theory for three coupled fields - density, velocity and spin - has a prohibitive degree of intricacy. A simplifying path consists in getting rid of density fluctuations by studying incompressible systems. This requires imposing a solenoidal constraint on the primary field, an unsolved problem even for equilibrium mode-coupling theories. Here, we perform an equilibrium dynamic renormalization group analysis of a mode-coupling field theory subject to a solenoidal constraint; using the classification of Halperin and Hohenberg, we can dub this case as a solenoidal Model G. We demonstrate that the constraint produces a new vertex that mixes static and dynamical coupling constants, and that this vertex is essential to grant the closure of the renormalization group structure and the consistency of dynamics with statics. Interestingly, although the solenoidal constraint leads to a modification of the static universality class, we find that it does not change the dynamical universality class, a result that seems to represent an exception to the general rule that dynamical universality classes are narrower than static ones. Our results constitute a solid stepping stone in the admittedly large chasm towards developing an off-equilibrium mode-coupling theory of biological groups.


title: Computational and Experimental Characterization of NF023, A Candidate Anticancer Compound Inhibiting cIAP2/TRAF2 Assembly

abstract: Protein-protein interactions are the basis of many important physiological processes and are currently promising, yet difficult, targets for drug discovery. In this context, inhibitor of apoptosis proteins (IAPs)-mediated interactions are pivotal for cancer cell survival; the interaction of the BIR1 domain of cIAP2 with TRAF2 was shown to lead the recruitment of cIAPs to the TNF receptor, promoting the activation of the NF-\kappa B survival pathway. In this work, using a combined in silico-in vitro approach, we identified a drug-like molecule, NF023, able to disrupt cIAP2 interaction with TRAF2. We demonstrated in vitro its ability to interfere with the assembly of the cIAP2-BIR1/TRAF2 complex and performed a thorough characterization of the compound's mode of action through 248 parallel unbiased molecular dynamics simulations of 300 ns (totaling almost 75 {\mu}s of all-atom sampling), which identified multiple binding modes to the BIR1 domain of cIAP2 via clustering and ensemble docking. NF023 is, thus, a promising protein-protein interaction disruptor, representing a starting point to develop modulators of NF-\kappa B-mediated cell survival in cancer. This study represents a model procedure that shows the use of large-scale molecular dynamics methods to typify promiscuous interactors.


title: Predicting Drug-Drug Interactions from Heterogeneous Data: An Embedding Approach

abstract: Predicting and discovering drug-drug interactions (DDIs) using machine learning has been studied extensively. However, most of the approaches have focused on text data or textual representation of the drug structures. We present the first work that uses multiple data sources such as drug structure images, drug structure string representation and relational representation of drug relationships as the input. To this effect, we exploit the recent advances in deep networks to integrate these varied sources of inputs in predicting DDIs. Our empirical evaluation against several state-of-the-art methods using standalone different data types for drugs clearly demonstrate the efficacy of combining heterogeneous data in predicting DDIs.


title: Spin-thermal shear coupling in a relativistic fluid

abstract: We show that spin polarization of a fermion in a relativistic fluid at local thermodynamic equilibrium can be generated by the symmetric derivative of the four-temperature vector, defined as thermal shear. As a consequence, besides vorticity, acceleration and temperature gradient, also the shear tensor contributes to the polarization of particles in a fluid. This contribution to the spin polarization vector, which is entirely non-dissipative, adds to the well known term proportional to thermal vorticity and may thus have important consequences for the solution of the local polarization puzzles observed in relativistic heavy ion collisions.


title: Play the Shannon Game With Language Models: A Human-Free Approach to Summary Evaluation

abstract: The goal of a summary is to concisely state the most important information in a document. With this principle in mind, we introduce new reference-free summary evaluation metrics that use a pretrained language model to estimate the information content shared between a document and its summary. These metrics are a modern take on the Shannon Game, a method for summary quality scoring proposed decades ago, where we replace human annotators with language models. We also view these metrics as an extension of BLANC, a recently proposed approach to summary quality measurement based on the performance of a language model with and without the help of a summary. Using transformer based language models, we empirically verify that our metrics achieve state-of-the-art correlation with human judgement of the summary quality dimensions of both coherence and relevance, as well as competitive correlation with human judgement of consistency and fluency.


title: Robustness via Cross-Domain Ensembles

abstract: We present a method for making neural network predictions robust to shifts from the training data distribution. The proposed method is based on making predictions via a diverse set of cues (called 'middle domains') and ensembling them into one strong prediction. The premise of the idea is that predictions made via different cues respond differently to a distribution shift, hence one should be able to merge them into one robust final prediction. We perform the merging in a straightforward but principled manner based on the uncertainty associated with each prediction. The evaluations are performed using multiple tasks and datasets (Taskonomy, Replica, ImageNet, CIFAR) under a wide range of adversarial and non-adversarial distribution shifts which demonstrate the proposed method is considerably more robust than its standard learning counterpart, conventional deep ensembles, and several other baselines.


title: Quasinormal modes in the field of a dyon-like dilatonic black hole

abstract: Quasinormal modes of massless test scalar field in the background of gravitational field for a non-extremal dilatonic dyonic black hole are explored. The dyon-like black hole solution is considered in the gravitational $4d$ model involving two scalar fields and two 2-forms. It is governed by two 2-dimensional dilatonic coupling vectors $\vec{\lambda}_i$ obeying $\vec{\lambda}_i (\vec{\lambda}_1 + \vec{\lambda}_2) > 0$, $i =1,2$. The first law of black hole thermodynamics is given and the Smarr relation is verified. Quasinormal modes for a massless scalar (test) field in the eikonal approximation are obtained and analysed. These modes depend upon a dimensionless parameter $a$ ($0 < a \leq 2$) which is a function of $\vec{\lambda}_i$. For limiting strong ($a = +0$) and weak ($a = 2$) coupling cases, they coincide with the well-known results for the Schwarzschild and Reissner-Nordstr\"om solutions. It is shown that the Hod conjecture, connecting the damping rate and the Hawking temperature, is satisfied for $0 < a \leq 1$ and all allowed values of parameters.


title: Tailored high-contrast attosecond electron pulses for coherent excitation and scattering

abstract: Temporally shaping the density of electron beams using light forms the basis for a wide range of established and emerging technologies, including free-electron lasers and attosecond electron microscopy. The modulation depth of compressed electron pulses is a key figure of merit limiting applications. In this work, we present an approach for generating background-free attosecond electron pulse trains by sequential inelastic electron-light scattering. Harnessing quantum interference in the fractional Talbot effect, we suppress unwanted background density in electron compression by several orders of magnitude. Our results will greatly enhance applications of coherent electron-light scattering, such as stimulated cathodoluminescence and streaking.


title: Landscape analysis for shallow neural networks: complete classification of critical points for affine target functions

abstract: In this paper, we analyze the landscape of the true loss of neural networks with one hidden layer and ReLU, leaky ReLU, or quadratic activation. In all three cases, we provide a complete classification of the critical points in the case where the target function is affine and one-dimensional. In particular, we show that there exist no local maxima and clarify the structure of saddle points. Moreover, we prove that non-global local minima can only be caused by `dead' ReLU neurons. In particular, they do not appear in the case of leaky ReLU or quadratic activation. Our approach is of a combinatorial nature and builds on a careful analysis of the different types of hidden neurons that can occur.


title: Thermodynamics in topological Josephson junctions

abstract: We study the thermodynamic properties of topological Josephson junctions using a quantum spin Hall (QSH) insulator-based junction as an example. In particular, we propose that phase-dependent measurements of the heat capacity offer an alternative to Josephson-current measurements to demonstrate key topological features. Even in an equilibrium situation, where the fermion parity is not conserved, the heat capacity exhibits a pronounced double peak in its phase dependence as a signature of the protected zero-energy crossing in the Andreev spectrum. This double-peak feature is robust against changes of the tunneling barrier and thus allows one to distinguish between topological and trivial junctions. At short time scales fermion parity is conserved and the heat capacity is $4\pi$-periodic in the superconducting phase difference. We propose a dispersive setup coupling the Josephson junction to a tank LC circuit to measure the heat capacity of the QSH-based Josephson junction sufficiently fast to detect the $4\pi$-periodicity. Although explicitly calculated for a short QSH-based Josephson junction, our results are also applicable to long as well as nanowire-based topological Josephson junctions.


title: Probing the topology of the two-photon bands via time-dependent quantum walks

abstract: Topological protection of quantum correlations opens new horizons and opportunities in quantum technologies. A variety of topological effects has recently been observed in qubit networks. However, the experimental identification of the topological phase still remains challenging, especially in the entangled many-body case. Here, we propose an approach to independently probe single- and two-photon topological invariants from the time evolution of the two-photon state in a one-dimensional array of qubits. Extending the bulk-boundary correspondence to the two-photon scenario, we show that an appropriate choice of the initial state enables the retrieval of the topological invariant for the different types of the two-photon states in the interacting Su-Schrieffer-Heeger model. Our analysis of the Zak phase reveals additional facets of topological protection in the case of collapse of bound photon pairs.


title: Functional portfolio optimization in stochastic portfolio theory

abstract: In this paper we develop a concrete and fully implementable approach to the optimization of functionally generated portfolios in stochastic portfolio theory. The main idea is to optimize over a family of rank-based portfolios parameterized by an exponentially concave function on the unit interval. This choice can be motivated by the long term stability of the capital distribution observed in large equity markets, and allows us to circumvent the curse of dimensionality. The resulting optimization problem, which is convex, allows for various regularizations and constraints to be imposed on the generating function. We prove an existence and uniqueness result for our optimization problem and provide a stability estimate in terms of a Wasserstein metric of the input measure. Then, we formulate a discretization which can be implemented numerically using available software packages and analyze its approximation error. Finally, we present empirical examples using CRSP data from the US stock market, including the performance of the portfolios allowing for dividends, defaults, and transaction costs.


title: Photon-emitter dressed states in a closed waveguide

abstract: We study a system made up of one or two two-level quantum emitters, coupled to a single transverse mode of a closed waveguide, in which photon wavenumbers and frequencies are discretized, and characterize the stable states in which one excitation is steadily shared between the field and the emitters. We unearth finite-size effects in the field-emitter interactions and identify a family of dressed states, that represent the forerunners of bound states in the continuum in the limit of an infinite waveguide. We finally consider the potential interest of such states for applications in the field of quantum information.


title: Exact edge landing probability for the bouncing coin toss and the three-sided die problem

abstract: Have you ever taken a disputed decision by tossing a coin and checking its landing side? This ancestral "heads or tails" practice is still widely used when facing undecided alternatives since it relies on the intuitive fairness of equiprobability. However, it critically disregards an interesting third outcome: the possibility of the coin coming at rest on its edge. Provided this evident yet elusive possibility, previous works have already focused on capturing all three landing probabilities of thick coins, but have only succeeded computationally as an exact and fully grounded analytical solution still evades the physics community. In this work we combine the classical equations of collisions with a statistical-mechanics approach to derive the exact analytical outcome probabilities of the partially inelastic coin toss. We validate the theoretical prediction by comparing it to previously reported simulations and experimental data; we discuss the moderate discrepancies arising at the highly inelastic regime; we describe the differences with previous, approximate models; and we finally propose the optimal geometry for the fair, cylindrical three-sided die based on the current results.


title: BERTSurv: BERT-Based Survival Models for Predicting Outcomes of Trauma Patients

abstract: Survival analysis is a technique to predict the times of specific outcomes, and is widely used in predicting the outcomes for intensive care unit (ICU) trauma patients. Recently, deep learning models have drawn increasing attention in healthcare. However, there is a lack of deep learning methods that can model the relationship between measurements, clinical notes and mortality outcomes. In this paper we introduce BERTSurv, a deep learning survival framework which applies Bidirectional Encoder Representations from Transformers (BERT) as a language representation model on unstructured clinical notes, for mortality prediction and survival analysis. We also incorporate clinical measurements in BERTSurv. With binary cross-entropy (BCE) loss, BERTSurv can predict mortality as a binary outcome (mortality prediction). With partial log-likelihood (PLL) loss, BERTSurv predicts the probability of mortality as a time-to-event outcome (survival analysis). We apply BERTSurv on Medical Information Mart for Intensive Care III (MIMIC III) trauma patient data. For mortality prediction, BERTSurv obtained an area under the curve of receiver operating characteristic curve (AUC-ROC) of 0.86, which is an improvement of 3.6% over baseline of multilayer perceptron (MLP) without notes. For survival analysis, BERTSurv achieved a concordance index (C-index) of 0.7. In addition, visualizations of BERT's attention heads help to extract patterns in clinical notes and improve model interpretability by showing how the model assigns weights to different inputs.


title: Empirical Analysis of Machine Learning Configurations for Prediction of Multiple Organ Failure in Trauma Patients

abstract: Multiple organ failure (MOF) is a life-threatening condition. Due to its urgency and high mortality rate, early detection is critical for clinicians to provide appropriate treatment. In this paper, we perform quantitative analysis on early MOF prediction with comprehensive machine learning (ML) configurations, including data preprocessing (missing value treatment, label balancing, feature scaling), feature selection, classifier choice, and hyperparameter tuning. Results show that classifier choice impacts both the performance improvement and variation most among all the configurations. In general, complex classifiers including ensemble methods can provide better performance than simple classifiers. However, blindly pursuing complex classifiers is unwise as it also brings the risk of greater performance variation.


title: Prediction of Hydraulic Blockage at Cross Drainage Structures using Regression Analysis

abstract: Hydraulic blockage of cross-drainage structures such as culverts is considered one of main contributor in triggering urban flash floods. However, due to lack of during floods data and highly non-linear nature of debris interaction, conventional modelling for hydraulic blockage is not possible. This paper proposes to use machine learning regression analysis for the prediction of hydraulic blockage. Relevant data has been collected by performing a scaled in-lab study and replicating different blockage scenarios. From the regression analysis, Artificial Neural Network (ANN) was reported best in hydraulic blockage prediction with $R^2$ of 0.89. With deployment of hydraulic sensors in smart cities, and availability of Big Data, regression analysis may prove helpful in addressing the blockage detection problem which is difficult to counter using conventional experimental and hydrological approaches.


title: Reduced order modeling of dynamical systems using artificial neural networks applied to water circulation

abstract: General circulation models are essential tools in weather and hydrodynamic simulation. They solve discretized, complex physical equations in order to compute evolutionary states of dynamical systems, such as the hydrodynamics of a lake. However, high-resolution numerical solutions using such models are extremely computational and time consuming, often requiring a high performance computing architecture to be executed satisfactorily. Machine learning (ML)-based low-dimensional surrogate models are a promising alternative to speed up these simulations without undermining the quality of predictions. In this work, we develop two examples of fast, reliable, low-dimensional surrogate models to produce a 36 hour forecast of the depth-averaged hydrodynamics at Lake George NY, USA. Our ML approach uses two widespread artificial neural network (ANN) architectures: fully connected neural networks and long short-term memory. These ANN architectures are first validated in the deterministic and chaotic regimes of the Lorenz system and then combined with proper orthogonal decomposition (to reduce the dimensionality of the incoming input data) to emulate the depth-averaged hydrodynamics of a flow simulator called SUNTANS. Results show the ANN-based reduced order models have promising accuracy levels (within 6% of the prediction range) and advocate for further investigation into hydrodynamic applications.


title: Multi-Time-Scale Input Approaches for Hourly-Scale Rainfall-Runoff Modeling based on Recurrent Neural Networks

abstract: This study proposes two straightforward yet effective approaches to reduce the required computational time of the training process for time-series modeling through a recurrent neural network (RNN) using multi-time-scale time-series data as input. One approach provides coarse and fine temporal resolutions of the input time-series to RNN in parallel. The other concatenates the coarse and fine temporal resolutions of the input time-series data over time before considering them as the input to RNN. In both approaches, first, finer temporal resolution data are utilized to learn the fine temporal scale behavior of the target data. Next, coarser temporal resolution data are expected to capture long-duration dependencies between the input and target variables. The proposed approaches were implemented for hourly rainfall-runoff modeling at a snow-dominated watershed by employing a long and short-term memory (LSTM) network, which is a newer type of RNN. Subsequently, the daily and hourly meteorological data were utilized as the input, and hourly flow discharge was considered as the target data. The results confirm that both of the proposed approaches can reduce the computational time for the training of RNN significantly (up to 32.4 times). Furthermore, one of the proposed approaches improves the estimation accuracy.


title: Hypothesis Tests on Rayleigh Wave Radiation Pattern Shapes: A Theoretical Assessment of Idealized Source Screening

abstract: Shallow seismic sources excite Rayleigh wave ground motion with azimuthally dependent radiation patterns. We place binary hypothesis tests on theoretical models of such radiation patterns to screen cylindrically symmetric sources (like explosions) from non-symmetric sources (like non-vertical dip-slip, or non-VDS faults). These models for data include sources with several unknown parameters, contaminated by Gaussian noise and embedded in a layered half-space. The generalized maximum likelihood ratio tests that we derive from these data models produce screening statistics and decision rules that depend on measured, noisy ground motion at discrete sensor locations. We explicitly quantify how the screening power of these statistics increase with the size of any dip-slip and strike-slip components of the source, relative to noise (faulting signal strength), and how they vary with network geometry. As applications of our theory, we apply these tests to (1) find optimal sensor locations that maximize the probability of screening non-circular radiation patterns, and (2) invert for the largest non-VDS faulting signal that could be mistakenly attributed to an explosion with damage, at a particular attribution probability. Lastly, we quantify how certain errors that are sourced by opening cracks increase screening rate errors. While such theoretical solutions are ideal and require future validation, they remain important in underground explosion monitoring scenarios because they provide fundamental physical limits on the discrimination power of tests that screen explosive from non-VDS faulting sources.


title: Sensitivity and robustness of Lagrangian coherent structures in coastal water systems

abstract: In coastal water systems, horizontal chaotic dispersion plays a significant role in the distribution and fate of pollutants. Lagrangian Coherent Structures (LCSs) provide useful tools to approach the problem of the transport of pollutants and have only recently been applied to coastal waters. While the fundamentals of the LCS approach using idealised analytical flow fields are well established in the literature, there are limited studies on their practical implementations in coastal waters where effects of boundaries and bathymetry frequently become significant. Due to their complex bathymetry and boundaries, unstructured grid systems are commonly used in modelling of coastal waters. For convenient derivation of LCS diagnostics, structured grids are commonly used. Here we examine the effect of mesh resolution, interpolation schemes and additive random noise on the LCS diagnostics in relation to coastal waters. Two kinematic model flows, the double gyre and the meandering jet, as well as validated outputs of a hydrodynamic model of Moreton Bay, Australia, on unstructured grids are used. The results show that LCSs are quite robust to the errors from interpolation schemes used in the data conversion from unstructured to structured grid. Attributed to the divergence in the underlying flow field, the results show that random errors in the order of 1-10 % cause a break down in the continuity of ridges of maximum finite time Lyapunov exponents and closed orbit elliptic LCSs. The result has significant implications on the suitability of applying LCS formulations based on a deterministic flow field to diffusive coastal waters.


title: Data-driven geophysical forecasting: Simple, low-cost, and accurate baselines with kernel methods

abstract: Modeling geophysical processes as low-dimensional dynamical systems and regressing their vector field from data is a promising approach for learning emulators of such systems. We show that when the kernel of these emulators is also learned from data (using kernel flows, a variant of cross-validation), then the resulting data-driven models are not only faster than equation-based models but are easier to train than neural networks such as the long short-term memory neural network. In addition, they are also more accurate and predictive than the latter. When trained on geophysical observational data, for example, the weekly averaged global sea-surface temperature, considerable gains are also observed by the proposed technique in comparison to classical partial differential equation-based models in terms of forecast computational cost and accuracy. When trained on publicly available re-analysis data for the daily temperature of the North-American continent, we see significant improvements over classical baselines such as climatology and persistence-based forecast techniques. Although our experiments concern specific examples, the proposed approach is general, and our results support the viability of kernel methods (with learned kernels) for interpretable and computationally efficient geophysical forecasting for a large diversity of processes.

categories: physics.ao-ph math.DS physics.flu-dyn physics.geo-ph stat.ML


title: Forecasting of Meteorological variables using statistical methods and tools

abstract: The need to understand the role of statistical methods for the forecasting of climatological parameters cannot be trivialized. This study gives an in depth review on the different variations of the Mann-Kendall (M-K) trend test and how they can be applied, regression techniques (Simple and Multiple), the Angstrom-Prescott model for solar radiation, etc. The study then goes ahead to apply some of them with data obtained from the Nigerian Meteorological Agency (NiMet), and applying tools like the python programming language and Wolfram Mathematica. Results show that the maximum ambient temperature for Calabar is increasing (Z=2.52) significantly after the calculated p-value < 0.05 (significant level). The seasonal M-K test was also applied for the dry and wet seasons and both were found to be increasing (Z=3.23 and Z=4.04 respectively) after their calculated p-values < 0.05. The relationship between refractivity and other meteorological parameters relating to it was discerned using partial differential equations giving the gradient of each with refractivity; this was compared with results from the correlation matrix to show that the water vapour contents of the atmosphere contributes significantly to the variation of refractivity. Multiple linear regression has also been adopted to give an accurate model for the prediction of refractivity in the region after the residual error between the calculated refractivity and predicted refractivity was minimal.


title: Locational Marginal Pricing: Towards a Free Market in Power

abstract: Nothing has done more to empower the free market, enterprise, and meritocracy than the spread of electricity and power to everyone. The power system has been the precursor to the greatest period of innovation in our history and has meant that visionaries with revolutionary ideas can compete with those with capital, political power, and means. Electricity, therefore, has been the great equalising force of the last 150 years, enhancing the productivity of the masses and granting prosperity to whole swathes of our nation. Whilst electricity has been one of the single largest innovations in enhancing the power of free markets, it is somewhat ironic that the way power is sold to consumers is largely unfree. The market is highly regulated, centralised, and is often used for political football by cynical politicians on both sides of the political spectrum. Introducing Locational Marginal Pricing into the UK grid system will increase economic freedom in the consumer markets for power, reduce prices for the poorest in the UK, decrease transmission losses, increase the permeation of low carbon generation in the grid, and incentivise investment in the UK's Northern Powerhouse initiative.


title: Quantum propensity in economics

abstract: This paper describes an approach to economics that is inspired by quantum computing, and is motivated by the need to develop a consistent quantum mathematical framework for economics. The traditional neoclassical approach assumes that rational utility-optimisers drive market prices to a stable equilibrium, subject to external perturbations. While this approach has been highly influential, it has come under increasing criticism following the financial crisis of 2007/8. The quantum approach, in contrast, is inherently probabilistic and dynamic. Decision-makers are described, not by a utility function, but by a propensity function which specifies the probability of transacting. We show how a number of cognitive phenomena such as preference reversal and the disjunction effect can be modelled by using a simple quantum circuit to generate an appropriate propensity function. Conversely, a general propensity function can be quantized to incorporate effects such as interference and entanglement that characterise human decision-making. Applications to some common problems in economics and finance are discussed.


title: The Landslide Velocity

abstract: Proper knowledge of velocity is required in accurately determining enormous destructive energy of a landslide. We present the first physics-based general analytical landslide velocity model that incorporates internal deformation and external forces: net driving force and viscous resistant. The model stands as a novel non-linear advective-dissipative system where classical Voellmy and inviscid Burgers' equation are specifications. Non-linear advection and external forcing fundamentally regulate the motion which substantially enhances our understanding of a coherently deforming landslide. Since analytical solutions provide fastest, cost-effective and best rigorous answer to the problem, we construct several new/general exact analytical solutions covering wider spectrum of landslide velocity. New solutions bridge existing gap between negligibly and massively deforming landslides. This provides a novel, rapid and consistent method for efficient coupling of different types of mass transports. Mechanism of landslide advection, stretching and approaching to steady-state has been explained. Shifting, up-lifting and stretching of velocity field stem from forcing and advection. Our solution describes fascinating breaking wave and emergence of landslide folding. This happens as the solution simultaneously introduces domain propagation, velocity up-lift and non-linear advection. Domain translation and stretching depends on net driving force and viscous drag controls shock wave generation, wave breaking, folding and velocity. Landslide dynamics are architectured by advection and reigned by system forcing. Analytically obtained velocities are close to observed values, constituting a new foundation of landslide velocity. This provides the practitioners with key information in instantly/accurately estimating impact force that is important in delineating hazard zones and mitigation of landslide hazards.


title: Fingerprints of the Cosmic Ray Driven Mechanism of the Ozone Hole

abstract: There is long research interest in electron-induced reactions of halogenated molecules. It has been two decades since the cosmic-ray (CR) driven electron-induced reaction (CRE) mechanism for the ozone hole formation was proposed. The derived CRE equation with stratospheric equivalent chlorine level and CR intensity as only two variables has well reproduced the observed data of stratospheric O3 and temperatures over the past 40 years. The CRE predictions of 11-year cyclic variations of the Antarctic O3 hole and associated stratospheric cooling have also been well confirmed. Measured altitude profiles of ozone and temperatures in Antarctic ozone holes provide convincing fingerprints of the CRE mechanism. A quantitative estimate indicates that the CRE-produced Cl atoms could completely deplete or even over-kill ozone in the CR-peak polar stratospheric region, consistent with observed altitude profiles of severest Antarctic ozone holes. After removing the natural CR effect, the hidden recovery in the Antarctic O3 hole since around 1995 is clearly discovered, while the recovery of O3 loss at mid-latitudes is being delayed by >=10 years. These results have provided strong evidence of the CRE mechanism. If the CR intensity keeps the current rising trend, the Antarctic O3 hole will return to the 1980 level by around 2060, while the returning of the O3 layer at mid-latitudes to the 1980 level will largely be delayed or will not even occur by the end of this century. The results strongly indicate that the CRE mechanism must be considered as a key factor in evaluating the O3 hole.


title: Portuguese eyewitness accounts of the great space weather event of 1582

abstract: Newly discovered descriptions about the great aurora observed in March 1582 are presented in this work. These records were made by Portuguese observers from Lisbon. Both records described the aurora like a great fire in the northern part of the sky. It was observed during three consecutive nights, according to one of the sources. Thus, we present a discussion of these auroral records in order to complement other works that studied the aurora sighted in March 1582.


title: A Survey of Stability Results for Redundancy Systems

abstract: Redundancy mechanisms consist in sending several copies of a same job to a subset of servers. It constitutes one of the most promising ways to exploit diversity in multiservers applications. However, its pros and cons are still not sufficiently understood in the context of realistic models with generic statistical properties of service-times distributions and correlation structures of copies. We aim at giving a survey of recent results concerning the stability-arguably the first benchmark of performance-of systems with cancel-oncompletion redundancy. We also point out open questions and conjectures.


title: NEO: Non Equilibrium Sampling on the Orbit of a Deterministic Transform

abstract: Sampling from a complex distribution $\pi$ and approximating its intractable normalizing constant Z are challenging problems. In this paper, a novel family of importance samplers (IS) and Markov chain Monte Carlo (MCMC) samplers is derived. Given an invertible map T, these schemes combine (with weights) elements from the forward and backward Orbits through points sampled from a proposal distribution $\rho$. The map T does not leave the target $\pi$ invariant, hence the name NEO, standing for Non-Equilibrium Orbits. NEO-IS provides unbiased estimators of the normalizing constant and self-normalized IS estimators of expectations under $\pi$ while NEO-MCMC combines multiple NEO-IS estimates of the normalizing constant and an iterated sampling-importance resampling mechanism to sample from $\pi$. For T chosen as a discrete-time integrator of a conformal Hamiltonian system, NEO-IS achieves state-of-the art performance on difficult benchmarks and NEO-MCMC is able to explore highly multimodal targets. Additionally, we provide detailed theoretical results for both methods. In particular, we show that NEO-MCMC is uniformly geometrically ergodic and establish explicit mixing time estimates under mild conditions.


title: Emergence of Structural Inequalities in Scientific Citation Networks

abstract: Structural inequalities persist in society, conferring systematic advantages to some people at the expense of others, for example, by giving them substantially more influence and opportunities. Using bibliometric data about authors of scientific publications, we identify two types of structural inequalities in scientific citations. First, female authors, who represent a minority of researchers, receive less recognition for their work (through citations) relative to male authors; second, authors affiliated with top-ranked institutions, who are also a minority, receive substantially more recognition compared to other authors. We present a model for the growth of directed citation networks and show that citations disparities arise from individual preferences to cite authors from the same group (homophily), highly cited or active authors (preferential attachment), as well as the size of the group and how frequently new authors join. We analyze the model and show that its predictions align well with real-world observations. Our theoretical and empirical analysis also suggests potential strategies to mitigate structural inequalities in science. In particular, we find that merely increasing the minority group size does little to narrow the disparities. Instead, reducing the homophily of each group, frequently adding new authors to a research field while providing them an accessible platform among existing, established authors, together with balanced group sizes can have the largest impact on reducing inequality. Our work highlights additional complexities of mitigating structural disparities stemming from asymmetric relations (e.g., directed citations) compared to symmetric relations (e.g., collaborations).

categories: physics.soc-ph cs.CY cs.DL cs.SI stat.AP


title: Relating boundary and interior solutions of the cohomological equation for cocycles by isometries of negatively curved spaces. The Livsic case

abstract: We consider the reducibility problem of cocycles by isometries of Gromov hyperbolic metric spaces in the Livsic setting. We show that provided that the boundary cocycle (that acts on a compact space) is reducible in a suitable H\"older class, then the original cocycle by isometries (that acts on an unbounded space) is also reducible.


title: From many-body quantum dynamics to the Hartree-Fock and Vlasov equations with singular potentials

abstract: We obtain the combined mean-field and semiclassical limit from the $N$-body Schr\"{o}dinger equation for fermions interacting via singular potentials. To obtain the result, we first prove the uniformity in Planck's constant $h$ propagation of regularity for solutions to the Hartree$\unicode{x2013}$Fock equation with singular pair interaction potentials of the form $\pm |x-y|^{-a}$, including the Coulomb and gravitational interactions. In the context of mixed states, we use these regularity properties to obtain quantitative estimates on the distance between solutions to the Schr\"{o}dinger equation and solutions to the Hartree$\unicode{x2013}$Fock and Vlasov equations in Schatten norms. For $a\in(0,1/2)$, we obtain local-in-time results when $N^{-1/2} \ll h \leq N^{-1/3}$. In particular, it leads to the derivation of the Vlasov equation with singular potentials. For $a\in[1/2,1]$, our results hold only on a small time scale, or with an $N$-dependent cutoff.


title: DNA Supercoiling Drives a Transition between Collective Modes of Gene Synthesis

abstract: Recent experiments showed that multiple copies of the molecular machine RNA polymerase (RNAP) can efficiently synthesize mRNA collectively in the active state of the promoter. However, environmentally-induced promoter repression results in long-distance antagonistic interactions that drastically reduce the speed of RNAPs and cause a quick arrest of mRNA synthesis. The mechanism underlying this transition between cooperative and antagonistic dynamics remains poorly understood. In this Letter, we introduce a continuum deterministic model for the translocation of RNAPs, where the speed of an RNAP is coupled to the local DNA supercoiling as well as the density of RNAPs on the gene. We assume that torsional stress experienced by individual RNAPs is exacerbated by high RNAP density on the gene and that transcription factors act as physical barriers to the diffusion of DNA supercoils. We show that this minimal model exhibits two transcription modes mediated by the torsional stress: a fluid mode when the promoter is active and a torsionally stressed mode when the promoter is repressed, in quantitative agreement with experimentally observed dynamics of co-transcribing RNAPs. Our work provides an important step towards understanding the collective dynamics of molecular machines involved in gene expression.


title: The Shape of Learning Curves: a Review

abstract: Learning curves provide insight into the dependence of a learner's generalization performance on the training set size. This important tool can be used for model selection, to predict the effect of more training data, and to reduce the computational complexity of model training and hyperparameter tuning. This review recounts the origins of the term, provides a formal definition of the learning curve, and briefly covers basics such as its estimation. Our main contribution is a comprehensive overview of the literature regarding the shape of learning curves. We discuss empirical and theoretical evidence that supports well-behaved curves that often have the shape of a power law or an exponential. We consider the learning curves of Gaussian processes, the complex shapes they can display, and the factors influencing them. We draw specific attention to examples of learning curves that are ill-behaved, showing worse learning performance with more training data. To wrap up, we point out various open problems that warrant deeper empirical and theoretical investigation. All in all, our review underscores that learning curves are surprisingly diverse and no universal model can be identified.


title: Refined Least Squares for Support Recovery

abstract: We study the problem of exact support recovery based on noisy observations and present Refined Least Squares (RLS). Given a set of noisy measurement $$ \myvec{y} = \myvec{X}\myvec{\theta}^* + \myvec{\omega},$$ and $\myvec{X} \in \mathbb{R}^{N \times D}$ which is a (known) Gaussian matrix and $\myvec{\omega} \in \mathbb{R}^N$ is an (unknown) Gaussian noise vector, our goal is to recover the support of the (unknown) sparse vector $\myvec{\theta}^* \in \left\{-1,0,1\right\}^D$. To recover the support of the $\myvec{\theta}^*$ we use an average of multiple least squares solutions, each computed based on a subset of the full set of equations. The support is estimated by identifying the most significant coefficients of the average least squares solution. We demonstrate that in a wide variety of settings our method outperforms state-of-the-art support recovery algorithms.


title: Longitudinal dynamics and chiral symmetry breaking in holographic light-front QCD

abstract: The breaking of chiral symmetry in holographic light-front QCD is encoded in its longitudinal dynamics with its chiral limit protected by the superconformal algebraic structure which governs its transverse dynamics. The scale in the longitudinal light-front Hamiltonian determines the confinement strength in this direction: It is also responsible for most of the light meson ground state mass consistent with the Gell-Mann-Oakes-Renner constraint. Longitudinal confinement and the breaking of chiral symmetry are found to be different manifestations of the same underlying dynamics as found in 't Hooft large $N_C$ QCD(1 + 1) model.


title: Paint by Word

abstract: We investigate the problem of zero-shot semantic image painting. Instead of painting modifications into an image using only concrete colors or a finite set of semantic concepts, we ask how to create semantic paint based on open full-text descriptions: our goal is to be able to point to a location in a synthesized image and apply an arbitrary new concept such as "rustic" or "opulent" or "happy dog." To do this, our method combines a state-of-the art generative model of realistic images with a state-of-the-art text-image semantic similarity network. We find that, to make large changes, it is important to use non-gradient methods to explore latent space, and it is important to relax the computations of the GAN to target changes to a specific region. We conduct user studies to compare our methods to several baselines.


title: Asymmetry underlies stability in power grids

abstract: Behavioral homogeneity is often critical for the functioning of network systems of interacting entities. In power grids, whose stable operation requires generator frequencies to be synchronized--and thus homogeneous--across the network, previous work suggests that the stability of synchronous states can be improved by making the generators homogeneous. Here, we show that a substantial additional improvement is possible by instead making the generators suitably heterogeneous. We develop a general method for attributing this counterintuitive effect to converse symmetry breaking, a recently established phenomenon in which the system must be asymmetric to maintain a stable symmetric state. These findings constitute the first demonstration of converse symmetry breaking in real-world systems, and our method promises to enable identification of this phenomenon in other networks whose functions rely on behavioral homogeneity.

categories: nlin.AO cond-mat.dis-nn cs.SY eess.SY physics.soc-ph


title: The topological modular forms of $\mathbb{R}P^2$ and $\mathbb{R}P^2 \wedge \mathbb{C}P^2$

abstract: In this paper, we study the elliptic spectral sequence computing $tmf_*(\mathbb{R} P^2)$ and $tmf_* (\mathbb{R} P^2 \wedge \mathbb{C} P^2)$. Specifically, we compute all differentials and resolve exotic extensions by 2, $\eta$, and $\nu$. For $tmf_* (\mathbb{R} P^2 \wedge \mathbb{C} P^2)$, we also compute the effect of the $v_1$-self maps of $\mathbb{R} P^2 \wedge \mathbb{C} P^2$ on $tmf$-homology.


title: From black hole to white hole via the intermediate static state

abstract: We discuss the macroscopic quantum tunneling from the black hole to the white hole of the same mass. Previous calculations in Ref.[1] demonstrated that the probability of the tunneling is $p \propto \exp(-2S_\text{BH})$, where $S_\text{BH}$ is the entropy of the Schwarzschild black hole. This in particular suggests that the entropy of the white hole is with minus sign the entropy of the black hole, $S_\text{WH}(M)=- S_\text{BH}(M)= - A/(4G)$. Here we use a different way of calculations. We consider three different types of the hole objects: black hole, white hole and the fully static intermediate state. The probability of tunneling transitions between these three states is found using singularities in the coordinate transformations between these objects. The black and white holes are described by the Painleve-Gullstrand coordinates with opposite shift vectors, while the intermediate state is described by the static Schwarzschild coordinates. The singularities in the coordinate transformations lead to the imaginary part in the action, which determines the tunneling exponent. For the white hole the negative entropy is obtained, while the intermediate state -- the fully static hole -- has zero entropy. This procedure is extended to the Reissner-Nordstr\"om black hole and to its white and static partners, and also to the entropy and temperature of the de Sitter Universe.