title: Non-Bloch PT symmetry breaking: Universal threshold and dimensional surprise

abstract: In the presence of non-Hermitian skin effect, non-Hermitian lattices generally have complex-valued eigenenergies under periodic boundary condition, but they can have non-Bloch PT symmetry and therefore completely real eigenenergies under open boundary condition. This novel PT symmetry and its breaking have been experimentally observed in one dimension. Here, we find that non-Bloch PT symmetry in two and higher dimensions exhibits drastically different behaviors compared to its one-dimensional counterpart. Whereas Bloch PT breaking and one-dimensional non-Bloch PT breaking generally have nonzero thresholds in the large-size limit, the threshold of two and higher-dimensional non-Bloch PT breaking universally approaches zero as the system size increases. A product measure, namely the product of bare non-Hermiticity and system size, is introduced to quantify the PT breaking tendency. This product being small is required for the perturbation theory to be valid, thus its growth with system size causes the breakdown of perturbation theory, which underlies the universal threshold. That the universal behaviors emerge only in two and higher dimensions indicates an unexpected interplay among PT symmetry, non-Hermitian skin effect, and spatial dimensionality. Our predictions can be confirmed on experimentally accessible platforms.

categories: cond-mat.mes-hall cond-mat.quant-gas physics.optics quant-ph


title: Mode Delocalization in Disordered Photonic Chern Insulator

abstract: In disordered two dimensional Chern insulators, a single bulk extended mode is predicted to exist per band, up to a critical disorder strength; all the other bulk modes are localized. This behavior contrasts strongly with topologically trivial two-dimensional phases, whose modes all become localized in the presence of disorder. Using a tight-binding model of a realistic photonic Chern insulator, we show that delocalized bulk eigenstates can be observed in an experimentally realistic setting. This requires the selective use of resonator losses to suppress topological edge states, and acquiring sufficiently large ensemble sizes using variable resonator detunings.


title: Testing No-Hair Theorem by Quasi-Periodic Oscillations: the quadrupole of GRO J1655$-$40

abstract: We perform an observational test of no-hair theorem using quasi-periodic oscillations within the relativistic precession model. Two well motivated metrics we apply are Kerr-Q and Hartle-Thorne metrics in which the quadrupole is the parameter that possibly encodes deviations from the Kerr black hole. The expressions for the quasi-periodic frequencies are derived before comparing the models with the observation. We encounter a degeneracy in constraining spin and quadrupole parameters that makes it difficult to measure their values. In particular, we here propose a novel test of no-hair theorem by adapting the Hartle-Thorne metric. It turns out that a Kerr black hole is a good description of the central object in GRO J1655$-$40 given the present observational precisions.


title: The Habitable-zone Planet Finder Detects a Terrestrial-mass Planet Candidate Closely Orbiting Gliese 1151: The Likely Source of Coherent Low-frequency Radio Emission from an Inactive Star

abstract: The coherent low-frequency radio emission detected by LOFAR from Gliese 1151, a quiescent M4.5 dwarf star, has radio emission properties consistent with theoretical expectations of star-planet interactions for an Earth-sized planet on a 1-5 day orbit. New near-infrared radial velocities from the Habitable-zone Planet Finder (HPF) spectrometer on the 10m Hobby-Eberly Telescope at McDonald Observatory, combined with previous velocities from HARPS-N, reveal a periodic Doppler signature consistent with an $m\sin i = 2.5 \pm 0.5 M_\oplus$ exoplanet on a 2.02-day orbit. Precise photometry from the Transiting Exoplanet Survey Satellite (TESS) shows no flares or activity signature, consistent with a quiescent M dwarf. While no planetary transit is detected in the TESS data, a weak photometric modulation is detectable in the photometry at a $\sim2$ day period. This independent detection of a candidate planet signal with the Doppler radial-velocity technique adds further weight to the claim of the first detection of star-exoplanet interactions at radio wavelengths, and helps validate this emerging technique for the detection of exoplanets.


title: Holography and black holes in asymptotically flat FLRW

abstract: In this paper, we extend the treatment of asymptotically decelerating spatially flat FLRW spacetimes initiated in [1]. We show that a certain class of those metrics is ruled by the asymptotic algebra $\mathfrak{b}\mathfrak{m}\mathfrak{s}_s$, which belongs to a one-parameter family of deformations of $\mathfrak{b}\mathfrak{m}\mathfrak{s}$. Furthermore, we enlarge our ansatz to include $\text{Diff}(S^2)$ transformations whose asymptotic algebra $\mathfrak{g}\mathfrak{b}\mathfrak{m}\mathfrak{s}_s$ is a one parameter deformation of $\mathfrak{g}\mathfrak{b}\mathfrak{m}\mathfrak{s}$. Therefore, the holographic algebras $\mathfrak{b}\mathfrak{m}\mathfrak{s}_s$ and $\mathfrak{g}\mathfrak{b}\mathfrak{m}\mathfrak{s}_s$ in FLRW can be related to their flat counterparts through a cosmological holographic flow. Finally, we introduce a logarithmic ansatz in order to account for cosmological black holes, which does not generally satisfy the peeling property but preserves the asymptotic algebra.


title: Characterizing the dynamical phase diagram of the Dicke model via classical and quantum probes

abstract: We theoretically study the dynamical phase diagram of the Dicke model in both classical and quantum limits using large, experimentally relevant system sizes. Our analysis elucidates that the model features dynamical critical points that are distinct from previously investigated excited-state equilibrium transitions. Moreover, our numerical calculations demonstrate that mean-field features of the dynamics remain valid in the exact quantum dynamics, but we also find that in regimes where quantum effects dominate signatures of the dynamical phases and chaos can persist in purely quantum metrics such as entanglement and correlations. Our predictions can be verified in current quantum simulators of the Dicke model including arrays of trapped ions.


title: Codimension one Ricci soliton subgroups of solvable Iwasawa groups

abstract: Recently, Jablonski proved that, to a large extent, a simply connected solvable Lie group endowed with a left-invariant Ricci soliton metric can be isometrically embedded into the solvable Iwasawa group of a non-compact symmetric space. Motivated by this result, we classify codimension one subgroups of the solvable Iwasawa groups of irreducible symmetric spaces of non-compact type whose induced metrics are Ricci solitons. We also obtain the classifications of codimension one Ricci soliton subgroups of Damek-Ricci spaces and generalized Heisenberg groups.


title: Generation of attached Langmuir circulations by a suspended macroalgal farm

abstract: In this study, we focus on Langmuir turbulence in the deep ocean with the presence of a large macroalgal farm using a Large Eddy Simulation method. The wave-current interactions are modelled by solving the wave-averaged equations. The hydrodynamic process over the farm is found to drive a persistent flow pattern similar to Langmuir circulations but is locked in space across the farm. These secondary circulations are generated because the cross-stream shear produced by the rows of canopy elements leads to a steady vertical vorticity field, which is then rotated to the downstream direction under the effect of vortex force. Since the driving mechanism is similar to the CraikLeibovich type 2 instability theory, these secondary circulations are also termed as attached Langmuir circulations. We then apply a triple decomposition on the flow field to unveil the underlying kinematics and energy transfer between the mean flow, the secondary flow resulting from the farm drag, and the transient eddies. Flow visualizations and statistics suggest that the attached Langmuir circulations result from the adjustment of the upper ocean mixed layer to the macroalgal farm, and they will weaken (if not disappear) when the flow reaches an equilibrium state within the farm. The tripledecomposed energy budgets reveal that the energy of the secondary flow is transferred from the mean flow under the action of canopy drag, while the transient eddies feed on wave energy transferred by the Stokes drift and energy conversion from the secondary flow.


title: Measurements of the differential cross sections of the production of Z + jets and $\gamma$ + jets and of Z boson emission collinear with a jet in pp collisions at $\sqrt{s} =$ 13 TeV

abstract: Measurements of the differential cross sections of Z + jets and $\gamma$ + jets production, and their ratio, are presented as a function of the boson transverse momentum. Measurements are also presented of the angular distribution between the Z boson and the closest jet. The analysis is based on pp collisions at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 35.9 fb$^{-1}$ recorded by the CMS experiment at the LHC. The results, corrected for detector effects, are compared with various theoretical predictions. In general, the predictions at higher orders in perturbation theory show better agreement with the measurements. This work provides the first measurement of the ratio of the differential cross sections of Z + jets and $\gamma$ + jets production at 13 TeV, as well as the first direct measurement of Z bosons emitted collinearly with a jet.


title: Late-Time Cosmology of Scalar-Coupled $f(R, \mathcal{G})$ Gravity

abstract: In this work by using a numerical analysis, we investigate in a quantitative way the late-time dynamics of scalar coupled $f(R,\mathcal{G})$ gravity. Particularly, we consider a Gauss-Bonnet term coupled to the scalar field coupling function $\xi(\phi)$, and we study three types of models, one with $f(R)$ terms that are known to provide a viable late-time phenomenology, and two Einstein-Gauss-Bonnet types of models. Our aim is to write the Friedmann equation in terms of appropriate statefinder quantities frequently used in the literature, and we numerically solve it by using physically motivated initial conditions. In the case that $f(R)$ gravity terms are present, the contribution of the Gauss-Bonnet related terms is minor, as we actually expected. This result is robust against changes in the initial conditions of the scalar field, and the reason is the dominating parts of the $f(R)$ gravity sector at late times. In the Einstein-Gauss-Bonnet type of models, we examine two distinct scenarios, firstly by choosing freely the scalar potential and the scalar Gauss-Bonnet coupling $\xi(\phi)$, in which case the resulting phenomenology is compatible with the latest Planck data and mimics the $\Lambda$-Cold-Dark-Matter model. In the second case, since there is no fundamental particle physics reason for the graviton to change its mass, we assume that primordially the tensor perturbations propagate with the speed equal to that of light's, and thus this constraint restricts the functional form of the scalar coupling function $\xi(\phi)$, which must satisfy the differential equation $\ddot{\xi}=H\dot{\xi}$.


