title: Climbing the Diagonal Clifford Hierarchy

abstract: Magic state distillation and the Shor factoring algorithm make essential use of logical diagonal gates. We introduce a method of synthesizing CSS codes that realize a target logical diagonal gate at some level $l$ in the Clifford hierarchy. The method combines three basic operations: concatenation, removal of $Z$-stabilizers, and addition of $X$-stabilizers. It explicitly tracks the logical gate induced by a diagonal physical gate that preserves a CSS code. The first step is concatenation, where the input is a CSS code and a physical diagonal gate at level $l$ inducing a logical diagonal gate at the same level. The output is a new code for which a physical diagonal gate at level $l+1$ induces the original logical gate. The next step is judicious removal of $Z$-stabilizers to increase the level of the induced logical operator. We identify three ways of climbing the logical Clifford hierarchy from level $l$ to level $l+1$, each built on a recursive relation on the Pauli coefficients of the induced logical operators. Removal of $Z$-stabilizers may reduce distance, and the purpose of the third basic operation, addition of $X$-stabilizers, is to compensate for such losses. For the coherent noise model, we describe how to switch between computation and storage of intermediate results in a decoherence-free subspace by simply applying Pauli $X$ matrices. The approach to logical gate synthesis taken in prior work focuses on the code states, and results in sufficient conditions for a CSS code to be fixed by a transversal $Z$-rotation. In contrast, we derive necessary and sufficient conditions by analyzing the action of a transversal diagonal gate on the stabilizer group that determines the code. The power of our approach is demonstrated by two proofs of concept: the $[[2^{l+1}-2,2,2]]$ triorthogonal code family, and the $[[2^m,\binom{m}{r},2^{\min\{r,m-r\}}]]$ quantum Reed-Muller code family.


title: Gapoera: Application Programming Interface for AI Environment of Indonesian Board Game

abstract: Currently, the development of computer games has shown a tremendous surge. The ease and speed of internet access today have also influenced the development of computer games, especially computer games that are played online. Internet technology has allowed computer games to be played in multiplayer mode. Interaction between players in a computer game can be built in several ways, one of which is by providing balanced opponents. Opponents can be developed using intelligent agents. On the other hand, research on developing intelligent agents is also growing rapidly. In computer game development, one of the easiest ways to measure the performance of an intelligent agent is to develop a virtual environment that allows the intelligent agent to interact with other players. In this research, we try to develop an intelligent agent and virtual environment for the board game. To be easily accessible, the intelligent agent and virtual environment are then developed into an Application Programming Interface (API) service called Gapoera API. The Gapoera API service that is built is expected to help game developers develop a game without having to think much about the artificial intelligence that will be embedded in the game. This service provides a basic multilevel intelligent agent that can provide users with playing board games commonly played in Indonesia. Although the Gapoera API can be used for various types of games, in this paper, we will focus on the discussion on a popular traditional board game in Indonesia, namely Mancala. The test results conclude that the multilevel agent concept developed has worked as expected. On the other hand, the development of the Gapoera API service has also been successfully accessed on several game platforms.


title: A Path-based Approach to Analyzing the Global Liner Shipping Network

abstract: The maritime shipping network is the backbone of global trade. Data about the movement of cargo through this network comes in various forms, from ship-level Automatic Identification System (AIS) data, to aggregated bilateral trade volume statistics. Multiple network representations of the shipping system can be derived from any one data source, each of which has advantages and disadvantages. In this work, we examine data in the form of liner shipping service routes, a list of walks through the port-to-port network aggregated from individual shipping companies by a large shipping logistics database. This data is inherently sequential, in that each route represents a sequence of ports called upon by a cargo ship. Previous work has analyzed this data without taking full advantage of the sequential information. Our contribution is to develop a path-based methodology for analyzing liner shipping service route data, computing navigational trajectories through the network that respect the routes and comparing these paths with those computed using other network representations of the same data. We further use these trajectories to re-analyze the role of a previously-identified structural core through the network, as well as to define and analyze a measure of betweenness centrality for nodes and edges.


title: Gauge and Scalar Fields on $\mathbb{CP}^2$: A Gauge-invariant Analysis I. The effective action from chiral scalars

abstract: A parametrization of gauge fields on complex projective spaces of arbitrary dimension is given as a generalization of the two-dimensional case. Gauge transformations act homogeneously on the fields, facilitating a manifestly gauge-invariant analysis. Specializing to four dimensions, we consider the nature of the effective action due to chiral scalars interacting with the gauge fields. The key qualitatively significant terms include a possible gauge-invariant mass term and a finite four-dimensional Wess-Zumino-Witten (WZW) action. We comment on relating the mass term to lattice simulations as well as Schwinger-Dyson analyses, and also on relating the WZW action to the instanton liquid picture of QCD.


title: Solving Large-Scale Granular Resource Allocation Problems Efficiently with POP

abstract: Resource allocation problems in many computer systems can be formulated as mathematical optimization problems. However, finding exact solutions to these problems using off-the-shelf solvers is often intractable for large problem sizes with tight SLAs, leading system designers to rely on cheap, heuristic algorithms. We observe, however, that many allocation problems are granular: they consist of a large number of clients and resources, each client requests a small fraction of the total number of resources, and clients can interchangeably use different resources. For these problems, we propose an alternative approach that reuses the original optimization problem formulation and leads to better allocations than domain-specific heuristics. Our technique, Partitioned Optimization Problems (POP), randomly splits the problem into smaller problems (with a subset of the clients and resources in the system) and coalesces the resulting sub-allocations into a global allocation for all clients. We provide theoretical and empirical evidence as to why random partitioning works well. In our experiments, POP achieves allocations within 1.5% of the optimal with orders-of-magnitude improvements in runtime compared to existing systems for cluster scheduling, traffic engineering, and load balancing.


title: Polyhedral expansions of compacta associated to finite approximations

abstract: This paper introduces some inverse sequences of different polyhedra all based on finite approximations of a compact metric space so they can be used to capture the shape type of the original space. It is shown that they are HPol-expansions, proving the so-called General Principle. We use these sequences to compute explicitly some inverse persistent homology groups of a space and measure its errors in the approximation process.


title: Double Trouble: How to not explain a text classifier's decisions using counterfactuals synthesized by masked language models?

abstract: A principle behind dozens of attribution methods is to take the prediction difference between before-and-after an input feature (here, a token) is removed as its attribution. A popular Input Marginalization (IM) method (Kim et al., 2020) uses BERT to replace a token, yielding more plausible counterfactuals. While Kim et al. (2020) reported that IM is effective, we find this conclusion not convincing as the DeletionBERT metric used in their paper is biased towards IM. Importantly, this bias exists in Deletion-based metrics, including Insertion, Sufficiency, and Comprehensiveness. Furthermore, our rigorous evaluation using 6 metrics and 3 datasets finds no evidence that IM is better than a Leave-One-Out (LOO) baseline. We find two reasons why IM is not better than LOO: (1) deleting a single word from the input only marginally reduces a classifier's accuracy; and (2) a highly predictable word is always given near-zero attribution, regardless of its true importance to the classifier. In contrast, making LIME samples more natural via BERT consistently improves LIME accuracy under several ROAR metrics.


title: A new perspective on cosmology through Supernovae Ia and Gamma Ray Bursts

abstract: The actual knowledge of the structure and future evolution of our universe is based on the use of cosmological models, which can be tested through the so-called 'probes', namely astrophysical phenomena, objects or structures with peculiar properties that can help to discriminate among different cosmological models. Among all the existing probes, of particular importance are the Supernovae Ia (SNe Ia) and the Gamma Ray Bursts (GRBs): the former are considered among the best standard candles so far discovered but suffer from the fact that can be observed until redshift $z=2.26$, while the latter are promising standardizable candles which have been observed up to $z=9.4$, surpassing even the farthest quasar known to date, which is at $z=7.64$. The standard candles can be used to test the cosmological models and to give the expected values of cosmological parameters, in particular the Hubble constant value. The Hubble constant is affected by the so-called \say{Hubble constant tension}, a discrepancy in more than 4 $\sigma$ between its value measured with local probes and its value measured through the cosmological probes. The increase in the number of observed SNe Ia, as well as the future standardization of GRBs through their correlations, will surely be of help in alleviating the Hubble constant tension and in explaining the structure of the universe at higher redshifts. A promising class of GRBs for future standardization is represented by the GRBs associated with Supernovae Ib/c, since these present features similar to the SNe Ia class and obey a tight correlation between their luminosity at the end of the plateau emission in X-rays and the time at the end of the plateau in the rest-frame.


title: Evolution of accretion disc reflection spectra due to a Type I X-ray burst

abstract: Irradiation of the accretion disc causes reflection signatures in the observed X-ray spectrum, encoding important information about the disc structure and density. A Type I X-ray burst will strongly irradiate the accretion disc and alter its properties. Previous numerical simulations predicted the evolution of the accretion disc due to an X-ray burst. Here, we process time-averaged simulation data of six time intervals to track changes in the reflection spectrum from the burst onset to just past its peak. We divide the reflecting region of the disc within $r\lesssim50$ km into 6-7 radial zones for every time interval and compute the reflection spectra for each zone. We integrate these reflection spectra to obtain a total reflection spectrum per time interval. The burst ionizes and heats the disc, which gradually weakens all emission lines. Compton scattering and bremsstrahlung rates increase in the disc during the burst rise, and the soft excess at $<$3 keV rises from $\approx4$% to $\approx38$% of the total emission at the burst peak. A soft excess is expected to be ubiquitous in the reflection spectra of X-ray bursts. Structural disc changes such as inflation because of heating or drainage of the inner disc due to Poynting-Robertson drag affect the strength of the soft excess. Further studies on the dependence of the reflection spectrum characteristics to changes in the accretion disc during an X-ray burst may lead to probes of the disc geometry.


title: Efficient Asymmetric photothermal source: Designing a heating Janus-Nanojet

