title: Bright Purcell enhanced single-photon source in the telecom O-band based on a quantum dot in a circular Bragg grating

abstract: The combination of semiconductor quantum dots (QDs) with photonic cavities is a promising way to realize non-classical light sources with state-of-the-art performances in terms of brightness, indistinguishability and repetition rate. In the present work we demonstrate the coupling of an InGaAs/GaAs QDs emitting in the telecom O-band to a circular Bragg grating cavity. We demonstrate a broadband geometric extraction efficiency enhancement by investigating two emission lines under above-band excitation, inside and detuned from the cavity mode, respectively. In the first case, a Purcell enhancement of 4 is attained. For the latter case, an end-to-end brightness of 1.4% with a brightness at the first lens of 23% is achieved. Using p-shell pumping, a combination of high count rate with pure single-photon emission (g(2)(0) = 0.01 in saturation) is achieved. Finally a good single-photon purity (g(2)(0) = 0.13) together with a high detector count rate of 191kcps is demonstrated for a temperature of up to 77K.


title: Probabilistic semi-nonnegative matrix factorization: a Skellam-based framework

abstract: We present a new probabilistic model to address semi-nonnegative matrix factorization (SNMF), called Skellam-SNMF. It is a hierarchical generative model consisting of prior components, Skellam-distributed hidden variables and observed data. Two inference algorithms are derived: Expectation-Maximization (EM) algorithm for maximum \emph{a posteriori} estimation and Variational Bayes EM (VBEM) for full Bayesian inference, including the estimation of parameters prior distribution. From this Skellam-based model, we also introduce a new divergence $\mathcal{D}$ between a real-valued target data $x$ and two nonnegative parameters $\lambda_{0}$ and $\lambda_{1}$ such that $\mathcal{D}\left(x\mid\lambda_{0},\lambda_{1}\right)=0\Leftrightarrow x=\lambda_{0}-\lambda_{1}$, which is a generalization of the Kullback-Leibler (KL) divergence. Finally, we conduct experimental studies on those new algorithms in order to understand their behavior and prove that they can outperform the classic SNMF approach on real data in a task of automatic clustering.


title: Climate Change Conspiracy Theories on Social Media

abstract: One of the critical emerging challenges in climate change communication is the prevalence of conspiracy theories. This paper discusses some of the major conspiracy theories related to climate change found in a large Twitter corpus. We use a state-of-the-art stance detection method to find whether conspiracy theories are more popular among Disbelievers or Believers of climate change. We then analyze which conspiracy theory is more popular than the others and how popularity changes with climate change belief. We find that Disbelievers of climate change are overwhelmingly responsible for sharing messages with conspiracy theory-related keywords, and not all conspiracy theories are equally shared. Lastly, we discuss the implications of our findings for climate change communication.


title: On the quantization of some polynomial minimal surfaces

abstract: A class of exact membrane solutions is quantized.


title: Interaction of Active Janus Particles with Passive Tracers

abstract: In this study, we investigated the motion of active SiO2-Pt Janus particles in the 2D bath of smaller silica tracers dispersed in varying areal densities. The effect on the organization of the tracer particles around the active JPs was also explored. Our experiments indicate that the interaction between the tracers and the active JPs mainly depend on the nature of collision marked by the duration of contact. For all the concentration regimes, we have shown that the short time collisions do not have significant impact on the motion of active JPs, however, during moderate/long-time collisions tracer(s) can lead to a significant change in active JPs motion and even cause them to rotate. In the concentrated regime, our experiments reveal the emergence of a novel organizational behavior of the passive tracers on the trailing Pt and the leading SiO2 with a strong dependence on the nature of collision.


title: Anomalous Hall antiferromagnets

abstract: The Hall effect, in which current flows perpendicular to an applied electrical bias, has played a prominent role in modern condensed matter physics over much of the subject's history. Appearing variously in classical, relativistic and quantum guises, it has among other roles contributed to establishing the band theory of solids, to research on new phases of strongly interacting electrons, and to the phenomenology of topological condensed matter. The dissipationless Hall current requires time-reversal symmetry breaking. For over a century it has either been ascribed to externally applied magnetic field and referred to as the ordinary Hall effect, or ascribed to spontaneous non-zero total internal magnetization (ferromagnetism) and referred to as the anomalous Hall effect. It has not commonly been associated with antiferromagnetic order. More recently, however, theoretical predictions and experimental observations have identified large Hall effects in some compensated magnetic crystals, governed by neither of the global magnetic-dipole symmetry breaking mechanisms mentioned above. The goals of this article are to systematically organize the present understanding of anomalous antiferromagnetic materials that generate a non-zero Hall effect, which we will call anomalous Hall antiferromagnets, and to discuss this class of materials in a broader fundamental and applied research context. Our motivation in drawing attention to anomalous Hall antiferromagnets is two-fold. First, since Hall effects that are not governed by magnetic dipole symmetry breaking are at odds with the traditional understanding of the phenomenon, the topic deserves attention on its own. Second, this new reincarnation has again placed the Hall effect in the middle of an emerging field of physics at the intersection of multipole magnetism, topological condensed matter, and spintronics.


title: An algorithmic view of $\ell_2$ regularization and some path-following algorithms

abstract: We establish an equivalence between the $\ell_2$-regularized solution path for a convex loss function, and the solution of an ordinary differentiable equation (ODE). Importantly, this equivalence reveals that the solution path can be viewed as the flow of a hybrid of gradient descent and Newton method applying to the empirical loss, which is similar to a widely used optimization technique called trust region method. This provides an interesting algorithmic view of $\ell_2$ regularization, and is in contrast to the conventional view that the $\ell_2$ regularization solution path is similar to the gradient flow of the empirical loss.New path-following algorithms based on homotopy methods and numerical ODE solvers are proposed to numerically approximate the solution path. In particular, we consider respectively Newton method and gradient descent method as the basis algorithm for the homotopy method, and establish their approximation error rates over the solution path. Importantly, our theory suggests novel schemes to choose grid points that guarantee an arbitrarily small suboptimality for the solution path. In terms of computational cost, we prove that in order to achieve an $\epsilon$-suboptimality for the entire solution path, the number of Newton steps required for the Newton method is $\mathcal O(\epsilon^{-1/2})$, while the number of gradient steps required for the gradient descent method is $\mathcal O\left(\epsilon^{-1} \ln(\epsilon^{-1})\right)$. Finally, we use $\ell_2$-regularized logistic regression as an illustrating example to demonstrate the effectiveness of the proposed path-following algorithms.


title: AGD-Autoencoder: Attention Gated Deep Convolutional Autoencoder for Brain Tumor Segmentation

abstract: Brain tumor segmentation is a challenging problem in medical image analysis. The endpoint is to generate the salient masks that accurately identify brain tumor regions in an fMRI screening. In this paper, we propose a novel attention gate (AG model) for brain tumor segmentation that utilizes both the edge detecting unit and the attention gated network to highlight and segment the salient regions from fMRI images. This feature enables us to eliminate the necessity of having to explicitly point towards the damaged area(external tissue localization) and classify(classification) as per classical computer vision techniques. AGs can easily be integrated within the deep convolutional neural networks(CNNs). Minimal computional overhead is required while the AGs increase the sensitivity scores significantly. We show that the edge detector along with an attention gated mechanism provide a sufficient enough method for brain segmentation reaching an IOU of 0.78


title: Enhancing an Intelligent Digital Twin with a Self-organized Reconfiguration Management based on Adaptive Process Models

abstract: Shorter product life cycles and increasing individualization of production leads to an increased reconfiguration demand in the domain of industrial automation systems, which will be dominated by cyber-physical production systems in the future. In constantly changing systems, however, not all configuration alternatives of the almost infinite state space are fully understood. Thus, certain configurations can lead to process instability, a reduction in quality or machine failures. Therefore, this paper presents an approach that enhances an intelligent Digital Twin with a self-organized reconfiguration management based on adaptive process models in order to find optimized configurations more comprehensively.


title: Robust Variable Selection and Estimation Via Adaptive Elastic Net S-Estimators for Linear Regression

abstract: Heavy-tailed error distributions and predictors with anomalous values are ubiquitous in high-dimensional regression problems and can seriously jeopardize the validity of statistical analyses if not properly addressed. For more reliable estimation under these adverse conditions, we propose a new robust regularized estimator for simultaneous variable selection and coefficient estimation. This estimator, called adaptive PENSE, possesses the oracle property without prior knowledge of the scale of the residuals and without any moment conditions on the error distribution. The proposed estimator gives reliable results even under very heavy-tailed error distributions and aberrant contamination in the predictors or residuals. Importantly, even in these challenging settings variable selection by adaptive PENSE remains stable. Numerical studies on simulated and real data sets highlight superior finite-sample performance in a vast range of settings compared to other robust regularized estimators in the case of contaminated samples and competitiveness compared to classical regularized estimators in clean samples.


title: Tate-Hochschild cohomology rings for eventually periodic Gorenstein algebras

abstract: Tate-Hochschild cohomology of an algebra is a generalization of ordinary Hochschild cohomology, which is defined on positive and negative degrees and has a ring structure. Our purpose of this paper is to study the eventual periodicity of an algebra by using the Tate-Hochschild cohomology ring. First, we deal with eventually periodic algebras and show that they are not necessarily Gorenstein algebras. Secondly, we characterize the eventual periodicity of a Gorenstein algebra as the existence of an invertible homogeneous element of the Tate-Hochschild cohomology ring of the algebra, which is our main result. Finally, we use tensor algebras to establish a way of constructing eventually periodic Gorenstein algebras.


title: Design an IT Policy Implementation Plan