title: Harvest -- An Open Source Toolkit for Extracting Posts and Post Metadata from Web Forums

abstract: Automatic extraction of forum posts and metadata is a crucial but challenging task since forums do not expose their content in a standardized structure. Content extraction methods, therefore, often need customizations such as adaptations to page templates and improvements of their extraction code before they can be deployed to new forums. Most of the current solutions are also built for the more general case of content extraction from web pages and lack key features important for understanding forum content such as the identification of author metadata and information on the thread structure. This paper, therefore, presents a method that determines the XPath of forum posts, eliminating incorrect mergers and splits of the extracted posts that were common in systems from the previous generation. Based on the individual posts further metadata such as authors, forum URL and structure are extracted. We also introduce Harvest, a new open source toolkit that implements the presented methods and create a gold standard extracted from 52 different Web forums for evaluating our approach. A comprehensive evaluation reveals that Harvest clearly outperforms competing systems.


title: The nature and origins of the low surface brightness outskirts of massive, central galaxies in Subaru HSC

abstract: We explore the stellar mass density and colour profiles of 118 low redshift, massive, central galaxies, selected to have assembled 90 percent of their stellar mass 6 Gyr ago, finding evidence of the minor merger activity expected to be the driver behind the size growth of quiescent galaxies. We use imaging data in the $g, r, i, z, y$ bands from the Subaru Hyper Suprime-Cam survey and perform SED fitting to construct spatially well-resolved radial profiles in colour and stellar mass surface density. Our visual morphological classification reveals that $\sim 42$ percent of our sample displays tidal features, similar to previous studies, $\sim 43$ percent of the remaining sample display a diffuse stellar halo and only $\sim 14$ percent display no features, down to a limiting $\mu_{r\mathrm{-band}}$ $\sim$ 28 mag arcsec$^{-2}$. We find good agreement between the stacked colour profiles of our sample to those derived from previous studies and an expected smooth, declining stellar mass surface density profile in the central regions (< 3 R$_{\mathrm{e}}$). However, we also see a flattening of the profile ($\Sigma_* \sim 10^{7.5}$ M$_\odot$ kpc$^{-2}$) in the outskirts (up to 10 R$_{\mathrm{e}}$), which is revealed by our method of specifically targeting tidal/accretion features. We find similar levels of tidal features and behaviour in the stellar mass surface density profiles in a younger comparison sample, however a lack of diffuse haloes. We also apply stacking techniques, similar to those in previous studies, finding such procedures wash out tidal features and thereby produces smooth declining profiles. The stellar material in the outskirts contributes on average $\sim 10^{10}$ M$_\odot$ or a few percent of the total stellar mass and has similar colours to SDSS satellites of similar stellar mass.


title: Chiral phase-coexistence in compressed double-twist elastomers

abstract: We adapt the theory of anisotropic rubber elasticity to model cross-linked double-twist liquid crystal cylinders such as exhibited in biological systems. In mechanical extension we recover strain-straightening, but with an exact expression in the small twist-angle limit. In compression, we observe coexistence between high and low twist phases. Coexistence begins at small compressive strains and is robustly observed for any anisotropic cross-links and for general double-twist functions -- but disappears at large twist angles. Within the coexistence region, significant compression of double-twist cylinders is allowed at constant stress. Our results are qualitatively consistent with previous observations of swollen or compressed collagen fibrils, indicating that this phenomenon may be readily accessible experimentally.


title: Information-theoretic Key Encapsulation and its Applications

abstract: A hybrid encryption scheme is a public-key encryption system that consists of a public-key part called the key encapsulation mechanism (KEM), and a (symmetric) secret-key part called data encapsulation mechanism (DEM): the public-key part is used to generate a shared secret key between two parties, and the symmetric key part is used to encrypt the message using the generated key. Hybrid encryption schemes are widely used for secure communication over the Internet. In this paper, we initiate the study of hybrid encryption in preprocessing model which assumes access to initial correlated variables by all parties (including the eavesdropper). We define information-theoretic KEM (iKEM) that, together with a (computationally) secure DEM, results in a hybrid encryption scheme in preprocessing model. We define the security of each building block, and prove a composition theorem that guarantees (computational) qe-chosen plaintext (CPA) security of the hybrid encryption system if the iKEM and the DEM satisfy qe-chosen encapculation attack and one-time security, respectively. We show that iKEM can be realized by a one-way SKA (OW-SKA) protocol with a revised security definition. Using an OW-SKA that satisfies this revised definition of security effectively allows the secret key that is generated by the OW-SKA to be used with a one-time symmetric key encryption system such as XORing a pseudorandom string with the message, and provide qe-CPA security for the hybrid encryption system.We discuss our results and directions for future work.


title: Bounds and Genericity of Sum-Rank-Metric Codes

abstract: We derive simplified sphere-packing and Gilbert--Varshamov bounds for codes in the sum-rank metric, which can be computed more efficiently than previous ones. They give rise to asymptotic bounds that cover the asymptotic setting that has not yet been considered in the literature: families of sum-rank-metric codes whose block size grows in the code length. We also provide two genericity results: we show that random linear codes achieve almost the sum-rank-metric Gilbert--Varshamov bound with high probability. Furthermore, we derive bounds on the probability that a random linear code attains the sum-rank-metric Singleton bound, showing that for large enough extension fields, almost all linear codes achieve it.


title: Siegel modular forms of degree two and three and invariant theory

abstract: This is a survey based on the construction of Siegel modular forms of degree 2 and 3 using invariant theory in joint work with Fabien Cl\'ery and Carel Faber.


title: The Forgotten Document-Oriented Database Management Systems: An Overview and Benchmark of Native XML DODBMSes in Comparison with JSON DODBMSes

abstract: In the current context of Big Data, a multitude of new NoSQL solutions for storing, managing, and extracting information and patterns from semi-structured data have been proposed and implemented. These solutions were developed to relieve the issue of rigid data structures present in relational databases, by introducing semi-structured and flexible schema design. As current data generated by different sources and devices, especially from IoT sensors and actuators, use either XML or JSON format, depending on the application, database technologies that store and query semi-structured data in XML format are needed. Thus, Native XML Databases, which were initially designed to manipulate XML data using standardized querying languages, i.e., XQuery and XPath, were rebranded as NoSQL Document-Oriented Databases Systems. Currently, the majority of these solutions have been replaced with the more modern JSON based Database Management Systems. However, we believe that XML-based solutions can still deliver performance in executing complex queries on heterogeneous collections. Unfortunately nowadays, research lacks a clear comparison of the scalability and performance for database technologies that store and query documents in XML versus the more modern JSON format. Moreover, to the best of our knowledge, there are no Big Data-compliant benchmarks for such database technologies. In this paper, we present a comparison for selected Document-Oriented Database Systems that either use the XML format to encode documents, i.e., BaseX, eXist-db, and Sedna, or the JSON format, i.e., MongoDB, CouchDB, and Couchbase. To underline the performance differences we also propose a benchmark that uses a heterogeneous complex schema on a large DBLP corpus.


title: Low-cost attacks on Ethereum 2.0 by sub-1/3 stakeholders

abstract: We outline two dishonest strategies that can be cheaply executed on the Ethereum 2.0 beacon chain, even by validators holding less than one-third of the total stake: malicious chain reorganizations ("reorgs") and finality delays. In a malicious reorg, an attacker withholds their blocks and attestations before releasing them at an opportune time in order to force a chain reorganization, which they can take advantage of by double-spending or front-running transactions. To execute a finality delay an attacker uses delayed block releases and withholding of attestations to increase the mean and variance of the time it takes blocks to become finalized. This impacts the efficiency and predictability of the system. We provide a probabilistic and cost analysis for each of these attacks, considering a validator with 30% of the total stake.


title: Reduced-order Model for Fluid Flows via Neural Ordinary Differential Equations

abstract: Reduced order models play an important role in the design, optimization and control of dynamical systems. In recent years, there has been an increasing interest in the application of data-driven techniques for model reduction that can decrease the computational burden of numerical solutions, while preserving the most important features of complex physical problems. In this paper, we use the proper orthogonal decomposition to reduce the dimensionality of the model and introduce a novel generative neural ODE (NODE) architecture to forecast the behavior of the temporal coefficients. With this methodology, we replace the classical Galerkin projection with an architecture characterized by the use of a continuous latent space. We exemplify the methodology on the dynamics of the Von Karman vortex street of the flow past a cylinder generated by a Large-eddy Simulation (LES)-based code. We compare the NODE methodology with an LSTM baseline to assess the extrapolation capabilities of the generative model and present some qualitative evaluations of the flow reconstructions.


title: Bubbling transition as a mechanism of destruction of synchronous oscillations of identical microbubble contrast agents

abstract: We study the process of destruction of synchronous oscillations in a model of two interacting microbubble contrast agents exposed to an external ultrasound field. Completely synchronous oscillations in this model are possible in case of identical bubbles, when the governing system of equations possess a symmetry leading to the existence of a synchronization manifold. Such synchronous oscillations can be destructed without breaking the corresponding symmetry of the governing dynamical system. Here we describe the phenomenological mechanism responsible for such destruction of synchronization and demonstrate its implementation in the studied model. We show that the appearance and expansion of transversally unstable areas in the synchronization manifold leads to transformation of a synchronous chaotic attractor into a hyperchaotic one. We also demonstrate that this bifurcation sequence is stable with respect to symmetry breaking perturbations.


title: Dietary Supplements and Nutraceuticals Under Investigation for COVID-19 Prevention and Treatment