abstract: Janus particles have flourished as subject of intensive research due to their synergetic properties and their promising use in different fields, especially in biomedicine. The combination of materials with radically different physical properties in the same nanostructure gives rise to the so-called Janus effects, allowing phenomena of a contrasting nature to occur in the same architecture. In particular, interesting advantages can be taken from a thermal Janus effect for photoinduced hyperthermia cancer therapies. Such therapies still have limitations associated to the heating control in terms of temperature stability and energy management. While previous studies have shown that some plasmonic single-material nanoparticles are somehow effective at killing cancer cells, it is necessary to investigate alternative plasmonic Janus nanoheaters to improve the heating efficiency and thermal control, mainly because the widespread single-material nanoheaters are highly homogeneous sources of heat, which implies that the surrounding biological medium is isotropically heated, equally affecting cancerous and healthy cells. A detailed thermoplasmonic study of the thermal Janus effect is still missing. Here we perform such study and demonstrate that doughnut-based Janus nanoparticles exhibit an outstanding photothermal control under practical illumination conditions, i.e., unpolarized light. Furthermore, we present novel and effective Janus nanoparticle designs that possess superior photothermal conversion features and unique directional heating capacity, being able to channel up to 91% of the total thermal energy onto a target. We discuss the implications of these innovative nanoparticles with regards to thermoplasmonics hyperthermia cancer therapy.


title: Mathematical Model of Volume Kinetics and Renal Function after Burn Injury and Resuscitation

abstract: This paper presents a mathematical model of blood volume kinetics and renal function in response to burn injury and resuscitation, which is applicable to the development and non-clinical testing of burn resuscitation protocols and algorithms. Prior mathematical models of burn injury and resuscitation are not ideally suited to such applications due to their limited credibility in predicting blood volume and urinary output observed in wide-ranging burn patients as well as in incorporating contemporary knowledge of burn pathophysiology. Our mathematical model consists of an established multi-compartmental model of blood volume kinetics, a hybrid mechanistic-phenomenological model of renal function, and novel lumped-parameter models of burn-induced perturbations in volume kinetics and renal function equipped with contemporary knowledge on burn-related physiology and pathophysiology. Using the dataset collected from 16 sheep, we showed that our mathematical model can be characterized with physiologically plausible parameter values to accurately predict blood volume kinetic and renal function responses to burn injury and resuscitation on an individual basis against a wide range of pathophysiological variability. Pending validation in humans, our mathematical model may serve as an effective basis for an in-depth understanding of complex burn-induced volume kinetic and renal function responses as well as development and non-clinical testing of burn resuscitation protocols and algorithms. Keywords: Mathematical Model, Burn Injury, Burn Resuscitation, Non-Clinical Testing, Burn Resuscitation Protocols, and Algorithms


title: Cleaning Dirty Books: Post-OCR Processing for Previously Scanned Texts

abstract: Substantial amounts of work are required to clean large collections of digitized books for NLP analysis, both because of the presence of errors in the scanned text and the presence of duplicate volumes in the corpora. In this paper, we consider the issue of deduplication in the presence of optical character recognition (OCR) errors. We present methods to handle these errors, evaluated on a collection of 19,347 texts from the Project Gutenberg dataset and 96,635 texts from the HathiTrust Library. We demonstrate that improvements in language models now enable the detection and correction of OCR errors without consideration of the scanning image itself. The inconsistencies found by aligning pairs of scans of the same underlying work provides training data to build models for detecting and correcting errors. We identify the canonical version for each of 17,136 repeatedly-scanned books from 58,808 scans. Finally, we investigate methods to detect and correct errors in single-copy texts. We show that on average, our method corrects over six times as many errors as it introduces. We also provide interesting analysis on the relation between scanning quality and other factors such as location and publication year.


title: Incommensurate magnetic orders and topological Hall effect in the square-net centrosymmetric EuGa$_2$Al$_2$ system

abstract: Neutron diffraction on the centrosymmetric square-net magnet EuGa$_2$Al$_2$ reveals multiple incommensurate magnetic states (AFM1,2,3) in zero field. In applied field, a new magnetic phase (A) is identified from magnetization and transport measurements, bounded by two of the $\mu_0H$~=~0 incommensurate magnetic phases (AFM1,helical and AFM3, cycloidal) with different moment orientations. Moreover, magneto-transport measurements indicate the presence of a topological Hall effect, with maximum values centered in the A phase. Together, these results render EuGa$_2$Al$_2$ a material with non-coplanar or topological spin texture in applied field. X-ray diffraction reveals an out-of-plane (OOP) charge density wave (CDW) below $T_{CDW} \sim$ 50 K while the magnetic propagation vector lies in plane below $T_N$ = 19.5 K. Together these data point to a new route to realizing in-plane non-collinear spin textures through an OOP CDW. In turn, these non-collinear spin textures may be unstable against the formation of topological spin textures in an applied field.


title: A fine property of Whitehead's algorithm

abstract: We develop a refinement of Whitehead's algorithm for primitive words in a free group. We generalize to subgroups, establishing a strengthened version of Whitehead's algorithm for free factors. We make use of these refinements in proving new results about primitive elements and free factors in a free group. These include a relative version of Whitehead's algorithm, and a criterion that tests whether a subgroup is a free factor just by looking at its primitive elements. We develop an algorithm to determine whether or not two vertices in the free factor complex have distance $d$ for $d=1,2,3$, as well as $d=4$ in a special case.


title: Non-expansive matrix number systems with bases similar to $J_n(1)$

abstract: We study representations of integral vectors in a number system with a matrix base $M$ and vector digits. We focus on the case when $M$ is similar to $J_n$, the Jordan block of $1$ of size $n$. If $M=J_2$, we classify digit sets of size 2 allowing representation of the whole $\mathbb{Z}^2$. For $J_n$ with $n\geq 3$, it is shown that three digits suffice to represent all of $\mathbb{Z}^n$. For bases similar to $J_n$, at most $n$ digits are required, with the exception of $n=1$. Moreover, the language of strings representing the zero vector with $M=J_2$ and the digits $(0,\pm 1)^T$ is shown not to be context-free, but to be recognizable by a Turing machine with logarithmic memory.


title: A Framework for Learning Assessment through Multimodal Analysis of Reading Behaviour and Language Comprehension

abstract: Reading comprehension, which has been defined as gaining an understanding of written text through a process of translating grapheme into meaning, is an important academic skill. Other language learning skills - writing, speaking and listening, all are connected to reading comprehension. There have been several measures proposed by researchers to automate the assessment of comprehension skills for second language (L2) learners, especially English as Second Language (ESL) and English as Foreign Language (EFL) learners. However, current methods measure particular skills without analysing the impact of reading frequency on comprehension skills. In this dissertation, we show how different skills could be measured and scored automatically. We also demonstrate, using example experiments on multiple forms of learners' responses, how frequent reading practices could impact on the variables of multimodal skills (reading pattern, writing, and oral fluency). This thesis comprises of five studies. The first and second studies are based on eye-tracking data collected from EFL readers in repeated reading (RR) sessions. The third and fourth studies are to evaluate free-text summary written by EFL readers in repeated reading sessions. The fifth and last study, described in the sixth chapter of the thesis, is to evaluate recorded oral summaries recited by EFL readers in repeated reading sessions. In a nutshell, through this dissertation, we show that multimodal skills of learners could be assessed to measure their comprehension skills as well as to measure the effect of repeated readings on these skills in the course of time, by finding significant features and by applying machine learning techniques with a combination of statistical models such as LMER.


title: Adaptive FEM for Helmholtz Equation with Large Wave Number

abstract: A posteriori upper and lower bounds are derived for the linear finite element method (FEM) for the Helmholtz equation with large wave number. It is proved rigorously that the standard residual type error estimator seriously underestimates the true error of the FE solution for the mesh size $h$ in the preasymptotic regime, which is first observed by Babu\v{s}ka, et al. for an one dimensional problem. By establishing an equivalence relationship between the error estimators for the FE solution and the corresponding elliptic projection of the exact solution, an adaptive algorithm is proposed and its convergence and quasi-optimality are proved under condition that $k^3h_0^{1+\al}$ is sufficiently small, where $h_0$ is the initial mesh size and $\frac12<\al\le 1$ is a regularity constant depending on the maximum reentrant angle of the domain. Numerical tests are given to verify the theoretical findings and to show that the adaptive continuous interior penalty finite element method (CIP-FEM) with appropriately selected penalty parameters can greatly reduce the pollution error and hence the residual type error estimator for this CIP-FEM is reliable and efficient even in the preasymptotic regime.


title: Logical Activation Functions: Logit-space equivalents of Probabilistic Boolean Operators

abstract: The choice of activation functions and their motivation is a long-standing issue within the neural network community. Neuronal representations within artificial neural networks are commonly understood as logits, representing the log-odds score of presence of features within the stimulus. We derive logit-space operators equivalent to probabilistic Boolean logic-gates AND, OR, and XNOR for independent probabilities. Such theories are important to formalize more complex dendritic operations in real neurons, and these operations can be used as activation functions within a neural network, introducing probabilistic Boolean-logic as the core operation of the neural network. Since these functions involve taking multiple exponents and logarithms, they are computationally expensive and not well suited to be directly used within neural networks. Consequently, we construct efficient approximations named $\text{AND}_\text{AIL}$ (the AND operator Approximate for Independent Logits), $\text{OR}_\text{AIL}$, and $\text{XNOR}_\text{AIL}$, which utilize only comparison and addition operations, have well-behaved gradients, and can be deployed as activation functions in neural networks. Like MaxOut, $\text{AND}_\text{AIL}$ and $\text{OR}_\text{AIL}$ are generalizations of ReLU to two-dimensions. While our primary aim is to formalize dendritic computations within a logit-space probabilistic-Boolean framework, we deploy these new activation functions, both in isolation and in conjunction to demonstrate their effectiveness on a variety of tasks including image classification, transfer learning, abstract reasoning, and compositional zero-shot learning.


title: Entropy-Stable Schemes in the Low-Mach-Number Regime: Flux-Preconditioning, Entropy Breakdowns, and Entropy Transfers