abstract: Information technology (IT) companies implement multi-dimensional policy plans that include procedures, sub-plans, and instructions to outline their business scopes, targets, and communications. This work outlined the IT policy implementation plan designed by an imaginary company with a random name called Northcentral Cloud Consulting Firm (NCCF), containing proposed IT policies, milestones and roadmaps, control framework, stakeholder responsibilities, knowledge transfer plan, and leadership roles. As NCCF's major customers seek data-driven solutions in cloud computing, the NCCF IT policy plan provides various data policies, including security and proper usage of machine learning services. The plan offers a detailed roadmap of its financial, geographical, and reputational expansion within three years. The IT policy plan also compromises an IT risk management, contingency, and emergency communication plan, mainly for protecting data and business continuity. Stakeholder responsibilities are incorporated into the IT policy plan, as NCCF considers any engagement with its customers as a collaborative effort in which both parties have and share several responsibilities.


title: An Analysis of X-Ray Hardness Ratios Between Asynchronous and Non-Asynchronous Polars

abstract: The subclass of magnetic Cataclysmic Variables (CV), known as asynchronous polars, are still relatively poorly understood. An asynchronous polar is a polar in which the spin period of the white dwarf is either shorter or longer than the binary orbital period (typically within a few percent). The asynchronous polars have been disproportionately detected in soft gamma-ray observations, leading us to consider the possibility that they have intrinsically harder X-ray spectra. We compared standard and asynchronous polars in order to examine the relationship between a CV's synchronization status and its spectral shape. Using the entire sample of asynchronous polars, we find that the asynchronous polars may, indeed, have harder spectra, but that the result is not statistically significant.


title: Radioactive Planet Formation

abstract: Young stellar objects are observed to have large X-ray fluxes and are thought to produce commensurate luminosities in energetic particles (cosmic rays). This particle radiation, in turn, can synthesize short-lived radioactive nuclei through spallation. With a focus on $^{26}$Al, this paper estimates the expected abundances of radioactive nulcei produced by spallation during the epoch of planet formation. In this model, cosmic rays are accelerated near the inner truncation radii of circumstellar disks, $r_{\scriptstyle X}\approx0.1$ AU, where intense magnetic activity takes place. For planets forming in this region, radioactive abundances can be enhanced over the values inferred for the early solar system (from meteoritic measurements) by factors of $\sim10-20$. These short-lived radioactive nuclei influence the process of planet formation and the properties of planets in several ways. The minimum size required for planetesimals to become fully molten decreases with increasing levels of radioactive enrichment, and such melting leads to loss of volatile components including water. Planets produced with an enhanced radioactive inventory have significant internal luminosity which can be comparable to that provided by the host star; this additional heating affects both atmospheric mass loss and chemical composition. Finally, the habitable zone of red dwarf stars is coincident with the magnetic reconnection region, so that planets forming at those locations will experience maximum exposure to particle radiation, and subsequent depletion of volatiles.


title: Machine-learning approach to finite-size effects in systems with strongly interacting fermions

abstract: We investigate the applicability of machine learning techniques in studying the finite-size effects associated with many-body physics. These techniques have an emerging presence in many-body theory as they have been used for interpolations, extrapolations, and in modeling wavefunctions. We will resolve several issues associated with machine learning and many-body calculations such as small datasets, outliers, and discontinuities, for the purpose of extrapolating finite calculations to macroscopic scales. We carry out a systematic investigation of two related systems by developing metrics that aim to avoid spurious effects and capture desired features. This work uses neural networks to extrapolate the Unitary Gas to the thermodynamic limit at zero-range, which is otherwise difficult to reach. The effective mass of strongly interacting neutron matter is also studied and makes use of the non-interacting problem to resolve discontinuous predictions. For this investigation, we also carried out new Auxiliary Field Diffusion Monte Carlo (AFDMC) calculations for a variety of densities and particle numbers. Ultimately, we demonstrate an effective utility for neural networks in this context.


title: KOALA: A Kalman Optimization Algorithm with Loss Adaptivity

abstract: Optimization is often cast as a deterministic problem, where the solution is found through some iterative procedure such as gradient descent. However, when training neural networks the loss function changes over (iteration) time due to the randomized selection of a subset of the samples. This randomization turns the optimization problem into a stochastic one. We propose to consider the loss as a noisy observation with respect to some reference optimum. This interpretation of the loss allows us to adopt Kalman filtering as an optimizer, as its recursive formulation is designed to estimate unknown parameters from noisy measurements. Moreover, we show that the Kalman Filter dynamical model for the evolution of the unknown parameters can be used to capture the gradient dynamics of advanced methods such as Momentum and Adam. We call this stochastic optimization method KOALA, which is short for Kalman Optimization Algorithm with Loss Adaptivity. KOALA is an easy to implement, scalable, and efficient method to train neural networks. We provide convergence analysis and show experimentally that it yields parameter estimates that are on par with or better than existing state of the art optimization algorithms across several neural network architectures and machine learning tasks, such as computer vision and language modeling.


title: SimCC: a Simple Coordinate Classification Perspective for Human Pose Estimation

abstract: The 2D heatmap-based approaches have dominated Human Pose Estimation (HPE) for years due to high performance. However, the long-standing quantization error problem in the 2D heatmap-based methods leads to several well-known drawbacks: 1) The performance for the low-resolution inputs is limited; 2) To improve the feature map resolution for higher localization precision, multiple costly upsampling layers are required; 3) Extra post-processing is adopted to reduce the quantization error. To address these issues, we aim to explore a brand new scheme, called \textit{SimCC}, which reformulates HPE as two classification tasks for horizontal and vertical coordinates. The proposed SimCC uniformly divides each pixel into several bins, thus achieving \emph{sub-pixel} localization precision and low quantization error. Benefiting from that, SimCC can omit additional refinement post-processing and exclude upsampling layers under certain settings, resulting in a more simple and effective pipeline for HPE. Extensive experiments conducted over COCO, CrowdPose, and MPII datasets show that SimCC outperforms heatmap-based counterparts, especially in low-resolution settings by a large margin.


title: Learning quantum many-body systems from a few copies

abstract: Estimating physical properties of quantum states from measurements is one of the most fundamental tasks in quantum science. In this work, we identify conditions on states under which it is possible to infer the expectation values of all quasi-local observables of a state from a number of copies that scales polylogarithmically with the system's size and polynomially on the locality of the target observables. We show that this constitutes a provable exponential improvement in the number of copies over state-of-the-art tomography protocols. We achieve our results by combining the maximum entropy method with tools from the emerging fields of classical shadows and quantum optimal transport. The latter allows us to fine-tune the error made in estimating the expectation value of an observable in terms of how local it is and how well we approximate the expectation value of a fixed set of few-body observables. We conjecture that our condition holds for all states exhibiting some form of decay of correlations and establish it for several subsets thereof. These include widely studied classes of states such as one-dimensional thermal and high-temperature Gibbs states of local commuting Hamiltonians on arbitrary hypergraphs or outputs of shallow circuits. Moreover, we show improvements of the maximum entropy method beyond the sample complexity that are of independent interest. These include identifying regimes in which it is possible to perform the postprocessing efficiently as well as novel bounds on the condition number of covariance matrices of many-body states.


title: Predicting the diversity of early epidemic spread on networks

abstract: The interplay of biological, social, structural and random factors makes disease forecasting extraordinarily complex. The course of an epidemic exhibits average growth dynamics determined by features of the pathogen and the population, yet also features significant variability reflecting the stochastic nature of disease spread. In this work, we reframe a stochastic branching process analysis in terms of probability generating functions and compare it to continuous time epidemic simulations on networks. In doing so, we predict the diversity of emerging epidemic courses on both homogeneous and heterogeneous networks. We show how the challenge of inferring the early course of an epidemic falls on the randomness of disease spread more so than on the heterogeneity of contact patterns. We provide an analysis which helps quantify, in real time, the probability that an epidemic goes supercritical or conversely, dies stochastically. These probabilities are often assumed to be one and zero, respectively, if the basic reproduction number, or R0, is greater than 1, ignoring the heterogeneity and randomness inherent to disease spread. This framework can give more insight into early epidemic spread by weighting standard deterministic models with likelihood to inform pandemic preparedness with probabilistic forecasts.


title: Vortex clustering, polarisation and circulation intermittency in classical and quantum turbulence

abstract: The understanding of turbulent flows is one of the biggest current challenges in physics, as no first-principles theory exists to explain their observed spatio-temporal intermittency. Turbulent flows may be regarded as an intricate collection of mutually-interacting vortices. This picture becomes accurate in quantum turbulence, which is built on tangles of discrete vortex filaments. Here, we study the statistics of velocity circulation in quantum and classical turbulence. We show that, in quantum flows, Kolmogorov turbulence emerges from the correlation of vortex orientations, while deviations -- associated with intermittency -- originate from their non-trivial spatial arrangement. We then link the spatial distribution of vortices in quantum turbulence to the coarse-grained energy dissipation in classical turbulence, enabling the application of existent models of classical turbulence intermittency to the quantum case. Our results provide a connection between the intermittency of quantum and classical turbulence and initiate a promising path to a better understanding of the latter.


title: Regularization-based Continual Learning for Fault Prediction in Lithium-Ion Batteries

abstract: In recent years, the use of lithium-ion batteries has greatly expanded into products from many industrial sectors, e.g. cars, power tools or medical devices. An early prediction and robust understanding of battery faults could therefore greatly increase product quality in those fields. While current approaches for data-driven fault prediction provide good results on the exact processes they were trained on, they often lack the ability to flexibly adapt to changes, e.g. in operational or environmental parameters. Continual learning promises such flexibility, allowing for an automatic adaption of previously learnt knowledge to new tasks. Therefore, this article discusses different continual learning approaches from the group of regularization strategies, which are implemented, evaluated and compared based on a real battery wear dataset. Online elastic weight consolidation delivers the best results, but, as with all examined approaches, its performance appears to be strongly dependent on task characteristics and task sequence.


title: Samplets: A new paradigm for data compression