abstract: Coronavirus disease 2019 (COVID-19) has caused global disruption and a significant loss of life. Existing treatments that can be repurposed as prophylactic and therapeutic agents could reduce the pandemic's devastation. Emerging evidence of potential applications in other therapeutic contexts has led to the investigation of dietary supplements and nutraceuticals for COVID-19. Such products include vitamin C, vitamin D, omega 3 polyunsaturated fatty acids, probiotics, and zinc, all of which are currently under clinical investigation. In this review, we critically appraise the evidence surrounding dietary supplements and nutraceuticals for the prophylaxis and treatment of COVID-19. Overall, further study is required before evidence-based recommendations can be formulated, but nutritional status plays a significant role in patient outcomes, and these products could help alleviate deficiencies. For example, evidence indicates that vitamin D deficiency may be associated with greater incidence of infection and severity of COVID-19, suggesting that vitamin D supplementation may hold prophylactic or therapeutic value. A growing number of scientific organizations are now considering recommending vitamin D supplementation to those at high risk of COVID-19. Because research in vitamin D and other nutraceuticals and supplements is preliminary, here we evaluate the extent to which these nutraceutical and dietary supplements hold potential in the COVID-19 crisis.


title: Gamma-ray detection toward the Coma cluster with Fermi-LAT: Implications for the cosmic ray content in the hadronic scenario

abstract: The presence of relativistic electrons within the diffuse gas phase of galaxy clusters is now well established, but their detailed origin remains unclear. Cosmic ray protons are also expected to accumulate during the formation of clusters and would lead to gamma-ray emission through hadronic interactions within the thermal gas. Recently, the detection of gamma-ray emission has been reported toward the Coma cluster with Fermi-LAT. Assuming that this gamma-ray emission arises from hadronic interactions in the ICM, we aim at exploring the implication of this signal on the cosmic ray content of the Coma cluster. We use the MINOT software to build a physical model of the cluster and apply it to the Fermi-LAT data. We also consider contamination from compact sources and the impact of various systematic effects. We confirm that a significant gamma-ray signal is observed within the characteristic radius $\theta_{500}$ of the Coma cluster, with a test statistic TS~27 for our baseline model. The presence of a possible point source may account for most of the observed signal. However, this source could also correspond to the peak of the diffuse emission of the cluster itself and extended models match the data better. We constrain the cosmic ray to thermal energy ratio within $R_{500}$ to $X_{\rm CRp}=1.79^{+1.11}_{-0.30}$\% and the slope of the energy spectrum of cosmic rays to $\alpha=2.80^{+0.67}_{-0.13}$. Finally, we compute the synchrotron emission associated with the secondary electrons produced in hadronic interactions assuming steady state. This emission is about four times lower than the overall observed radio signal, so that primary cosmic ray electrons or reacceleration of secondary electrons is necessary to explain the total emission. Assuming an hadronic origin of the signal, our results provide the first quantitative measurement of the cosmic ray proton content in a cluster.[Abridged]


title: Source Independent Quantum Walk Random Number Generation

abstract: Source independent quantum random number generators (SI-QRNG) are cryptographic protocols which attempt to extract random strings from quantum sources where the source is under the control of an adversary while the measurement devices are fully characterized. This represents a middle-ground between fully-trusted and full-device independence, allowing for fast bit-generation rates with current-day technology, while also providing a strong security guarantee. In this paper we analyze an SI-QRNG protocol based on quantum walks and develop a new proof technique to show security.


title: Higher Spin Gravities and Presymplectic AKSZ Models

abstract: As a step towards quantization of Higher Spin Gravities we construct the presymplectic AKSZ sigma-model for $4d$ Higher Spin Gravity which is AdS/CFT dual of Chern-Simons vector models. It is shown that the presymplectic structure leads to the correct quantum commutator of higher spin fields and to the correct algebra of the global higher spin symmetry currents. The presymplectic AKSZ model is proved to be unique, it depends on two coupling constants in accordance with the AdS/CFT duality, and it passes some simple checks of interactions.


title: A self-calibrating superconducting pair-breaking detector

abstract: We propose and experimentally demonstrate a self-calibrating detector of Cooper pair depairing in a superconductor based on a mesoscopic superconducting island coupled to normal metal leads. On average, exactly one electron passes through the device per broken Cooper pair, independent of the absorber volume, device or material parameters. The device operation is explained by a simple analytical model and verified with numerical simulations in quantitative agreement with experiment. In a proof-of-concept experiment, we use such a detector to measure the high-frequency phonons generated by another, electrically decoupled superconducting island, with a measurable signal resulting from less than 10 fW of dissipated power.


title: Probabilistic Interval Analysis for the Analytic Prediction of the Pattern Tolerance Distribution in Linear Phased Arrays With Random Excitation Errors

abstract: A statistical approach based on the interval analysis (IA) is proposed for the analysis of the effects, on the radiation patterns radiated by phased arrays, of random errors and tolerances in the amplitudes and phases of the array-elements excitations. Starting from the efficient, reliable, and inclusive computation of the bounds of the complex-valued interval array pattern function by means of IA, an analytic method is presented to yield closed-form expressions for the probability of occurrence of a user-chosen value of the power pattern or of its features within the corresponding IA-derived bounds. A set of numerical examples is reported and discussed to assess the reliability of the proposed probabilistic interval analysis (PIA) method with the results from Monte Carlo simulations as well as to point out its effectiveness and potentialities/advantages/efficiency in real applications of great industrial interest.


title: Exact Diagonalization for Magic-Angle Twisted Bilayer Graphene

abstract: We report on finite-size exact-diagonalization calculations in a Hilbert space defined by the continuum-model flat moir\'e bands of magic angle twisted bilayer graphene (MATBG). For moir\'e band filling $3>|\nu|>2$, where superconductivity is strongest, we obtain evidence that the ground state is a spin ferromagnet. Near $|\nu|=3$, we find Chern insulator ground states that have spontaneous spin, valley, and sublattice polarization, and demonstrate that the anisotropy energy in this order-parameter space is strongly band-filling-factor dependent. We emphasize that inclusion of the remote band self-energy is necessary for a reliable description of MATBG flat band correlations.


title: Cuprate superconductors as viewed through a striped lens

abstract: Understanding the electron pairing in hole-doped cuprate superconductors has been a challenge, in particular because the "normal" state from which it evolves is unprecedented. Now, after three and a half decades of research, involving a wide range of experimental characterizations, it is possible to delineate a clear and consistent cuprate story. It starts with doping holes into a charge-transfer insulator, resulting in in-gap states. These states exhibit a pseudogap resulting from the competition between antiferromagnetic superexchange $J$ between nearest-neighbor Cu atoms (a real-space interaction) and the kinetic energy of the doped holes, which, in the absence of interactions, would lead to extended Bloch-wave states whose occupancy is characterized in reciprocal space. To develop some degree of coherence on cooling, the spin and charge correlations must self-organize in a cooperative fashion. A specific example of resulting emergent order is that of spin and charge stripes, as observed in La$_{2-x}$Ba$_x$CuO$_4$. While stripe order frustrates bulk superconductivity, it nevertheless develops pairing and superconducting order of an unusual character. The antiphase order of the spin stripes decouples them from the charge stripes, which can be viewed as hole-doped, two-leg, spin-$\frac12$ ladders. To achieve superconducting order, the pair correlations in neighboring ladders must develop phase order. In the presence of spin stripe order, antiphase Josephson coupling can lead to pair-density-wave superconductivity. Alternatively, in-phase superconductivity requires that the spin stripes have an energy gap, which empirically limits the coherent superconducting gap. Hence, superconducting order in the cuprates involves a compromise between the pairing scale, which is maximized at $x\sim\frac18$, and phase coherence, which is optimized at $x\sim0.2$.


title: Optimal Excitation Matching Strategy for Sub-Arrayed Phased Linear Arrays Generating Arbitrary-Shaped Beams

abstract: The design of phased arrays able to generate arbitrary-shaped beams through a sub-arrayed architecture is addressed here. The synthesis problem is cast in the excitation matching framework, so as to yield clustered phased arrays providing optimal trade-offs between the complexity of the array architecture (i.e., the minimum number of control points at the sub-array level) and the matching of a reference pattern. A synthesis tool based on the k-means algorithm is proposed for jointly optimizing the sub-array configuration and the complex sub-array coefficients. Selected numerical results, including pencil beams with sidelobe notches and asymmetric lobes as well as shaped main lobes, are reported and discussed to highlight the peculiarities of the proposed approach also in comparison with some extensions to complex excitations of state-of-the-art sub-array design methods.


title: Interspecific competition shapes the structural stability of mutualistic networks

abstract: Mutualistic networks have attracted increasing attention in the ecological literature in the last decades as they play a key role in the maintenance of biodiversity. Here, we develop an analytical framework to study the structural stability of these networks including both mutualistic and competitive interactions. Analytical and numerical analyses show that the structure of the competitive network fundamentally alters the necessary conditions for species coexistence in communities. Using 50 real mutualistic networks, we show that when the relative importance of shared partners is incorporated via weighted competition, the feasibility area in the parameter space is highly correlated with May's stability criteria and can be predicted by a functional relationship between the number of species, the network connectance and the average interaction strength in the community. Our work reopens a decade-long debate about the complexity-stability relationship in ecological communities, and highlights the role of the relative structures of different interaction types.


title: Proper Scoring Rules and Domination

abstract: I generalize a theorem of Predd, et al.~(2009) on domination and strictly proper scoring rules to the case of non-additive scoring rules.


title: An Irregular Two-Sizes Square Tiling Method for the Design of Isophoric Phased Arrays

abstract: The design of isophoric phased arrays composed of two-sized square-shaped tiles that fully cover rectangular apertures is dealt with. The number and the positions of the tiles within the array aperture are optimized to fit desired specifications on the power pattern features. Toward this end, starting from the derivation of theoretical conditions for the complete tileability of the aperture, an ad hoc coding of the admissible arrangements, which implies a drastic reduction of the cardinality of the solution space, and their compact representation with a graph are exploited to profitably apply an effective optimizer based on an integer-coded genetic algorithm. A set of representative numerical examples, concerned with state-of-the-art benchmark problems, is reported and discussed to give some insights on the effectiveness of both the proposed tiled architectures and the synthesis strategy.


title: Modular Design of Hexagonal Phased Arrays Through Diamond Tiles

