title: Topological reconstruction theorems for varieties

abstract: We study Torelli-type theorems in the Zariski topology for varieties of dimension at least 2, over arbitrary fields. In place of the Hodge structure, we use the linear equivalence relation on Weil divisors. Using this setup, we prove a universal Torelli theorem in the sense of Bogomolov and Tschinkel. The proofs rely heavily on new variants of the classical Fundamental Theorem of Projective Geometry of Veblen and Young. For proper normal varieties over uncountable algebraically closed fields of characteristic 0, we show that the Zariski topological space can be used to recover the linear equivalence relation on divisors. As a consequence, we show that the underlying scheme of any such variety is uniquely determined by its Zariski topological space. We use this to prove a topological version of Gabriel's theorem, stating that a proper normal variety over an uncountable algebraically closed field of characteristic 0 is determined by its category of constructible abelian \'etale sheaves. We also discuss a conjecture in arbitrary characteristic, relating the Zariski topological space to the perfection of a proper normal variety.


title: On the FRB luminosity function -- II. Event rate density

abstract: The luminosity function of Fast Radio Bursts (FRBs), defined as the event rate per unit cosmic co-moving volume per unit luminosity, may help to reveal the possible origins of FRBs and design the optimal searching strategy. With the Bayesian modelling, we measure the FRB luminosity function using 46 known FRBs. Our Bayesian framework self-consistently models the selection effects, including the survey sensitivity, the telescope beam response, and the electron distributions from Milky Way / the host galaxy / local environment of FRBs. Different from the previous companion paper, we pay attention to the FRB event rate density and model the event counts of FRB surveys based on the Poisson statistics. Assuming a Schechter luminosity function form, we infer (at the 95% confidence level) that the characteristic FRB event rate density at the upper cut-off luminosity $L^*=2.9_{-1.7}^{+11.9}\times10^{44}\,\rm erg\, s^{-1}$ is $\phi^*=339_{-313}^{+1074}\,\rm Gpc^{-3}\, yr^{-1}$, the power-law index is $\alpha=-1.79_{-0.35}^{+0.31}$, and the lower cut-off luminosity is $L_0\le9.1\times10^{41}\,\rm erg\, s^{-1}$. The event rate density of FRBs is found to be $3.5_{-2.4}^{+5.7}\times10^4\,\rm Gpc^{-3}\, yr^{-1}$ above $10^{42}\,\rm erg\, s^{-1}$, $5.0_{-2.3}^{+3.2}\times10^3\,\rm Gpc^{-3}\, yr^{-1}$ above $10^{43}\,\rm erg\, s^{-1}$, and $3.7_{-2.0}^{+3.5}\times10^2\,\rm Gpc^{-3}\, yr^{-1}$ above $10^{44}\,\rm erg\, s^{-1}$. As a result, we find that, for searches conducted at 1.4 GHz, the optimal diameter of single-dish radio telescopes to detect FRBs is 30-40 m. The possible astrophysical implications of the measured event rate density are also discussed in the current paper.


title: Evidence for unbounded growth of the number entropy in many-body localized phases

abstract: We investigate the number entropy $S_N$---which characterizes particle-number fluctuations between subsystems---following a quench in one-dimensional interacting many-body systems with potential disorder. We find evidence that in the regime which is expected to show many-body localization (MBL) and where the entanglement entropy grows as $S\sim \ln t$ as function of time $t$, the number entropy grows as $S_N\sim\ln\ln t$, indicating continuing particle transport at a very slow rate. We demonstrate that this growth is consistent with a relation between entanglement and number entropy recently established for non-interacting systems.


title: Covariant model for the Dalitz decay of the $N(1535)$ resonance

abstract: We develop a covariant model for the $\gamma^\ast N \to N(1535)$ transition in the timelike kinematical region, the region where the square momentum transfer $q^2$ is positive. Our starting point is the covariant spectator quark model constrained by data in the spacelike kinematical region ($Q^2 = -q^2 >0$). The model is used to estimate the contributions of valence quarks to the transition form factors, and one obtains a fair description of the Dirac form factor at intermediate and large $Q^2$. For the Pauli form factor there is evidence that beyond the quark-core contributions there are also significant contributions of meson cloud effects. Combining the quark-core model with an effective description of the meson cloud effects, we derive a parametrization of the spacelike data that can be extended covariantly to the timelike region. This extension enabled us to estimate the Dalitz decay widths of the $N(1535)$ resonance, among other observables. Our calculations can help in the interpretation of the present experiments at HADES ($pp$ collisions and others).


title: Revisiting the distributions of Jupiter's irregular moons: II. orbital characteristics

abstract: This paper statistically describes the orbital distribution laws of Jupiter's irregular moons, most of which are members of the Ananke, Carme and Pasiphae groups. By comparing 19 known continuous distributions, it is verified that suitable distribution functions exist to describe the orbital distributions of these natural satellites. For each distribution type, interval estimation is used to estimate the corresponding parameter values. At a given significance level, a one-sample Kolmogorov-Smirnov non-parametric test is applied to verify the specified distribution, and we often select the one with the largest $p$-value. The results show that the semi-major axis, mean inclination and orbital period of the moons in the Ananke group and Carme group obey Stable distributions. In addition, according to Kepler's third law of planetary motion and by comparing the theoretically calculated best-fitting cumulative distribution function (CDF) with the observed CDF, we demonstrate that the theoretical distribution is in good agreement with the empirical distribution. Therefore, these characteristics of Jupiter's irregular moons are indeed very likely to follow some specific distribution laws, and it will be possible to use these laws to help study certain features of poorly investigated moons or even predict undiscovered ones.


title: PANDA: A Gigapixel-level Human-centric Video Dataset

abstract: We present PANDA, the first gigaPixel-level humAN-centric viDeo dAtaset, for large-scale, long-term, and multi-object visual analysis. The videos in PANDA were captured by a gigapixel camera and cover real-world scenes with both wide field-of-view (~1 square kilometer area) and high-resolution details (~gigapixel-level/frame). The scenes may contain 4k head counts with over 100x scale variation. PANDA provides enriched and hierarchical ground-truth annotations, including 15,974.6k bounding boxes, 111.8k fine-grained attribute labels, 12.7k trajectories, 2.2k groups and 2.9k interactions. We benchmark the human detection and tracking tasks. Due to the vast variance of pedestrian pose, scale, occlusion and trajectory, existing approaches are challenged by both accuracy and efficiency. Given the uniqueness of PANDA with both wide FoV and high resolution, a new task of interaction-aware group detection is introduced. We design a 'global-to-local zoom-in' framework, where global trajectories and local interactions are simultaneously encoded, yielding promising results. We believe PANDA will contribute to the community of artificial intelligence and praxeology by understanding human behaviors and interactions in large-scale real-world scenes. PANDA Website: http://www.panda-dataset.com.


title: Torsional Anomalies and Bulk-Dislocation Correspondence in Weyl Systems

abstract: Based on the supersymmetric quantum mechanical approach, we have systematically studied both the $U\left(1\right)$ gauge anomaly and the diffeomorphism anomaly in Weyl systems with torsion, curvature and external electromagnetic fields. These anomalies relate to the chiral current (or current) non-conservation and chiral energy-momentum (or energy-momentum) non-conservation, respectively, which can be applied to the $^{3}\text{He-A}$ phase, the chiral superconductors and the Weyl semimetals with dislocations and disclinations. In sharp difference with other anomalies, there exist torsional anomalies depending on the position of Weyl nodes in the energy-momentum space. These anomalies originate from particles pumped up through the Weyl nodes and they are thus insensitive to the ultra-violet physics, while the Nieh-Yan anomaly is from the particle inflow through the ultra-violet cut-off. The current non-conservation as well as the energy-momentum non-conservation are found, which stem from the zero modes trapped in the dislocations and they can be understood from the Callan-Harvey mechanism. Finally, by comparing our results with the well-established momentum anomaly in the $^{3}\text{He-A}$ phase, the Nieh-Yan term as well as other cut-off dependent terms are shown to be negligible, because the ratio between the Lorentz symmetry breaking scale and the chemical potential is of order $10^{-5}$.


title: The primary importance of the research question: Implications for understanding natural versus controlled direct effects and the 'cross-world independence assumption'

abstract: When developing new interventions to minimize the harmful effects of an exposure, investigators usually target the mechanisms that mediate the causal effect of the exposure on the outcome. Predicting the causal effect of these new interventions is generally done through identifying either (1) the controlled direct effect, or (2) the pure (natural) direct effect. In this opinion piece, we use the interventionist approach to discuss how these two approaches answer different questions, and the additional underlying assumptions of each compared to the other. We use a specific example for the development of a new intervention that might reduce the harmful effects of smoking on chronic obstructive pulmonary disease by removing the inhalation of harmful chemicals.


title: Modeling Multiscale Variable Renewable Energy and Inflow Scenarios in Very Large Regions with Nonparametric Bayesian Networks

abstract: In this paper, we propose a non-parametric Bayesian network method to generate synthetic scenarios of hourly generation for variable renewable energy(VRE) plants. The methodology consists of a non-parametric estimation of the probability distribution of VRE generation, followed by an inverse probability integral transform, in order to obtain normally distributed variables of VRE generation. Then, we build a Bayesian network based on the evaluation of the spatial correlation between variables (VRE generation and hydro inflows, but load forecast, temperature, and other types of random variables could also be used with the proposed framework), to generate future synthetic scenarios while keeping the historical spatial correlation structure. Finally, we present a real-life case study, that uses real data from the Brazilian power system, to show the improvements that the present methodology allows for real-life studies.


title: Searching for Thermal Inversion Agents in the Transmission Spectrum of MASCARA-2b/KELT-20b: Detection of Neutral Iron and Ionised Calcium H$\&$K Lines