abstract: Entropy-Stable (ES) schemes, specifically those built from [Tadmor \textit{Math. Comput.} 49 (1987) 91], have been gaining interest over the past decade, especially in the context of under-resolved simulations of compressible turbulent flows using high-order methods. These schemes are attractive because they can provide stability in a global and nonlinear sense (consistency with thermodynamics). However, fully realizing the potential of ES schemes requires a better grasp of their local behavior. Entropy-stability itself does not imply good local behavior [Gouasmi \textit{et al.} \textit{J. Sci. Comp.} 78 (2019) 971, Gouasmi \textit{et al.} \textit{Comput. Methd. Appl. M.} 363 (2020) 112912]. In this spirit, we studied ES schemes in problems where \textit{global stability is not the core issue}. In the present work, we consider the accuracy degradation issues typically encountered by upwind-type schemes in the low-Mach-number regime [Turkel \textit{Annu. Rev. Fluid Mech.} 31 (1999) 285] and their treatment using \textit{Flux-Preconditioning} [Turkel \textit{J. Comput. Phys.} 72 (1987) 277, Miczek \textit{et al.} \textit{A \& A} 576 (2015) A50]. ES schemes suffer from the same issues and Flux-Preconditioning can improve their behavior without interfering with entropy-stability. This is first demonstrated analytically: using similarity and congruence transforms we were able to establish conditions for a preconditioned flux to be ES, and introduce the ES variants of the Miczek's and Turkel's preconditioned fluxes. This is then demonstrated numerically through first-order simulations of two simple test problems representative of the incompressible and acoustic limits, the Gresho Vortex and a right-moving acoustic wave. The results are overall consistent with previous studies [...]


title: On semi-classical spectral series for an atom in a periodic polarized electric field

abstract: In this report we present preliminary results about the tunneling problem for a magnetic Schr\"odinger operator. As a motivation we consider the 3-D time-dependent Schr\"odinger operator $H(t)=-h^2\Delta+V+E(t)\cdot x$ where $V$ is a radial potential and $E(t)$ a circularly polarized field with uniform frequency $\omega$. The quantum monodromy operator (QMO) that takes the system through a complete period $T=2\pi/\omega$, turns out to be unitarily equivalent to $e^{iTP_A(x,hD_x)/h}$, where $P_A(x,hD_x))$ identifies with a magnetic Schr\"odinger operator. When $V$ is sufficiently confining, $P_A(x,hD_x))$ presents a double magnetic well. Then we construct its semi-classical ground state and examine the splitting between its two first eigenvalues.


title: Solving N-player dynamic routing games with congestion: a mean field approach

abstract: The recent emergence of navigational tools has changed traffic patterns and has now enabled new types of congestion-aware routing control like dynamic road pricing. Using the fundamental diagram of traffic flows - applied in macroscopic and mesoscopic traffic modeling - the article introduces a new N-player dynamic routing game with explicit congestion dynamics. The model is well-posed and can reproduce heterogeneous departure times and congestion spill back phenomena. However, as Nash equilibrium computations are PPAD-complete, solving the game becomes intractable for large but realistic numbers of vehicles N. Therefore, the corresponding mean field game is also introduced. Experiments were performed on several classical benchmark networks of the traffic community: the Pigou, Braess, and Sioux Falls networks with heterogeneous origin, destination and departure time tuples. The Pigou and the Braess examples reveal that the mean field approximation is generally very accurate and computationally efficient as soon as the number of vehicles exceeds a few dozen. On the Sioux Falls network (76 links, 100 time steps), this approach enables learning traffic dynamics with more than 14,000 vehicles.

categories: math.DS cs.MA cs.NI cs.SY eess.SY math.OC


title: Solutions to the MiniBooNE Anomaly from New Physics in Charged Meson Decays

abstract: We point out that production of new bosons by charged meson decays can greatly enhance the sensitivity of beam-focused accelerator-based experiments to new physics signals. This enhancement arises since the charged mesons are focused and their three-body decays do not suffer from helicity suppression in the same way as their usual two-body decays. As a realistic application, we attempt to explain the MiniBooNE low energy excess utilizing this overlooked mechanism, uniquely realizing dark-sector interpretations as plausible solutions to the excess. As proof of the principle, we consider two well-motivated classes of dark-sector models, models of vector-portal dark matter and models of long-lived (pseudo)scalar. We argue that the model parameter values to accommodate the excess are consistent with existing limits and that they can be tested at current and future accelerator-based neutrino experiments.


title: SOFT: Softmax-free Transformer with Linear Complexity

abstract: Vision transformers (ViTs) have pushed the state-of-the-art for various visual recognition tasks by patch-wise image tokenization followed by self-attention. However, the employment of self-attention modules results in a quadratic complexity in both computation and memory usage. Various attempts on approximating the self-attention computation with linear complexity have been made in Natural Language Processing. However, an in-depth analysis in this work shows that they are either theoretically flawed or empirically ineffective for visual recognition. We further identify that their limitations are rooted in keeping the softmax self-attention during approximations. Specifically, conventional self-attention is computed by normalizing the scaled dot-product between token feature vectors. Keeping this softmax operation challenges any subsequent linearization efforts. Based on this insight, for the first time, a softmax-free transformer or SOFT is proposed. To remove softmax in self-attention, Gaussian kernel function is used to replace the dot-product similarity without further normalization. This enables a full self-attention matrix to be approximated via a low-rank matrix decomposition. The robustness of the approximation is achieved by calculating its Moore-Penrose inverse using a Newton-Raphson method. Extensive experiments on ImageNet show that our SOFT significantly improves the computational efficiency of existing ViT variants. Crucially, with a linear complexity, much longer token sequences are permitted in SOFT, resulting in superior trade-off between accuracy and complexity.


title: Detection of OH in the ultra-hot Jupiter WASP-76b

abstract: Ultra-hot Jupiters have dayside temperatures at which most molecules are expected to thermally dissociate. The dissociation of water vapour results in the production of the hydroxyl radical (OH). While OH absorption is easily observed in near-infrared spectra of M dwarfs, which have similar effective temperatures as ultra-hot Jupiters, it is often not considered when studying the atmospheres of ultra-hot Jupiters. We use high-resolution spectroscopic observations of a transit of WASP-76b obtained using CARMENES to study the presence of OH. After validating the OH line list, we generate model transit spectra of WASP-76b with petitRADTRANS. The data are corrected for stellar and telluric contamination and cross-correlated with the model spectra. After combining all cross-correlation functions from the transit, a detection map is constructed. Constraints on the planet properties from the OH absorption are obtained from a Markov chain Monte Carlo analysis. OH is detected in the atmosphere of WASP-76b with a peak signal-to-noise ratio of 6.1. From the retrieval we obtain $K_p=232 \pm 12$ km/s and a blueshift of $-13.2 \pm 1.6$ km/s, which are offset from the expected velocities. Considering the fast spin rotation of the planet, the blueshift is best explained with the signal predominantly originating from the evening terminator and the presence of a strong dayside-to-nightside wind. The increased $K_p$ over its expected value (196.5 km/s) is, however, a bit puzzling. The signal is found to be broad, with a full width at half maximum of $16.8^{+4.6}_{-4.0}$ km/s. The retrieval results in a weak constraint on the mean temperature of 2700-3700 K at the pressure range of the OH signal. We show that OH is readily observable in the transit spectra of ultra-hot Jupiters. Studying this molecule can provide insights into the molecular dissociation processes in the atmospheres of such planets.


title: Replica Wormholes and Holographic Entanglement Negativity

abstract: Recent work has shown how to understand the Page curve of an evaporating black hole from replica wormholes. However, more detailed information about the structure of its quantum state is needed to fully understand the dynamics of black hole evaporation. Here we study entanglement negativity, an important measure of quantum entanglement in mixed states, in a couple of toy models of evaporating black holes. We find four phases dominated by different types of geometries: the disconnected, cyclically connected, anti-cyclically connected, and pairwise connected geometries. The last of these geometries are new replica wormholes that break the replica symmetry spontaneously. We also analyzed the transitions between these four phases by summing more generic replica geometries using a Schwinger-Dyson equation. In particular, we find enhanced corrections to various negativity measures near the transition between the cyclic and pairwise phase.


title: Learning Proposals for Practical Energy-Based Regression

abstract: Energy-based models (EBMs) have experienced a resurgence within machine learning in recent years, including as a promising alternative for probabilistic regression. However, energy-based regression requires a proposal distribution to be manually designed for training, and an initial estimate has to be provided at test-time. We address both of these issues by introducing a conceptually simple method to automatically learn an effective proposal distribution, which is parameterized by a separate network head. To this end, we derive a surprising result, leading to a unified training objective that jointly minimizes the KL divergence from the proposal to the EBM, and the negative log-likelihood of the EBM. At test-time, we can then employ importance sampling with the trained proposal to efficiently evaluate the learned EBM and produce stand-alone predictions. Furthermore, we utilize our derived training objective to learn mixture density networks (MDNs) with a jointly trained energy-based teacher, consistently outperforming conventional MDN training on four real-world regression tasks within computer vision. Code is available at https://github.com/fregu856/ebms_proposals.


title: Hydrodynamic superradiance in wave-mediated cooperative tunneling

abstract: Superradiance and subradiance occur in quantum optics when the emission rate of photons from multiple atoms is enhanced and diminished, respectively, owing to interaction between neighboring atoms. We here demonstrate a classical analog thereof in a theoretical model of droplets walking on a vibrating bath. Two droplets are confined to identical two-level systems, a pair of wells between which the drops may tunnel, joined by an intervening coupling cavity. The resulting classical superradiance is rationalized in terms of the system's non-Markovian, pilot-wave dynamics.


title: Adversarial robustness for latent models: Revisiting the robust-standard accuracies tradeoff

abstract: Over the past few years, several adversarial training methods have been proposed to improve the robustness of machine learning models against adversarial perturbations in the input. Despite remarkable progress in this regard, adversarial training is often observed to drop the standard test accuracy. This phenomenon has intrigued the research community to investigate the potential tradeoff between standard accuracy (a.k.a generalization) and robust accuracy (a.k.a robust generalization) as two performance measures. In this paper, we revisit this tradeoff for latent models and argue that this tradeoff is mitigated when the data enjoys a low-dimensional structure. In particular, we consider binary classification under two data generative models, namely Gaussian mixture model and generalized linear model, where the features data lie on a low-dimensional manifold. We develop a theory to show that the low-dimensional manifold structure allows one to obtain models that are nearly optimal with respect to both, the standard accuracy and the robust accuracy measures. We further corroborate our theory with several numerical experiments, including Mixture of Factor Analyzers (MFA) model trained on the MNIST dataset.


title: Missing the Point: Non-Convergence in Iterative Imputation Algorithms