abstract: In this article, we introduce the concept of samplets by transferring the construction of Tausch-White wavelets to the realm of data. This way we obtain a multilevel representation of discrete data which directly enables data compression, detection of singularities and adaptivity. Applying samplets to represent kernel matrices, as they arise in kernel based learning or Gaussian process regression, we end up with quasi-sparse matrices. By thresholding small entries, these matrices are compressible to O(N log N) relevant entries, where N is the number of data points. This feature allows for the use of fill-in reducing reorderings to obtain a sparse factorization of the compressed matrices. Besides the comprehensive introduction to samplets and their properties, we present extensive numerical studies to benchmark the approach. Our results demonstrate that samplets mark a considerable step in the direction of making large data sets accessible for analysis.


title: Discontinuous phase transitions in the q-voter model with generalized anticonformity on random graphs

abstract: We study the binary $q$-voter model with generalized anticonformity on random Erd\H{o}s-R\'enyi graphs. In such a model, two types of social responses, conformity and anticonformity, occur with complementary probabilities and the size of the source of influence $q_c$ in case of conformity is independent from the size of the source of influence $q_a$ in case of anticonformity. For $q_c=q_a=q$ the model reduces to the original $q$-voter model with anticonformity. Previously, such a generalized model was studied only on the complete graph, which corresponds to the mean-field approach. It was shown that it can display discontinuous phase transitions for $q_c \ge q_a + \Delta q$, where $\Delta q=4$ for $q_a \le 3$ and $\Delta q=3$ for $q_a>3$. In this paper, we pose the question if discontinuous phase transitions survive on random graphs with an average node degree $\langle k\rangle \le 150$ observed empirically in social networks. Using the pair approximation, as well as Monte Carlo simulations, we show that discontinuous phase transitions indeed can survive, even for relatively small values of $\langle k\rangle$. Moreover, we show that for $q_a < q_c - 1$ pair approximation results overlap the Monte Carlo ones. On the other hand, for $q_a \ge q_c - 1$ pair approximation gives qualitatively wrong results indicating discontinuous phase transitions neither observed in the simulations nor within the mean-field approach. Finally, we report an intriguing result showing that the difference between the spinodals obtained within the pair approximation and the mean-field approach follows a power law with respect to $\langle k\rangle$, as long as the pair approximation indicates correctly the type of the phase transition.


title: Optimization of Electrolyte Rebalancing in Vanadium Redox Flow Batteries

abstract: This paper presents a novel algorithm to optimize energy capacity restoration of vanadium redox flow batteries (VRFBs). VRFB technologies can have their lives prolonged through a partially restoration of the lost capacity by electrolyte rebalancing. Our algorithm finds the optimal number and time of these rebalancing services to minimize the service cost, while maximizing the revenues from energy arbitrage. We show that the linearized form of this problem can be analytically solved, and that the objective function is convex. To solve the complete problem, we develop a two-step mixed integer linear programming (MILP) algorithm, which first finds the bounds for optimal number of services and then optimizes the number, and time of the services. We then present a theoretical analysis and optimization results for a case study of energy arbitrage in New York ISO.


title: Pseudo-Model-Free Hedging for Variable Annuities via Deep Reinforcement Learning

abstract: This paper proposes a two-phase deep reinforcement learning approach, for hedging variable annuity contracts with both GMMB and GMDB riders, which can address model miscalibration in Black-Scholes financial and constant force of mortality actuarial market environments. In the training phase, an infant reinforcement learning agent interacts with a pre-designed training environment, collects sequential anchor-hedging reward signals, and gradually learns how to hedge the contracts. As expected, after a sufficient number of training steps, the trained reinforcement learning agent hedges, in the training environment, equally well as the correct Delta while outperforms misspecified Deltas. In the online learning phase, the trained reinforcement learning agent interacts with the market environment in real time, collects single terminal reward signals, and self-revises its hedging strategy. The hedging performance of the further trained reinforcement learning agent is demonstrated via an illustrative example on a rolling basis to reveal the self-revision capability on the hedging strategy by online learning.


title: Burrows Wheeler Transform on a Large Scale: Algorithms Implemented in Apache Spark

abstract: With the rapid growth of Next Generation Sequencing (NGS) technologies, large amounts of "omics" data are daily collected and need to be processed. Indexing and compressing large sequences datasets are some of the most important tasks in this context. Here we propose algorithms for the computation of Burrows Wheeler transform relying on Big Data technologies, i.e., Apache Spark and Hadoop. Our algorithms are the first ones that distribute the index computation and not only the input dataset, allowing to fully benefit of the available cloud resources.


title: A Survey of Uncertainty in Deep Neural Networks

abstract: Due to their increasing spread, confidence in neural network predictions became more and more important. However, basic neural networks do not deliver certainty estimates or suffer from over or under confidence. Many researchers have been working on understanding and quantifying uncertainty in a neural network's prediction. As a result, different types and sources of uncertainty have been identified and a variety of approaches to measure and quantify uncertainty in neural networks have been proposed. This work gives a comprehensive overview of uncertainty estimation in neural networks, reviews recent advances in the field, highlights current challenges, and identifies potential research opportunities. It is intended to give anyone interested in uncertainty estimation in neural networks a broad overview and introduction, without presupposing prior knowledge in this field. A comprehensive introduction to the most crucial sources of uncertainty is given and their separation into reducible model uncertainty and not reducible data uncertainty is presented. The modeling of these uncertainties based on deterministic neural networks, Bayesian neural networks, ensemble of neural networks, and test-time data augmentation approaches is introduced and different branches of these fields as well as the latest developments are discussed. For a practical application, we discuss different measures of uncertainty, approaches for the calibration of neural networks and give an overview of existing baselines and implementations. Different examples from the wide spectrum of challenges in different fields give an idea of the needs and challenges regarding uncertainties in practical applications. Additionally, the practical limitations of current methods for mission- and safety-critical real world applications are discussed and an outlook on the next steps towards a broader usage of such methods is given.


title: Splitting of the superfluid transition in $^3$He by nematic aerogels

abstract: Role of inhomogeneous perturbations of nematic aerogel on the form of the order parameter emerging at the transition of liquid $^3$He in the superfluid state is considered. It is shown that a region of stability of the polar distorted ABM phase can begin right from the transition temperature. The symmetry argument is given, which selects the most favorable aerogel for stabilization of pure polar phase.


title: Time Encoding of Finite-Rate-of-Innovation Signals

abstract: Time-encoding of continuous-time signals is an alternative sampling paradigm to conventional methods such as Shannon's sampling. In time-encoding, the signal is encoded using a sequence of time instants where an event occurs, and hence fall under event-driven sampling methods. Time-encoding can be designed agnostic to the global clock of the sampling hardware, which makes sampling asynchronous. Moreover, the encoding is sparse. This makes time-encoding energy efficient. However, the signal representation is nonstandard and in general, nonuniform. In this paper, we consider time-encoding of finite-rate-of-innovation signals, and in particular, periodic signals composed of weighted and time-shifted versions of a known pulse. We consider encoding using both crossing-time-encoding machine (C-TEM) and integrate-and-fire time-encoding machine (IF-TEM). We analyze how time-encoding manifests in the Fourier domain and arrive at the familiar sum-of-sinusoids structure of the Fourier coefficients that can be obtained starting from the time-encoded measurements via a suitable linear transformation. Thereafter, standard FRI techniques become applicable. Further, we extend the theory to multichannel time-encoding such that each channel operates with a lower sampling requirement. We also study the effect of measurement noise, where the temporal measurements are perturbed by additive noise. To combat the effect of noise, we propose a robust optimization framework to simultaneously denoise the Fourier coefficients and estimate the annihilating filter accurately. We provide sufficient conditions for time-encoding and perfect reconstruction using C-TEM and IF-TEM, and furnish extensive simulations to substantiate our findings.


title: Post-Metering Value-Added Services for Low Voltage Electricity Users: Lessons Learned From the Italian Experience of CHAIN 2

abstract: Electrical energy consumption data accessibility for low voltage end users is one of the pillars of smart grids. In some countries, despite the presence of smart meters, a fragmentary data availability and/or the lack of standardization hinders the creation of post-metering value-added services and confines such innovative solutions to the prototypal and experimental level. We take inspiration from the technology adopted in Italy, where the national regulatory authority actively supported the definition of a solution agreed upon by all the involved stakeholders. In this context, smart meters are enabled to convey data to low voltage end users through a power line communication channel (CHAIN 2) in near real-time. The aim of this paper is twofold. On the one hand, it describes the proof of concept that the channel underwent and its subsequent validation (with performances nearing 99% success rate). On the other hand, it defines a classification framework (I2MA) for post-metering value-added services, in order to categorize each use case based on both level of service and expected benefits, and understand its maturity level. As an example, we apply the methodology to the 16 use cases defined in Italy. The lessons learned from the regulatory, technological, and functional approach of the Italian experience bring us to the provision of recommendations for researchers and industry experts. In particular, we argue that a well-functioning post-metering value-added services' market can flourish when: i) distribution system operators certify the measurements coming from smart meters; ii) national regulatory authorities support the technological innovation needed for setting up this market; and iii) service providers create customer-oriented solutions based on smart meters' data.


title: Using spin alignment of inelastically-excited fast beams to make spin assignments: the spectroscopy of 13O as a test case

abstract: Excited states in 13O were investigated using inelastic scattering of an E/A=69.5-MeV 13O beam off of a 9Be target. The excited states were identified in the invariant-mass spectra of the decay products. Both single proton and sequential two-proton decays of the excited states were examined. For a number of the excited states, the protons were emitted with strong anisotropy where emissions transverse to the beam axis are favored. The measured proton-decay angular distributions were compared to predictions from distorted-wave born-approximation (DWBA) calculations of the spin alignment which was shown to be largely independent of the excitation mechanism. The deduced $^{13}$O level scheme is compared to ab initio no-core shell model with continuum (NCSMC) predictions. The lowest-energy excited states decay isotropically consistent with predictions of strong proton 1s1/2 structure. Above these states in the level scheme, we observed a number of higher-spin states not predicted within the model. Possibly these are associated with rotational bands built on deformed cluster configurations predicted by antisymmetrized molecular dynamics (AMD) calculations. The spin alignment mechanism is shown to be useful for making spin assignments and may have widespread use.