abstract: The modular design of planar phased array antennas with hexagonal apertures is addressed by means of innovative diamond-shaped tiling techniques. Both tiling configuration and subarray coefficients are optimized to fit user-defined power-mask constraints on the radiation pattern. Toward this end, suitable surface-tiling mathematical theorems are customized to the problem at hand to guarantee optimal performance in case of low/medium-size arrays, while the computationally hard tiling of large arrays is yielded thanks to an effective integer-coded GA-based exploration of the arising high-cardinality solution spaces. By considering ideal as well as real array models, a set of representative benchmark problems is dealt with to assess the effectiveness of the proposed architectures and tiling strategies. Moreover, comparisons with alternative tiling architectures are also performed to show to the interested readers the advantages and the potentialities of the diamond subarraying of hexagonal apertures.


title: Robust model benchmarking and bias-imbalance in data-driven materials science: a case study on MODNet

abstract: As the number of novel data-driven approaches to material science continues to grow, it is crucial to perform consistent quality, reliability and applicability assessments of model performance. In this paper, we benchmark the Materials Optimal Descriptor Network (MODNet) method and architecture against the recently released MatBench v0.1, a curated test suite of materials datasets. MODNet is shown to outperform current leaders on 6 of the 13 tasks, whilst closely matching the current leaders on a further 2 tasks; MODNet performs particularly well when the number of samples is below 10,000. Attention is paid to two topics of concern when benchmarking models. First, we encourage the reporting of a more diverse set of metrics as it leads to a more comprehensive and holistic comparison of model performance. Second, an equally important task is the uncertainty assessment of a model towards a target domain. Significant variations in validation errors can be observed, depending on the imbalance and bias in the training set (i.e., similarity between training and application space). By using an ensemble MODNet model, confidence intervals can be built and the uncertainty on individual predictions can be quantified. Imbalance and bias issues are often overlooked, and yet are important for successful real-world applications of machine learning in materials science and condensed matter.


title: Perverse Leray filtration and specialization with applications to the Hitchin morphism

abstract: We initiate and develop a framework to handle the specialization morphism as a filtered morphism for the perverse, and for the perverse Leray filtration, on the cohomology with constructible coefficients of varieties and morphisms parameterized by a curve. As an application, we use this framework to carry out a detailed study of filtered specialization for the Hitchin morphisms associated with the compactification of Dolbeault moduli spaces in [de-2018].


title: Power laws and phase transitions in heterogenous car following with reaction times

abstract: We study the effect of reaction times on the kinetics of relaxation to stationary states and on congestion transitions in heterogeneous traffic. Heterogeneity is modeled as quenched disorders in the parameters of the car following model and in the reaction times of the drivers. We observed that at low densities, the relaxation to stationary state from a homogeneous initial state is governed by the same power laws as derived by E. Ben-Naim et al., Kinetics of clustering in traffic flow, Phys. Rev. E 50, 822 (1994). The stationary state, at low densities, is a single giant platoon of vehicles with the slowest vehicle being the leader. We observed formation of spontaneous jams inside the giant platoon which move upstream as stop-go waves and dissipate at its tail. The transition happens when the head of the giant platoon interacts with its tail, stable stop-go waves form, which circulate in the ring without dissipating. We observed that the system behaves differently when the transition density is approached from above that it does when approached from below. When the transition density is approached from below, the gap distribution behind the leader has a double peak and is fat-tailed but has a bounded support and thus the maximum gap in the system and the variance of the gap distribution tend to size-independent values. When the transition density is approached from above, the gap distribution becomes a power law and, consequently, the maximum gap in the system and the variance in the gaps diverge as a power law, thereby creating a discontinuity at the transition. Thus, we observe a phase transition of unusual kind in which both a discontinuity and a power law are observed at the transition density. These unusual features vanish in the absence of reaction time (e.g., automated driving).


title: Reconstruction of Bloch wavefunctions of holes in a semiconductor

abstract: A central goal of condensed-matter physics is to understand how the diverse electronic and optical properties of crystalline materials emerge from the wavelike motion of electrons through periodically arranged atoms. However, more than 90 years after Bloch derived the functional forms of electronic waves in crystals [1] (now known as Bloch wavefunctions), rapid scattering processes have so far prevented their direct experimental reconstruction. In high-order sideband generation [2-9], electrons and holes generated in semiconductors by a near-infrared laser are accelerated to a high kinetic energy by a strong terahertz field, and recollide to emit near-infrared sidebands before they are scattered. Here we reconstruct the Bloch wavefunctions of two types of hole in gallium arsenide at wavelengths much longer than the spacing between atoms by experimentally measuring sideband polarizations and introducing an elegant theory that ties those polarizations to quantum interference between different recollision pathways. These Bloch wavefunctions are compactly visualized on the surface of a sphere. High-order sideband generation can, in principle, be observed from any direct-gap semiconductor or insulator. We thus expect that the method introduced here can be used to reconstruct low-energy Bloch wavefunctions in many of these materials, enabling important insights into the origin and engineering of the electronic and optical properties of condensed matter.


title: DEFT: Detection Embeddings for Tracking

abstract: Most modern multiple object tracking (MOT) systems follow the tracking-by-detection paradigm, consisting of a detector followed by a method for associating detections into tracks. There is a long history in tracking of combining motion and appearance features to provide robustness to occlusions and other challenges, but typically this comes with the trade-off of a more complex and slower implementation. Recent successes on popular 2D tracking benchmarks indicate that top-scores can be achieved using a state-of-the-art detector and relatively simple associations relying on single-frame spatial offsets -- notably outperforming contemporary methods that leverage learned appearance features to help re-identify lost tracks. In this paper, we propose an efficient joint detection and tracking model named DEFT, or "Detection Embeddings for Tracking." Our approach relies on an appearance-based object matching network jointly-learned with an underlying object detection network. An LSTM is also added to capture motion constraints. DEFT has comparable accuracy and speed to the top methods on 2D online tracking leaderboards while having significant advantages in robustness when applied to more challenging tracking data. DEFT raises the bar on the nuScenes monocular 3D tracking challenge, more than doubling the performance of the previous top method. Code is publicly available.


title: A Heuristic for Dynamic Output Predictive Control Design for Uncertain Nonlinear Systems

abstract: In this paper, a simple heuristic is proposed for the design of uncertainty aware predictive controllers for nonlinear models involving uncertain parameters. The method relies on Machine Learning-based approximation of ideal deterministic MPC solutions with perfectly known parameters. An efficient construction of the learning data set from these off-line solutions is proposed in which each solution provides many samples in the learning data. This enables a drastic reduction of the required number of Non Linear Programming problems to be solved off-line while explicitly exploiting the statistics of the parameters dispersion. The learning data is then used to design a fast on-line output dynamic feedback that explicitly incorporate information of the statistics of the parameters dispersion. An example is provided to illustrate the efficiency and the relevance of the proposed framework. It is in particular shown that the proposed solution recovers up to 78\% of the expected advantage of having a perfect knowledge of the parameters compared to nominal design.


title: Numerical Differentiation using local Chebyshev-Approximation

abstract: In applied mathematics, especially in optimization, functions are often only provided as so called "Black-Boxes" provided by software packages, or very complex algorithms, which make automatic differentation very complicated or even impossible. Hence one seeks the numerical approximation of the derivative. Unfortunately numerical differentation is a difficult task in itself, and it is well known that it is numerical instable. There are many works on this topic, including the usage of (global) Chebyshev approximations. Chebyshev approximations have the great property that they converge very fast, if the function is smooth. Nevertheless those approches have several drawbacks, since in practice functions are not smooth, and a global approximation needs many function evalutions. Nevertheless there is hope. Since functions in real world applications are most times smooth except for finite points, corners or edges. This motivates to use a local Chebyshev approach, where the function is only approximated locally, and hence the Chebyshev approximations still yields a fast approximation of the desired function. We will study such an approch in this work, and will provide a numerical example


title: Confusion2vec 2.0: Enriching Ambiguous Spoken Language Representations with Subwords

abstract: Word vector representations enable machines to encode human language for spoken language understanding and processing. Confusion2vec, motivated from human speech production and perception, is a word vector representation which encodes ambiguities present in human spoken language in addition to semantics and syntactic information. Confusion2vec provides a robust spoken language representation by considering inherent human language ambiguities. In this paper, we propose a novel word vector space estimation by unsupervised learning on lattices output by an automatic speech recognition (ASR) system. We encode each word in confusion2vec vector space by its constituent subword character n-grams. We show the subword encoding helps better represent the acoustic perceptual ambiguities in human spoken language via information modeled on lattice structured ASR output. The usefulness of the proposed Confusion2vec representation is evaluated using semantic, syntactic and acoustic analogy and word similarity tasks. We also show the benefits of subword modeling for acoustic ambiguity representation on the task of spoken language intent detection. The results significantly outperform existing word vector representations when evaluated on erroneous ASR outputs. We demonstrate that Confusion2vec subword modeling eliminates the need for retraining/adapting the natural language understanding models on ASR transcripts.


title: Can a single PDE govern well the propagation of the electric wave field in a heterogeneous medium in 3D?

abstract: It is demonstrated in this paper that the propagation of the electric wave field in a heterogeneous medium in 3D can sometimes be governed well by a single PDE, which is derived from the Maxwell's equations. The corresponding component of the electric field dominates two other components. This justifies some past results of the second author with coauthors about numerical solutions of coefficient inverse problems with experimental electromagnetic data. In addition, since it is simpler to work in applications with a single PDE rather than with the complete Maxwell's system, then the result of this paper might be useful to researchers working on applied issues of the propagation of electromagnetic waves in inhomogeneous media.


title: Coupled Hypergraph Maps and Chaotic Cluster Synchronization