abstract: We analyse the transmission spectra of KELT-20b/MASCARA-2b to search for possible thermal inversion agents. The data consist of three transits obtained using HARPSN and one using CARMENES. We removed stellar and telluric lines before cross-correlating the residuals with spectroscopic templates produced using a 1D plane-parallel model assuming an isothermal atmosphere and chemical equilibrium at solar metallicity. Using a likelihood-mapping method, we detect Fe\,{\sc i} at $>$ 13-$\sigma$, Ca\,{\sc ii} H$\&$K at $>$ 6-$\sigma$ and confirm the previous detections of Fe\,{\sc ii}, Ca\,{\sc ii} IRT and Na\,{\sc i} D. The detected signal of Fe\,{\sc i} is shifted by -3.4$\pm$0.4 km s$^{-1}$ from the planetary rest frame, which indicates a strong day-night wind. Our likelihood-mapping technique also reveals that the absorption features of the detected species extend to different altitudes in the planet's atmosphere. Assuming that the line lists are accurate, we do not detect other potential thermal inversion agents (NaH, MgH, AlO, SH, CaO, VO, FeH and TiO) suggesting that non-chemical equilibrium mechanisms (e.g. a cold-trap) might have removed Ti- and V-bearing species from the upper atmosphere. Our results, therefore, shows that KELT-20b/MASCARA-2b cannot possess an inversion layer caused by a TiO/VO-related mechanism. The presence of an inversion layer would therefore likely be caused by metal atoms such as Fe\,{\sc i} and Fe\,{\sc ii}. Finally, we report a double-peak structure in the Fe\,{\sc i} signal in all of our data-sets that could be a signature of atmospheric dynamics. However, further investigation is needed to robustly determine the origin of the signal.


title: Image Restoration for Under-Display Camera

abstract: The new trend of full-screen devices encourages us to position a camera behind a screen. Removing the bezel and centralizing the camera under the screen brings larger display-to-body ratio and enhances eye contact in video chat, but also causes image degradation. In this paper, we focus on a newly-defined Under-Display Camera (UDC), as a novel real-world single image restoration problem. First, we take a 4k Transparent OLED (T-OLED) and a phone Pentile OLED (P-OLED) and analyze their optical systems to understand the degradation. Second, we design a Monitor-Camera Imaging System (MCIS) for easier real pair data acquisition, and a model-based data synthesizing pipeline to generate Point Spread Function (PSF) and UDC data only from display pattern and camera measurements. Finally, we resolve the complicated degradation using deconvolution-based pipeline and learning-based methods. Our model demonstrates a real-time high-quality restoration. The presented methods and results reveal the promising research values and directions of UDC.


title: Unpaired Image-to-Image Translation using Adversarial Consistency Loss

abstract: Unpaired image-to-image translation is a class of vision problems whose goal is to find the mapping between different image domains using unpaired training data. Cycle-consistency loss is a widely used constraint for such problems. However, due to the strict pixel-level constraint, it cannot perform geometric changes, remove large objects, or ignore irrelevant texture. In this paper, we propose a novel adversarial-consistency loss for image-to-image translation. This loss does not require the translated image to be translated back to be a specific source image but can encourage the translated images to retain important features of the source images and overcome the drawbacks of cycle-consistency loss noted above. Our method achieves state-of-the-art results on three challenging tasks: glasses removal, male-to-female translation, and selfie-to-anime translation.


title: Scattering matrix of arbitrary tight-binding Hamiltonians

abstract: A novel efficient method to calculate the scattering matrix (SM) of arbitrary tight-binding Hamiltonians is proposed, including cases with multiterminal structures. In particular, the SM of two kind of fundamental structures are given, which can be used to obtain the SM of bigger systems iteratively. Also, a procedure to obtain the SM of layer-composed periodic leads is described. This method allows renormalization approaches, which permits computations over macroscopic length systems without introducing additional approximations. Finally, the transmission coefficient of a ring-shaped multiterminal system and the transmission function of a square-lattice nanoribbon with a reduced width region are calculated.


title: Quantifying Defects in Thin Films using Machine Vision

abstract: The sensitivity of thin-film materials and devices to defects motivates extensive research into the optimization of film morphology. This research could be accelerated by automated experiments that characterize the response of film morphology to synthesis conditions. Optical imaging can resolve morphological defects in thin films and is readily integrated into automated experiments but the large volumes of images produced by such systems require automated analysis. Existing approaches to automatically analyzing film morphologies in optical images require application-specific customization by software experts and are not robust to changes in image content or imaging conditions. Here we present a versatile convolutional neural network (CNN) for thin-film image analysis which can identify and quantify the extent of a variety of defects and is applicable to multiple materials and imaging conditions. This CNN is readily adapted to new thin-film image analysis tasks and will facilitate the use of imaging in automated thin-film research systems.


title: Fast, Convexified Stochastic Optimal Open-Loop Control For Linear Systems Using Empirical Characteristic Functions

abstract: We consider the problem of stochastic optimal control in the presence of an unknown disturbance. We characterize the disturbance via empirical characteristic functions, and employ a chance constrained approach. By exploiting properties of characteristic functions and underapproximating cumulative distribution functions, we can reformulate a nonconvex problem by a conic, convex under-approximation. This results in extremely fast solutions that are assured to maintain probabilistic constraints. We construct algorithms for optimal open-loop control using piecewise linear approximations of the empirical characteristic function, and demonstrate our approach on two examples.


title: Compensation for undefined behaviors during robot task execution by switching controllers depending on embedded dynamics in RNN

abstract: Robotic applications require both correct task performance and compensation for undefined behaviors. Although deep learning is a promising approach to perform complex tasks, the response to undefined behaviors that are not reflected in the training dataset remains challenging. In a human-robot collaborative task, the robot may adopt an unexpected posture due to collisions and other unexpected events. Therefore, robots should be able to recover from disturbances for completing the execution of the intended task. We propose a compensation method for undefined behaviors by switching between two controllers. Specifically, the proposed method switches between learning-based and model-based controllers depending on the internal representation of a recurrent neural network that learns task dynamics. We applied the proposed method to a pick-and-place task and evaluated the compensation for undefined behaviors. Experimental results from simulations and on a real robot demonstrate the effectiveness and high performance of the proposed method.


title: Circulation Control for Faster Minimum Cost Flow in Unit-Capacity Graphs

abstract: We present an $m^{4/3+o(1)}\log W$-time algorithm for solving the minimum cost flow problem in graphs with unit capacity, where $W$ is the maximum absolute value of any edge weight. For sparse graphs, this improves over the best known running time for this problem and, by well-known reductions, also implies improved running times for the shortest path problem with negative weights, minimum cost bipartite $\boldsymbol{\mathit{b}}$-matching when $\|\boldsymbol{\mathit{b}}\|_1 = O(m)$, and recovers the running time of the currently fastest algorithm for maximum flow in graphs with unit capacities (Liu-Sidford, 2020). Our algorithm relies on developing an interior point method-based framework which acts on the space of circulations in the underlying graph. From the combinatorial point of view, this framework can be viewed as iteratively improving the cost of a suboptimal solution by pushing flow around circulations. These circulations are derived by computing a regularized version of the standard Newton step, which is partially inspired by previous work on the unit-capacity maximum flow problem (Liu-Sidford, 2019), and subsequently refined based on the very recent progress on this problem (Liu-Sidford, 2020). The resulting step problem can then be computed efficiently using the recent work on $\ell_p$-norm minimizing flows (Kyng-Peng-Sachdeva-Wang, 2019). We obtain our faster algorithm by combining this new step primitive with a customized preconditioning method, which aims to ensure that the graph on which these circulations are computed has sufficiently large conductance.


title: Combining neutrino experimental light-curves for pointing to the next Galactic Core-Collapse Supernova

abstract: The multi-messenger observation of the next galactic core-collapse supernova will shed light on the different physical processes involved in these energetic explosions. Good timing and pointing capabilities of neutrino detectors would help in the search for an electromagnetic or gravitational-wave counterparts. An approach for the determination of the arrival time delay of the neutrino signal at different experiments using a direct detected neutrino light-curve matching is discussed. A simplified supernova model and detector simulation are used for its application. The arrival time delay and its uncertainty between two neutrino detectors are estimated with chi-square and cross-correlation methods. The direct comparison of the detected light-curves offers the advantage to be model-independent. Millisecond time resolution on the arrival time delay at two different detectors is needed. Using the computed time delay between different combinations of currently operational and future detectors, a triangulation method is used to infer the supernova localisation in the sky. The combination of IceCube, Hyper-Kamiokande, JUNO and KM3NeT/ARCA provides a 90% confidence area of 140$\pm$20 deg$^2$. These low-latency analysis methods can be implemented in the SNEWS alert system.


title: Video Caption Dataset for Describing Human Actions in Japanese

abstract: In recent years, automatic video caption generation has attracted considerable attention. This paper focuses on the generation of Japanese captions for describing human actions. While most currently available video caption datasets have been constructed for English, there is no equivalent Japanese dataset. To address this, we constructed a large-scale Japanese video caption dataset consisting of 79,822 videos and 399,233 captions. Each caption in our dataset describes a video in the form of "who does what and where." To describe human actions, it is important to identify the details of a person, place, and action. Indeed, when we describe human actions, we usually mention the scene, person, and action. In our experiments, we evaluated two caption generation methods to obtain benchmark results. Further, we investigated whether those generation methods could specify "who does what and where."


title: Multi-SimLex: A Large-Scale Evaluation of Multilingual and Cross-Lingual Lexical Semantic Similarity

abstract: We introduce Multi-SimLex, a large-scale lexical resource and evaluation benchmark covering datasets for 12 typologically diverse languages, including major languages (e.g., Mandarin Chinese, Spanish, Russian) as well as less-resourced ones (e.g., Welsh, Kiswahili). Each language dataset is annotated for the lexical relation of semantic similarity and contains 1,888 semantically aligned concept pairs, providing a representative coverage of word classes (nouns, verbs, adjectives, adverbs), frequency ranks, similarity intervals, lexical fields, and concreteness levels. Additionally, owing to the alignment of concepts across languages, we provide a suite of 66 cross-lingual semantic similarity datasets. Due to its extensive size and language coverage, Multi-SimLex provides entirely novel opportunities for experimental evaluation and analysis. On its monolingual and cross-lingual benchmarks, we evaluate and analyze a wide array of recent state-of-the-art monolingual and cross-lingual representation models, including static and contextualized word embeddings (such as fastText, M-BERT and XLM), externally informed lexical representations, as well as fully unsupervised and (weakly) supervised cross-lingual word embeddings. We also present a step-by-step dataset creation protocol for creating consistent, Multi-Simlex-style resources for additional languages. We make these contributions -- the public release of Multi-SimLex datasets, their creation protocol, strong baseline results, and in-depth analyses which can be be helpful in guiding future developments in multilingual lexical semantics and representation learning -- available via a website which will encourage community effort in further expansion of Multi-Simlex to many more languages. Such a large-scale semantic resource could inspire significant further advances in NLP across languages.