abstract: Iterative imputation is a popular tool to accommodate missing data. While it is widely accepted that valid inferences can be obtained with this technique, these inferences all rely on algorithmic convergence. There is no consensus on how to evaluate the convergence properties of the method. Our study provides insight into identifying non-convergence in iterative imputation algorithms. We found that--in the cases considered--inferential validity was achieved after five to ten iterations, much earlier than indicated by diagnostic methods. We conclude that it never hurts to iterate longer, but such calculations hardly bring added value.


title: Optimal randomized classification trees

abstract: Classification and Regression Trees (CARTs) are off-the-shelf techniques in modern Statistics and Machine Learning. CARTs are traditionally built by means of a greedy procedure, sequentially deciding the splitting predictor variable(s) and the associated threshold. This greedy approach trains trees very fast, but, by its nature, their classification accuracy may not be competitive against other state-of-the-art procedures. Moreover, controlling critical issues, such as the misclassification rates in each of the classes, is difficult. To address these shortcomings, optimal decision trees have been recently proposed in the literature, which use discrete decision variables to model the path each observation will follow in the tree. Instead, we propose a new approach based on continuous optimization. Our classifier can be seen as a randomized tree, since at each node of the decision tree a random decision is made. The computational experience reported demonstrates the good performance of our procedure.


title: Wind Farm Icing Loss Forecast Pertinent to Winter Extremes

abstract: The 2021 Texas power crisis has highlighted the vulnerability of the power system under wind extremes, particularly with the increasing penetration of energy resources that depend on weather conditions (e.g., wind energy). The current wind power forecast models do not effectively consider the impact of such extreme weather events. In the present study, we provide a fast and robust statistical model developed using ten years of utility-scale turbine data at the Eolos Wind Energy Research Station to forecast the icing losses under such weather conditions. This model covers different cold climate impacts, including precipitation icing, frost contamination, and low-temperature effect. This model has been assessed using three large-scale (larger than 100 MW) wind farm data involving turbines with different capacities and from different manufacturers across multiple geographic regions. Notably, the model has been used to predict the wind power losses in entire Texas (larger than 91% of total wind installation in Texas) during 2021 Texas power crisis. The proposed model can be easily integrated into the existing wind farm and power grid operations, allowing the power system operators can develop more appropriate and pinpointed plans to balance the severe and sudden energy deficits and increase the system integrity under winter extremes.


title: Variational Probabilistic Multi-Hypothesis Tracking

abstract: This paper proposes a novel multi-target tracking (MTT) algorithm for scenarios with arbitrary numbers of measurements per target. We propose the variational probabilistic multi-hypothesis tracking (VPMHT) algorithm based on the variational Bayesian expectation-maximisation (VBEM) algorithm to resolve the MTT problem in the classic PMHT algorithm. With the introduction of variational inference, the proposed VPMHT handles track-loss much better than the conventional probabilistic multi-hypothesis tracking (PMHT) while preserving a similar or even better tracking accuracy. Extensive numerical simulations are conducted to demonstrate the effectiveness of the proposed algorithm.


title: Imprecise Subset Simulation

abstract: The objective of this work is to quantify the uncertainty in probability of failure estimates resulting from incomplete knowledge of the probability distributions for the input random variables. We propose a framework that couples the widely used Subset simulation (SuS) with Bayesian/information theoretic multi-model inference. The process starts with data used to infer probability distributions for the model inputs. Often such data sets are small. Multi-model inference is used to assess uncertainty associated with the model-form and parameters of these random variables in the form of model probabilities and the associated joint parameter probability densities. A sampling procedure is used to construct a set of equally probable candidate probability distributions and an optimal importance sampling distribution is determined analytically from this set. Subset simulation is then performed using this optimal sampling density and the resulting conditional probabilities are re-weighted using importance sampling. The result of this process are empirical probability distributions of failure probabilities that provide direct estimates of the uncertainty in failure probability estimates that result from inference on small data sets. The method is demonstrated to be both computationally efficient -- requiring only a single subset simulation and nominal cost of sample re-weighting -- and to provide reasonable estimates of the uncertainty in failure probabilities.


title: Probing the initial mass function of the first stars with transients

abstract: The emergence of the first, so-called Population III (Pop III), stars shaped early cosmic history in ways that crucially depends on their initial mass function (IMF). However, because of the absence of direct observational constraints, the detailed IMF remains elusive. Nevertheless, numerical simulations agree in broad terms that the first stars were typically massive and should often end their lives in violent, explosive deaths. These fates include extremely luminous pair-instability supernovae (PISNe) and bright gamma-ray bursts (GRBs), the latter arising from the collapse of rapidly rotating progenitor stars into black holes. These high-redshift transients are expected to be within the detection limits of upcoming space telescope missions, allowing to place effective constraints on the shape of the primordial IMF that is not easily accessible with other probes. This paper presents a framework to probe the Pop III IMF, utilizing the cosmological source densities of high-redshift PISNe and GRBs. Considering these transients separately could provide useful constraints on the Pop III IMF, but tighter bounds are obtainable by combining PISN and GRB counts. This combined diagnostic is more robust as it is independent of the underlying Pop III star formation rate density, an unknown prior. Future surveys promise to capture most high-redshift GRBs across the entire sky, but high-redshift PISN searches with future telescopes, e.g. Roman Space Telescope, will likely be substantially incomplete. Nevertheless, we demonstrate that even such lower bounds on the PISN count will be able to provide key constraints on the primordial IMF, in particular, if it is top-heavy or not.

categories: astro-ph.HE astro-ph.CO astro-ph.GA astro-ph.SR


title: Cosmological "constant" in a universe born in the metastable false vacuum state

abstract: The cosmological constant $\Lambda$ is a measure of the energy density of the vacuum. Therefore properties of the energy of the system in the metastable vacuum state reflect properties of $\Lambda = \Lambda(t)$. We analyze properties of the energy, $E(t)$, of a general quantum system in the metastable state in various phases of the decay process: In the exponential phase, in the transition phase between the exponential decay and the later phase, where decay law as a function of time $t$ is in the form of powers of $1/t$, and also in this last phase. We found that this energy having an approximate value resulting from the Weisskopf--Wigner theory in the exponential decay phase is reduced very fast in the transition phase to its asymptotic value $E(t) \simeq E_{min} + \alpha_{2}/t^{2}+\ldots$ in the late last phase of the decay process. (Here $E_{min}$ is the minimal energy of the system). This quantum mechanism reduces the energy of the system in the unstable state by a dozen or even several dozen orders or more. We show that if to assume that a universe was born in metastable false vacuum state then according to this quantum mechanism the cosmological constant $\Lambda$ can have a very great value resulting from the quantum field theory calculations in the early universe in the inflationary era, $\Lambda \simeq \Lambda_{qft}$, and then it can later be quickly reduced to the very, very small values.


title: Ultimate capacity limit of a multi-span link with phase-insensitive amplification

abstract: The Shannon capacity of a point-to-point link with an optimised configuration of optical amplifiers is compared with general detection strategies permitted by quantum mechanics. Results suggest that the primary application area of receivers based on such strategies may be unamplified short-distance links or free-space optical communication


title: QLP Data Release Notes 001: K2 + TESS Analysis

abstract: TESS will target the ecliptic plane in Sectors 42 - 46. These sectors overlap with campaigns from the K2 mission, providing a unique opportunity for multi-mission light curve analysis. This data release note describes the combined analysis of K2 and TESS light curves as part of the Quick-Look Pipeline (QLP) procedure, which processes light curves for all targets in TESS Full-Frame Images (FFIs) down to TESS magnitude T = 13.5. We describe updates to our codebase, and the planet transit search, candidate triage, and report generation that are affected by this combined analysis.


title: ReLAX: Reinforcement Learning Agent eXplainer for Arbitrary Predictive Models

abstract: Counterfactual examples (CFs) are one of the most popular methods for attaching post-hoc explanations to machine learning (ML) models. However, existing CF generation methods either exploit the internals of specific models or depend on each sample's neighborhood, thus they are hard to generalize for complex models and inefficient for large datasets. This work aims to overcome these limitations and introduces ReLAX, a model-agnostic algorithm to generate optimal counterfactual explanations. Specifically, we formulate the problem of crafting CFs as a sequential decision-making task and then find the optimal CFs via deep reinforcement learning (DRL) with discrete-continuous hybrid action space. Extensive experiments conducted on several tabular datasets have shown that ReLAX outperforms existing CF generation baselines, as it produces sparser counterfactuals, is more scalable to complex target models to explain, and generalizes to both classification and regression tasks. Finally, to demonstrate the usefulness of our method in a real-world use case, we leverage CFs generated by ReLAX to suggest actions that a country should take to reduce the risk of mortality due to COVID-19. Interestingly enough, the actions recommended by our method correspond to the strategies that many countries have actually implemented to counter the COVID-19 pandemic.


title: The Synchrony of Production & Escape: Half the Bright Ly$\alpha$ Emitters at $z\approx2$ have Lyman Continuum Escape Fractions $\approx50\%$

abstract: The ionizing photon escape fraction (LyC $f_{\rm{esc}}$) of star-forming galaxies is the single greatest unknown in the reionization budget. Stochastic sightline effects prohibit the direct separation of LyC leakers from non-leakers at significant redshift. Here we circumvent this uncertainty by inferring $f_{\rm{esc}}$ with resolved (R>4000) LyA profiles from the X-SHOOTER LyA survey at z=2 (XLS-z2). We select leakers ($f_{\rm{esc}}>20$%) and non-leakers ($f_{\rm{esc}}<5$%) from a representative sample of $>0.2 L^{*}$ LyA emitters (LAEs). With median stacked spectra of these subsets covering 1000-8000 {\AA} (rest-frame) we investigate the conditions for LyC $f_{\rm{esc}}$. We find the following differences between leakers vs. non-leakers: (i) strong nebular CIV and HeII emission vs. non-detections, (ii) O32~8.5 vs. ~3, (iii) Ha/Hb indicating no dust vs. E(B-V)~0.3, (iv) MgII emission close to the systemic velocity vs. redshifted, optically thick MgII, (v) LyA $f_{\rm{esc}}$ of ~50% vs. ~10%. The extreme EWs in leakers (O3+Hb~1100 {\AA}) constrain the characteristic timescale of LyC escape to ~3-10 Myr bursts when short-lived stars with the hardest ionizing spectra shine. The defining traits of leakers -- extremely ionizing stellar populations, low column densities, a dust-free, high ionization state ISM -- occur simultaneously in the $f_{\rm{esc}}>20\%$ stack, suggesting they are causally connected, and motivating why indicators like O32 may suffice to constrain $f_{\rm{esc}}$ at z>6 with JWST. The leakers comprise half our sample, have a median LyC $f_{\rm{esc}}$~50%, and an ionising production efficiency $\log({\xi_{\rm{ion}}/\rm{Hz\ erg^{-1}}})$~25.9. These results show LAEs -- the type of galaxies rare at z=2, but that become the norm at higher redshift -- are highly efficient ionizers, with extreme $\xi_{\rm{ion}}$ and prolific $f_{\rm{esc}}$ occurring in sync. (ABRIDGED)