abstract: Coupled map lattices (CMLs) are prototypical dynamical systems on networks/graphs. They exhibit complex patterns generated via the interplay of diffusive/Laplacian coupling and nonlinear reactions modelled by a single iterated map at each node; the maps are often taken as unimodal, e.g., logistic or tent maps. In this letter, we propose a class of higher-order coupled dynamical systems involving the hypergraph Laplacian, which we call coupled hypergraph maps (CHMs). By combining linearized (in-)stability analysis of synchronized states, hypergraph spectral theory, and numerical methods, we detect robust regions of chaotic cluster synchronization occurring in parameter space upon varying coupling strength and the main bifurcation parameter of the unimodal map. Furthermore, we find key differences between Laplacian and hypergraph Laplacian coupling and detect various other classes of periodic and quasi-periodic patterns. The results show the high complexity of coupled graph maps and indicate that they might be an excellent universal model class to understand the similarities and differences between dynamics on classical graphs and dynamics on hypergraphs.


title: Stability and performance verification of dynamical systems controlled by neural networks: algorithms and complexity

abstract: This work makes several contributions on stability and performance verification of nonlinear dynamical systems controlled by neural networks. First, we show that the stability and performance of a polynomial dynamical system controlled by a neural network with semialgebraically representable activation functions (e.g., ReLU) can be certified by convex semidefinite programming. The result is based on the fact that the semialgebraic representation of the activation functions and polynomial dynamics allows one to search for a Lyapunov function using polynomial sum-of-squares methods. Second, we remark that even in the case of a linear system controlled by a neural network with ReLU activation functions, the problem of verifying asymptotic stability is undecidable. Finally, under additional assumptions, we establish a converse result on the existence of a polynomial Lyapunov function for this class of systems. Numerical results with code available online on examples of state-space dimension up to 50 and neural networks with several hundred neurons and up to 30 layers demonstrate the method.


title: Neural Recursive Belief States in Multi-Agent Reinforcement Learning

abstract: In multi-agent reinforcement learning, the problem of learning to act is particularly difficult because the policies of co-players may be heavily conditioned on information only observed by them. On the other hand, humans readily form beliefs about the knowledge possessed by their peers and leverage beliefs to inform decision-making. Such abilities underlie individual success in a wide range of Markov games, from bluffing in Poker to conditional cooperation in the Prisoner's Dilemma, to convention-building in Bridge. Classical methods are usually not applicable to complex domains due to the intractable nature of hierarchical beliefs (i.e. beliefs of other agents' beliefs). We propose a scalable method to approximate these belief structures using recursive deep generative models, and to use the belief models to obtain representations useful to acting in complex tasks. Our agents trained with belief models outperform model-free baselines with equivalent representational capacity using common training paradigms. We also show that higher-order belief models outperform agents with lower-order models.


title: Modelling and simulation of a wave energy converter

abstract: In this work we present the mathematical model and simulations of a particular wave energy converter, the so-called oscillating water column. In this device, waves governed by the one-dimensional nonlinear shallow water equations arrive from offshore, encounter a step in the bottom and then arrive into a chamber to change the volume of the air to activate the turbine. The system is reformulated as two transmission problems: one is related to the wave motion over the stepped topography and the other one is related to the wave-structure interaction at the entrance of the chamber. We finally use the characteristic equations of Riemann invariants to obtain the discretized transmission conditions and we implement the Lax-Friedrichs scheme to get numerical solutions.


title: A Privacy-Preserving Distributed Control of Optimal Power Flow

abstract: We consider a distributed optimal power flow formulated as an optimization problem that maximizes a nondifferentiable concave function. Solving such a problem by the existing distributed algorithms can lead to data privacy issues because the solution information exchanged within the algorithms can be utilized by an adversary to infer the data. To preserve data privacy, in this paper we propose a differentially private projected subgradient (DP-PS) algorithm that includes a solution encryption step. We show that a sequence generated by DP-PS converges in expectation, in probability, and with probability 1. Moreover, we show that the rate of convergence in expectation is affected by a target privacy level of DP-PS chosen by the user. We conduct numerical experiments that demonstrate the convergence and data privacy preservation of DP-PS.


title: Sub-grid scale characterization and asymptotic behavior of multi-dimensional upwind schemes for the vorticity transport equations

abstract: We study the sub-grid scale characteristics of a vorticity-transport-based approach for large-eddy simulations. In particular, we consider a multi-dimensional upwind scheme for the vorticity transport equations and establish its properties in the under-resolved regime. The asymptotic behavior of key turbulence statistics of velocity gradients, vorticity, and invariants is studied in detail. Modified equation analysis indicates that dissipation can be controlled locally via non-linear limiting of the gradient employed for the vorticity reconstruction on the cell face such that low numerical diffusion is obtained in well-resolved regimes and high numerical diffusion is realized in under-resolved regions. The enstrophy budget highlights the remarkable ability of the truncation terms to mimic the true sub-grid scale dissipation and diffusion. The modified equation also reveals diffusive terms that are similar to several commonly employed sub-grid scale models including tensor-gradient and hyper-viscosity models. Investigations on several canonical turbulence flow cases show that large-scale features are adequately represented and remain consistent in terms of spectral energy over a range of grid resolutions. Numerical dissipation in under-resolved simulations is consistent and can be characterized by diffusion terms discovered in the modified equation analysis. A minimum state of scale separation necessary to obtain asymptotic behavior is characterized using metrics such as effective Reynolds number and effective grid spacing. Temporally-evolving jet simulations, characterized by large-scale vortical structures, demonstrate that high Reynolds number vortex-dominated flows are captured when criteria is met and necessitate diffusive non-linear limiting of vorticity reconstruction be employed to realize accuracy in under-resolved simulations.


title: Learning to Isolate Muons

abstract: Distinguishing between prompt muons produced in heavy boson decay and muons produced in association with heavy-flavor jet production is an important task in analysis of collider physics data. We explore whether there is information available in calorimeter deposits that is not captured by the standard approach of isolation cones. We find that convolutional networks and particle-flow networks accessing the calorimeter cells surpass the performance of isolation cones, suggesting that the radial energy distribution and the angular structure of the calorimeter deposits surrounding the muon contain unused discrimination power. We assemble a small set of high-level observables which summarize the calorimeter information and close the performance gap with networks which analyze the calorimeter cells directly. These observables are theoretically well-defined and can be studied with collider data.


title: Insiders and Outsiders in Research on Machine Learning and Society

abstract: A subset of machine learning research intersects with societal issues, including fairness, accountability and transparency, as well as the use of machine learning for social good. In this work, we analyze the scholars contributing to this research at the intersection of machine learning and society through the lens of the sociology of science. By analyzing the authorship of all machine learning papers posted to arXiv, we show that compared to researchers from overrepresented backgrounds (defined by gender and race/ethnicity), researchers from underrepresented backgrounds are more likely to conduct research at this intersection than other kinds of machine learning research. This state of affairs leads to contention between two perspectives on insiders and outsiders in the scientific enterprise: outsiders being those outside the group being studied, and outsiders being those who have not participated as researchers in an area historically. This contention manifests as an epistemic question on the validity of knowledge derived from lived experience in machine learning research, and predicts boundary work that we see in a real-world example.


title: Prime zeta function statistics and Riemann zero-difference repulsion

abstract: We present a derivation of the numerical phenomenon that differences between the Riemann zeta function's nontrivial zeros tend to avoid being equal to the imaginary parts of the zeros themselves, a property called statistical "repulsion" between the zeros and their differences. Our derivation relies on the statistical properties of the prime zeta function, whose singularity structure specifies the positions of the Riemann zeros. We show that the prime zeta function on the critical line is asymptotically normally distributed with a covariance function that is closely approximated by the logarithm of the Riemann zeta function's magnitude on the 1-line. This creates notable negative covariance at separations approximately equal to the imaginary parts of the Riemann zeros. This covariance function and the singularity structure of the prime zeta function combine to create a conditional statistical bias at the locations of the Riemann zeros that predicts the zero-difference repulsion effect. Our method readily generalizes to describe similar effects in the zeros of related L-functions.

categories: math.NT cond-mat.stat-mech math-ph math.MP math.PR


title: Hypergraph States in SU(N)1, N odd prime, Chern-Simons Theory

abstract: Graph states and hypergraph states can be constructed from products of basic operations that appear in SU(N)1. The level-rank dual of a theorem of Salton, Swingle, and Walter implies that these operations can be prepared topologically in the n-torus Hilbert space of Chern-Simons theory for N neq 5 mod 4. For SU(N)1, N = 5 mod 4, only stabilizer states can be prepared on the n-torus Hilbert space, which restricts the construction to graph states.


title: Downbeat Tracking with Tempo-Invariant Convolutional Neural Networks

abstract: The human ability to track musical downbeats is robust to changes in tempo, and it extends to tempi never previously encountered. We propose a deterministic time-warping operation that enables this skill in a convolutional neural network (CNN) by allowing the network to learn rhythmic patterns independently of tempo. Unlike conventional deep learning approaches, which learn rhythmic patterns at the tempi present in the training dataset, the patterns learned in our model are tempo-invariant, leading to better tempo generalisation and more efficient usage of the network capacity. We test the generalisation property on a synthetic dataset created by rendering the Groove MIDI Dataset using FluidSynth, split into a training set containing the original performances and a test set containing tempo-scaled versions rendered with different SoundFonts (test-time augmentation). The proposed model generalises nearly perfectly to unseen tempi (F-measure of 0.89 on both training and test sets), whereas a comparable conventional CNN achieves similar accuracy only for the training set (0.89) and drops to 0.54 on the test set. The generalisation advantage of the proposed model extends to real music, as shown by results on the GTZAN and Ballroom datasets.


title: Measurement of the W$\gamma$ production cross section in proton-proton collisions at $\sqrt{s} =$ 13 TeV and constraints on effective field theory coefficients

abstract: A fiducial cross section for W$\gamma$ production in proton-proton collisions is measured at a center-of-mass energy of 13 TeV in 137 fb$^{-1}$ of data collected using the CMS detector at the LHC. The W $\to$ e$\nu$ and $\mu\nu$ decay modes are used in a maximum-likelihood fit to the lepton-photon invariant mass distribution to extract the combined cross section. The measured cross section is compared with theoretical expectations at next-to-leading order in quantum chromodynamics. In addition, 95% confidence level intervals are reported for anomalous triple-gauge couplings within the framework of effective field theory.