title: Quantum sensing networks for the estimation of linear functions

abstract: The theoretical framework for networked quantum sensing has been developed to a great extent in the past few years, but there are still a number of open questions. Among these, a problem of great significance, both fundamentally and for constructing efficient sensing networks, is that of the role of inter-sensor correlations in the simultaneous estimation of multiple linear functions, where the latter are taken over a collection local parameters and can thus be seen as global properties. In this work we provide a solution to this when each node is a qubit and the state of the network is sensor-symmetric. First we derive a general expression linking the amount of inter-sensor correlations and the geometry of the vectors associated with the functions, such that the asymptotic error is optimal. Using this we show that if the vectors are clustered around two special subspaces, then the optimum is achieved when the correlation strength approaches its extreme values, while there is a monotonic transition between such extremes for any other geometry. Furthermore, we demonstrate that entanglement can be detrimental for estimating non-trivial global properties, and that sometimes it is in fact irrelevant. Finally, we perform a non-asymptotic analysis of these results using a Bayesian approach, finding that the amount of correlations needed to enhance the precision crucially depends on the number of measurement data. Our results will serve as a basis to investigate how to harness correlations in networks of quantum sensors operating both in and out of the asymptotic regime.


title: This PIN Can Be Easily Guessed: Analyzing the Security of Smartphone Unlock PINs

abstract: In this paper, we provide the first comprehensive study of user-chosen 4- and 6-digit PINs (n=1220) collected on smartphones with participants being explicitly primed for device unlocking. We find that against a throttled attacker (with 10, 30, or 100 guesses, matching the smartphone unlock setting), using 6-digit PINs instead of 4-digit PINs provides little to no increase in security, and surprisingly may even decrease security. We also study the effects of blocklists, where a set of "easy to guess" PINs is disallowed during selection. Two such blocklists are in use today by iOS, for 4-digits (274 PINs) as well as 6-digits (2910 PINs). We extracted both blocklists compared them with four other blocklists, including a small 4-digit (27 PINs), a large 4-digit (2740 PINs), and two placebo blocklists for 4- and 6-digit PINs that always excluded the first-choice PIN. We find that relatively small blocklists in use today by iOS offer little or no benefit against a throttled guessing attack. Security gains are only observed when the blocklists are much larger, which in turn comes at the cost of increased user frustration. Our analysis suggests that a blocklist at about 10% of the PIN space may provide the best balance between usability and security.


title: Resonance peak shift in the photo-current of ultrahigh-mobility two-dimensional electron systems

abstract: We report on a theoretical study on the rise of strong peaks at the harmonics of the cyclotron resonance in the irradiated magnetoresistance in ultraclean two-dimensional electron systems. The motivation is the experimental observation of a totally unexpected strong resistance peak showing up at the second harmonic. We extend the radiation-driven electron orbit model (previously developed to study photocurrent oscillations and zero resistance states) to a ultraclean scenario that implies longer scattering time and longer mean free path. Thus, when the mean free path is equivalent, in terms of energy, to twice the cyclotron energy ($2\hbar w_{c}$), the electron behaves as under an effective magnetic field twice the one really applied. Then, at high radiation power and/or low temperature, a resistance spike can be observed {\it at the second harmonic}. For even cleaner samples the energy distance could increase to three or four times the cyclotron energy giving rise to resistance peaks at higher harmonics (third, fourth, etc.), i.e., a resonance peak shift to lower magnetic fields as the quality of the sample increases. Thus, by selecting the sample mobility one automatically would select the radiation resonance response without altering the radiation frequency.


title: Koopman Operator Methods for Global Phase Space Exploration of Equivariant Dynamical Systems

abstract: In this paper, we develop the Koopman operator theory for dynamical systems with symmetry. In particular, we investigate how the Koopman operator and eigenfunctions behave under the action of the symmetry group of the underlying dynamical system. Further, exploring the underlying symmetry, we propose an algorithm to construct a global Koopman operator from local Koopman operators. In particular, we show, by exploiting the symmetry, data from all the invariant sets are not required for constructing the global Koopman operator; that is, local knowledge of the system is enough to infer the global dynamics.


title: Through a Jet Speed Darkly: The Emergence of Robust Euro-Atlantic Regimes in the Absence of Jet Speed Variability

abstract: Euro-Atlantic regimes are typically identified using either the latitude of the eddy-driven jet, or clustering algorithms in the phase space of 500hPa geopotential height (Z500). However, while robust trimodality is visibly apparent in jet latitude indices, Z500 clusters require highly sensitive significance tests to distinguish them from autocorrelated noise. As a result, even small shifts in the time-period considered can notably alter the diagnosed regimes. Fixing the optimal regime number is also hard to justify. We argue that the jet speed, a near-Gaussian distribution projecting strongly onto the Z500 field, is the source of this lack of robustness. Once its influence is removed, the Z500 phase space becomes visibly non-Gaussian, and clustering algorithms easily recover three extremely stable regimes, corresponding to the jet latitude regimes. Further analysis supports the existence of two additional regimes, corresponding to a tilted and split jet. This framework therefore naturally unifies the two regime perspectives.


title: Modelling the He I triplet absorption at 10830 Angstroms in the atmosphere of HD 209458 b

abstract: HD 209458b is an exoplanet with an upper atmosphere undergoing blow-off escape that has mainly been studied using measurements of the Ly-alpha absorption. Recently, high-resolution measurements of absorption in the He I triplet line at 10830 angstroms of several exoplanets (including HD 209458b) have been reported, creating a new opportunity to probe escaping atmospheres. We aim to better understand the atmospheric regions of HD 209458b from where the escape originates. We developed a 1D hydrodynamic model with spherical symmetry for the HD 209458 b thermosphere coupled with a non-local thermodynamic model for the population of the He triplet state. In addition, we performed high-resolution radiative transfer calculations of synthetic spectra for the He triplet lines and compared them with the measured absorption spectrum in order to retrieve information about the atmospheric parameters. We find that the measured spectrum constrains the [H]/[H$^{+}$] transition altitude occurring in the range of 1.2 to 1.9Rp. H is almost fully ionised at altitudes above 2.9Rp. We also find that the X-ray and EUV absorption takes place at effective radii from 1.16 to 1.30Rp, and that the He triplet peak density occurs at altitudes from 1.04 to 1.60Rp. Additionally, the averaged mmw is confined to the 0.61-0.73 g/mole interval, and the thermospheric H/He ratio should be larger than 90/10, and most likely approximately 98/2. We also provide a one-to-one relationship between mass-loss rate and temperature. Based on the energy-limited escape approach and assuming heating efficiencies of 0.1-0.2, we find a mass-loss rate in the range of (0.42-1.00)$\times 10^{11}$ g/s and a corresponding temperature range of 7125 to 8125K. The analysis of the measured He triplet absorption spectrum significantly constrains the thermospheric structure of HD 209458b and advances our knowledge of its escaping atmosphere.


title: Moving Target Monte Carlo

abstract: The Markov Chain Monte Carlo (MCMC) methods are popular when considering sampling from a high-dimensional random variable $\mathbf{x}$ with possibly unnormalised probability density $p$ and observed data $\mathbf{d}$. However, MCMC requires evaluating the posterior distribution $p(\mathbf{x}|\mathbf{d})$ of the proposed candidate $\mathbf{x}$ at each iteration when constructing the acceptance rate. This is costly when such evaluations are intractable. In this paper, we introduce a new non-Markovian sampling algorithm called Moving Target Monte Carlo (MTMC). The acceptance rate at $n$-th iteration is constructed using an iteratively updated approximation of the posterior distribution $a_n(\mathbf{x})$ instead of $p(\mathbf{x}|\mathbf{d})$. The true value of the posterior $p(\mathbf{x}|\mathbf{d})$ is only calculated if the candidate $\mathbf{x}$ is accepted. The approximation $a_n$ utilises these evaluations and converges to $p$ as $n \rightarrow \infty$. A proof of convergence and estimation of convergence rate in different situations are given.


title: Wasserstein Distributionally Robust Look-Ahead Economic Dispatch

abstract: We consider the problem of look-ahead economic dispatch (LAED) with uncertain renewable energy generation. The goal of this problem is to minimize the cost of conventional energy generation subject to uncertain operational constraints. The risk of violating these constraints must be below a given threshold for a family of probability distributions with characteristics similar to observed past data or predictions. We present two data-driven approaches based on two novel mathematical reformulations of this distributionally robust decision problem. The first one is a tractable convex program in which the uncertain constraints are defined via the distributionally robust conditional-value-at-risk. The second one is a scalable robust optimization program that yields an approximate distributionally robust chance-constrained LAED. Numerical experiments on the IEEE 39-bus system with real solar production data and forecasts illustrate the effectiveness of these approaches. We discuss how system operators should tune these techniques in order to seek the desired robustness-performance trade-off and we compare their computational scalability.


title: Gravity field modelling for the Hannover 10m atom interferometer