title: Optimal Edge Weight Perturbations to Attack Shortest Paths

abstract: Finding shortest paths in a given network (e.g., a computer network or a road network) is a well-studied task with many applications. We consider this task under the presence of an adversary, who can manipulate the network by perturbing its edge weights to gain an advantage over others. Specifically, we introduce the Force Path Problem as follows. Given a network, the adversary's goal is to make a specific path the shortest by adding weights to edges in the network. The version of this problem in which the adversary can cut edges is NP-complete. However, we show that Force Path can be solved to within arbitrary numerical precision in polynomial time. We propose the PATHPERTURB algorithm, which uses constraint generation to build a set of constraints that require paths other than the adversary's target to be sufficiently long. Across a highly varied set of synthetic and real networks, we show that the optimal solution often reduces the required perturbation budget by about half when compared to a greedy baseline method.


title: Learning nonequilibrium control forces to characterize dynamical phase transitions

abstract: Sampling the collective, dynamical fluctuations that lead to nonequilibrium pattern formation requires probing rare regions of trajectory space. Recent approaches to this problem based on importance sampling, cloning, and spectral approximations, have yielded significant insight into nonequilibrium systems, but tend to scale poorly with the size of the system, especially near dynamical phase transitions. Here we propose a machine learning algorithm that samples rare trajectories and estimates the associated large deviation functions using a many-body control force by leveraging the flexible function representation provided by deep neural networks, importance sampling in trajectory space, and stochastic optimal control theory. We show that this approach scales to hundreds of interacting particles and remains robust at dynamical phase transitions.


title: Spitzer phase curve observations and circulation models of the inflated ultra-hot Jupiter WASP-76b

abstract: The large radii of many hot Jupiters can only be matched by models that have hot interior adiabats, and recent theoretical work has shown that the interior evolution of hot Jupiters has a significant impact on their atmospheric structure. Due to its inflated radius, low gravity, and ultra-hot equilibrium temperature, WASP-76b is an ideal case study for the impact of internal evolution on observable properties. Hot interiors should most strongly affect the non-irradiated side of the planet, and thus full phase curve observations are critical to ascertain the effect of the interior on the atmospheres of hot Jupiters. In this work, we present the first Spitzer phase curve observations of WASP-76b. We find that WASP-76b has an ultra-hot day side and relatively cold nightside with brightness temperatures of $2471 \pm 27~\mathrm{K}$/$1518 \pm 61~\mathrm{K}$ at $3.6~\micron$ and $2699 \pm 32~\mathrm{K}$/$1259 \pm 44~\mathrm{K}$ at $4.5~\micron$, respectively. These results provide evidence for a dayside thermal inversion. Both channels exhibit small phase offsets of $0.68 \pm 0.48^{\circ}$ at $3.6~\micron$ and $0.67 \pm 0.2^{\circ}$ at $4.5~\mu\mathrm{m}$. We compare our observations to a suite of general circulation models that consider two end-members of interior temperature along with a broad range of frictional drag strengths. Strong frictional drag is necessary to match the small phase offsets and cold nightside temperatures observed. From our suite of cloud-free GCMs, we find that only cases with a cold interior can reproduce the cold nightsides and large phase curve amplitude at $4.5~\micron$, hinting that the hot interior adiabat of WASP-76b does not significantly impact its atmospheric dynamics or that clouds blanket its nightside.


title: Ice inheritance in dynamical disk models

abstract: The compositions of planet-forming disks are set by a combination of material inherited from the interstellar medium and material reprocessed during disk formation and evolution. Indeed, comets and primitive meteorites exhibit interstellar-like isotopic ratios and/or volatile compositions, supporting that some pristine material was incorporated intact into icy planetesimals in the Solar Nebula. To date, the survival of volatile interstellar material in the disk stage has not been modeled using realistic disk physics. Here, we present a modeling framework to track the destruction of interstellar ices on dust grains undergoing transport processes within a disk, with a particular focus on explaining the incorporation of pristine material into icy planetesimals. We find it is difficult to explain inheritance through the local assembly of comets, as ice destruction is rapid for small (<10um) grains in the inner few tens of au. Instead, a plausible pathway to inheritance is to form pebbles at larger disk radii, which then drift inwards to the comet-forming zone with their ices mostly preserved. Small grains beyond ~100 au can experience ice photodissociation at the tens of percent level, however little of the ice is actually lost from the grain, likely making this a robust site for in situ ice chemistry. Our models also indicate that many complex organic species should survive passage through the disk intact. This raises the possibility that organics synthesized in the interstellar medium can be delivered to terrestrial planets by icy body impact and thus potentially participate in origins of life chemistry.


title: Quasi-positive curvature on Bazaikin spaces

abstract: We completely characterize the sectional curvature of all of the $13$-dimensional Bazaikin spaces. In particular, we show that all Bazaikin spaces admit a quasi-positively curved Riemannian metric, and that, up to isometry, there is a unique Bazaikin space which is almost positively curved but not positively curved.


title: IntraLoss: Further Margin via Gradient-Enhancing Term for Deep Face Recognition