title: Gate-tunable heavy fermion quantum criticality in a moir\'e Kondo lattice

abstract: We propose a realization of Kondo-lattice physics in moir\'e superlattices at the interface between a WX$_2$ homobilayer and MoX$_2$ monolayer (where X=S,Se). Under appropriate gating conditions, the interface-WX$_2$-layer forms a triangular lattice of local moments that couple to itinerant electrons in the other WX$_2$-layer via a gate-tunable Kondo exchange interaction. Using a parton mean-field approach we identify a range of twist-angles which support a gate-tuned quantum phase transition between a heavy-fermion liquid with large anomalous Hall conductance and a fractionalized chiral spin-liquid coexisting with a light Fermi liquid, and describe experimental signatures to distinguish among competing theoretical scenarios.


title: The impact of natal kicks on galactic r-process enrichment by neutron star mergers

abstract: We study galactic enrichment with rapid neutron capture (r-process) elements in cosmological, magnetohydrodynamical simulations of a Milky Way-mass galaxy. We include a variety of enrichment models, based on either neutron star mergers or a rare class of core-collapse supernova as sole r-process sources. For the first time in cosmological simulations, we implement neutron star natal kicks on-the-fly to study their impact. With kicks, neutron star mergers are more likely to occur outside the galaxy disc, but how far the binaries travel before merging also depends on the kick velocity distribution and shape of the delay time distribution for neutron star mergers. In our fiducial model, the median r-process abundance ratio is somewhat lower and the trend with metallicity is slightly steeper when kicks are included. In a model 'optimized' to better match observations, with a higher rate of early neutron star mergers, the median r-process abundances are fairly unaffected by kicks. In both models, the scatter in r-process abundances is much larger with natal kicks, especially at low metallicity, giving rise to more r-process enhanced stars. We experimented with a range of kick velocities and find that with lower velocities, the scatter is reduced, but still larger than without natal kicks. We discuss the possibility that the observed scatter in r-process abundances is predominantly caused by natal kicks removing the r-process sources far from their birth sites, making enrichment more inhomogeneous, rather than the usual interpretation that the scatter is set by the rarity of its production source.


title: Molecular Chemistry for Dark Matter II: Recombination, Molecule Formation, and Halo Mass Function in Atomic Dark Matter