abstract: Absolute gravimeters (AG) are used in geodesy, geophysics, and physics for a wide spectrum of applications. Stable gravimetric measurements over timescales from several days to decades are required to provide relevant insight into geophysical processes. Users of AGs participate in comparisons with a metrological reference in order to monitor the temporal stability of the instruments and determine the bias to that reference. However, since no measurement standard of higher-order accuracy currently exists, users of AGs participate in key comparisons led by the CIPM. These comparisons provide the reference values of highest accuracy compared to the calibration against a single AG. The construction of stationary, large scale atom interferometers paves the way towards a new measurement standard in absolute gravimetry used as a reference with a potential stability up to 1 nm/s$^2$ at 1 s integration time. At the Leibniz University Hannover, we are currently building such a very long baseline atom interferometer with a 10 m long interaction zone. The knowledge of local gravity and its gradient along and around the baseline is required to establish the instrument's uncertainty budget and enable transfers of gravimetric measurements to nearby devices for comparison and calibration purposes. We therefore established a control network for relative gravimeters and repeatedly measured its connections during the construction of the atom interferometer. We additionally developed a 3D model of the host building to investigate the self-attraction effect and studied the impact of mass changes due to groundwater hydrology on the gravity field around the reference instrument. The gravitational effect from the building 3D model is in excellent agreement with the latest gravimetric measurement campaign which opens the possibility to transfer gravity values with an uncertainty below the 10 nm/s$^2$ level.


title: Generating functions and topological complexity

abstract: We examine the rationality conjecture which states that (a) the formal power series $\sum_{r\ge 1} \tc_{r+1}(X)\cdot x^r$ represents a rational function of $x$ with a single pole of order 2 at $x=1$ and (b) the leading coefficient of the pole equals $\cat(X)$. Here $X$ is a finite CW-complex and for $r\ge 2$ the symbol $\tc_r(X)$ denotes its $r$-th sequential topological complexity. We analyse an example (violating the Ganea conjecture) and conclude that part (b) of the rationality conjecture is false in general. Besides, we establish a cohomological version of the rationality conjecture.


title: Re-inflation of warm and hot Jupiters

abstract: Understanding the anomalous radii of many transiting hot gas giant planets is a fundamental problem of planetary science. Recent detections of re-inflated warm Jupiters orbiting post-main-sequence stars and the re-inflation of hot Jupiters while their host stars evolve on the main-sequence may help constrain models for the anomalous radii of hot Jupiters. In this work, we present evolution models studying the re-inflation of gas giants to determine how varying the depth and intensity of deposited heating affects both main-sequence re-inflation of hot Jupiters and post-main-sequence re-inflation of warm Jupiters. We find that deeper heating is required to re-inflate hot Jupiters than is needed to suppress their cooling, and that the timescale of re-inflation decreases with increasing heating rate and depth. We find a strong degeneracy between heating rate and depth, with either strong shallow heating or weak deep heating providing an explanation for main-sequence re-inflation of hot Jupiters. This degeneracy between heating rate and depth can be broken in the case of post-main-sequence re-inflation of warm Jupiters, as the inflation must be rapid to occur within post-main-sequence evolution timescales. We also show that the dependence of heating rate on incident stellar flux inferred from the sample of hot Jupiters can explain re-inflation of both warm and hot Jupiters. TESS will obtain a large sample of warm Jupiters orbiting post-main-sequence stars, which will help to constrain the mechanism(s) causing the anomalous radii of gas giant planets.


title: Calculating Fluctuations and Self-Correlations Numerically for Causal Charge Diffusion in Relativistic Heavy-Ion Collisions

abstract: We study the propagation and diffusion of electric charge fluctuations in the Bjorken hydrodynamic model with both white and Catteneo noise using purely numerical methods. We show that a global lattice of noise fluctuations is required to fully calculate the two-point correlators of charge. We solve the stochastic differential equations that arise from the charge conservation equation on the lattice. We explicitly identify the self-correlation term in the case of Catteneo noise and provide a physical interpretation. We provide a numerical recipe to remove this contribution from the full two-point correlators. Finally, we calculate the balance functions for charged hadrons. By limiting the speed of signal propagation, we observe the expected narrowing of the balance functions after removing the self-correlations.


title: Implementation of a Walsh-Hadamard gate in a superconducting qutrit

abstract: We have implemented a Walsh-Hadamard gate, which performs a quantum Fourier transform, in a superconducting qutrit. The qutrit is encoded in the lowest three energy levels of a capacitively shunted flux device, operated at the optimal flux-symmetry point. We use an efficient decomposition of the Walsh-Hadamard gate into two unitaries, generated by off-diagonal and diagonal Hamiltonians respectively. The gate implementation utilizes simultaneous driving of all three transitions between the three pairs of energy levels of the qutrit, one of which is implemented with a two-photon process. The gate has a duration of 35 ns and an average fidelity over a representative set of states, including preparation and tomography errors, of 99.2%, characterized with quantum state tomography. Compensation of ac-Stark and Bloch-Siegert shifts is essential for reaching high gate fidelities.


title: Nonparametric constraints on neutron star matter with existing and upcoming gravitational wave and pulsar observations

abstract: Observations of neutron stars, whether in binaries or in isolation, provide information about the internal structure of the most extreme material objects in the Universe. In this work, we combine information from recent observations to place joint constraints on the properties of neutron star matter. We use (i) lower limits on the maximum mass of neutron stars obtained through radio observations of heavy pulsars, (ii) constraints on tidal properties inferred through the gravitational waves neutron star binaries emit as they coalesce, and (iii) information about neutron stars' masses and radii obtained through X-ray emission from surface hot spots. In order to combine information from such distinct messengers while avoiding the kind of modeling systematics intrinsic to parametric inference schemes, we employ a nonparametric representation of the neutron-star equation of state based on Gaussian processes conditioned on nuclear theory models. We find that existing astronomical observations imply $R_{1.4}=12.32^{+1.09}_{-1.47}\,\mathrm{km}$ for the radius of a $1.4\,M_{\odot}$ neutron star and $p(2\rho_\mathrm{nuc})=3.8^{+2.7}_{-2.9}\times10^{34}\,\mathrm{dyn}/\mathrm{cm}^2$ for the pressure at twice nuclear saturation density at the 90% credible level. The upper bounds are driven by the gravitational wave observations, while X-ray and heavy pulsar observations drive the lower bounds. Additionally, we compute expected constraints from potential future astronomical observations and find that they can jointly determine $R_{1.4}$ to ${\cal{O}}(1)\,\mathrm{km}$ and $p(2\rho_\mathrm{nuc})$ to 80% relative uncertainty in the next five years.


title: Pruned Neural Networks are Surprisingly Modular

abstract: The learned weights of a neural network are often considered devoid of scrutable internal structure. To discern structure in these weights, we introduce a measurable notion of modularity for multi-layer perceptrons (MLPs), and investigate the modular structure of MLPs trained on datasets of small images. Our notion of modularity comes from the graph clustering literature: a "module" is a set of neurons with strong internal connectivity but weak external connectivity. We find that training and weight pruning produces MLPs that are more modular than randomly initialized ones, and often significantly more modular than random MLPs with the same (sparse) distribution of weights. Interestingly, they are much more modular when trained with dropout. We also present exploratory analyses of the importance of different modules for performance and how modules depend on each other. Understanding the modular structure of neural networks, when such structure exists, will hopefully render their inner workings more interpretable to engineers. Note that this paper has been superceded by "Clusterability in Neural Networks", arxiv:2103.03386 and "Quantifying Local Specialization in Deep Neural Networks", arxiv:2110.08058!


title: Optimization-Based Hierarchical Motion Planning for Autonomous Racing

abstract: In this paper we propose a hierarchical controller for autonomous racing where the same vehicle model is used in a two level optimization framework for motion planning. The high-level controller computes a trajectory that minimizes the lap time, and the low-level nonlinear model predictive path following controller tracks the computed trajectory online. Following a computed optimal trajectory avoids online planning and enables fast computational times. The efficiency is further enhanced by the coupling of the two levels through a terminal constraint, computed in the high-level controller. Including this constraint in the real-time optimization level ensures that the prediction horizon can be shortened, while safety is guaranteed. This proves crucial for the experimental validation of the approach on a full size driverless race car. The vehicle in question won two international student racing competitions using the proposed framework; moreover, our hierarchical controller achieved an improvement of 20% in the lap time compared to the state of the art result achieved using a very similar car and track.


title: Maximum Likelihood Speed Estimation of Moving Objects in Video Signals

abstract: Video processing solutions for motion analysis are key tasks in many computer vision applications, ranging from human activity recognition to object detection. In particular, speed estimation algorithms may be relevant in contexts such as street monitoring and environment surveillance. In most realistic scenarios, the projection of a framed object of interest onto the image plane is likely to be affected by dynamic changes mainly related to perspectival transformations or periodic behaviours. Therefore, advanced speed estimation techniques need to rely on robust algorithms for object detection that are able to deal with potential geometrical modifications. The proposed method is composed of a sequence of pre-processing operations, that aim to reduce or neglect perspetival effects affecting the objects of interest, followed by the estimation phase based on the Maximum Likelihood (ML) principle, where the speed of the foreground objects is estimated. The ML estimation method represents, indeed, a consolidated statistical tool that may be exploited to obtain reliable results. The performance of the proposed algorithm is evaluated on a set of real video recordings and compared with a block-matching motion estimation algorithm. The obtained results indicate that the proposed method shows good and robust performance.


title: Cryptanalytic Extraction of Neural Network Models

abstract: We argue that the machine learning problem of model extraction is actually a cryptanalytic problem in disguise, and should be studied as such. Given oracle access to a neural network, we introduce a differential attack that can efficiently steal the parameters of the remote model up to floating point precision. Our attack relies on the fact that ReLU neural networks are piecewise linear functions, and thus queries at the critical points reveal information about the model parameters. We evaluate our attack on multiple neural network models and extract models that are 2^20 times more precise and require 100x fewer queries than prior work. For example, we extract a 100,000 parameter neural network trained on the MNIST digit recognition task with 2^21.5 queries in under an hour, such that the extracted model agrees with the oracle on all inputs up to a worst-case error of 2^-25, or a model with 4,000 parameters in 2^18.5 queries with worst-case error of 2^-40.4. Code is available at https://github.com/google-research/cryptanalytic-model-extraction.


title: Magnetic flux noise in superconducting qubits and the gap states continuum