abstract: Existing classification-based face recognition methods have achieved remarkable progress, introducing large margin into hypersphere manifold to learn discriminative facial representations. However, the feature distribution is ignored. Poor feature distribution will wipe out the performance improvement brought about by margin scheme. Recent studies focus on the unbalanced inter-class distribution and form a equidistributed feature representations by penalizing the angle between identity and its nearest neighbor. But the problem is more than that, we also found the anisotropy of intra-class distribution. In this paper, we propose the `gradient-enhancing term' that concentrates on the distribution characteristics within the class. This method, named IntraLoss, explicitly performs gradient enhancement in the anisotropic region so that the intra-class distribution continues to shrink, resulting in isotropic and more compact intra-class distribution and further margin between identities. The experimental results on LFW, YTF and CFP-FP show that our outperforms state-of-the-art methods by gradient enhancement, demonstrating the superiority of our method. In addition, our method has intuitive geometric interpretation and can be easily combined with existing methods to solve the previously ignored problems.


title: Tidal Evolution of the Earth-Moon System with a High Initial Obliquity

abstract: A giant impact origin for the Moon is generally accepted, but many aspects of lunar formation remain poorly understood and debated. \'Cuk et al. (2016) proposed that an impact that left the Earth-Moon system with high obliquity and angular momentum could explain the Moon's orbital inclination and isotopic similarity to Earth. In this scenario, instability during the Laplace Plane transition, when the Moon's orbit transitions from the gravitational influence of Earth's figure to that of the Sun, would both lower the system's angular momentum to its present-day value and generate the Moon's orbital inclination. Recently, Tian and Wisdom (2020) discovered new dynamical constraints on the Laplace Plane transition and concluded that the Earth-Moon system could not have evolved from an initial state with high obliquity. Here we demonstrate that the Earth-Moon system with an initially high obliquity can evolve into the present state, and we identify a spin-orbit secular resonance as a key dynamical mechanism in the later stages of the Laplace Plane transition. Some of the simulations by Tian and Wisdom (2020) did not encounter this late secular resonance, as their model suppressed obliquity tides and the resulting inclination damping. Our results demonstrate that a giant impact that left Earth with high angular momentum and high obliquity ($\theta > 61^{\circ}$) is a promising scenario for explaining many properties of the Earth-Moon system, including its angular momentum and obliquity, the geochemistry of Earth and the Moon, and the lunar inclination.


title: Mitigating Performance Saturation in Neural Marked Point Processes: Architectures and Loss Functions

abstract: Attributed event sequences are commonly encountered in practice. A recent research line focuses on incorporating neural networks with the statistical model -- marked point processes, which is the conventional tool for dealing with attributed event sequences. Neural marked point processes possess good interpretability of probabilistic models as well as the representational power of neural networks. However, we find that performance of neural marked point processes is not always increasing as the network architecture becomes more complicated and larger, which is what we call the performance saturation phenomenon. This is due to the fact that the generalization error of neural marked point processes is determined by both the network representational ability and the model specification at the same time. Therefore we can draw two major conclusions: first, simple network structures can perform no worse than complicated ones for some cases; second, using a proper probabilistic assumption is as equally, if not more, important as improving the complexity of the network. Based on this observation, we propose a simple graph-based network structure called GCHP, which utilizes only graph convolutional layers, thus it can be easily accelerated by the parallel mechanism. We directly consider the distribution of interarrival times instead of imposing a specific assumption on the conditional intensity function, and propose to use a likelihood ratio loss with a moment matching mechanism for optimization and model selection. Experimental results show that GCHP can significantly reduce training time and the likelihood ratio loss with interarrival time probability assumptions can greatly improve the model performance.


title: Quasi-isometry invariance of relative filling functions

abstract: For a finitely generated group $G$ and collection of subgroups $\mathcal{P}$ we prove that the relative Dehn function of a pair $(G,\mathcal{P})$ is invariant under quasi-isometry of pairs. Along the way we show quasi-isometries of pairs preserve almost malnormality of the collection and fineness of the associated coned off Cayley graphs. We also prove that for a cocompact simply connected combinatorial $G$-$2$-complex $X$ with finite edge stabilisers, the combinatorial Dehn function is well-defined if and only if the $1$-skeleton of $X$ is fine. We also show that if $H$ is a hyperbolically embedded subgroup of a finitely presented group $G$, then the relative Dehn function of the pair $(G, H)$ is well-defined. In the appendix, it is shown that show that the Baumslag-Solitar group $\mathrm{BS}(k,l)$ has a well-defined Dehn function with respect to the cyclic subgroup generated by the stable letter if and only if neither $k$ divides $l$ nor $l$ divides $k$.


title: M-FAC: Efficient Matrix-Free Approximations of Second-Order Information

abstract: Efficiently approximating local curvature information of the loss function is a key tool for optimization and compression of deep neural networks. Yet, most existing methods to approximate second-order information have high computational or storage costs, which can limit their practicality. In this work, we investigate matrix-free, linear-time approaches for estimating Inverse-Hessian Vector Products (IHVPs) for the case when the Hessian can be approximated as a sum of rank-one matrices, as in the classic approximation of the Hessian by the empirical Fisher matrix. We propose two new algorithms as part of a framework called M-FAC: the first algorithm is tailored towards network compression and can compute the IHVP for dimension $d$, if the Hessian is given as a sum of $m$ rank-one matrices, using $O(dm^2)$ precomputation, $O(dm)$ cost for computing the IHVP, and query cost $O(m)$ for any single element of the inverse Hessian. The second algorithm targets an optimization setting, where we wish to compute the product between the inverse Hessian, estimated over a sliding window of optimization steps, and a given gradient direction, as required for preconditioned SGD. We give an algorithm with cost $O(dm + m^2)$ for computing the IHVP and $O(dm + m^3)$ for adding or removing any gradient from the sliding window. These two algorithms yield state-of-the-art results for network pruning and optimization with lower computational overhead relative to existing second-order methods. Implementations are available at [9] and [17].


title: Performance Evaluation of Mixed-Precision Runge-Kutta Methods

abstract: Additive Runge-Kutta methods designed for preserving highly accurate solutions in mixed-precision computation were proposed and analyzed in [8]. These specially designed methods use reduced precision or the implicit computations and full precision for the explicit computations. We develop a FORTRAN code to solve a nonlinear system of ordinary differential equations using the mixed precision additive Runge-Kutta (MP-ARK) methods on IBM POWER9 and Intel x86\_64 chips. The convergence, accuracy, runtime, and energy consumption of these methods is explored. We show that these MP-ARK methods efficiently produce accurate solutions with significant reductions in runtime (and by extension energy consumption).


title: Novel Visual Category Discovery with Dual Ranking Statistics and Mutual Knowledge Distillation

abstract: In this paper, we tackle the problem of novel visual category discovery, i.e., grouping unlabelled images from new classes into different semantic partitions by leveraging a labelled dataset that contains images from other different but relevant categories. This is a more realistic and challenging setting than conventional semi-supervised learning. We propose a two-branch learning framework for this problem, with one branch focusing on local part-level information and the other branch focusing on overall characteristics. To transfer knowledge from the labelled data to the unlabelled, we propose using dual ranking statistics on both branches to generate pseudo labels for training on the unlabelled data. We further introduce a mutual knowledge distillation method to allow information exchange and encourage agreement between the two branches for discovering new categories, allowing our model to enjoy the benefits of global and local features. We comprehensively evaluate our method on public benchmarks for generic object classification, as well as the more challenging datasets for fine-grained visual recognition, achieving state-of-the-art performance.


title: Coupling charge and topological reconstructions at polar oxide interfaces

abstract: In oxide heterostructures, different materials are integrated into a single artificial crystal, resulting in a breaking of inversion-symmetry across the heterointerfaces. A notable example is the interface between polar and non-polar materials, where valence discontinuities lead to otherwise inaccessible charge and spin states. This approach paved the way to the discovery of numerous unconventional properties absent in the bulk constituents. However, control of the geometric structure of the electronic wavefunctions in correlated oxides remains an open challenge. Here, we create heterostructures consisting of ultrathin SrRuO$_3$, an itinerant ferromagnet hosting momentum-space sources of Berry curvature, and LaAlO$_3$, a polar wide-bandgap insulator. Transmission electron microscopy reveals an atomically sharp LaO/RuO$_2$/SrO interface configuration, leading to excess charge being pinned near the LaAlO$_3$/SrRuO$_3$ interface. We demonstrate through magneto-optical characterization, theoretical calculations and transport measurements that the real-space charge reconstruction modifies the momentum-space Berry curvature in SrRuO$_3$, driving a reorganization of the topological charges in the band structure. Our results illustrate how the topological and magnetic features of oxides can be manipulated by engineering charge discontinuities at oxide interfaces.


title: Prospects of a Pb$^{2+}$ ion clock

abstract: We propose a high-performance atomic clock based on the 1.81 PHz transition between the ground and first-excited state of doubly ionized lead. Utilizing an even isotope of lead, both clock states have $I=J=F=0$, where $I$, $J$, and $F$ are the conventional quantum numbers specifying nuclear, electronic, and total angular momentum, respectively. The clock states are nondegenerate and completely immune to nonscalar perturbations, including first order Zeeman and electric quadrupole shifts. Additionally, the proposed clock is relatively insusceptible to other frequency shifts (blackbody radiation, second order Zeeman, Doppler), accommodates "magic" rf trapping, and is robust against decoherence mechanisms that can otherwise limit clock stability. By driving the transition as a two-photon $E1$+$M1$ process, the accompanying probe Stark shift is appreciable yet manageable for practical Rabi frequencies.


title: A comparative study of various Deep Learning techniques for spatio-temporal Super-Resolution reconstruction of Forced Isotropic Turbulent flows

abstract: Super-resolution is an innovative technique that upscales the resolution of an image or a video and thus enables us to reconstruct high-fidelity images from low-resolution data. This study performs super-resolution analysis on turbulent flow fields spatially and temporally using various state-of-the-art machine learning techniques like ESPCN, ESRGAN and TecoGAN to reconstruct high-resolution flow fields from low-resolution flow field data, especially keeping in mind the need for low resource consumption and rapid results production/verification. The dataset used for this study is extracted from the 'isotropic 1024 coarse' dataset which is a part of Johns Hopkins Turbulence Databases (JHTDB). We have utilized pre-trained models and fine tuned them to our needs, so as to minimize the computational resources and the time required for the implementation of the super-resolution models. The advantages presented by this method far exceed the expectations and the outcomes of regular single structure models. The results obtained through these models are then compared using MSE, PSNR, SAM, VIF and SCC metrics in order to evaluate the upscaled results, find the balance between computational power and output quality, and then identify the most accurate and efficient model for spatial and temporal super-resolution of turbulent flow fields.


title: The automorphism group of a valued field of generalised formal power series

abstract: Let $ k $ be a field, $ G $ a totally ordered abelian group and $ \mathbb K = k((G)) $ the maximal field of generalised power series, endowed with the canonical valuation $ v $. We study the group $ v \mathrm{-Aut} K $ of valuation preserving automorphisms of a subfield $ k(G)\subseteq K\subseteq \mathbb K $, where $ k(G) $ is the fraction field of the group ring $ k[G] $. Under the assumption that $ K $ satisfies two lifting properties we are able to generalise and refine Hofberger's decomposition of $ v \mathrm{-Aut}\mathbb K $ and prove a structure theorem decomposing $ v\mathrm{-Aut} K $ into a 4-factor semi-direct product of notable subgroups. We then identify a large class of Hahn fields satisfying the two aforementioned lifting properties. Next we focus on the group of strongly additive automorphisms of $ K $. We give an explicit description of the group of strongly additive internal automorphisms in terms of the groups of homomorphisms $ \mathrm{Hom}(G,k^\times) $ of $ G $ into $ k^\times $ and $ \mathrm{Hom}(G,1+I_K) $ of $ G $ into the group of 1-units of the valuation ring of $ K $. Finally, we specialise our results to some relevant special cases. In particular, we extend the work of Schilling on the field of Laurent series and that of Deschamps on the field of Puiseux series.


title: Critical point asymptotics for Gaussian random waves with densities of any Sobolev regularity

abstract: We consider Gaussian random monochromatic waves $u$ on the plane depending on a real parameter $s$ that is directly related to the regularity of its Fourier transform. Specifically, the Fourier transform of $u$ is $f\,d\sigma$, where $d\sigma$ is the Hausdorff measure on the unit circle and the density $f$ is a function on the circle that, roughly speaking, has exactly $s-\frac12$ derivatives in $L^2$ almost surely. When $s=0$, one recovers the classical setting for random waves with a translation-invariant covariance-kernel. The main thrust of this paper is to explore the connection between the regularity parameter $s$ and the asymptotic behavior of the number $N(\nabla u,R)$ of critical points that are contained in the disk of radius $R\gg1$. More precisely, we show that the expectation $\mathbb{E}N(\nabla u,R)$ grows like the area of the disk when the regularity is low enough ($s<\frac32$) and like the diameter when the regularity is high enough ($s>\frac52$), and that the corresponding exponent changes according to a linear interpolation law in the intermediate regime. The transitions occurring at the endpoint cases involve the square root of the logarithm of the radius. Interestingly, the highest asymptotic growth rate occurs only in the classical translation-invariant setting, $s=0$. A key step of the proof of this result is the obtention of precise asymptotic expansions for certain Neumann series of Bessel functions. When the regularity parameter is $s>5$, we show that in fact $N(\nabla u,R)$ grows like the diameter with probability 1, albeit the ratio is not a universal constant but a random variable.


title: Barrier billiard and random matrices

abstract: The barrier billiard is the simplest example of pseudo-integrable models with interesting and intricate classical and quantum properties. Using the Wiener-Hopf method it is demonstrated that quantum mechanics of a rectangular billiard with a barrier in the centre can be reduced to the investigation of a certain unitary matrix. Under heuristic assumptions this matrix is substituted by a special low-complexity random unitary matrix of independent interest. The main results of the paper are (i) spectral statistics of such billiards is insensitive to the barrier height and (ii) it is well described by the semi-Poisson distributions.


title: Regularity of the SLE$_4$ uniformizing map and the SLE$_8$ trace

abstract: We show that the modulus of continuity of the SLE$_4$ uniformizing map is given by $(\log \delta^{-1})^{-1/3+o(1)}$ as $\delta \to 0$. As a consequence of our analysis, we show that the Jones-Smirnov condition for conformal removability (with quasihyperbolic geodesics) does not hold for SLE$_4$. We also show that the modulus of continuity for SLE$_8$ with the capacity time parameterization is given by $(\log \delta^{-1})^{-1/4+o(1)}$ as $\delta \to 0$, proving a conjecture of Alvisio and Lawler.


title: Estimation and Inference in Factor Copula Models with Exogenous Covariates

abstract: A factor copula model is proposed in which factors are either simulable or estimable from exogenous information. Point estimation and inference are based on a simulated methods of moments (SMM) approach with non-overlapping simulation draws. Consistency and limiting normality of the estimator is established and the validity of bootstrap standard errors is shown. Doing so, previous results from the literature are verified under low-level conditions imposed on the individual components of the factor structure. Monte Carlo evidence confirms the accuracy of the asymptotic theory in finite samples and an empirical application illustrates the usefulness of the model to explain the cross-sectional dependence between stock returns.


title: Building Stable Off-chain Payment Networks

abstract: Payment channel is a protocol which allows cryptocurrency users to route multiple transactions through network without committing them to the main blockchain network (mainnet). This ability makes them the most prominent solution to blockchains' scalability problem. Each modification of payment channels requires a transaction on the mainnet and therefore, big transaction fees. In this paper, we assume that a set of payment transactions are given (batch or online) and we study the problem of scheduling modificiations on payment channels to route all of the transactions with minimum modification cost. We investigate two cost models for aforementioned problem: the step cost function in which every channel modification has a constant cost and the linear cost function in which modification costs are proportional to the amount of change. For the step cost function model, we prove impossibility results for both batch and online case. Moreover, some heuristic methods for the batch case are presented and compared. For the linear cost we propose a polynomial time algorithm using linear programming for the batch case.


title: Partial 3D Object Retrieval using Local Binary QUICCI Descriptors and Dissimilarity Tree Indexing

abstract: A complete pipeline is presented for accurate and efficient partial 3D object retrieval based on Quick Intersection Count Change Image (QUICCI) binary local descriptors and a novel indexing tree. It is shown how a modification to the QUICCI query descriptor makes it ideal for partial retrieval. An indexing structure called Dissimilarity Tree is proposed which can significantly accelerate searching the large space of local descriptors; this is applicable to QUICCI and other binary descriptors. The index exploits the distribution of bits within descriptors for efficient retrieval. The retrieval pipeline is tested on the artificial part of SHREC'16 dataset with near-ideal retrieval results.


title: A comment on a possible inadequacy of new redefinitions of heat and work in quantum thermodynamics

abstract: We analyze the new redefinitions of heat Q and work W recently presented in [arXiv: 1912.01939; arXiv:1912.01983v5] in the quantum thermodynamics domain. According to these redefinitions, heat must be associated with the variation of entropy, while work must be associated with variation of state vectors. Analyzing the behavior of two specific examples, we show some peculiarities of these new redefinitions which, based on the counterexample presented, seems to point to a possible inadequacy of these redefinitions.


title: Nodal count for Dirichlet-to-Neumann operators with potential

abstract: We consider Dirichlet-to-Neumann operators associated to $\Delta+q$ on a Lipschitz domain in a smooth manifold, where $q$ is an $L^{\infty}$ potential. We prove a Courant-type bound for the nodal count of the extensions $u_k$ of the $k$th Dirichlet-to-Neumann eigenfunctions $\phi_k$ to the interior satisfying $(\Delta+q)u_k=0$. The classical Courant nodal domain theorem is known to hold for Steklov eigenfunctions, which are the harmonic extension of the Dirichlet-to-Neumann eigenfunctions associated to $\Delta$. Our result extends it to a larger family of Dirichlet-to-Neumann operators. Our proof makes use of the duality between the Steklov and Robin problems.


title: Joint Measurement of the $^{235}$U Antineutrino Spectrum by Prospect and Stereo

abstract: The PROSPECT and STEREO collaborations present a combined measurement of the pure $^{235}$U antineutrino spectrum, without site specific corrections or detector-dependent effects. The spectral measurements of the two highest precision experiments at research reactors are found to be compatible with $\chi^2/\mathrm{ndf} = 24.1/21$, allowing a joint unfolding of the prompt energy measurements into antineutrino energy. This $\bar{\nu}_e$ energy spectrum is provided to the community, and an excess of events relative to the Huber model is found in the 5-6 MeV region. When a Gaussian bump is fitted to the excess, the data-model $\chi^2$ value is improved, corresponding to a $2.4\sigma$ significance.


title: Theoretical modeling of dendrite growth from conductive wire electropolymerization

abstract: Electropolymerization is a bottom-up materials engineering process of micro and nano-scale that utilizes electrical signals to deposit conducting dendrites' morphologies by a redox reaction in the liquid phase. It resembles synaptogenesis in the brain, in which electrical stimulation in the brain causes the formation of synapses from the cellular neural composites. The strategy has been recently explored for neuromorphic engineering by establishing link between the electrical signals and the dendrites' shapes. Since the geometry of these structures determines their electrochemical properties, understanding the mechanisms that regulate the polymer assembly under electrically programmed conditions is an important aspect. In this manuscript, we simulate this phenomenon using mesoscale simulations, taking into account the important features of spatial-temporal potential mapping based on the time-varying signal, the motion of charged particles in the liquid due to the electric field, and the attachment of particles on the electrode. The study helps in visualizing the motion of particles in different electrical conditions, which is not possible to probe experimentally. Consistent with the experiments, the higher AC frequency of electrical activities favors linear wire-like growth, while lower frequency leads to more dense and fractal dendrites growth, and voltage offset leads to asymmetrical growth. We find that dendrites' shape and growth process systematically depend on particle concentration and random scattering. We discover that the different dendrites' architectures are associated with different Laplace and diffusion fields, which govern the monomers trajectory and subsequent dendrites' growth. Such unconventional engineering routes could have a variety of applications from neuromorphic engineering to bottom-up computing strategies.

categories: cond-mat.dis-nn cond-mat.stat-mech physics.app-ph q-bio.NC


title: Non-linear ballistic response of quantum spin-Hall edge states

abstract: Topological edge states exhibit dissipationless transport and electrically-driven topological phase transitions, making them ideal for next-generation transistors that are not constrained by Moore's law. Nevertheless, their dispersion has never been probed and is often assumed to be simply linear, without any rigorous justification. Here we determine the non-linear electrical response of topological edge states in the ballistic regime and demonstrate the way this response ascertains the presence of symmetry breaking terms in the edge dispersion, such as deviations from non-linearity and tilted spin quantization axes. The non-linear response stems from discontinuities in the band occupation on either side of a Zeeman gap, and its direction is set by the spin orientation with respect to the Zeeman field. We determine the edge dispersion for several classes of topological materials and discuss experimental measurement.


title: Evaluating Large Language Models Trained on Code

abstract: We introduce Codex, a GPT language model fine-tuned on publicly available code from GitHub, and study its Python code-writing capabilities. A distinct production version of Codex powers GitHub Copilot. On HumanEval, a new evaluation set we release to measure functional correctness for synthesizing programs from docstrings, our model solves 28.8% of the problems, while GPT-3 solves 0% and GPT-J solves 11.4%. Furthermore, we find that repeated sampling from the model is a surprisingly effective strategy for producing working solutions to difficult prompts. Using this method, we solve 70.2% of our problems with 100 samples per problem. Careful investigation of our model reveals its limitations, including difficulty with docstrings describing long chains of operations and with binding operations to variables. Finally, we discuss the potential broader impacts of deploying powerful code generation technologies, covering safety, security, and economics.


title: Differentiable Architecture Pruning for Transfer Learning

abstract: We propose a new gradient-based approach for extracting sub-architectures from a given large model. Contrarily to existing pruning methods, which are unable to disentangle the network architecture and the corresponding weights, our architecture-pruning scheme produces transferable new structures that can be successfully retrained to solve different tasks. We focus on a transfer-learning setup where architectures can be trained on a large data set but very few data points are available for fine-tuning them on new tasks. We define a new gradient-based algorithm that trains architectures of arbitrarily low complexity independently from the attached weights. Given a search space defined by an existing large neural model, we reformulate the architecture search task as a complexity-penalized subset-selection problem and solve it through a two-temperature relaxation scheme. We provide theoretical convergence guarantees and validate the proposed transfer-learning strategy on real data.


title: Asymptotic theory of hydrodynamic interactions between slender filaments

abstract: Hydrodynamic interactions (HIs) are important in biophysics research because they influence both the collective and the individual behaviour of microorganisms and self-propelled particles. For instance, HIs at the micro-swimmer level determine the attraction or repulsion between individuals, and hence their collective behaviour. Meanwhile, HIs between swimming appendages (e.g. cilia and flagella) influence the emergence of swimming gaits, synchronised bundles and metachronal waves. In this study, we address the issue of HIs between slender filaments separated by a distance larger than their contour length (d>L) by means of asymptotic calculations and numerical simulations. We first derive analytical expressions for the extended resistance matrix of two arbitrarily-shaped rigid filaments as a series expansion in inverse powers of d/L>1. The coefficients in our asymptotic series expansion are then evaluated using two well-established methods for slender filaments, resistive-force theory (RFT) and slender-body theory (SBT), and our asymptotic theory is verified using numerical simulations based on SBT for the case of two parallel helices. The theory captures the qualitative features of the interactions in the regime d/L>1, which opens the path to a deeper physical understanding of hydrodynamically governed phenomena such as inter-filament synchronisation and multiflagellar propulsion. To demonstrate the usefulness of our results, we next apply our theory to the case of two helices rotating side-by-side, where we quantify the dependence of all forces and torques on the distance and phase difference between them. Using our understanding of pairwise HIs, we then provide physical intuition for the case of a circular array of rotating helices. Our theoretical results will be useful for the study of HIs between bacterial flagella, nodal cilia, and slender microswimmers.


title: Long Short-Term Transformer for Online Action Detection

abstract: We present Long Short-term TRansformer (LSTR), a temporal modeling algorithm for online action detection, which employs a long- and short-term memory mechanism to model prolonged sequence data. It consists of an LSTR encoder that dynamically leverages coarse-scale historical information from an extended temporal window (e.g., 2048 frames spanning of up to 8 minutes), together with an LSTR decoder that focuses on a short time window (e.g., 32 frames spanning 8 seconds) to model the fine-scale characteristics of the data. Compared to prior work, LSTR provides an effective and efficient method to model long videos with fewer heuristics, which is validated by extensive empirical analysis. LSTR achieves state-of-the-art performance on three standard online action detection benchmarks, THUMOS'14, TVSeries, and HACS Segment. Code has been made available at: https://xumingze0308.github.io/projects/lstr


title: Emergence of the Classical from within the Quantum Universe

abstract: Decoherence shows how the openness of quantum systems -- interaction with their environment -- suppresses flagrant manifestations of quantumness. Einselection accounts for the emergence of preferred quasi-classical pointer states. Quantum Darwinism goes beyond decoherence. It posits that the information acquired by the monitoring environment responsible for decoherence is disseminated, in many copies, in the environment, and thus becomes accessible to observers. This indirect nature of the acquisition of information by observers who use the environment as a communication channel is the mechanism through which objective classical reality emerges from the quantum substrate: States of the systems of interest are not subjected to direct measurements (hence, not perturbed) by the agents acquiring information about them. Thus, they can exist unaffected by the information gained by observers.


title: Prospects for probing gravitational waves from primordial black hole binaries

abstract: We study the prospects of future gravitational wave (GW) detectors in probing primordial black hole (PBH) binaries. We show that across a broad mass range from $10^{-5}M_\odot$ to $10^7M_\odot$, future GW interferometers provide a potential probe of the PBH abundance that is more sensitive than any currently existing experiment. In particular, we find that galactic PBH binaries with masses as low as $10^{-5}M_\odot$ may be probed with ET, AEDGE and LISA by searching for nearly monochromatic continuous GW signals. Such searches could independently test the PBH interpretation of the ultrashort microlensing events observed by OGLE. We also consider the possibility of observing GWs from asteroid mass PBH binaries through graviton-photon conversion.


abstract: The ability to autonomously learn behaviors via direct interactions in uninstrumented environments can lead to generalist robots capable of enhancing productivity or providing care in unstructured settings like homes. Such uninstrumented settings warrant operations only using the robot's proprioceptive sensor such as onboard cameras, joint encoders, etc which can be challenging for policy learning owing to the high dimensionality and partial observability issues. We propose RRL: Resnet as representation for Reinforcement Learning -- a straightforward yet effective approach that can learn complex behaviors directly from proprioceptive inputs. RRL fuses features extracted from pre-trained Resnet into the standard reinforcement learning pipeline and delivers results comparable to learning directly from the state. In a simulated dexterous manipulation benchmark, where the state of the art methods fail to make significant progress, RRL delivers contact rich behaviors. The appeal of RRL lies in its simplicity in bringing together progress from the fields of Representation Learning, Imitation Learning, and Reinforcement Learning. Its effectiveness in learning behaviors directly from visual inputs with performance and sample efficiency matching learning directly from the state, even in complex high dimensional domains, is far from obvious.


title: Peaked bulk crystal nucleation in charged sphere melts from salt concentration dependent crystallization experiments at very low metastability

abstract: We determined bulk crystal nucleation rates in aqueous suspensions of charged spheres at low metastability. Experiments were performed in dependence on electrolyte concen-tration and for two different particle number densities. The time-dependent nucleation rate shows a pronounced initial peak, while post-solidification crystal size distributions are skewed towards larger crystallite sizes. At each concentration, the nucleation rate density initially drops exponentially with increasing salt concentration. The complete data set, however, shows an unexpected scaling of the nucleation rate densities with met-astability times the number density of particles. Parameterization of our results in terms of Classical Nucleation Theory reveals unusually low interfacial free energies of the nu-cleus surfaces and nucleation barriers well below the thermal energy. We tentatively attribute our observations to the presence of doublets introduced by the employed con-ditioning technique and acting as nucleation seeds.


title: BEDE: Bayesian Estimates of Dust Evolution For Nearby Galaxies

abstract: We build a rigorous statistical framework to provide constraints on the chemical and dust evolution parameters for nearby late-type galaxies with a wide range of gas fractions ($3\%<f_g<94\%$). A Bayesian Monte Carlo Markov Chain framework provides statistical constraints on the parameters used in chemical evolution models. Nearly a million one-zone chemical and dust evolution models were compared to 340 galaxies. Relative probabilities were calculated from the $\chi^2$ between data and models, marginalised over the different time steps, galaxy masses and star formation histories. We applied this method to find `best fitting' model parameters related to metallicity, and subsequently fix these metal parameters to study the dust parameters. For the metal parameters, a degeneracy was found between the choice of initial mass function, supernova metal yield tables and outflow prescription. For the dust parameters, the uncertainties on the best fit values are often large except for the fraction of metals available for grain growth, which is well constrained. We find a number of degeneracies between the dust parameters, limiting our ability to discriminate between chemical models using observations only. For example, we show that the low dust content of low-metallicity galaxies can be resolved by either reducing the supernova dust yields and/or including photo-fragmentation. We also show that supernova dust dominates the dust mass for low metallicity galaxies and grain growth dominates for high metallicity galaxies. The transition occurs around $12+\log({\rm O/H})=7.75$, which is lower than found in most studies in the literature.