title: Consequences of constant elevation scans for instrumental systematics in Cosmic Microwave Background Experiments

abstract: Instrumental systematics need to be controlled to high precision for upcoming Cosmic Microwave Background (CMB) experiments. The level of contamination caused by these systematics is often linked to the scan strategy, and scan strategies for satellite experiments can significantly mitigate these systematics. However, no detailed study has been performed for ground-based experiments. Here we show that under the assumption of constant elevation scans (CESs), the ability of the scan strategy to mitigate these systematics is strongly limited, irrespective of the detailed structure of the scan strategy. We calculate typical values and maps of the quantities coupling the scan to the systematics, and show how these quantities vary with the choice of observing elevations. These values and maps can be used to calculate and forecast the magnitude of different instrumental systematics without requiring detailed scan strategy simulations. As a reference point, we show that inclusion of even a single boresight rotation angle significantly improves over sky rotation alone for mitigating these systematics. A standard metric for evaluating cross-linking is related to one of the parameters studied in this work, so a corollary of our work is that the cross-linking will suffer from the same CES limitations and therefore upcoming CMB surveys will unavoidably have poorly cross-linked regions if they use CESs, regardless of detailed scheduling choices. Our results are also relevant for non-CMB surveys that perform constant elevation scans and may have scan-coupled systematics, such as intensity mapping surveys.


title: Covariate adjustment in subgroup analyses of randomized clinical trials: A propensity score approach

abstract: Background: Subgroup analyses are frequently conducted in randomized clinical trials to assess evidence of heterogeneous treatment effect across patient subpopulations. Although randomization balances covariates within subgroups in expectation, chance imbalance may be amplified in small subgroups and harm the precision. Covariate adjustment in overall analysis of RCT is often conducted, via either ANCOVA or propensity score weighting, but for subgroup analysis has been rarely discussed. In this article, we develop propensity score weighting methodology for covariate adjustment to improve the precision and power of subgroup analyses in RCTs. Methods: We extend the propensity score weighting methodology to subgroup analyses by fitting a logistic propensity model with pre-specified covariate-subgroup interactions. We show that, by construction, overlap weighting exactly balances the covariates with interaction terms in each subgroup. Extensive simulations were performed to compare the operating characteristics of unadjusted, different propensity score weighting and the ANCOVA estimator. We apply these methods to the HF-ACTION trial to evaluate the effect of exercise training on 6-minute walk test in pre-specified subgroups. Results: Standard errors of the adjusted estimators are smaller than those of the unadjusted estimator. The propensity score weighting estimator is as efficient as ANCOVA, and is often more efficient when subgroup sample size is small (e.g.<125), and/or when outcome model is misspecified. The weighting estimators with full-interaction propensity model consistently outperform the standard main-effect propensity model. Conclusion: Propensity score weighting is a transparent and objective method to adjust chance imbalance of important covariates in subgroup analyses of RCTs. It is crucial to include the full covariate-subgroup interactions in the propensity score model.


title: HiCOPS: High Performance Computing Framework for Tera-Scale Database Search of Mass Spectrometry based Omics Data

abstract: Database-search algorithms, that deduce peptides from Mass Spectrometry (MS) data, have tried to improve the computational efficiency to accomplish larger, and more complex systems biology studies. Existing serial, and high-performance computing (HPC) search engines, otherwise highly successful, are known to exhibit poor-scalability with increasing size of theoretical search-space needed for increased complexity of modern non-model, multi-species MS-based omics analysis. Consequently, the bottleneck for computational techniques is the communication costs of moving the data between hierarchy of memory, or processing units, and not the arithmetic operations. This post-Moore change in architecture, and demands of modern systems biology experiments have dampened the overall effectiveness of the existing HPC workflows. We present a novel efficient parallel computational method, and its implementation on memory-distributed architectures for peptide identification tool called HiCOPS, that enables more than 100-fold improvement in speed over most existing HPC proteome database search tools. HiCOPS empowers the supercomputing database search concept for comprehensive identification of peptides, and all their modified forms within a reasonable time-frame. We demonstrate this by searching Gigabytes of experimental MS data against Terabytes of databases where HiCOPS completes peptide identification in few minutes using 72 parallel nodes (1728 cores) compared to several weeks required by existing state-of-the-art tools using 1 node (24 cores); 100 minutes vs 5 weeks; 500x speedup. Finally, we formulate a theoretical framework for our overhead-avoiding strategy, and report superior performance evaluation results for key metrics including execution time, CPU utilization, speedups, and I/O efficiency. The software will be made available at: hicops.github.io


title: Echo-SyncNet: Self-supervised Cardiac View Synchronization in Echocardiography

abstract: In echocardiography (echo), an electrocardiogram (ECG) is conventionally used to temporally align different cardiac views for assessing critical measurements. However, in emergencies or point-of-care situations, acquiring an ECG is often not an option, hence motivating the need for alternative temporal synchronization methods. Here, we propose Echo-SyncNet, a self-supervised learning framework to synchronize various cross-sectional 2D echo series without any external input. The proposed framework takes advantage of both intra-view and inter-view self supervisions. The former relies on spatiotemporal patterns found between the frames of a single echo cine and the latter on the interdependencies between multiple cines. The combined supervisions are used to learn a feature-rich embedding space where multiple echo cines can be temporally synchronized. We evaluate the framework with multiple experiments: 1) Using data from 998 patients, Echo-SyncNet shows promising results for synchronizing Apical 2 chamber and Apical 4 chamber cardiac views; 2) Using data from 3070 patients, our experiments reveal that the learned representations of Echo-SyncNet outperform a supervised deep learning method that is optimized for automatic detection of fine-grained cardiac phase; 3) We show the usefulness of the learned representations in a one-shot learning scenario of cardiac keyframe detection. Without any fine-tuning, keyframes in 1188 validation patient studies are identified by synchronizing them with only one labeled reference study. We do not make any prior assumption about what specific cardiac views are used for training and show that Echo-SyncNet can accurately generalize to views not present in its training set. Project repository: github.com/fatemehtd/Echo-SyncNet.


title: $^{229m}$Th isomer from a nuclear model perspective

abstract: The physical conditions for the emergence of the extremely low-lying nuclear isomer $^{229m}$Th at approximately 8 eV are investigated in the framework of our recently proposed nuclear structure model. Our theoretical approach explains the $^{229m}$Th-isomer phenomenon as the result of a very fine interplay between collective quadrupole-octupole and single-particle dynamics in the nucleus. We find that the isomeric state can only appear in a rather limited model space of quadrupole-octupole deformations in the single-particle potential, with the octupole deformation being of a crucial importance for its formation. Within this deformation space the model-described quantities exhibit a rather smooth behaviour close to the line of isomer-ground state quasi-degeneracy determined by the crossing of the corresponding single-particle orbitals. Our comprehensive analysis confirms the previous model predictions for reduced transition probabilities and the isomer magnetic moment, while showing a possibility for limited variation in the ground-state magnetic moment theoretical value. These findings prove the reliability of the model and suggest that the same dynamical mechanism could manifest in other actinide nuclei giving a general prescription for the search and exploration of similar isomer phenomena.


title: Equilibration and Typicality in Quantum Processes

abstract: Forgetfulness is a common feature of nature. Moreover, without forgetfulness, repeatability would be impossible. Despite this, small systems constantly leak information about their state to their surroundings, and quantum mechanics tells us that this information cannot be deleted, invariably returning to influence their future behaviour. How can physical nature be forgetful if it is forbidden to forget? The results within this thesis bridge this gap between what we see in the real world and what idealised physical theories say, explaining the emergence of forgetful processes from closed quantum dynamics and allowing to quantify the rate at which memory effects become relevant.


title: Discovery and exclusion prospects for staus produced by heavy Higgs bosons decays at the LHC

abstract: In a previous work we developed a search strategy for staus produced by the decay of the heavy CP-even Higgs boson $H$ within the context of the large $\tan\beta$ regime of the minimal supersymmetric standard model (MSSM) in an scenario of large stau mixing. Here we study the performance of such search strategy by confronting it with the complementary mixing pattern in which decays of both the CP-even and CP-odd heavy Higgs bosons contribute to the production of $\widetilde{\tau}_1\widetilde{\tau}_2^{*} \;+\; c.c$ pairs. Again, we focus on final states with two opposite-sign tau leptons and large missing transverse energy. We find that our proposed search strategy, although optimized for the large stau mixing scenario, is still quite sensitive to the complementary mixing pattern. For instance, with a total integrated luminosity of only 100 fb$^{-1}$ we are able to exclude heavy Higgs masses above 850 GeV for average stau masses higher than 290 GeV. We also extend the results reported in the preceding work for the large mixing scenario by including now the exclusion limits at 100 fb$^{-1}$ and the prospects both for exclusion and discovery in a potential high luminosity phase of the LHC (1000 fb$^{-1}$). Finally, we discuss the possibility to distinguish the two mixing scenarios when they share the same relevant mass spectrum and both reach the discovery level with our search strategy.


title: Nearest Neighbor-based Importance Weighting

abstract: Importance weighting is widely applicable in machine learning in general and in techniques dealing with data covariate shift problems in particular. A novel, direct approach to determine such importance weighting is presented. It relies on a nearest neighbor classification scheme and is relatively straightforward to implement. Comparative experiments on various classification tasks demonstrate the effectiveness of our so-called nearest neighbor weighting (NNeW) scheme. Considering its performance, our procedure can act as a simple and effective baseline method for importance weighting.


title: Predicting the probability distribution of bus travel time to move towards reliable planning of public transport services