abstract: In the present study we investigate the selected local aspects of the metal-induced gap states (MIGSs) at the disordered metal-insulator interface, that were previously proposed to produce magnetic moments responsible for the magnetic flux noise in some of the superconducting qubit modalities. Our analysis attempts to supplement the available studies and provide new theoretical contribution toward their validation. In particular, we explicitly discuss the behavior of the MIGSs in the momentum space as a function of the local onsite energy deviation, that mimics random potential disorder at the interface. It is found, that when the difference between the characteristic electronic potentials in the insulator increases, the corresponding MIGSs become more localized. This effect is associated with the increasing degree of the potential disorder that was earlier observed to produce highly localized MIGSs in the superconducting qubits. At the same time, the presented findings show that the disorder-induced localization of the MIGSs can be related directly to the decay characteristics of these states as well as to the bulk electronic properties of the insulator. As a result, our study reinforces plausibility of the previous corresponding investigations on the origin of the flux noise, but also allows to draw future directions toward their better verification.


title: Hadronic vacuum polarization: $(g-2)_\mu$ versus global electroweak fits

abstract: Hadronic vacuum polarization (HVP) is not only a critical part of the Standard Model (SM) prediction for the anomalous magnetic moment of the muon $(g-2)_\mu$, but also a crucial ingredient for global fits to electroweak (EW) precision observables due to its contribution to the running of the fine-structure constant encoded in $\Delta\alpha^{(5)}_\text{had}$. We find that with modern EW precision data, including the measurement of the Higgs mass, the global fit alone provides a competitive, independent determination of $\Delta \alpha^{(5)}_\text{had}\big|_\text{EW}=270.2(3.0)\times 10^{-4}$. This value actually lies below the range derived from $e^+e^-\to\text{hadrons}$ cross-section data, and thus goes into the opposite direction as would be required if a change in HVP were to bring the SM prediction for $(g-2)_\mu$ into agreement with the Brookhaven measurement. Depending on the energy where the bulk of the changes in the cross section occurs, reconciling experiment and SM prediction for $(g-2)_\mu$ by adjusting HVP would thus not necessarily weaken the case for physics beyond the SM (BSM), but to some extent shift it from $(g-2)_\mu$ to the EW fit. We briefly explore some options of BSM scenarios that could conceivably explain the ensuing tension.


title: ReZero is All You Need: Fast Convergence at Large Depth

abstract: Deep networks often suffer from vanishing or exploding gradients due to inefficient signal propagation, leading to long training times or convergence difficulties. Various architecture designs, sophisticated residual-style networks, and initialization schemes have been shown to improve deep signal propagation. Recently, Pennington et al. used free probability theory to show that dynamical isometry plays an integral role in efficient deep learning. We show that the simplest architecture change of gating each residual connection using a single zero-initialized parameter satisfies initial dynamical isometry and outperforms more complex approaches. Although much simpler than its predecessors, this gate enables training thousands of fully connected layers with fast convergence and better test performance for ResNets trained on CIFAR-10. We apply this technique to language modeling and find that we can easily train 120-layer Transformers. When applied to 12 layer Transformers, it converges 56% faster on enwiki8.


title: Learning Diverse Fashion Collocation by Neural Graph Filtering

abstract: Fashion recommendation systems are highly desired by customers to find visually-collocated fashion items, such as clothes, shoes, bags, etc. While existing methods demonstrate promising results, they remain lacking in flexibility and diversity, e.g. assuming a fixed number of items or favoring safe but boring recommendations. In this paper, we propose a novel fashion collocation framework, Neural Graph Filtering, that models a flexible set of fashion items via a graph neural network. Specifically, we consider the visual embeddings of each garment as a node in the graph, and describe the inter-garment relationship as the edge between nodes. By applying symmetric operations on the edge vectors, this framework allows varying numbers of inputs/outputs and is invariant to their ordering. We further include a style classifier augmented with focal loss to enable the collocation of significantly diverse styles, which are inherently imbalanced in the training set. To facilitate a comprehensive study on diverse fashion collocation, we reorganize Amazon Fashion dataset with carefully designed evaluation protocols. We evaluate the proposed approach on three popular benchmarks, the Polyvore dataset, the Polyvore-D dataset, and our reorganized Amazon Fashion dataset. Extensive experimental results show that our approach significantly outperforms the state-of-the-art methods with over 10% improvements on the standard AUC metric on the established tasks. More importantly, 82.5% of the users prefer our diverse-style recommendations over other alternatives in a real-world perception study.


title: On the Localization of Unmanned Aerial Vehicles with Cellular Networks

abstract: Localization plays a key role for safe operation of UAVs enabling beyond visual line of sight applications. Compared to GPS based localization, cellular networks can reduce the positioning error and cost since cellular connectivity is becoming a prominent solution as a communication system for UAVs. As a first step towards localization, UAV needs to receive sufficient number of localization signals each having a signal to interference plus noise ratio (SINR) greater than a threshold, which is defined as localizability. On the other hand, three-dimensional mobility of UAVs, altitude dependent channel characteristics between base stations (BSs) and UAVs, its line of sight and non-line of sight conditions, and resulting interference from the neighboring BSs pose challenges to receive usable signals from the required number of BSs. In this paper, we utilize a tractable approach to calculate localizability, probability of successfully receiving usable signals from a certain number of BSs, in a two-tier hexagonal cellular network. Localizability has an impact on overall localization performance regardless of the localization technique to be used. In our simulation study, we investigate the relation between the probability of localizability with respect to number of participating BSs, post-processing SINR requirement, air-to-ground channel characteristics, and network coordination, which are shown to be the most important factors for the localizability performance of UAVs.


title: Effects of fixatives on myelin molecular order probed with RP-CARS microscopy

abstract: When live imaging is not feasible, sample fixation allows preserving the ultrastructure of biological samples for subsequent microscopy analysis. This process could be performed with various methods, each one affecting differently the biological structure of the sample. While these alterations were well-characterized using traditional microscopy, little information is available about the effects of the fixatives on the spatial molecular orientation of the biological tissue. We tackled this issue by employing Rotating-Polarization Coherent Anti-Stokes Raman Scattering (RP-CARS) microscopy to study the effects of different fixatives on the myelin sub-micrometric molecular order and micrometric morphology. RP-CARS is a novel technique derived from CARS microscopy that allows probing spatial orientation of molecular bonds while maintaining the intrinsic chemical selectivity of CARS microscopy. By characterizing the effects of the fixation procedures, the present work represents a useful guide for the choice of the best fixation technique(s), in particular for polarisation-resolved CARS microscopy. Finally, we show that the combination of paraformaldehyde and glutaraldehyde can be effectively employed as a fixative for RP-CARS microscopy, as long as the effects on the molecular spatial distribution, here characterized, are taken into account.


title: Data-mining for Fault-Zone Detection of Distance Relay in FACTS-Based Transmission

abstract: In this study, the problem of fault zone detection of distance relaying in FACTS-based transmission lines is analyzed. Existence of FACTS devices on the transmission line, when they are included in the fault zone, from the distance relay point of view, causes different problems in determining the exact location of the fault by changing the impedance seen by the relay. The extent of these changes depends on the parameters that have been set in FACTS devices. To solve the problem associated with these compensators, two instruments for separation and analysis of three-line currents, from the relay point of view at fault instance, have been utilized. The wavelet transform was used to separate the high-frequency components of the three-line currents, and the support vector machine (using methods for multi-class usage) was used for classification of fault location into three protection regions of distance relay. Besides, to investigate the effects of TCSC location on fault zone detection of distance relay, two places, one in fifty percent of line length and the other in seventy-five percent of line length, have been considered as two scenarios for confirmation of the proposed method. Simulations indicate that this method is effective in the protection of FACTS-based transmission lines.


title: RealityCheck: Bringing Modularity, Hierarchy, and Abstraction to Automated Microarchitectural Memory Consistency Verification

abstract: Modern SoCs are heterogeneous parallel systems comprised of components developed by distinct teams and possibly even different vendors. The memory consistency model (MCM) of processors in such SoCs specifies the ordering rules which constrain the values that can be read by load instructions in parallel programs running on such systems. The implementation of required MCM orderings can span components which may be designed and implemented by many different teams. Ideally, each team would be able to specify the orderings enforced by their components independently and then connect them together when conducting MCM verification. However, no prior automated approach for formal hardware MCM verification provided this. To bring automated hardware MCM verification in line with the realities of the design process, we present RealityCheck, a methodology and tool for automated formal MCM verification of modular microarchitectural ordering specifications. RealityCheck allows users to specify their designs as a hierarchy of distinct modules connected to each other rather than a single flat specification. It can then automatically verify litmus test programs against these modular specifications. RealityCheck also provides support for abstraction, which enables scalable verification by breaking up the verification of the entire design into smaller verification problems. We present results for verifying litmus tests on 7 different designs using RealityCheck. These include in-order and out-of-order pipelines, a non-blocking cache, and a heterogeneous processor. Our case studies cover the TSO and RISC-V (RVWMO) weak memory models. RealityCheck is capable of verifying 98 RVWMO litmus tests in under 4 minutes each, and its capability for abstraction enables up to a 32.1% reduction in litmus test verification time for RVWMO.


title: Stability of Gravastars with Exterior Regular Black Holes

abstract: This paper examines the stability of thin-shell gravastars in the context of regular spacetimes (Bardeen and Bardeen-de Sitter black holes). We apply cut and paste approach to construct gravastars through the matching of interior non-singular de Sitter geometry with exterior regular black hole. This model contains three regions, i.e., interior, thin-shell and exterior. The interior and exterior regions are connected at thin-shell. We investigate physical viability of the developed model by the energy conditions and explore its stability by using radial perturbation about the equilibrium shell radius. It is found that thin-shell gravastars show large stable regions for the Bardeen-de Sitter black hole as compared to the Bardeen black hole. It is concluded that stable regions exist near the formation of expected event horizon.


title: HEMlets PoSh: Learning Part-Centric Heatmap Triplets for 3D Human Pose and Shape Estimation