title: Assessing putative bias in prediction of anti-microbial resistance from real-world genotyping data under explicit causal assumptions

abstract: Whole genome sequencing (WGS) is quickly becoming the customary means for identification of antimicrobial resistance (AMR) due to its ability to obtain high resolution information about the genes and mechanisms that are causing resistance and driving pathogen mobility. By contrast, traditional phenotypic (antibiogram) testing cannot easily elucidate such information. Yet development of AMR prediction tools from genotype-phenotype data can be biased, since sampling is non-randomized. Sample provenience, period of collection, and species representation can confound the association of genetic traits with AMR. Thus, prediction models can perform poorly on new data with sampling distribution shifts. In this work -- under an explicit set of causal assumptions -- we evaluate the effectiveness of propensity-based rebalancing and confounding adjustment on AMR prediction using genotype-phenotype AMR data from the Pathosystems Resource Integration Center (PATRIC). We select bacterial genotypes (encoded as k-mer signatures, i.e. DNA fragments of length k), country, year, species, and AMR phenotypes for the tetracycline drug class, preparing test data with recent genomes coming from a single country. We test boosted logistic regression (BLR) and random forests (RF) with/without bias-handling. On 10,936 instances, we find evidence of species, location and year imbalance with respect to the AMR phenotype. The crude versus bias-adjusted change in effect of genetic signatures on AMR varies but only moderately (selecting the top 20,000 out of 40+ million k-mers). The area under the receiver operating characteristic (AUROC) of the RF (0.95) is comparable to that of BLR (0.94) on both out-of-bag samples from bootstrap and the external test (n=1,085), where AUROCs do not decrease. We observe a 1%-5% gain in AUROC with bias-handling compared to the sole use of genetic signatures. ...