abstract: An important aspect of the quality of a public transport service is its reliability, which is defined as the invariability of the service attributes. Preventive measures taken during planning can reduce risks of unreliability throughout operations. In order to tackle reliability during the service planning phase, a key piece of information is the long-term prediction of the density of the travel time, which conveys the uncertainty of travel times. We introduce a reliable approach to one of the problems of service planning in public transport, namely the Multiple Depot Vehicle Scheduling Problem (MDVSP), which takes as input a set of trips and the probability density function (p.d.f.) of the travel time of each trip in order to output delay-tolerant vehicle schedules. This work empirically compares probabilistic models for the prediction of the conditional p.d.f. of the travel time, as a first step towards reliable MDVSP solutions. Two types of probabilistic models, namely similarity-based density estimation models and a smoothed Logistic Regression for probabilistic classification model, are compared on a dataset of more than 41,000 trips and 50 bus routes of the city of Montr\'eal. The result of a vast majority of probabilistic models outperforms that of a Random Forests model, which is not inherently probabilistic, thus highlighting the added value of modeling the conditional p.d.f. of the travel time with probabilistic models. A similarity-based density estimation model using a $k$ Nearest Neighbors method and a Kernel Density Estimation predicted the best estimate of the true conditional p.d.f. on this dataset.


title: Enhanced Convergence of Quantum Typicality using a Randomized Low-Rank Approximation

abstract: We present a method to reduce the variance of stochastic trace estimators used in quantum typicality (QT) methods via a randomized low-rank approximation of the finite-temperature density matrix $e^{-\beta H}$. The trace can be evaluated with higher accuracy in the low-rank subspace while using the QT estimator to approximate the trace in the complementary subspace. We present two variants of the trace estimator and demonstrate their efficacy using numerical experiments. The experiments show that the low-rank approximation outperforms the standard QT trace estimator for moderate- to low-temperature. We argue this is due to the low-rank approximation accurately represent the density matrix at low temperatures, allowing for accurate results for the trace.

categories: quant-ph cond-mat.quant-gas cond-mat.stat-mech cond-mat.str-el physics.comp-ph


title: COVID-19 spreading under containment actions

abstract: We propose an epidemiological model that includes the mobility patterns of the individuals, in the spirit to those considered in (Barmak, 2011, 2016) and (Medus, 2011). We assume that people move around in a city of 120x120 blocks with 300 inhabitants in each block. The mobility pattern is associated to a complex network in which nodes represent blocks while the links represent the traveling path of the individuals. We implemented three confinement strategies in order to mitigate the disease spreading: 1) global confinement, 2) partial restriction to mobility, and 3) localized confinement. In the first case, it was observed that a global isolation policy prevents the massive outbreak of the disease. In the second case, a partial restriction to mobility could lead to a massive contagion if this was not complemented with sanitary measures such as the use of masks and social distancing. Finally, a local isolation policy was proposed, conditioned to the health status of each block. It was observed that this mitigation strategy was able to contain and even reduce the outbreak of the disease by intervening in specific regions of the city according to their level of contagion. It was also observed that this strategy is capable of controlling the epidemic in the case that a certain proportion of those infected are asymptomatic.


title: Variational Bayes survival analysis for unemployment modelling

abstract: Mathematical modelling of unemployment dynamics attempts to predict the probability of a job seeker finding a job as a function of time. This is typically achieved by using information in unemployment records. These records are right censored, making survival analysis a suitable approach for parameter estimation. The proposed model uses a deep artificial neural network (ANN) as a non-linear hazard function. Through embedding, high-cardinality categorical features are analysed efficiently. The posterior distribution of the ANN parameters are estimated using a variational Bayes method. The model is evaluated on a time-to-employment data set spanning from 2011 to 2020 provided by the Slovenian public employment service. It is used to determine the employment probability over time for each individual on the record. Similar models could be applied to other questions with multi-dimensional, high-cardinality categorical data including censored records. Such data is often encountered in personal records, for example in medical records.


title: The fractional nonlinear PT dimer

abstract: We examine a fractional Discrete Nonlinear Schrodinger dimer, where the usual first-order derivative of the time evolution is replaced by a non integer-order derivative. The dimer is nonlinear (Kerr) and PT -symmetric, and we examine the exchange dynamics between both sites. By means of the Laplace transformation technique, the linear PT dimer is solved in closed form in terms of Mittag-Leffler functions, while for the nonlinear regime, we resort to numerical computations using the direct explicit Grunwald algorithm. In general, the main effect of the fractional derivative is the onset of a monotonically decreasing time envelope for the amplitude of the oscillatory exchange. In the presence of PT symmetry, the dynamics shows damped oscillations for small gain/loss in both sites, while at higher gain/loss parameter values, the amplitudes of both sites grows unbounded. In the presence of nonlinearity, selftrapping is still possible although the trapped fraction decreases as the nonlinearity is increased past threshold, in marked contrast with the standard case.


title: pcoxtime: Penalized Cox Proportional Hazard Model for Time-dependent Covariates

abstract: The penalized Cox proportional hazard model is a popular analytical approach for survival data with a large number of covariates. Such problems are especially challenging when covariates vary over follow-up time (i.e., the covariates are time-dependent). The standard R packages for fully penalized Cox models cannot currently incorporate time-dependent covariates. To address this gap, we implement a variant of gradient descent algorithm (proximal gradient descent) for fitting penalized Cox models. We apply our implementation to real and simulated data sets.


title: Super-replication with transaction costs under model uncertainty for continuous processes

abstract: We formulate a superhedging theorem in the presence of transaction costs and model uncertainty. Asset prices are assumed continuous and uncertainty is modelled in a parametric setting. Our proof relies on a new topological framework in which no Krein-Smulian type theorem is available.


title: Asymptotically linear iterated function systems on the real line

abstract: Given a sequence of i.i.d. random functions $\Psi_{n}:\mathbb{R}\to\mathbb{R}$, $n\in\mathbb{N}$, we consider the iterated function system and Markov chain which is recursively defined by $X_{0}^{x}:=x$ and $X_{n}^{x}:=\Psi_{n-1}(X_{n-1}^{x})$ for $x\in\mathbb{R}$ and $n\in\mathbb{N}$. Under the two basic assumptions that the $\Psi_{n}$ are a.s. continuous at any point in $\mathbb{R}$ and asymptotically linear at the "endpoints" $\pm\infty$, we study the tail behavior of the stationary laws of such Markov chains by means of Markov renewal theory. Our approach provides an extension of Goldie's implicit renewal theory and can also be viewed as an adaptation of Kesten's work on products of random matrices to one-dimensional function systems as described. Our results have applications in quite different areas of applied probability like queuing theory, econometrics, mathematical finance and population dynamics. Our results have applications in quite different areas of applied probability like queuing theory, econometrics, mathematical finance and population dynamics, e.g. ARCH models and random logistic transforms.


title: A Massive, Clumpy Molecular Gas Distribution and Displaced AGN in Zw 3146

abstract: We present a recent ALMA observation of the CO(1-0) line emission in the central galaxy of the Zw 3146 galaxy cluster ($z=0.2906$). We also present updated X-ray cavity measurements from archival Chandra observations. The $5\times 10^{10}\,M_{\odot}$ supply of molecular gas, which is confined to the central 4 kpc, is marginally resolved into three extensions that are reminiscent of the filaments observed in similar systems. No velocity structure that would be indicative of ordered motion is observed. The three molecular extensions all trail X-ray cavities, and are potentially formed from the condensation of intracluster gas lifted in the wakes of the rising bubbles. Many cycles of feedback would be require to account for the entire molecular gas reservoir. The molecular gas and continuum source are mutually offset by 2.6 kpc, with no detected line emission coincident with the continuum source. It is the molecular gas, not the continuum source, that lies at the gravitational center of the brightest cluster galaxy. As the brightest cluster galaxy contains possible tidal features, the displaced continuum source may correspond to the nucleus of a merging galaxy. We also discuss the possibility that a gravitational wave recoil following a black hole merger may account for the displacement.


title: Parallax estimation for push-frame satellite imagery: application to super-resolution and 3D surface modeling from Skysat products

abstract: Recent constellations of satellites, including the Skysat constellation, are able to acquire bursts of images. This new acquisition mode allows for modern image restoration techniques, including multi-frame super-resolution. As the satellite moves during the acquisition of the burst, elevation changes in the scene translate into noticeable parallax. This parallax hinders the results of the restoration. To cope with this issue, we propose a novel parallax estimation method. The method is composed of a linear Plane+Parallax decomposition of the apparent motion and a multi-frame optical flow algorithm that exploits all frames simultaneously. Using SkySat L1A images, we show that the estimated per-pixel displacements are important for applying multi-frame super-resolution on scenes containing elevation changes and that can also be used to estimate a coarse 3D surface model.


title: Cleora: A Simple, Strong and Scalable Graph Embedding Scheme

abstract: The area of graph embeddings is currently dominated by contrastive learning methods, which demand formulation of an explicit objective function and sampling of positive and negative examples. This creates a conceptual and computational overhead. Simple, classic unsupervised approaches like Multidimensional Scaling (MSD) or the Laplacian eigenmap skip the necessity of tedious objective optimization, directly exploiting data geometry. Unfortunately, their reliance on very costly operations such as matrix eigendecomposition make them unable to scale to large graphs that are common in today's digital world. In this paper we present Cleora: an algorithm which gets the best of two worlds, being both unsupervised and highly scalable. We show that high quality embeddings can be produced without the popular step-wise learning framework with example sampling. An intuitive learning objective of our algorithm is that a node should be similar to its neighbors, without explicitly pushing disconnected nodes apart. The objective is achieved by iterative weighted averaging of node neigbors' embeddings, followed by normalization across dimensions. Thanks to the averaging operation the algorithm makes rapid strides across the embedding space and usually reaches optimal embeddings in just a few iterations. Cleora runs faster than other state-of-the-art CPU algorithms and produces embeddings of competitive quality as measured on downstream tasks: link prediction and node classification. We show that Cleora learns a data abstraction that is similar to contrastive methods, yet at much lower computational cost. We open-source Cleora under the MIT license allowing commercial use under https://github.com/Synerise/cleora.


title: Jetting during oblique impacts of spherical impactors