abstract: Estimating 3D human pose from a single image is a challenging task. This work attempts to address the uncertainty of lifting the detected 2D joints to the 3D space by introducing an intermediate state-Part-Centric Heatmap Triplets (HEMlets), which shortens the gap between the 2D observation and the 3D interpretation. The HEMlets utilize three joint-heatmaps to represent the relative depth information of the end-joints for each skeletal body part. In our approach, a Convolutional Network (ConvNet) is first trained to predict HEMlets from the input image, followed by a volumetric joint-heatmap regression. We leverage on the integral operation to extract the joint locations from the volumetric heatmaps, guaranteeing end-to-end learning. Despite the simplicity of the network design, the quantitative comparisons show a significant performance improvement over the best-of-grade methods (e.g. $20\%$ on Human3.6M). The proposed method naturally supports training with "in-the-wild" images, where only weakly-annotated relative depth information of skeletal joints is available. This further improves the generalization ability of our model, as validated by qualitative comparisons on outdoor images. Leveraging the strength of the HEMlets pose estimation, we further design and append a shallow yet effective network module to regress the SMPL parameters of the body pose and shape. We term the entire HEMlets-based human pose and shape recovery pipeline HEMlets PoSh. Extensive quantitative and qualitative experiments on the existing human body recovery benchmarks justify the state-of-the-art results obtained with our HEMlets PoSh approach.


title: Gravastars in $f(R,T,R_{\mu\nu}T^{\mu\nu})$ Gravity

abstract: This work is devoted to study the analytical and regular solutions of a particular self-gravitating object (i.e., gravastar) in a particular theory of gravity. We derive the corresponding field equations in the presence of effective energy momentum tensor associated with the perfect fluid configuration of a spherical system. We then describe the mathematical formulations of the three respective regions i.e., inner, shell and exterior of a gravastar separately. Additionally, the significance and physical characteristics along with the graphical representation of gravastars are discussed in detail. It is seen that under some specific constraints, $f(R,T,R_{\mu\nu}T^{\mu\nu})$ gravity is likely to host gravastars.


title: Unbiased Estimation of the Gradient of the Log-Likelihood in Inverse Problems

abstract: We consider the problem of estimating a parameter associated to a Bayesian inverse problem. Treating the unknown initial condition as a nuisance parameter, typically one must resort to a numerical approximation of gradient of the log-likelihood and also adopt a discretization of the problem in space and/or time. We develop a new methodology to unbiasedly estimate the gradient of the log-likelihood with respect to the unknown parameter, i.e. the expectation of the estimate has no discretization bias. Such a property is not only useful for estimation in terms of the original stochastic model of interest, but can be used in stochastic gradient algorithms which benefit from unbiased estimates. Under appropriate assumptions, we prove that our estimator is not only unbiased but of finite variance. In addition, when implemented on a single processor, we show that the cost to achieve a given level of error is comparable to multilevel Monte Carlo methods, both practically and theoretically. However, the new algorithm provides the possibility for parallel computation on arbitrarily many processors without any loss of efficiency, asymptotically. In practice, this means any precision can be achieved in a fixed, finite constant time, provided that enough processors are available.


title: Distribution of Water Vapor in Molecular Clouds. II

abstract: The depth-dependent abundance of both gas-phase and solid-state water within dense, quiescent, molecular clouds is important to both the cloud chemistry and gas cooling. Where water is in the gas phase, it's free to participate in the network of ion-neutral reactions that lead to a host of oxygen-bearing molecules, and its many ortho and para energy levels make it an effective coolant for gas temperatures greater than 20K. Where water is abundant as ice on grain surfaces, and unavailable to cool the gas, significant amounts of oxygen are removed from the gas phase, suppressing the gas-phase chemical reactions that lead to a number of oxygen-bearing species, including O2. Models of FUV-illuminated clouds predict that the gas-phase water abundance peaks in the range Av ~3 and 8mag of the cloud surface, depending on the gas density and FUV field strength. Deeper within such clouds, water is predicted to exist mainly as ice on grain surfaces. More broadly, these models are used to analyze a variety of other regions, including outflow cavities associated with young stellar objects and the surface layers of protoplanetary disks. In this paper, we report the results of observational tests of FUV-illuminated cloud models toward the Orion Molecular Ridge and Cepheus B using data obtained from the Herschel Space Observatory and the Five College Radio Astronomy Observatory. Toward Orion, 2220 spatial positions were observed along the face-on Orion Ridge in the H2O 110-101 557GHz and NH3 J,K=1,0-0,0 572GHz lines. Toward Cepheus B, two strip scans were made in the same lines across the edge-on ionization front. These new observations demonstrate that gas-phase water exists primarily within a few magnitudes of dense cloud surfaces, strengthening the conclusions of an earlier study based on a much smaller data set, and indirectly supports the prediction that water ice is quite abundant in dense clouds.


title: The pyrochlore S=1/2 Heisenberg antiferromagnet at finite temperature

abstract: Frustrated three dimensional quantum magnets are notoriously impervious to theoretical analysis. Here we use a combination of three computational methods to investigate the three dimensional pyrochlore $S=1/2$ quantum antiferromagnet, an archetypical frustrated magnet, at finite temperature, $T$: canonical typicality for a finite cluster of $2\times 2 \times 2$ unit cells (i.e. $32$ sites), a finite-$T$ matrix product state method on a larger cluster with $48$ sites, and the numerical linked cluster expansion (NLCE) using clusters up to $25$ lattice sites, which include non-trivial hexagonal and octagonal loops. We focus on thermodynamic properties (energy, specific heat capacity, entropy, susceptibility, magnetisation) next to the static structure factor. We find a pronounced maximum in the specific heat at $T = 0.57 J$, which is stable across finite size clusters and converged in the series expansion. This is well-separated from a residual amount of spectral weight of $0.47 k_B \ln 2$ per spin which has not been released even at $T\approx0.25 J$, the limit of convergence of our results. This is a large value compared to a number of highly frustrated models and materials, such as spin ice or the kagome $S=1/2$ Heisenberg antiferromagnet. We also find a non-monotonic dependence on $T$ of the magnetisation at low magnetic fields, reflecting the dominantly non-magnetic character of the low-energy spectral weight. A detailed comparison of our results to measurements for the $S=1$ material NaCaNi$_2$F$_7$ yields rough agreement of the functional form of the specific heat maximum, which in turn differs from the sharper maximum of the heat capacity of the spin ice material Dy$_2$Ti$_2$O$_7$, all of which are yet qualitatively distinct from conventional, unfrustrated magnets.


title: Avoiding gauge ambiguities in cavity QED

abstract: Systems of interacting charges and fields are ubiquitous in physics. Recently, it has been shown that Hamiltonians derived using different gauges can yield different physical results when matter degrees of freedom are truncated to a few low-lying energy eigenstates. This effect is particularly prominent in the ultra-strong coupling regime. Such ambiguities arise because transformations reshuffle the partition between light and matter degrees of freedom and so level truncation is a gauge dependent approximation. To avoid this gauge ambiguity, we redefine the electromagnetic fields in terms of potentials for which the resulting canonical momenta and Hamiltonian are explicitly unchanged by the gauge choice of this theory. Instead the light/matter partition is assigned by the intuitive choice of separating an electric field between displacement and polarisation contributions. This approach is an attractive choice in typical cavity quantum electrodynamics situations.


title: Heavy Thermal Relics from Zombie Collisions

abstract: We propose a new thermal freezeout mechanism which results in dark matter masses exceeding the unitarity bound by many orders of magnitude, without violating perturbative unitarity or modifying the standard cosmology. The process determining the relic abundance is $\chi \zeta^\dagger \to \zeta \zeta$, where $\chi$ is the dark matter candidate. For $ m_\zeta < m_\chi < 3 m_\zeta$, $\chi$ is cosmologically long-lived and scatters against the exponentially more abundant $\zeta$. Therefore, such a process allows for exponentially heavier dark matter for the same interaction strength as a particle undergoing ordinary $2 \to 2$ freezeout; or equivalently, exponentially weaker interactions for the same mass. We demonstrate this mechanism in a leptophilic dark matter model, which allows for dark matter masses up to $10^{9}$ GeV.


title: Disorder-free localization in an interacting two-dimensional lattice gauge theory

abstract: Disorder-free localization has been recently introduced as a mechanism for ergodicity breaking in low-dimensional homogeneous lattice gauge theories caused by local constraints imposed by gauge invariance. We show that also genuinely interacting systems in two spatial dimensions can become nonergodic as a consequence of this mechanism. Specifically, we prove nonergodic behavior in the quantum link model by obtaining a rigorous bound on the localization-delocalization transition through a classical correlated percolation problem implying a fragmentation of Hilbert space on the nonergodic side of the transition. We study the quantum dynamics in this system by means of an efficient and perturbatively controlled representation of the wavefunction in terms of a variational network of classical spins akin to artificial neural networks. We identify a distinguishing dynamical signature by studying the propagation of line defects, yielding different light cone structures in the localized and ergodic phases, respectively. The methods we introduce in this work can be applied to any lattice gauge theory with finite-dimensional local Hilbert spaces irrespective of spatial dimensionality.


title: A landscape of orientifold vacua

abstract: We present a vast landscape of O3/O7 orientifolds that descends from the famous set of complete intersection Calabi-Yau threefolds (CICY). We give distributions of topological data relevant for phenomenology such as the orientifold-odd Hodge numbers, the D3-tadpole, and multiplicities of O3 and O7-planes. Somewhat surprisingly, almost all of these orientifolds have conifold singularities whose deformation branches are projected out by the orientifolding. However, they can be resolved, so most of the orientifolds actually descend from a much larger and possibly new set of CY threefolds that can be reached from the CICYs via conifold transitions. We observe an interesting class of $\mathcal{N}=1$ geometric transitions involving colliding O-planes. Finally, as an application, we use our dataset to produce examples of orientifolds that satisfy the topological requirements for the existence of ultra-light throat axions (\textit{thraxions}) as proposed in \cite{Hebecker:2018yxs}. The database can be accessed at https://www.desy.de/~westphal/orientifold_webpage/cicy_orientifolds.html