title: Differential Equations of Genus Four Hyperelliptic $\wp$ Functions

abstract: In order to find higher dimensional integrable models, we study differential equations of hyperelliptic $\wp$ functions up to genus four. For genus two, differential equations of hyperelliptic $\wp$ functions can be written in the Hirota form. If the genus is more than one, we have KdV equation. If the genus is more than two, we have KdV and another KdV equations. If the genus becomes more than three, there appear differential equations which cannot be written in the Hirota form, which means that the Hirota form is not enough to characterize the integrable differential equations. We have shown that some differential equations are satisfied for general genus. We can obtain differential equations for general genus step by step.


title: Rating and aspect-based opinion graph embeddings for explainable recommendations

abstract: The success of neural network embeddings has entailed a renewed interest in using knowledge graphs for a wide variety of machine learning and information retrieval tasks. In particular, recent recommendation methods based on graph embeddings have shown state-of-the-art performance. In general, these methods encode latent rating patterns and content features. Differently from previous work, in this paper, we propose to exploit embeddings extracted from graphs that combine information from ratings and aspect-based opinions expressed in textual reviews. We then adapt and evaluate state-of-the-art graph embedding techniques over graphs generated from Amazon and Yelp reviews on six domains, outperforming baseline recommenders. Additionally, our method has the advantage of providing explanations that involve the coverage of aspect-based opinions given by users about recommended items.


title: Completely compact Herz-Schur multipliers of dynamical systems

abstract: We prove that if $G$ is a discrete group and $(A,G,\alpha)$ is a C*-dynamical system such that the reduced crossed product $A\rtimes_{r,\alpha} G$ possesses property (SOAP) then every completely compact Herz-Schur $(A,G,\alpha)$-multiplier can be approximated in the completely bounded norm by Herz-Schur $(A,G,\alpha)$-multipliers of finite rank. As a consequence, if $G$ has the approximation property (AP) then the completely compact Herz-Schur multipliers of $A(G)$ coincide with the closure of $A(G)$ in the completely bounded multiplier norm. We study the class of invariant completely compact Herz-Schur multipliers of $A\rtimes_{r,\alpha} G$ and provide a description of this class in the case of the irrational rotation algebra.


title: Sleep syndromes onset detection based on automatic sleep staging algorithm