abstract: Dissipative dark matter predicts rich observable phenomena that can be tested with future large-scale structure surveys. As a specific example, we study atomic dark matter, consisting of a heavy particle and a light particle charged under a dark electromagnetism. In particular, we calculate the cosmological evolution of atomic dark matter focusing on dark recombination and dark-molecule formation. We have obtained the relevant interaction-rate coefficients by re-scaling the rates for normal hydrogen, and evolved the abundances for ionized, atomic, and molecular states using a modified version of Recfast++ (which we have released publicly at https://github.com/jamesgurian/RecfastJulia). We also provide an analytical approximation for the final abundances. We then calculate the effects of the atomic dark matter on the linear power spectrum, which enter through a dark-photon diffusion and dark acoustic oscillations. At the formation time, the atomic dark matter model suppresses halo abundances on scales smaller than the diffusion scale, just like the warm dark matter models suppress the abundance below the free-streaming scale. The subsequent evolution with radiative cooling, however, will alter the halo mass function further.


title: A universal tripartite entanglement signature of ungappable edge states

abstract: Gapped two-dimensional topological phases can feature ungappable edge states which are robust even in the absence of protecting symmetries. In this work we show that a multipartite entanglement measure recently proposed in the context of holography, the Markov gap, provides a universal diagnostic of ungappable edge states. Defined as a difference of the reflected entropy and mutual information $h(A:B) = S_R(A:B) - I(A:B)$ between two parties, we argue that for $A,B$ being adjacent subregions in the bulk, $h=\frac{c_+}{3}\log 2$, where $c_+$ is the minimal total central charge of the boundary theory. As evidence, we prove that $h=0$ for string-net models, and numerically verify that $h=\frac{|C|}{3}\log 2$ for a Chern-$C$ insulator. Our work establishes a unique bulk entanglement criteria for the presence of a conformal field theory on the boundary.


title: The dark matter halo masses of elliptical galaxies as a function of observationally robust quantities

abstract: Context. The assembly history of the stellar component of a massive elliptical galaxy is closely related to that of its dark matter halo. Measuring how the properties of galaxies correlate with their halo mass can help understand their evolution. Aims. We investigate how the dark matter halo mass of elliptical galaxies varies as a function of their properties, using weak gravitational lensing observations. To minimise the chances of biases, we focus on galaxy properties that can be determined robustly: the surface brightness profile and the colour. Methods. We selected 2409 central massive elliptical galaxies from the SDSS spectroscopic sample. We first measured their surface brightness profile and colours by fitting Sersic models to photometric data from the Kilo-Degree Survey (KiDS). We fitted their halo mass distribution as a function of redshift, rest-frame $r-$band luminosity, half-light radius and rest-frame $u-g$ colour, using KiDS weak lensing data and a Bayesian hierarchical approach. For the sake of robustness to assumptions on the large-radii behaviour of the surface brightness, we repeated the analysis replacing total luminosity and half-light radius with the luminosity within a 10~kpc aperture, $L_{r,10}$, and the light-weighted surface brightness slope, $\Gamma_{10}$. Results. We did not detect any correlation between halo mass and either half-light radius or colour, at fixed redshift and luminosity. Conclusions. Our results indicate that the average star formation efficiency of massive elliptical galaxies has little dependence on their final size or colour. This suggests that the origin of the diversity in the size and colour distribution of these objects lies with properties other than halo mass.


title: (Re)Solving Reionization with Ly{\alpha}: How Bright Ly{\alpha} Emitters account for the $z\approx2-8$ Cosmic Ionizing Background

abstract: The cosmic ionizing emissivity from star-forming galaxies has long been anchored to UV luminosity functions. Here we introduce an emissivity framework based on Ly$\alpha$ emitters (LAEs), which naturally hones in on the subset of galaxies responsible for the ionizing background due to the intimate connection between the production and escape of Ly$\alpha$ and LyC photons. Using constraints on the escape fractions of bright LAEs ($L_{\rm{Ly\alpha}}>0.2 L^{*}$) at $z\approx2$ obtained from resolved Ly$\alpha$ profiles, and arguing for their redshift-invariance, we show that: (i) quasars and LAEs together reproduce the relatively flat emissivity at $z\approx2-6$, which is non-trivial given the strong evolution in both the star-formation density and quasar number density at these epochs and (ii) LAEs produce late and rapid reionization between $z\approx6-9$ under plausible assumptions. Within this framework, the $>10\times$ rise in the UV population-averaged $f_{\rm{esc}}$ between $z\approx3-7$ naturally arises due to the same phenomena that drive the growing Ly$\alpha$ emitter fraction with redshift. Generally, a LAE dominated emissivity yields a peak in the distribution of the ionizing budget with UV luminosity as reported in latest simulations. Using our adopted parameters ($f_{\rm{esc}}=50\%$, $\xi_{\rm{ion}}=10^{25.9}$ Hz erg$^{-1}$ for half the bright LAEs), a highly ionizing minority of galaxies with $M_{\rm UV}<-17$ accounts for the entire ionizing budget from star-forming galaxies. Rapid flashes of LyC from such rare galaxies produce a "disco" ionizing background. We conclude proposing tests to further develop our suggested Ly$\alpha$-anchored formalism.


title: Jet Launching from Merging Magnetized Binary Neutron Stars with Realistic Equations of State

abstract: We perform general relativistic, magnetohydrodynamic (GRMHD) simulations of binary neutron stars in quasi-circular orbit that merge and undergo delayed or prompt collapse to a black hole (BH). The stars are irrotational and modeled using an SLy or an H4 nuclear equation of state. To assess the impact of the initial magnetic field configuration on jet launching, we endow the stars with a purely poloidal magnetic field that is initially unimportant dynamically and is either confined to the stellar interior or extends from the interior into the exterior as in typical pulsars. Consistent with our previous results, we find that only the BH + disk remnants originating from binaries that form hypermassive neutron stars (HMNSs) and undergo delayed collapse can drive magnetically-powered jets. We find that the closer the total mass of the binary is to the threshold value for prompt collapse, the shorter is the time delay between the gravitational wave peak amplitude and jet launching. This time delay also strongly depends on the initial magnetic field configuration. We also find that seed magnetic fields confined to the stellar interior can launch a jet over $\sim 25\,\rm ms$ later than those with pulsar-like magnetic fields. The lifetime of the jet [$\Delta t\lesssim 150\,\rm ms$] and its outgoing Poynting luminosity [$L_{\rm EM}\sim 10^{52\pm 1}\rm erg/s$] are consistent with typical short gamma-ray burst central engine lifetimes, as well as with the Blandford--Znajek mechanism for launching jets and their associated Poynting luminosities. Our numerical results also suggest that the dynamical ejection of matter can be enhanced by the magnetic field. Therefore, GRMHD studies are required to fully understand kilonova signals from GW170818-like events.


title: Four flavours, triality and bimodular forms

abstract: We consider $\mathcal{N}=2$ supersymmetric $\text{SU}(2)$ gauge theory with $N_f=4$ massive hypermultiplets. The duality group of this theory contains transformations acting on the UV-coupling $\tau_{\text{UV}}$ as well as on the running coupling $\tau$. We establish that subgroups of the duality group act separately on $\tau_{\text{UV}}$ and $\tau$, while a larger group acts simultaneously on $\tau_{\text{UV}}$ and $\tau$. For special choices of the masses, we find that the duality groups can be identified with congruence subgroups of $\text{SL}(2,\mathbb Z)$. We demonstrate that in such cases, the order parameters are instances of bimodular forms with arguments $\tau$ and $\tau_{\text{UV}}$. Since the UV duality group of the theory contains the triality group of outer automorphisms of the flavour symmetry $\text{SO}(8)$, the duality action gives rise to an orbit of mass configurations. Consequently, the corresponding order parameters combine to vector-valued bimodular forms with $\text{SL}(2,\mathbb Z)$ acting simultaneously on the two couplings.


title: From EMBER to FIRE: predicting high resolution baryon fields from dark matter simulations with Deep Learning

abstract: Hydrodynamic simulations provide a powerful, but computationally expensive, approach to study the interplay of dark matter and baryons in cosmological structure formation. Here we introduce the EMulating Baryonic EnRichment (EMBER) Deep Learning framework to predict baryon fields based on dark-matter-only simulations thereby reducing computational cost. EMBER comprises two network architectures, U-Net and Wasserstein Generative Adversarial Networks (WGANs), to predict two-dimensional gas and HI densities from dark matter fields. We design the conditional WGANs as stochastic emulators, such that multiple target fields can be sampled from the same dark matter input. For training we combine cosmological volume and zoom-in hydrodynamical simulations from the Feedback in Realistic Environments (FIRE) project to represent a large range of scales. Our fiducial WGAN model reproduces the gas and HI power spectra within 10% accuracy down to ~10 kpc scales. Furthermore, we investigate the capability of EMBER to predict high resolution baryon fields from low resolution dark matter inputs through upsampling techniques. As a practical application, we use this methodology to emulate high-resolution HI maps for a dark matter simulation of a L=100 Mpc/h comoving cosmological box. The gas content of dark matter haloes and the HI column density distributions predicted by EMBER agree well with results of large volume cosmological simulations and abundance matching models. Our method provides a computationally efficient, stochastic emulator for augmenting dark matter only simulations with physically consistent maps of baryon fields.


title: Molecular Chemistry for Dark Matter III: DarkKROME

abstract: Dark matter that is dissipative may cool sufficiently to form compact objects, including black holes. Determining the abundance and mass spectrum of those objects requires an accurate model of the chemistry relevant for the cooling of the dark matter gas. Here we introduce a chemistry tool for dark matter, DarkKROME, an extension of the KROME software package. DarkKROME is designed to include all atomic and molecular processes relevant for dark matter with two unequal-mass fundamental fermions, interacting via a massless-photon mediated $U(1)$ force. We use DarkKROME to perform one-zone collapse simulations and study the evolution of temperature-density phase diagrams for various dark-sector parameters.


title: Strong lensing in UNIONS: Toward a pipeline from discovery to modeling

abstract: We present a search for galaxy-scale strong gravitational lenses in the initial 2 500 square degrees of the Canada-France Imaging Survey (CFIS). We designed a convolutional neural network (CNN) committee that we applied to a selection of 2 344 002 exquisite-seeing $r$-band images of color-selected luminous red galaxies (LRGs). Our classification uses a realistic training set where the lensing galaxies and the lensed sources are both taken from real data, namely the CFIS $r$-band images themselves and the Hubble Space Telescope (HST). A total of 9 460 candidates obtain a score above 0.5 with the CNN committee. After a visual inspection of the candidates, we find a total of 133 lens candidates, of which 104 are completely new. The set of false positives mainly contains ring, spiral, and merger galaxies, and to a lesser extent galaxies with nearby companions. We classify 32 of the lens candidates as secure lenses and 101 as maybe lenses. For the 32 highest quality lenses, we also fit a singular isothermal ellipsoid mass profile with external shear along with an elliptical Sersic profile for the lens and source light. This automated modeling step provides distributions of properties for both sources and lenses that have Einstein radii in the range $0.5\arcsec<\theta_E<2.5\arcsec$. Finally, we introduce a new lens and/or source single-band deblending algorithm based on auto-encoder representation of our candidates. This is the first time an end-to-end lens-finding and modeling pipeline is assembled together, in view of future lens searches in a single band, as will be possible with Euclid.


title: Modular invariance, misalignment and finiteness in non-supersymmetric strings

abstract: In this article we show that finite perturbative corrections in non-supersymmetric strings can be understood via an interplay between modular invariance and misaligned supersymmetry. While modular invariance is known to be crucial in closed-string models, its presence and role for open strings is more subtle. Nevertheless, we argue that it leads to cancellations in physical quantities such as the one-loop cosmological constant and prevents them from diverging. In particular, we show that if the sector-averaged number of states does not grow exponentially, as predicted by misaligned supersymmetry, all exponential divergences in the one-loop cosmological constant cancel out as well. To account for the absence of power-law divergences, instead, we need to resort to the modular structure of the partition function. We finally comment on the presence of misaligned supersymmetry in the known 10-dimensional tachyon-free non-supersymmetric string theories.


title: Discovery of a double detonation thermonuclear supernova progenitor

abstract: We present the discovery of a new double detonation progenitor system consisting of a hot subdwarf B (sdB) binary with a white dwarf companion with an P=76.34179(2) min orbital period. Spectroscopic observations are consistent with an sdB star during helium core burning residing on the extreme horizontal branch. Chimera light curves are dominated by ellipsoidal deformation of the sdB star and a weak eclipse of the companion white dwarf. Combining spectroscopic and light curve fits we find a low mass sdB star, $M_{\rm sdB}=0.383\pm0.028$ M$_\odot$ with a massive white dwarf companion, $M_{\rm WD}=0.725\pm0.026$ M$_\odot$. From the eclipses we find a blackbody temperature for the white dwarf of 26,800 K resulting in a cooling age of $\approx$25 Myrs whereas our MESA model predicts an sdB age of $\approx$170 Myrs. We conclude that the sdB formed first through stable mass transfer followed by a common envelope which led to the formation of the white dwarf companion $\approx$25 Myrs ago. Using the MESA stellar evolutionary code we find that the sdB star will start mass transfer in $\approx$6 Myrs and in $\approx$60 Myrs the white dwarf will reach a total mass of $0.92$ M$_\odot$ with a thick helium layer of $0.17$ M$_\odot$. This will lead to a detonation that will likely destroy the white dwarf in a peculiar thermonuclear supernova. PTF1 2238+7430 is only the second confirmed candidate for a double detonation thermonuclear supernova. Using both systems we estimate that at least $\approx$1% of white dwarf thermonuclear supernovae originate from sdB+WD binaries with thick helium layers, consistent with the small number of observed peculiar thermonuclear explosions.


title: Background of radio photons from primordial black holes

abstract: We compute the isotropic radiation background due to Hawking emission from primordial black holes (PBHs), and examine if this background is a viable option in explaining the excess radiowave background observed by the ARCADE2 and LWA1 experiments at $\lesssim 1\,$GHz. We find that even under the extreme assumption that all of the dark matter is in the form of PBHs, the radio brightness temperature induced by Hawking evaporation of PBHs is $\mathcal{O}(10^{-46})\,$K, highly subdominant compared to the cosmic microwave background. The main reason for this is that for PBHs in the mass range $\sim10^{12}$-$10^{14}\,$kg, which can be constrained by Hawking emission, the spectrum peaks at $10^7$ to $10^5\,$eV. As the Hawking spectrum is power law suppressed towards lower energies, negligible flux of $\mu$eV photons is obtained. The peak of the Hawking spectrum shifts to lower energies for higher masses, but the number density is low and so is the specific intensity. Because Hawking emission from PBHs is thus unable to explain the observed excess radio background, we also consider the alternative possibility of radio emission from gas accretion onto supermassive PBHs. These PBHs can readily produce strong radio emission that could easily explain the ARCADE2/LWA1 excess.


title: The Dark Side of Pluto

abstract: During its departure from Pluto, New Horizons used its LORRI camera to image a portion of Pluto's southern hemisphere that was in a decades-long seasonal winter darkness, but still very faintly illuminated by sunlight reflected by Charon. Recovery of this faint signal was technically challenging. The bright ring of sunlight forward-scattered by haze in the Plutonian atmosphere encircling the nightside hemisphere was severely overexposed, defeating the standard smeared-charge removal required for LORRI images. Reconstruction of the overexposed portions of the raw images, however, allowed adequate corrections to be accomplished. The small solar elongation of Pluto during the departure phase also generated a complex scattered-sunlight background in the images that was three orders of magnitude stronger than the estimated Charon-light flux (the Charon-light flux is similar to the flux of moonlight on Earth a few days before first quarter). A model background image was constructed for each Pluto image based on principal component analysis (PCA) applied to an ensemble of scattered-sunlight images taken at identical Sun-spacecraft geometry to the Pluto images. The recovered Charon-light image revealed a high-albedo region in the southern hemisphere. We argue that this may be a regional deposit of N_2 or CH_4 ice. The Charon-light image also shows that the south polar region currently has markedly lower albedo than the north polar region of Pluto, which may reflect the sublimation of N_2 ice or the deposition of haze particulates during the recent southern summer.


title: The SSA22 HI Tomography Survey (SSA22-HIT). I. Data Set and Compiled Redshift Catalog

abstract: We conducted a deep spectroscopic survey, named SSA22-HIT, in the SSA22 field with the DEep Imaging MultiObject Spectrograph (DEIMOS) on the Keck telescope, designed to tomographically map high-z HI gas through analysis of Lya absorption in background galaxies' spectra. In total, 198 galaxies were spectroscopically confirmed at 2.5 < z < 6 with a few low-z exceptions in the 26 x 15 arcmin^2 area, of which 148 were newly determined in this study. Our redshift measurements were merged with previously confirmed redshifts available in the 34 x 27 arcmin^2 area of the SSA22 field. This compiled catalog containing 730 galaxies of various types at z > 2 is useful for various applications, and it is made publicly available. Our SSA22-HIT survey has increased by approximately twice the number of spectroscopic redshifts of sources at z > 3.2 in the observed field. From a comparison with publicly available redshift catalogs, we show that our compiled redshift catalog in the SSA22 field is comparable to those among major extragalactic survey fields in terms of a combination of wide area and high surface number density of objects at z > 2. About 40 % of the spectroscopically confirmed objects in SSA22-HIT show reasonable quality of spectra in the wavelengths shorter than Lya when a sufficient amount of smoothing is adopted. Our data set enables us to make the HI tomographic map at z > 3, which we present in a parallel study.


title: Renormalization group in quantum critical theories with Harris-marginal disorder

abstract: We develop a renormalization group for weak Harris-marginal disorder in otherwise strongly interacting quantum critical theories, focusing on systems which have emergent conformal invariance. Using conformal perturbation theory, we argue that previously proposed random lines of fixed points with Lifshitz scaling in fact flow towards other universal fixed points, and this flow is captured by a "one-loop" analysis. Our approach appears best controlled in theories with only a few operators with low scaling dimension. In this regime, we compare our predictions for the flow of disorder to holographic models, and find complete agreement.


title: $\mathcal{O}(\alpha_s^3)$ corrections to semileptonic $b \to c$ decays in the heavy daughter approximation

abstract: We present our recent calculation of the third order corrections to the semileptonic $b \to c$ and the muon decays. The calculation has been performed in an expansion around the limit $m_c \sim m_b$, but shows decent convergence even for $m_c=0$ from which the contribution to the muon decay can be extracted. For the semileptonic $b \to c$ decay we find large perturbative corrections in the on-shell scheme which can be significantly reduced by changing to the kinetic scheme for the bottom quark mass. These results are important input for the inclusive determination of $|V_{cb}|$ and the Fermi coupling constant $G_F$.


title: Multipartitioning topological phases by vertex states and quantum entanglement

abstract: We discuss multipartitions of the gapped ground states of (2+1)-dimensional topological liquids into three (or more) spatial regions that are adjacent to each other and meet at points. By considering the reduced density matrix obtained by tracing over a subset of the regions, we compute various correlation measures, such as entanglement negativity, reflected entropy, and associated spectra. We utilize the bulk-boundary correspondence to show that such multipartitions can be achieved by using what we call vertex states in (1+1)-dimensional conformal field theory -- these are a type of state used to define an interaction vertex in string field theory and can be thought of as a proper generalization of conformal boundary states. This approach allows an explicit construction of the reduced density matrix near the entangling boundaries. We find the fingerprints of topological liquid in these quantities, such as (universal pieces in) the scaling of the entanglement negativity, and a non-trivial distribution of the spectrum of the partially transposed density matrix. For reflected entropy, we test the recent claim that states the difference between reflected entropy and mutual information is given, once short-range correlations are properly removed, by $(c/3)\ln 2$ where $c$ is the central charge of the topological liquid that measures ungappable edge degrees of freedom. As specific examples, we consider topological chiral $p$-wave superconductors and Chern insulators. We also study a specific lattice fermion model realizing Chern insulator phases and calculate the correlation measures numerically, both in its gapped phases and at critical points separating them.


title: How to Quantify Polarization in Models of Opinion Dynamics

abstract: It is widely believed that society is becoming increasingly polarized around important issues, a dynamic that does not align with common mathematical models of opinion formation in social networks. In particular, measures of polarization based on opinion variance are known to decrease over time in frameworks such as the popular DeGroot model. Complementing recent work that seeks to resolve this apparent inconsistency by modifying opinion models, we instead resolve the inconsistency by proposing changes to how polarization is quantified. We present a natural class of group-based polarization measures that capture the extent to which opinions are clustered into distinct groups. Using theoretical arguments and empirical evidence, we show that these group-based measures display interesting, non-monotonic dynamics, even in the simple DeGroot model. In particular, for many natural social networks, group-based metrics can increase over time, and thereby correctly capture perceptions of increasing polarization. Our results build on work by DeMarzo et al., who introduced a group-based polarization metric based on ideological alignment. We show that a central tool from that work, a limit analysis of individual opinions under the DeGroot model, can be extended to the dynamics of other group-based polarization measures, including established statistical measures like bimodality. We also consider local measures of polarization that operationalize how polarization is perceived in a network setting. In conjunction with evidence from prior work that group-based measures better align with real-world perceptions of polarization, our work provides formal support for the use of these measures in place of variance-based polarization in future studies of opinion dynamics.


title: New constraints on the minimum mass for thermonuclear lithium burning in brown dwarfs

abstract: The theory of substellar evolution predicts that there is a sharp mass boundary between lithium and non-lithium brown dwarfs, not far below the substellar-mass limit. The imprint of thermonuclear burning is carved on the surface lithium abundance of substellar-mass objects during the first few hundred million years of their evolution, leading to a sharp boundary between lithium and non-lithium brown dwarfs, so-called, the lithium test. New optical spectroscopic observations of the binaries DENIS+J063001.4-184014 and DENIS+J225210.7-173013 obtained using the 10.4-m Gran Telescopio de Canarias are reported here. They allow us to re-determine their combined optical spectral types (M9.5 and L6.5, respectively) and to search for the presence of the LiI resonance doublet. The non detection of the LiI feature in the combined spectrum of DENIS\,J063001.4$-$184014AB is converted into estimates for the depletion of lithium in the individual components of this binary system. In DENIS\,J225210.7$-$173013AB we report the detection of a weak LiI feature which we tentatively ascribe as arising from the contribution of the T3.5-type secondary. Combining our results with data for seven other brown dwarf binaries in the literature treated in a self-consistent way, we confirm that there is indeed a sharp transition in mass for lithium depletion in brown dwarfs, as expected from theoretical calculations. We estimate such mass boundary is observationally located at 51.48$^{+0.22}_{-4.00}$ $M_\mathrm{Jup}$, which is lower than the theoretical determinations.


title: Temperature-dependent periodicity of the persistent current in strongly interacting systems

abstract: The persistent current in small isolated rings enclosing magnetic flux is the current circulating in equilibrium in the absence of an external excitation. While initially studied in superconducting and normal metals, recently, atomic persistent currents have been generated in ultracold gases spurring a new wave of theoretical investigations. Nevertheless, our understanding of the persistent currents in interacting systems is far from complete, especially at finite temperatures. Here we consider the fermionic one-dimensional Hubbard model and show that in the strong-interacting limit, the current can change its flux period and sign (diamagnetic or paramagnetic) as a function of temperature, features that cannot be explained within the single-particle or Luttinger liquid techniques. Also, the magnitude of the current can counterintuitively increase with temperature, in addition to presenting different rates of decay depending on the polarization of the system. Our work highlights the properties of the strongly-interacting multi-component systems which are missed by conventional approximation techniques, but can be important for the interpretation of experiments on persistent currents in ultracold gases.


title: Law Smells: Defining and Detecting Problematic Patterns in Legal Drafting

abstract: Building on the computer science concept of code smells, we initiate the study of law smells, i.e., patterns in legal texts that pose threats to the comprehensibility and maintainability of the law. With five intuitive law smells as running examples - namely, duplicated phrase, long element, large reference tree, ambiguous syntax, and natural language obsession -, we develop a comprehensive law smell taxonomy. This taxonomy classifies law smells by when they can be detected, which aspects of law they relate to, and how they can be discovered. We introduce text-based and graph-based methods to identify instances of law smells, confirming their utility in practice using the United States Code as a test case. Our work demonstrates how ideas from software engineering can be leveraged to assess and improve the quality of legal code, thus drawing attention to an understudied area in the intersection of law and computer science and highlighting the potential of computational legal drafting.


title: The MASSIVE SURVEY XVI. The Stellar Initial Mass Function in the Center of MASSIVE Early-Type Galaxies

abstract: The stellar initial mass function (IMF) is a fundamental property in the measurement of stellar masses and galaxy star formation histories. In this work we focus on the most massive galaxies in the nearby universe $\log(M_{\star}/M_{\odot})>11.2$. We obtain high quality Magellan/LDSS-3 long slit spectroscopy with a wide wavelength coverage of $0.4\mu{\rm m}-1.01\mu{\rm m}$ for 41 early-type galaxies (ETGs) in the MASSIVE survey, and derive high S/N spectra within an aperture of $R_{\rm e}/8$. Using detailed stellar synthesis models, we constrain the elemental abundances and stellar IMF of each galaxy through full spectral modeling. All the ETGs in our sample have an IMF that is steeper than a Milky Way (Kroupa) IMF. The best-fit IMF mismatch parameter, $\alpha_{\rm IMF}=(M/L)/(M/L)_{\rm MW}$, ranges from 1.12 to 3.05, with an average of $\langle \alpha_{\rm IMF} \rangle=1.84$, suggesting that on average, the IMF is more bottom-heavy than Salpeter. Comparing the estimated stellar mass with the dynamical mass, we find that most galaxies have stellar masses smaller than their dynamical masses within the $1\sigma$ uncertainty. We complement our sample with lower-mass galaxies from the literature, and confirm that $\log(\alpha_{\rm IMF})$ is positively correlated with $\log(\sigma)$, $\log(M_{\star})$, and $\log(M_{\rm dyn})$. The IMF in the centers of more massive ETGs is more bottom-heavy. In addition, we find that $\log(\alpha_{\rm IMF})$ is positively correlated with both [Mg/Fe] and the estimated total metallicity [Z/H]. We find suggestive evidence that the effective stellar surface density $\Sigma_{\rm Kroupa}$ might be responsible for the variation of $\alpha_{\rm IMF}$. We conclude that $\sigma$, [Mg/Fe] and [Z/H] are the primary drivers of the global stellar IMF variation.


title: Local, Interactive, and Actionable: a Pandemic Behavioral Nudge

abstract: The informational environment surrounding the Covid-19 pandemic has been widely recognized as fragmented, politicized, and complex [1]. This has resulted in polarized public views regarding the veracity of scientific communication, the severity of the threat posed by the virus, and the necessity of nonpharmaceutical interventions (NPIs) which can slow the spread of infections [2]. This paper describes CovidCommitment.org, an effort toward enhancing NPI adoption through the combination of a social behavioral commitment device and interactive map-based visualizations of localized infection data as tabulated via a 1-hourdrive-time isochrone. This paper describes the system design and presents a preliminary analysis of user behavior within the system.


title: Improving Robustness of Malware Classifiers using Adversarial Strings Generated from Perturbed Latent Representations

abstract: In malware behavioral analysis, the list of accessed and created files very often indicates whether the examined file is malicious or benign. However, malware authors are trying to avoid detection by generating random filenames and/or modifying used filenames with new versions of the malware. These changes represent real-world adversarial examples. The goal of this work is to generate realistic adversarial examples and improve the classifier's robustness against these attacks. Our approach learns latent representations of input strings in an unsupervised fashion and uses gradient-based adversarial attack methods in the latent domain to generate adversarial examples in the input domain. We use these examples to improve the classifier's robustness by training on the generated adversarial set of strings. Compared to classifiers trained only on perturbed latent vectors, our approach produces classifiers that are significantly more robust without a large trade-off in standard accuracy.


title: Bicontinuous soft solids with a gradient in channel size designed for energy storage applications

abstract: We present examples of bicontinuous interfacially jammed emulsion gels ("bijels") with a designed gradient in the channel size along the sample. These samples are created by quenching binary fluids which have a gradient in particle concentration along the sample, since the channel size is determined by the local particle concentration. A gradient in local particle concentration is achieved using a two-stage loading process, with different particle volume fractions in each stage. Confocal microscopy and image analysis were used to quantitatively measure the channel size of the bijels. Bijels with a gradient in channel size of up to 2.8%/mm have been created. Such tailored soft materials could act as templates for energy materials optimised for both high ionic transport rates (high power) and high interfacial area (high energy density), potentially making them useful in novel energy applications.


title: The building up of observed stellar scaling relations of massive galaxies and the connection to black hole growth in the TNG50 simulation

abstract: [abridged] We study how mock-observed stellar morphological and structural properties of massive galaxies are built up between $z=0.5$ and $z=3$ in the TNG50 cosmological simulation. We generate mock images with the properties of the CANDELS survey and derive Sersic parameters and optical rest-frame morphologies as usually done in the observations. Overall, the simulation reproduces the observed evolution of the abundances of different galaxy morphological types of star-forming and quiescent galaxies. The $\log{M_*}-\log R_e$ and $\log{M_*}-\log\Sigma_1$ relations of the simulated star-forming and quenched galaxies also match the observed slopes and zeropoints to within 1-$\sigma$. In the simulation, galaxies increase their observed central stellar mass density ($\Sigma_1$) and transform in morphology from irregular/clumpy systems to normal Hubble-type systems in the Star Formation Main Sequence at a characteristic stellar mass of $\sim 10^{10.5}~M_\odot$. This morphological transformation is connected to the activity of the central Super Massive Black Holes (SMBHs). At low stellar masses ($10^9$ < $M_*/M_\odot$ < $10^{10}$) SMBHs grow rapidly, while at higher mass SMBHs switch into the kinetic feedback mode and grow more slowly. During this low-accretion phase, SMBH feedback leads to the quenching of star-formation, along with a simultaneous growth in $\Sigma_1$. More compact massive galaxies grow their SMBHs faster than extended ones of the same mass and end up quenching earlier. In the TNG50 simulation, SMBHs predominantly grow via gas accretion before galaxies quench, and $\Sigma_1$ increases substantially after SMBH growth slows down. The simulation predicts therefore that quiescent galaxies have higher $\Sigma_1$ values than star-forming galaxies for the same SMBH mass, which disagrees with alternative models, and may potentially be in tension with some observations.


title: Modeling the scattering polarization of the solar Ca i 4227 {\AA} line with angle-dependent partial frequency redistribution

abstract: Context. The correct modeling of the scattering polarization signals observed in several strong resonance lines requires taking partial frequency redistribution (PRD) phenomena into account. Aims. This work aims at assessing the impact and the range of validity of the angle-averaged AA approximation with respect to the general angle-dependent (AD) treatment of PRD effects in the modeling of scattering polarization in strong resonance lines, with focus on the solar Ca i 4227 {\AA} line. Methods. Spectral line polarization is modeled by solving the radiative transfer problem for polarized radiation, under nonlocal thermodynamic equilibrium conditions, taking PRD effects into account, in static one-dimensional semi-empirical atmospheric models presenting arbitrary magnetic fields. The problem is solved through a two-step approach. In step 1, the problem is solved for intensity only, considering a multi-level atom. In step 2, the problem is solved including polarization, considering a two-level atom with an unpolarized and infinitely sharp lower level, and fixing the lower level population calculated at step 1. Results. The results for the Ca i 4227 {\AA} line show a good agreement between the AA and AD calculations for the Q/I and U/I wings signals. However, AA calculations reveal an artificial trough in the line-core peak of the linear polarization profiles, whereas AD calculations show a sharper peak in agreement with observations. Conclusions. An AD treatment of PRD effects is essential to correctly model the line-core peak of the scattering polarization signal of the Ca i 4227 {\AA} line. By contrast, in the considered static case, the AA approximation seems to be suitable to model the wing scattering polarization lobes and their magnetic sensitivity through magneto-optical effects.


title: A Reinforcement Learning Approach to Parameter Selection for Distributed Optimal Power Flow

abstract: With the increasing penetration of distributed energy resources, distributed optimization algorithms have attracted significant attention for power systems applications due to their potential for superior scalability, privacy, and robustness to a single point-of-failure. The Alternating Direction Method of Multipliers (ADMM) is a popular distributed optimization algorithm; however, its convergence performance is highly dependent on the selection of penalty parameters, which are usually chosen heuristically. In this work, we use reinforcement learning (RL) to develop an adaptive penalty parameter selection policy for the AC optimal power flow (ACOPF) problem solved via ADMM with the goal of minimizing the number of iterations until convergence. We train our RL policy using deep Q-learning, and show that this policy can result in significantly accelerated convergence (up to a 59% reduction in the number of iterations compared to existing, curvature-informed penalty parameter selection methods). Furthermore, we show that our RL policy demonstrates promise for generalizability, performing well under unseen loading schemes as well as under unseen losses of lines and generators (up to a 50% reduction in iterations). This work thus provides a proof-of-concept for using RL for parameter selection in ADMM for power systems applications.


title: Volatile transport modeling on Triton with new observational constraints

abstract: Neptune's moon Triton shares many similarities with Pluto, including volatile cycles of N2, CH4 and CO, and represents a benchmark case for the study of surface-atmosphere interactions on volatile-rich KBOs. Within the context of New Horizons observations of Pluto as well as recent Earth-based observations of Triton, we adapt a Plutonian VTM to Triton, and test its ability to simulate its volatile cycles, thereby aiding our understanding of its climate. We present VTM simulations exploring the volatile cycles on Triton over long-term and seasonal timescales for varying model parameters. We explore what scenarios and model parameters allow for a best match of the available observations. In particular, our set of observational constraints include Voyager 2 observations, ground-based NIR (0.8 to 2.4 {\mu}m) disk-integrated spectra and the evolution of surface pressure as retrieved from stellar occultations. Our results show that Triton's poles act as cold traps for volatile ices and favor the formation of polar caps extending to lower latitudes through glacial flow. As previously evidenced by other VTMs, North-South asymmetries in surface properties can favor the development of one cap over the other. Our best-case simulations are obtained for a global reservoir of N2 ice thicker than 200 m and a bedrock thermal inertia larger than 500 SI. The large N2 ice reservoir implies a permanent N2 southern cap extending to the equator. Our results also suggest that a small permanent polar cap exists in the northern (currently winter) hemisphere if the internal heat flux remains radiogenic (< 3 mW m-2). Finally, we provide predictions for the evolution of ice distribution, surface pressure, CO and CH4 atmospheric mixing ratios in the next decades. We also model the thermal lightcurves of Triton in 2022, which serve as predictions for future JWST observations.


title: Measurements of heavy-flavor production as a function of multiplicity with ALICE at the LHC

abstract: In this contribution, the production of heavy-flavor hadrons as a function of multiplicity, via the study of the D-meson and heavy-flavor hadron decay leptons self-normalized yields in pp collisions at the center of mass energy $\sqrt{s} = 13$ TeV is discussed. Comparisons are made with similar measurements of J/$\psi$ at $\sqrt{s} = 13$ TeV and various model calculations. The $\mathrm{\Lambda_{c}^{+}/D^{0}}$ and $\mathrm{D^{+}_{s}/D^{0}}$ yield ratios in different multiplicity intervals in pp collisions at $\sqrt{s} = 13$ TeV are also reported. In addition, the ALICE measurement of $\mathrm{\Lambda_{c}^{+}}$ production in p-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV down to transverse momentum ($p_{\rm T}$) = 0 GeV/$c$ is presented. Finally, the nuclear modification factor is shown for open charm hadrons at $\sqrt{s} = 5.02$ TeV in p-Pb collisions. Finally, measurements of the elliptic flow of heavy-flavor hadron decay leptons in p-Pb systems are presented, which hint towards a possible collective behaviour in high multiplicity p-Pb collisions.


title: Precision Calibration of Radio Interferometers for 21 cm Cosmology with No Redundancy and Little Knowledge of Antenna Beams and the Radio Sky

abstract: We introduce CALAMITY, a precision bandpass calibration method for radio interferometry. CALAMITY can solve for direction independent gains with arbitrary frequency structure to the high precision required for 21 cm cosmology with minimal knowledge of foregrounds or antenna beams and does not require any degree of redundancy (repeated identical measurements of the same baseline). We have achieved this through two key innovations. Firstly, we model the foregrounds on each baseline independently using a flexible and highly efficient set of basis functions that have minimal overlap with 21 cm modes and enforce spectral smoothness in the calibrated foregrounds. Secondly, we use an off-the-shelf GPU accelerated API (TENSORFLOW) to solve for per-baseline foregrounds simultaneously with per-frequency antenna gains in a single optimization loop. GPU acceleration is critical for our technique to be able to solve for the large numbers of foreground and gain parameters simultaneously across all frequencies for an interferometer with $\gtrsim 10$ antennas in a reasonable amount of time. In this paper, we give an overview of our technique and using realistic simulations and demonstrate its performance in solving for and removing pathological gain structures to the level necessary to measure fluctuations in the 21 cm emission field from Hydrogen gas during the Cosmic Dawn and Reionization. If you want to start using CALAMITY now, you can find a tutorial notebook at https://github.com/aewallwi/calamity/blob/main/examples/Calamity_Tutorial.ipynb .