title: The QCD Axion at Finite Density

abstract: We show how the properties of the QCD axion change in systems at finite baryonic density, such as neutron stars. At nuclear saturation densities, where corrections can be reliably computed, we find a mild reduction of the axion mass and up to an order of magnitude enhancement in the model-independent axion coupling to neutrons. At moderately higher densities, if realized, meson (kaon) condensation can trigger axion condensation. We also study the axion potential at asymptotically large densities, where the color-superconducting phase of QCD potentially leads to axion condensation, and the mass of the axion is generically several orders of magnitude smaller than in vacuum due to the suppressed instantons. Several phenomenological consequences of the axion being sourced by neutron stars are discussed, such as its contribution to their total mass, the presence of an axionic brane, or axion-photon conversion in the magnetosphere.


title: Cometary Activity Discovered on a Distant Centaur: A Non-Aqueous Sublimation Mechanism

abstract: Centaurs are minor planets thought to have originated in the outer Solar System region known as the Kuiper Belt. Active Centaurs enigmatically display comet-like features (e.g., tails, comae) even though they orbit in the gas giant region where it is too cold for water to readily sublimate. Only 18 active Centaurs have been identified since 1927 and, consequently, the underlying activity mechanism(s) have remained largely unknown up to this point. Here we report the discovery of activity emanating from Centaur 2014 OG392, based on archival images we uncovered plus our own new observational evidence acquired with the Dark Energy Camera (Cerro Tololo Inter-American Observatory Blanco 4 m telescope), the Inamori-Magellan Areal Camera & Spectrograph (Las Campanas Observatory 6.5 m Walter Baade Telescope) and the Large Monolithic Imager (Lowell Observatory 4.3 m Discovery Channel Telescope). We detect a coma as far as 400,000 km from 2014 OG392, and our novel analysis of sublimation processes and dynamical lifetime suggest carbon dioxide and/or ammonia are the most likely candidates for causing activity on this and other active Centaurs. We find 2014 OG392 is optically red, but CO2 and NH3 are spectrally neutral in this wavelength regime so the reddening agent is as yet unidentified.


title: Constraining the reionization history using deep learning from 21cm tomography with the Square Kilometre Array

abstract: Upcoming 21cm surveys with the SKA1-LOW telescope will enable imaging of the neutral hydrogen distribution on cosmological scales in the early Universe. These surveys are expected to generate huge imaging datasets that will encode more information than the power spectrum. This provides an alternative unique way to constrain the reionization history, which might break the degeneracy in the power spectral analysis. Using Convolutional Neural Networks (CNN), we create a fast estimator of the neutral fraction from the 21cm maps that are produced by our large semi-numerical simulation. Our estimator is able to efficiently recover the neutral fraction ($x_{\rm HI}$) at several redshifts with a high accuracy of 99\% as quantified by the coefficient of determination $R^{2}$. Adding the instrumental effects from the SKA design slightly increases the loss function, but nevertheless we are still able to recover the neutral fraction with a similar high accuracy of 98\%, which is only 1 per cent less. While a weak dependence on redshift is observed, the accuracy increases rapidly with decreasing neutral fraction. This is due to the fact that the instrumental noise increases towards high redshift where the Universe is highly neutral. Our results show the promise of directly using 21cm-tomography to constrain the reionization history in a model independent way, complementing similar efforts, such as those of the optical depth measurements from the Cosmic Microwave Background (CMB) observations by {\it Planck}.


title: Multidimensional super- and subradiance in waveguide quantum electrodynamics

abstract: We study the collective decay rates of multi-dimensional quantum networks in which one-dimensional waveguides form an intersecting hyper-rectangular lattice, with qubits located at the lattice points. We introduce and motivate the \emph{dimensional reduction of poles} (DRoP) conjecture, which identifies all collective decay rates of such networks via a connection to waveguides with a one-dimensional topology (e.g. a linear chain of qubits). Using DRoP, we consider many-body effects such as superradiance, subradiance, and bound-states in continuum in multi-dimensional quantum networks. We find that, unlike one-dimensional linear chains, multi-dimensional quantum networks have superradiance in distinct levels, which we call multi-dimensional superradiance. Furthermore, we generalize the $N^{-3}$ scaling of subradiance in a linear chain to $d$-dimensional networks.


title: Highly connected 7-manifolds, the linking form and non-negative curvature

abstract: In a recent article, the authors constructed a six-parameter family of highly connected 7-manifolds which admit an SO(3)-invariant metric of non-negative sectional curvature. Each member of this family is the total space of a Seifert fibration with generic fibre $S^3$ and, in particular, has the cohomology ring of an $S^3$-bundle over $S^4$. In the present article, the linking form of these manifolds is computed and used to demonstrate that the family contains infinitely many manifolds which are not even homotopy equivalent to an $S^3$-bundle over $S^4$, the first time that any such spaces have been shown to admit non-negative sectional curvature.


title: Gravitational Interactions and Neutrino Masses

abstract: We describe a scenario where the smallness of neutrino masses is related to a global symmetry that is only violated by quantum gravitational effects. The coupling of neutrinos to gauge singlet right-handed fermions is attributed to symmetry preserving gravitational operators suppressed by the Planck mass, in this framework. The proposed scenario leads to axion particles that decay into neutrinos, which could be probed through cosmological measurements and may help explain the Hubble parameter tension. Depending on the details of the implementation, the scenario could provide axion dark matter candidates.


title: Topological phase transitions in four dimensions

abstract: We show that four-dimensional systems may exhibit a topological phase transition analogous to the well-known Berezinskii-Kosterlitz-Thouless vortex unbinding transition in two-dimensional systems. The realisation of an engineered quantum system, where the predicted phase transition shall occur, is also presented. We study a suitable generalization of the sine-Gordon model in four dimensions and the renormalization group flow equation of its couplings, showing that the critical value of the frequency is the square of the corresponding value in $2D$. The value of the anomalous dimension at the critical point is determined ($\eta=1/32$) and a conjecture for the universal jump of the superfluid stiffness ($4/\pi^2$) presented.


title: Taking EGSnrc to new lows: Development of egs++ lattice geometry and testing with microscopic

abstract: Purpose: This work introduces a new lattice geometry library, egs_lattice, into the EGSnrc Monte Carlo code, which can be used for both modeling very large (previously unfeasible) quantities of geometries and establishing recursive boundary conditions. The reliability of egs_lattice, as well as EGSnrc in general, is cross-validated and tested. Methods: New Bravais, cubic, and hexagonal lattice geometries are defined in egs_lattice and their transport algorithms are described. Simulations of cells and Gold NanoParticle (GNP) containing cavities are implemented to compare to published Geant4-DNA and PENELOPE results. Recursive boundary conditions, implemented through a cubic lattice, are used to perform electron Fano cavity tests. Results: Lattices are successfully implemented in EGSnrc. EGSnrc calculated doses to cell cytoplasm and nucleus when irradiated by an internal electron source with a median difference of 0.6% compared to Geant4-DNA. EGSnrc calculated the ratio of dose to a microscopic cavity containing GNPs over dose to a cavity containing a homogeneous mixture of gold, and results generally agree (within 1%) with PENELOPE. The Fano test is passed (sub-0.1%) for all energies/ cells considered. Additionally, the recursive boundary conditions used for the Fano test provided a factor of over a million increase in efficiency in some cases. Conclusions: The egs_lattice geometry library, currently available as a pull request on the EGSnrc GitHub develop branch, is now freely accessible as open-source code. Lattice geometry implementations cross-validated with independent simulations in other MC codes and verified with the electron Fano cavity test demonstrate not only the reliability of egs_lattice, but, by extension, EGSnrc's ability to simulate transport in nanometer geometries and score in microscopic cavities.


title: The Smallest Eigenvalue Distribution of the Jacobi Unitary Ensembles

abstract: In the hard edge scaling limit of the Jacobi unitary ensemble generated by the weight $x^{\alpha}(1-x)^{\beta},~x\in[0,1],~\alpha,\beta>0$, the probability that all eigenvalues of Hermitian matrices from this ensemble lie in the interval $[t,1]$ is given by the Fredholm determinant of the Bessel kernel. We derive the constant in the asymptotics of this Bessel-kernel determinant. A specialization of the results gives the constant in the asymptotics of the probability that the interval $(-a,a),a>0,$ is free of eigenvalues in the Jacobi unitary ensemble with the symmetric weight $(1-x^2)^{\beta}, x\in[-1,1]$.


title: Flip-sort and combinatorial aspects of pop-stack sorting

abstract: Flip-sort is a natural sorting procedure which raises fascinating combinatorial questions. It finds its roots in the seminal work of Knuth on stack-based sorting algorithms and leads to many links with permutation patterns. We present several structural, enumerative, and algorithmic results on permutations that need few (resp. many) iterations of this procedure to be sorted. In particular, we give the shape of the permutations after one iteration, and characterize several families of permutations related to the best and worst cases of flip-sort. En passant, we also give some links between pop-stack sorting, automata, and lattice paths, and introduce several tactics of bijective proofs which have their own interest.


title: Reconstructing the time evolution of wall-bounded turbulent flows from non-time resolved PIV measurements

abstract: Particle Image Velocimetry (PIV) systems are often limited in their ability to fully resolve the spatiotemporal fluctuations inherent in turbulent flows due to hardware constraints. In this study, we develop models based on Rapid Distortion Theory (RDT) and Taylor's Hypothesis (TH) to reconstruct the time evolution of a turbulent flow field in the intermediate period between consecutive PIV snapshots obtained using a non-time resolved system. The linear governing equations are evolved forwards and backwards in time using the PIV snapshots as initial conditions. The flow field in the intervening period is then reconstructed by taking a weighted sum of the forward and backward estimates. This spatiotemporal weighting function is designed to account for the advective nature of the RDT and TH equations. Reconstruction accuracy is evaluated as a function of spatial resolution and reconstruction time horizon using Direct Numerical Simulation data for turbulent channel flow from the Johns Hopkins Turbulence Database. This method reconstructs single-point turbulence statistics well and resolves velocity spectra at frequencies higher than the temporal Nyquist limit of the acquisition system. Reconstructions obtained using a characteristics-based evolution of the flow field under TH prove to be more accurate compared to reconstructions obtained from numerical integration of the discretized forms of RDT and TH. The effect of measurement noise on reconstruction error is also evaluated.