abstract: During the early stages of an impact a small amount material may be jetted and ejected at speeds exceeding the impact velocity. Jetting is an important process for producing melt during relatively low velocity impacts. How impact angle affects the jetting process has yet to be fully understood. Here, we simulate jetting during oblique impacts using the iSALE shock physics code. Assuming both the target and impactor have the same composition (dunite), we examine the jetted material which exceeds the impact velocity. Our results show that oblique impacts always produce more jetted ejecta than vertical impacts, except for grazing impacts with impact angles $< 15^{\circ}$. A 45$^{\circ}$ impact with an impact velocity of 3 km/s produces jetted material equal to $\sim$ 7 \% of the impactor mass. This is 6 times the jetted mass produced by a vertical impact with similar impact conditions. We also find that the origin of jetted ejecta depends on impact angle; for impact angles less than 45$^{\circ}$, most of the jet is composed of impactor material, while at higher impact angles the jet is dominated by target material. Our findings are consistent with previous experimental work. In all cases, jetted materials are preferentially distributed downrange of the impactor.


title: Improved Cooperation by Exploiting a Common Signal

abstract: Can artificial agents benefit from human conventions? Human societies manage to successfully self-organize and resolve the tragedy of the commons in common-pool resources, in spite of the bleak prediction of non-cooperative game theory. On top of that, real-world problems are inherently large-scale and of low observability. One key concept that facilitates human coordination in such settings is the use of conventions. Inspired by human behavior, we investigate the learning dynamics and emergence of temporal conventions, focusing on common-pool resources. Extra emphasis was given in designing a realistic evaluation setting: (a) environment dynamics are modeled on real-world fisheries, (b) we assume decentralized learning, where agents can observe only their own history, and (c) we run large-scale simulations (up to 64 agents). Uncoupled policies and low observability make cooperation hard to achieve; as the number of agents grow, the probability of taking a correct gradient direction decreases exponentially. By introducing an arbitrary common signal (e.g., date, time, or any periodic set of numbers) as a means to couple the learning process, we show that temporal conventions can emerge and agents reach sustainable harvesting strategies. The introduction of the signal consistently improves the social welfare (by 258% on average, up to 3306%), the range of environmental parameters where sustainability can be achieved (by 46% on average, up to 300%), and the convergence speed in low abundance settings (by 13% on average, up to 53%).


title: Correlation of the sunspot number and the waiting time distribution of solar flares, coronal mass ejections, and solar wind switchback events observed with the Parker Solar Probe

abstract: Waiting time distributions of solar flares and {\sl coronal mass ejections (CMEs)} exhibit power law-like distribution functions with slopes in the range of $\alpha_{\tau} \approx 1.4-3.2$, as observed in annual data sets during 4 solar cycles (1974-2012). We find a close correlation between the waiting time power law slope $\alpha_\tau$ and the {\sl sunspot number (SN)}, i.e., $\alpha_\tau$ = 1.38 + 0.01 $\times$ SN. The waiting time distribution can be fitted with a Pareto-type function of the form $N(\tau) = N_0$ $(\tau_0 + \tau)^{-\alpha_{\tau}}$, where the offset $\tau_0$ depends on the instrumental sensitivity, the detection threshold of events, and pulse pile-up effects. The time-dependent power law slope $\alpha_{\tau}(t)$ of waiting time distributions depends only on the global solar magnetic flux (quantified by the sunspot number) or flaring rate, independent of other physical parameters of {\sl self-organized criticality (SOC)} or {\sl magneto-hydrodynamic (MHD)} turbulence models. Power law slopes of $\alpha_{\tau}\approx 1.2-1.6$ were also found in solar wind switchback events, as observed with the {\sl Parker Solar Probe (PSP)}. We conclude that the annual variability of switchback events in the heliospheric solar wind is modulated by flare and CME rates originating in the photosphere and lower corona.


title: Uniform Distribution of Sequences and its interplay with Functional Analysis

abstract: In this paper we apply ideas from the theory of Uniform Distribution of sequences to Functional Analysis and then drawing inspiration from the consequent results, we study concepts and results in Uniform Distribution itself. So let $E$ be a Banach space. Then we prove:\\ (a) If $F$ is a bounded subset of $E$ and $x \in \overline{\co}(F)$ (= the closed convex hull of $F$), then there is a sequence $(x_n) \subseteq F$ which is Ces\`{a}ro summable to $x$.\\ (b) If $E$ is separable, $F \subseteq E^*$ bounded and $f \in \overline{\co}^{w^*}(F)$, then there is a sequence $(f_n) \subseteq F$ whose sequence of arithmetic means $\frac{f_1+\dots+f_N}{N}$, $N \ge 1$ weak$^*$-converges to $f$. By the aid of the Krein-Milman theorem, both (a) and (b) have interesting implications for closed, convex and bounded subsets $\Omega$ of $E$ such that $\Omega=\overline{\co}(\ex \Omega)$ and for weak$^*$ compact and convex subsets of $E^*$. Of particular interest is the case when $\Omega=B_{C(K)^*}$, where $K$ is a compact metric space. By further expanding the previous ideas and results, we are able to generalize a classical theorem of Uniform Distribution which is valid for increasing functions $\varphi:I=[0,1] \rightarrow \mathbb{R}$ with $\varphi(0)=0$ and $\varphi(1)=1$, for functions $\varphi$ of bounded variation on $I$ with $\varphi(0)=0$ and total variation $V_0^1 \varphi=1$.


title: Typing Errors in Factual Knowledge Graphs: Severity and Possible Ways Out

abstract: Factual knowledge graphs (KGs) such as DBpedia and Wikidata have served as part of various downstream tasks and are also widely adopted by artificial intelligence research communities as benchmark datasets. However, we found these KGs to be surprisingly noisy. In this study, we question the quality of these KGs, where the typing error rate is estimated to be 27% for coarse-grained types on average, and even 73% for certain fine-grained types. In pursuit of solutions, we propose an active typing error detection algorithm that maximizes the utilization of both gold and noisy labels. We also comprehensively discuss and compare unsupervised, semi-supervised, and supervised paradigms to deal with typing errors in factual KGs. The outcomes of this study provide guidelines for researchers to use noisy factual KGs. To help practitioners deploy the techniques and conduct further research, we published our code and data.


title: Fuzzing Hardware Like Software

abstract: Hardware flaws are permanent and potent: hardware cannot be patched once fabricated, and any flaws may undermine any software executing on top. Consequently, verification time dominates implementation time. The gold standard in hardware Design Verification (DV) is concentrated at two extremes: random dynamic verification and formal verification. Both struggle to root out the subtle flaws in complex hardware that often manifest as security vulnerabilities. The root problem with random verification is its undirected nature, making it inefficient, while formal verification is constrained by the state-space explosion problem, making it infeasible against complex designs. What is needed is a solution that is directed, yet under-constrained. Instead of making incremental improvements to existing DV approaches, we leverage the observation that existing software fuzzers already provide such a solution, and adapt them for hardware DV. Specifically, we translate RTL hardware to a software model and fuzz that model. The central challenge we address is how best to mitigate the differences between the hardware execution model and software execution model. This includes: 1) how to represent test cases, 2) what is the hardware equivalent of a crash, 3) what is an appropriate coverage metric, and 4) how to create a general-purpose fuzzing harness for hardware. To evaluate our approach, we fuzz four IP blocks from Google's OpenTitan SoC. Our experiments reveal a two orders-of-magnitude reduction in run time to achieve Finite State Machine (FSM) coverage over traditional dynamic verification schemes. Moreover, with our design-agnostic harness, we achieve over 88% HDL line coverage in three out of four of our designs -- even without any initial seeds.


title: Sensitivity to New Physics of Isotope Shift Studies using the Coronal Lines of Highly Charged Calcium Ions

abstract: Promising searches for new physics beyond the current Standard Model (SM) of particle physics are feasible through isotope-shift spectroscopy, which is sensitive to a hypothetical fifth force between the neutrons of the nucleus and the electrons of the shell. Such an interaction would be mediated by a new particle which could in principle be associated with dark matter. In so-called King plots, the mass-scaled frequency shifts of two optical transitions are plotted against each other for a series of isotopes. Subtle deviations from the expected linearity could reveal such a fifth force. Here, we study experimentally and theoretically six transitions in highly charged ions of Ca, an element with five stable isotopes of zero nuclear spin. Some of the transitions are suitable for upcoming high-precision coherent laser spectroscopy and optical clocks. Our results provide a sufficient number of clock transitions for -- in combination with those of singly charged Ca$^+$ -- application of the generalized King plot method. This will allow future high-precision measurements to remove higher-order SM-related nonlinearities and open a new door to yet more sensitive searches for unknown forces and particles.


title: Epoch of reionization parameter estimation with the 21-cm bispectrum

abstract: We present the first application of the isosceles bispectrum to MCMC parameter inference from the cosmic 21-cm signal. We extend the MCMC sampler 21cmMC to use the fast bispectrum code, BiFFT, when computing the likelihood. We create mock 1000h observations with SKA1-low, using PyObs21 to account for uv-sampling and thermal noise. Assuming the spin temperature is much higher than that of the CMB, we consider two different reionization histories for our mock observations: fiducial and late-reionization. For both models we find that bias on the inferred parameter means and 1-$\sigma$ credible intervals can be substantially reduced by using the isosceles bispectrum (calculated for a wide range of scales and triangle shapes) together with the power spectrum (as opposed to just using one of the statistics). We find that making the simplifying assumption of a Gaussian likelihood with a diagonal covariance matrix does not notably bias parameter constraints for the three-parameter reionization model and basic instrumental effects considered here. This is true even if we use extreme (unlikely) initial conditions which would be expected to amplify biases. We also find that using the cosmic variance error calculated with Monte-Carlo simulations using the fiducial model parameters whilst assuming the late-reionization model for the simulated data also does not strongly bias the inference. This implies we may be able to sparsely sample and interpolate the cosmic variance error over the parameter space, substantially reducing computational costs. All codes used in this work are publicly-available.