abstract: In this paper, we propose a novel method and a practical approach to predicting early onsets of sleep syndromes, including restless leg syndrome, insomnia, based on an algorithm that is comprised of two modules. A Fast Fourier Transform is applied to 30 seconds long epochs of EEG recordings to provide localized time-frequency information, and a deep convolutional LSTM neural network is trained for sleep stage classification. Automating sleep stages detection from EEG data offers great potential to tackling sleep irregularities on a daily basis. Thereby, a novel approach for sleep stage classification is proposed which combines the best of signal processing and statistics. In this study, we used the PhysioNet Sleep European Data Format (EDF) Database. The code evaluation showed impressive results, reaching an accuracy of 86.43, precision of 77.76, recall of 93,32, F1-score of 89.12 with the final mean false error loss of 0.09.


title: Downsizing revised: Star formation timescales for elliptical galaxies with an environment-dependent IMF and number of SNIa

abstract: Previous studies of the stellar mean metallicity and [Mg/Fe] values of massive elliptical (E)~galaxies suggest that their stars were formed in a very short timescale which cannot be reconciled with estimates from stellar population synthesis (SPS) studies and with hierarchical-assembly. Applying the previously developed chemical evolution code, GalIMF, which allows an environment-dependent stellar initial mass function (IMF) to be applied in the integrated galaxy initial mass function (IGIMF) theory instead of an invariant canonical IMF, the star formation timescales (SFT) of E galaxies are re-evaluated. The code's uniqueness lies in it allowing the galaxy-wide IMF and associated chemical enrichment to evolve as the physical conditions in the galaxy change. The calculated SFTs become consistent with the independent SPS results if the number of type Ia supernovae (SNIa) per unit stellar mass increases for more massive E~galaxies. This is a natural outcome of galaxies with higher star-formation rates producing more massive star clusters, spawning a larger number of SNIa progenitors per star. The calculations show E~galaxies with a stellar mass $\approx 10^{9.5} M_\odot$ to have had the longest mean SFTs of $\approx2\,$Gyr. The bulk of more massive E~galaxies were formed faster (SFT$\,\approx 1\,$Gyr) leading to domination by M~dwarf stars and larger dynamical mass-to-light ratios as observed, while lower-mass galaxies tend to lose their gas supply more easily due to their shallower potential and therefore also have similarly-short mean SFTs. This work achieves, for the first time, consistency of the SFTs for early-type galaxies between chemical-enrichment and SPS modelling and leads to an improved understanding of how the star formation environment may affect the total number of SNIa per unit stellar mass formed.


title: Natural-Scalaron Inflation

abstract: A pseudo Nambu-Goldstone boson (such as an axion-like particle) is a theoretically well-motivated inflaton as it features a naturally flat potential (natural inflation). This is because Goldstone's theorem protects its potential from sizable quantum corrections. Such corrections, however, generically generates an $R^2$ term in the action, which leads to another inflaton candidate because of the equivalence between the $R^2$ term and a scalar field, the scalaron, with a quasi flat potential (Starobinsky inflation). Here it is investigated a new multifield scenario in which both the scalaron and a pseudo Nambu-Goldstone boson are active (natural-scalaron inflation). For generality, also a non-minimal coupling is included, which is shown to emerge from microscopic theories. It is demonstrated that a robust inflationary attractor is present even when the masses of the two inflatons are comparable. Moreover, the presence of the scalaron allows to satisfy all observational bounds in a large region of the parameter space, unlike what happens in pure-natural inflation.


title: Inconsistency of Islands in Theories with Long-Range Gravity

abstract: In ordinary gravitational theories, any local bulk operator in an entanglement wedge is accompanied by a long-range gravitational dressing that extends to the asymptotic part of the wedge. Islands are the only known examples of entanglement wedges that are disconnected from the asymptotic region of spacetime. In this paper, we show that the lack of an asymptotic region in islands creates a potential puzzle that involves the gravitational Gauss law, independently of whether or not there is a non-gravitational bath. In a theory with long-range gravity, the energy of an excitation localized to the island can be detected from outside the island, in contradiction with the principle that operators in an entanglement wedge should commute with operators from its complement. In several known examples, we show that this tension is resolved because islands appear in conjunction with a massive graviton. We also derive some additional consistency conditions that must be obeyed by islands in decoupled systems. Our arguments suggest that islands might not constitute consistent entanglement wedges in standard theories of massless gravity where the Gauss law applies.


title: On the choice of entropy variables in multifield inflation

abstract: We discuss the usefulness and theoretical consistency of different entropy variables used in the literature to describe isocurvature perturbations in multifield inflationary models with a generic curved field space. We clarify which is the proper entropy variable to be used to match the evolution of isocurvature modes during inflation to the one after the reheating epoch in order to compare with observational constraints. In particular, we find that commonly used variables, as the relative entropy perturbation or the one associated to the decomposition in tangent and normal perturbations with respect to the inflationary trajectory, even if more useful to perform numerical studies, can lead to results which are wrong by several orders of magnitude, or even to apparent destabilisation effects which are unphysical for cases with light kinetically coupled spectator fields.


title: Spin-orbit pumping

abstract: We study theoretically the effect of a rotating electric field on a diffusive nanowire and find an effect that is analogous to spin pumping, which refers to the generation of spin through a rotating magnetic field. The electron spin couples to the electric field because the particle motion induces an effective magnetic field in its rest frame. In a diffusive system the velocity of the particle, and therefore also its effective magnetic field, rapidly and randomly changes direction. Nevertheless, we demonstrate analytically and via a physical argument why the combination of the two effects described above produces a finite magnetization along the axis of rotation. This manifests as a measurable spin-voltage in the range of tens of microvolts.


title: Operator scaling dimensions and multifractality at measurement-induced transitions

abstract: Repeated local measurements of quantum many body systems can induce a phase transition in their entanglement structure. These measurement-induced phase transitions (MIPTs) have been studied for various types of dynamics, yet most cases yield quantitatively similar values of the critical exponents, making it unclear if there is only one underlying universality class. Here, we directly probe the properties of the conformal field theories governing these MIPTs using a numerical transfer-matrix method, which allows us to extract the effective central charge, as well as the first few low-lying scaling dimensions of operators at these critical points. Our results provide convincing evidence that the generic and Clifford MIPTs for qubits lie in different universality classes and that both are distinct from the percolation transition for qudits in the limit of large onsite Hilbert space dimension. For the generic case, we find strong evidence of multifractal scaling of correlation functions at the critical point, reflected in a continuous spectrum of scaling dimensions.

categories: cond-mat.dis-nn cond-mat.stat-mech cond-mat.str-el quant-ph


title: Dynamical Formation Channels for Fast Radio Bursts in Globular Clusters

abstract: The repeating fast radio burst (FRB) localized to a globular cluster in M81 challenges our understanding of FRB models. In this Letter, we explore dynamical formation scenarios for objects in old globular clusters that may plausibly power FRBs. Using N-body simulations, we demonstrate that young neutron stars may form in globular clusters at a rate of up to $\sim50\,\rm{Gpc}^{-3}\,\rm{yr}^{-1}$ through a combination of binary white dwarf mergers, white dwarf--neutron star mergers, binary neutron star mergers, and accretion induced collapse of massive white dwarfs in binary systems. We consider two FRB emission mechanisms: First, we show that a magnetically-powered source (e.g., a magnetar with field strength $\gtrsim10^{14}\,$G) is viable for radio emission efficiencies $\gtrsim10^{-4}$. This would require magnetic activity lifetimes longer than the associated spin-down timescales and longer than empirically-constrained lifetimes of Galactic magnetars. Alternatively, if these dynamical formation channels produce young rotation-powered neutron stars with spin periods of $\sim10\,$ms and magnetic fields of $\sim10^{11}\,$G (corresponding to spin-down lifetimes of $\gtrsim10^5\,$yr), the inferred event rate and energetics can be reasonably reproduced for order unity duty cycles. Additionally, we show that recycled millisecond pulsars or low-mass X-ray binaries similar to those well-observed in Galactic globular clusters may also be plausible channels, but only if their duty cycle for producing bursts similar to the M81 FRB is small.


title: Thermal Instabilities and Shattering in the High-Redshift WHIM: Convergence Criteria and Implications for Low-Metallicity Strong HI Absorbers

abstract: Using a novel suite of cosmological simulations zooming in on a Mpc-scale intergalactic sheet or "pancake" at z~3-5, we conduct an in-depth study of the thermal properties and HI content of the warm-hot intergalactic medium (WHIM) at those redshifts. The simulations span nearly three orders of magnitude in gas-cell mass, from ~(7.7x10^6-1.5x10^4)Msun, one of the highest resolution simulations of such a large patch of the inter-galactic medium (IGM) to date. At z~5, a strong accretion shock develops around the main pancake following a collision between two smaller sheets. Gas in the post-shock region proceeds to cool rapidly, triggering thermal instabilities and the formation of a multiphase medium. We find neither the mass, nor the morphology, nor the distribution of HI in the WHIM to be converged at our highest resolution. Interestingly, the lack of convergence is more severe for the less dense, more metal-poor, intra-pancake medium (IPM) in between filaments and far from any star-forming galaxies. As the resolution increases, the IPM develops a shattered structure, with ~kpc scale clouds containing most of the HI. From our lowest to highest resolution, the covering fraction of metal-poor (Z<10^{-3}Zsun) Lyman-limit systems (NHI>10^{17.2}/cm^2) in the IPM at z~4 increases from (3-15)%, while that of Damped Lyman-alpha Absorbers (NHI>10^{20}/cm^2) with similar metallicity increases threefold, from (0.2-0.6)%, with no sign of convergence. We find that a necessary condition for the formation of a multiphase, shattered structure is resolving the cooling length, lcool=cs*tcool, at T~10^5K. If this scale is unresolved, gas "piles up" at these temperatures and cooling to lower temperatures becomes very inefficient. We conclude that state-of-the-art cosmological simulations are still unable to resolve the multi-phase structure of the low-density IGM, with potentially far-reaching implications.