title: A geometry-based model for spreading drops applied to drops on a silicon wafer and a swellable polymer brush film

abstract: We investigate the dynamics of spreading in a regime where the shape of the drop is close to a spherical cap. The latter simplification is applicable in the late (viscous) stage of spreading for highly viscous drops with a diameter below the capillary length. Moreover, it applies to the spreading of a drop on a swellable polymer brush, where the complex interaction with the substrate leads to a very slow spreading dynamics. The spherical cap geometry allows to derive a closed ordinary differential equation (ODE) for the spreading if the capillary number is a function of the contact angle as it is the case for empirical contact angle models. The latter approach has been introduced by de Gennes (Reviews of Modern Physics, 1985) for small contact angles. In the present work, we generalize the method to arbitrary contact angles. The method is applied to experimental data of spreading water-glycerol drops on a silicon wafer and spreading water drops on a PNIPAm coated silicon wafer. It is found that the ODE-model is able to describe the spreading kinetics in the case of partial wetting. Moreover, the model can predict the spreading dynamics of spherical cap-shaped droplets if the relationship between the contact angle and the capillary number is universal.


title: Managing Data Lineage of O&G Machine Learning Models: The Sweet Spot for Shale Use Case

abstract: Machine Learning (ML) has increased its role, becoming essential in several industries. However, questions around training data lineage, such as "where has the dataset used to train this model come from?"; the introduction of several new data protection legislation; and, the need for data governance requirements, have hindered the adoption of ML models in the real world. In this paper, we discuss how data lineage can be leveraged to benefit the ML lifecycle to build ML models to discover sweet-spots for shale oil and gas production, a major application in the Oil and Gas O&G Industry.


title: A Practical Approach to Navigating the Tradeoff Between Privacy and Precise Utility

abstract: Due to the recent popularity of online social networks, coupled with people's propensity to disclose personal information in an effort to achieve certain gratifications, the problem of navigating the tradeoff between privacy and utility attracted a lot of recent interest and generated a rich body of research. A critical prerequisite to solving the problem is to appropriately capture the privacy and the utility aspects in the problem formulation. Most of the existing works' focus is on the notion of privacy, while utility loss is often treated as the undesirable but necessary distortion of the true data, introduced by the privacy mechanism. By contrast, we are interested in modelling utility differently, by associating it with specific attributes of a user, just like privacy is associated with specific private attributes in the literature. Our model of utility facilitates a better and more precise privacy mechanism, and achieves better privacy-utility tradeoffs. We further incorporate into our problem formulation a practical constraint on acceptable loss in utility per unit gain in privacy, which allows users to customize the privacy mechanisms in order to account for the relative values that each user associates with their own privacy and utility. This paper discusses the intricacies of our utility model and the corresponding privacy-utility tradeoff, and introduces a heuristic greedy algorithm to solve the problem.


title: A Fractional-Order Normalized Bouc-Wen Model for Piezoelectric Hysteresis Nonlinearity

abstract: This paper presents a new fractional-order normalized Bouc-Wen (BW) (FONBW) model to describe the asymmetric and rate-dependent hysteresis nonlinearity of piezoelectric actuators (PEAs). In view of the fact that the classical BW (CBW) model is only efficient for the symmetric and rate-independent hysteresis description, the FONBW model is devoted to characterizing the asymmetric and rate-dependent behaviors of the hysteresis in PEAs by adopting an Nth-order polynomial input function and two fractional operators, respectively. Different from the traditional modified BW models, the proposed FONBW model also eliminates the redundancy of parameters in the CBW model via the normalization processing. By this way, the developed FONBW model has a relatively simple mathematic expression with fewer parameters to simultaneously characterize the asymmetric and rate-dependent hysteresis behaviors of PEAs. Model parameters are identified by the self-adaptive differential evolution algorithm. To validate the effectiveness of the proposed model, a series of model verification and inverse-multiplicative-structure-based feedforward control experiments are carried out on a PEA system. Results show that the proposed model is superior to the CBW model and traditional modified BW model in modeling accuracy and hysteresis compensation.


title: A Density version of Waring's problem

abstract: In this paper, we study a density version of Waring's problem. We prove that a positive density subset of $k$th-powers forms an asymptotic additive basis of order $O(k^2)$ provided that the relative lower density of the set is greater than $(1 - \mathcal{Z}_k^{-1}/2)^{1/k}$, where $\mathcal{Z}_k$ is certain constant depending on $k$ for which it holds that $\mathcal{Z}_k > 1$ for every $k$ and $\lim_{k \rightarrow \infty} \mathcal{Z}_k = 1$.


title: Integrating Scientific Knowledge with Machine Learning for Engineering and Environmental Systems

abstract: There is a growing consensus that solutions to complex science and engineering problems require novel methodologies that are able to integrate traditional physics-based modeling approaches with state-of-the-art machine learning (ML) techniques. This paper provides a structured overview of such techniques. Application-centric objective areas for which these approaches have been applied are summarized, and then classes of methodologies used to construct physics-guided ML models and hybrid physics-ML frameworks are described. We then provide a taxonomy of these existing techniques, which uncovers knowledge gaps and potential crossovers of methods between disciplines that can serve as ideas for future research.


title: GPU Parallelization of Policy Iteration RRT#

abstract: Sampling-based planning has become a de facto standard for complex robots given its superior ability to rapidly explore high-dimensional configuration spaces. Most existing optimal sampling-based planning algorithms are sequential in nature and cannot take advantage of wide parallelism available on modern computer hardware. Further, tight synchronization of exploration and exploitation phases in these algorithms limits sample throughput and planner performance. Policy Iteration RRT# (PI-RRT#) exposes fine-grained parallelism during the exploitation phase, but this parallelism has not yet been evaluated using a concrete implementation. We first present a novel GPU implementation of PI-RRT#'s exploitation phase and discuss data structure considerations to maximize parallel performance. Our implementation achieves 3-4x speedup over a serial PI-RRT# implementation for a 77.9% decrease in overall planning time on average. As a second contribution, we introduce the Batched-Extension RRT# algorithm, which loosens the synchronization present in PI-RRT# to realize independent 12.97x and 12.54x speedups under serial and parallel exploitation, respectively.


title: Direct measurement of the $^7$Be $L/K$ capture ratio in Ta-based superconducting tunnel junctions

abstract: We report a high-statistics measurement of the $L/K$ orbital electron capture (EC) ratio in $^7$Be embedded in cryogenic Ta. The thin Ta film formed part of a high-resolution superconducting tunnel junction (STJ) radiation detector that was used to identify the signals from different decay channels. The measured $L/K$ capture ratio of 0.070(7) is significantly larger than the only previous measurement of this quantity and the theoretical predictions that include in-medium effects. This value is a uniquely sensitive probe of the 1s and 2s orbital overlaps with the nucleus, and is of relevance to nuclear and atomic physics, as well as Li production in novae and other astrophysical scenarios. This is the first experiment that uses STJs for nuclear-recoil detection, opening a new experimental avenue for low-energy precision measurements with rare isotopes.


title: Machine-learning assisted cross-domain prediction of ionic conductivity in sodium and lithium-based superionic conductors

abstract: Solid state lithium- and sodium-ion batteries utilize solid ionicly conducting compounds as electrolytes. However, the ionic conductivity of such materials tends to be lower than their liquid counterparts, necessitating research efforts into finding suitable alternatives. The process of electrolyte screening is often based on a mixture of domain expertise and trial-and-error, both of which are time and resource-intensive. Data-driven and machine learning approaches have recently come to the fore to accelerate learnings towards discovery. In this work, we present a simple machine-learning based approach to predict the ionic conductivity of sodium and lithium-based SICON compounds. Using primarily theoretical elemental feature descriptors derivable from tabulated information on the unit cell and the atomic properties of the components of a target compound on a limited dataset of 70 NASICON-examples, we have designed a logistic regression-based model capable of distinguishing between poor and good superionic conductors with a cross-validation accuracy of over 82%. Moreover, we demonstrate how such a system is capable of cross-domain classification on lithium-based examples at the same accuracy, despite being introduced to zero lithium-based compounds during training. Through a systematic permutation-based evaluation process, we reduced the number of considered features from 47 to 7, reduction of over 83%, while simultaneously improving model performance. The contributions of different electronic and structural features to overall ionic conductivity is also discussed, and contrasted with accepted theories in literature. Our results demonstrate the utility of such a simple, yet interpretable tool provides opportunities for initial screening of potential candidates as solid-state electrolytes through the use of existing data.


title: A Review of Reduced-Order Models for Microgrids: Simplifications vs Accuracy

abstract: Inverter-based microgrids are an important technology for sustainable electrical power systems and typically use droop-controlled grid-forming inverters to interface distributed energy resources to the network and control the voltage and frequency. Ensuring stability of such microgrids is a key issue, which requires the use of appropriate models for analysis and control system design. Full-order detailed models can be more difficult to analyze and increase computational complexity, hence a number of reduced-order models have been proposed in the literature which present various trade-offs between accuracy and complexity. However, these simplifications present the risk of failing to adequately capture important dynamics of the microgrid. Therefore, there is a need for a comprehensive review and assessment of their relative quality, which is something that has not been systematically carried out thus far in the literature and we aim to address in this paper. In particular, we review various inverter-based microgrid reduced-order models and investigate the accuracy of their predictions for stability via a comparison with a corresponding detailed average model. Our study shows that the simplifications reduced order models rely upon can affect their accuracy in various regimes of the line R/X ratios, and that inappropriate model choices can result in substantially inaccurate stability results. Finally, we present recommendations on the use of reduced order models for the stability analysis of microgrids.