title: Topological Entropy for Arbitrary Subsets of Infinite Product Spaces

abstract: In this note a notion of generalized topological entropy for arbitrary subsets of the space of all sequences in a compact topological space is introduced. It is shown that for a continuous map on a compact space the generalized topological entropy of the set of all orbits of the map coincides with the classical topological entropy of the map. Some basic properties of this new notion of entropy are considered; among them are: the behavior of the entropy with respect to disjoint union, cartesian product, component restriction and dilation, shift mapping, and some continuity properties with respect to Vietoris topology. As an example, it is shown that any self-similar structure of a fractal given by a finite family of contractions gives rise to a notion of intrinsic topological entropy for subsets of the fractal. A generalized notion of Bowen's entropy associated to any increasing sequence of compatible semimetrics on a topological space is introduced and some of its basic properties are considered. As a special case for $1\leq p\leq\infty$ the Bowen $p$-entropy of sets of sequences of any metric space is introduced. It is shown that the notions of generalized topological entropy and Bowen $\infty$-entropy for compact metric spaces coincide.


title: GP-ETAS: Semiparametric Bayesian inference for the spatio-temporal Epidemic Type Aftershock Sequence model

abstract: The spatio-temporal Epidemic Type Aftershock Sequence (ETAS) model is widely used to describe the self-exciting nature of earthquake occurrences. While traditional inference methods provide only point estimates of the model parameters, we aim at a full Bayesian treatment of model inference, allowing naturally to incorporate prior knowledge and uncertainty quantification of the resulting estimates. Therefore, we introduce a highly flexible, non-parametric representation for the spatially varying ETAS background intensity through a Gaussian process (GP) prior. Combined with classical triggering functions this results in a new model formulation, namely the GP-ETAS model. We enable tractable and efficient Gibbs sampling by deriving an augmented form of the GP-ETAS inference problem. This novel sampling approach allows us to assess the posterior model variables conditioned on observed earthquake catalogues, i.e., the spatial background intensity and the parameters of the triggering function. Empirical results on two synthetic data sets indicate that GP-ETAS outperforms standard models and thus demonstrate the predictive power for observed earthquake catalogues including uncertainty quantification for the estimated parameters. Finally, a case study for the l'Aquila region, Italy, with the devastating event on 6 April 2009, is presented.


title: New multi-part collisional model of the main belt: The contribution to near-Earth asteroids

abstract: Aims. We developed a six-part collisional evolution model of the main asteroid belt (MB) and used it to study the contribution of the different regions of the MB to the near-earth asteroids (NEAs). Methods. We built a statistical code called ACDC that simulates the collisional evolution of the MB split into six regions (namely Inner, Middle, Pristine, Outer, Cybele and High-Inclination belts) according to the positions of the major resonances present there ($\nu_{6}$, 3:1J, 5:2J, 7:3J and 2:1J). We consider the Yarkovsky effect and the mentioned resonances as the main mechanism that removes asteroids from the different regions of the MB and delivers them to the NEA region. We calculated the evolution of the NEAs coming from the different source regions by considering the bodies delivered by the resonances and mean dynamical timescales in the NEA population. Results. Our model is in agreement with the major observational constraints associated with the MB, such as the size distributions of the different regions of the MB and the number of large asteroid families. It is also able to reproduce the observed NEAs with H < 16 and agrees with recent estimations for H < 20, but deviates for smaller sizes. We find that most sources make a significant contribution to the NEAs; however the Inner and Middle belts stand out as the most important source of NEAs followed by the Outer belt. The contributions of the Pristine and Cybele regions are minor. The High-Inclination belt is the source of only a fraction of the actual observed NEAs with high inclination, as there are dynamical processes in that region that enable asteroids to increase and decrease their inclinations.


title: Fast charging of quantum battery assisted by noise

abstract: We investigate the performance of a quantum battery exposed to local Markovian and non-Markovian dephasing noises. The battery is initially prepared as the ground state of a one-dimensional transverse $XY$ model with open boundary condition and is charged (discharged) via interactions with local bosonic reservoirs. We show that in the transient regime, quantum battery (QB) can store energy faster and has a higher maximum extractable work, quantified via ergotropy, when it is affected by local phase-flip or bit-flip Markovian noise compared to the case when there is no noise in the system. In both the charging and discharging processes, we report the enhancement in work-output as well as in ergotropy when all the spins are affected by non-Markovian Ohmic bath both in the transient and the steady-state regimes, thereby showing a counter-intuitive advantage of decoherence in QB. Both in Markovian and non-Markovian cases, we identify the system parameters and the corresponding noise models which lead to maximum enhancement of work-output and ergotropy. Moreover, we show that the benefit due to noise persists even with the initial state being prepared at a moderate temperature.


title: Recovery of surfaces and functions in high dimensions: sampling theory and links to neural networks

abstract: Several imaging algorithms including patch-based image denoising, image time series recovery, and convolutional neural networks can be thought of as methods that exploit the manifold structure of signals. While the empirical performance of these algorithms is impressive, the understanding of recovery of the signals and functions that live on manifold is less understood. In this paper, we focus on the recovery of signals that live on a union of surfaces. In particular, we consider signals living on a union of smooth band-limited surfaces in high dimensions. We show that an exponential mapping transforms the data to a union of low-dimensional subspaces. Using this relation, we introduce a sampling theoretical framework for the recovery of smooth surfaces from few samples and the learning of functions living on smooth surfaces. The low-rank property of the features is used to determine the number of measurements needed to recover the surface. Moreover, the low-rank property of the features also provides an efficient approach, which resembles a neural network, for the local representation of multidimensional functions on the surface.


title: Finding an induced path that is not a shortest path

abstract: We give a polynomial-time algorithm that, with input a graph $G$ and two vertices $u,v$ of $G$, decides whether there is an induced $uv$-path that is longer than the shortest $uv$-path.


title: Geometrical and Mechanical Characterisation of Hollow Thermoplastic Microspheres for Syntactic Foam Applications

abstract: Recently, hollow thermoplastic microspheres, such as Expancel made by Nouryon, have emerged as an innovative filler material for use in polymer-matrix composites. The resulting all-polymer syntactic foam takes on excellent damage tolerance properties, strong recoverability under large strains, and favourable energy dissipation characteristics. Despite finding increasing usage in various industries and applications, including in coatings, films, sealants, packaging, composites for microfluidics, medical ultrasonics and cementious composites, there is a near-complete absence of statistical geometrical information for Expancel microspheres. Further, their mechanical properties have not yet been reported. In this work we characterise the geometrical quantities of two classes of Expancel thermoplastic microspheres using X-ray computed tomography, focused ion beam and electron microscopy. We also observe the spatial distribution of microspheres within a polyurethane-matrix syntactic foam. We show that the volume-weighted polydisperse shell diameter in both classes of microsphere follows a normal distribution. Interestingly, polydispersity of the shell wall thickness is not observed and in particular the shell thickness is not correlated to the shell diameter. We employ the measured geometrical information in analytical micromechanical techniques in the small strain regime to determine, for the first time, estimates of the Young's modulus and Poisson's ratio of the microsphere shell material. Our results contribute to potential future improvements in the design and fabrication of syntactic foams that employ thermoplastic microspheres. Given the breadth of fields which utilise thermoplastic microspheres, we anticipate that our results, together with the methods used, will be of use in a much broader context in future materials research.


title: Instantons and Khovanov skein homology on $I\times T^2$

abstract: Asaeda, Przytycki and Sikora defined a generalization of Khovanov homology for links in $I$-bundles over compact surfaces. We prove that for a link $L\subset (-1,1)\times T^2$, the Asaeda-Przytycki-Sikora homology of $L$ has rank $2$ with $\mathbb{Z}/2$-coefficients if and only if $L$ is isotopic to an embedded knot in $\{0\}\times T^2$.


title: Time-Variant Variational Transfer for Value Functions

abstract: In most of the transfer learning approaches to reinforcement learning (RL) the distribution over the tasks is assumed to be stationary. Therefore, the target and source tasks are i.i.d. samples of the same distribution. In the context of this work, we consider the problem of transferring value functions through a variational method when the distribution that generates the tasks is time-variant, proposing a solution that leverages this temporal structure inherent in the task generating process. Furthermore, by means of a finite-sample analysis, the previously mentioned solution is theoretically compared to its time-invariant version. Finally, we will provide an experimental evaluation of the proposed technique with three distinct temporal dynamics in three different RL environments.


title: Ab initio electronic stopping power for protons in Ga$_{0.5}$In$_{0.5}$P/GaAs/Ge triple-junction solar cells for space applications

abstract: Motivated by the radiation damage of solar panels in space, firstly, the results of Monte Carlo particle transport simulations are presented for proton impact on triple-junction Ga$_{0.5}$In$_{0.5}$P/GaAs/Ge solar cells, showing the proton projectile penetration in the cells as a function of energy. It is followed by a systematic {\it ab initio} investigation of the electronic stopping power for protons in different layers of the cell at the relevant velocities via real-time time-dependent density functional theory (RT-TDDFT) calculations. The electronic stopping power is found to depend significantly on different channeling conditions, which should affect the low velocity damage predictions, and which are understood in terms of impact parameter and electron density along the path. Additionally, we explore the effect of the interface between the layers of the multilayer structure on the energy loss of a proton, along with the effect of strain in the lattice-matched solar cell. Both effects are found to be small compared with the main bulk effect. The interface energy loss has been found to increase with decreasing proton velocity, and in one case, there is an effective interface energy gain.


title: Background magnetic field and quantum correlations in the Schwinger effect

abstract: In this work we consider two complex scalar fields distinguished by their masses coupled to constant background electric and magnetic fields in the $(3+1)$-dimensional Minkowski spacetime and subsequently investigate a few measures quantifying the quantum correlations between the created particle-antiparticle Schwinger pairs. Since the background magnetic field itself cannot cause the decay of the Minkowski vacuum, our chief motivation here is to investigate the interplay between the effects due to the electric and magnetic fields. We start by computing the entanglement entropy for the vacuum state of a single scalar field. Second, we consider some maximally entangled states for the two-scalar field system and compute the logarithmic negativity and the mutual information. Qualitative differences of these results pertaining to the charge content of the states are emphasised. Based upon these results, we suggest some possible effects of a background magnetic field on the degradation of entanglement between states in an accelerated frame, for charged quantum fields.


title: A Systematic Study on the Absorption Features of Interstellar Ices in Presence of Impurities

abstract: Spectroscopic studies play a key role in the identification and analysis of interstellar ices and their structure. Some molecules have been identified within the interstellar ices either as pure, mixed, or even as layered structures. Absorption band features of water ice can significantly change with the presence of different types of impurities (CO, CO2, CH3OH, H2CO, etc.). In this work, we carried out a theoretical investigation to understand the behavior of water band frequency, and strength in the presence of impurities. The computational study has been supported and complemented by some infrared spectroscopy experiments aimed at verifying the effect of HCOOH, NH3 , and CH3 OH on the band profiles of pure H2O ice. Specifically, we explored the effect on the band strength of libration, bending, bulk stretching, and free-OH stretching modes. Computed band strength profiles have been compared with our new and existing experimental results, thus pointing out that vibrational modes of H2O and their intensities can change considerably in the presence of impurities at different concentrations. In most cases, the bulk stretching mode is the most affected vibration, while the bending is the least affected mode. HCOOH was found to have a strong influence on the libration, bending, and bulk stretching band profiles. In the case of NH3, the free-OH stretching band disappears when the impurity concentration becomes 50%. This work will ultimately aid a correct interpretation of future detailed spaceborne observations of interstellar ices by means of the upcoming JWST mission.


title: K$_2$CoS$_2$: A new two-dimensional in-plane antiferromagnetic insulator

abstract: The recent discovery of two-dimensional (2D) magnetic materials has brought magnetism to the flatland. This has opened up exciting opportunities for the exploration of fundamental physics as well as for novel device applications. Here, we predict a new thermodynamically stable 2D magnetic material, K$_2$CoS$_2$, which retains its in-plane anti-ferromagnetic order down to the monolayer and bilayer limits. We find that the magnetic moments ($2.5 \mu_B/$Co) are aligned along the intra-Co chains, from monolayer to bulk. The non-magnetic electronic spectrum of both the monolayer and bilayer films is found to host flat bands and van-Hove singularities, which in instrumental in giving rise to the the magnetic ground state. Based on classical Monte-Carlo simulations, we estimate the Neel temperature for the antiferromagnetic monolayer to be $\approx 15$K. Our study thus establishes that K$_2$CoS$_2$ hosts a robust antiferromagnetic state which persists from the monolayer limit to the bulk material.


title: Global X-Ray Properties of the Vela and Puppis A Supernova Remnants

abstract: The Vela and Puppis A supernova remnants (SNRs) comprise a large emission region of $\sim 8^{\circ}$ diameter in the soft X-ray sky. The HaloSat CubeSat mission provides the first soft X-ray ($0.4-7$ keV) observation of the entire Vela SNR and Puppis A SNR region with a single pointing and moderate spectral resolution. HaloSat observations of the Vela SNR are best fit with a two-temperature thermal plasma model consisting of a cooler component with $kT_{1} = 0.19^{+0.01}_{-0.01}$ keV in collisional ionization equilibrium and a hotter component with $kT_{2} = 1.06^{+0.45}_{-0.27}$ keV in non-equilibrium ionization. Observations of the Puppis A SNR are best fit with a single-component plane-parallel shocked plasma model with $kT = 0.86^{+0.06}_{-0.05}$ keV in non-equilibrium ionization. For the first time, we find the total X-ray luminosities of both components of the Vela SNR spectrum in the $0.5-7$ keV energy band to be $L_X = 4.4^{+1.4}_{-1.4} \times 10^{34}$ erg s$^{-1}$ for the cooler component and $L_X = 4.1^{+1.8}_{-1.5} \times 10^{34}$ erg s$^{-1}$ for the hotter component. We find the total X-ray luminosities of the Vela and Puppis A SNRs to be $L_{\text{X}} = 8.4 \times 10^{34}$ erg s$^{-1}$ and $L_X = 6.7^{+1.1}_{-0.9} \times 10^{36}$ erg s$^{-1}$.


title: Superfluid and supersolid phases of 4He on the second layer of graphite

abstract: We revisited the phase diagram of the second layer of 4He on top of graphite using quantum Monte Carlo methods. Our aim was to explore the existence of the novel phases suggested recently in experimental works, and determine their properties and stability limits. We found evidence of a superfluid quantum phase with hexatic correlations, induced by the corrugation of the first Helium layer, and a quasi-two-dimensional supersolid corresponding to a 7/12 registered phase. The 4/7 commensurate solid was found to be unstable, while the triangular incommensurate crystals, stable at large densities, were normal.


title: Vortices evolution in the solar atmosphere: A dynamical equation for the swirling strength

abstract: We study vortex dynamics in the solar atmosphere by employing and deriving the analytical evolution equations of two vortex identification criteria. The two criteria used are vorticity and the swirling strength. Vorticity can be biased in the presence of shear flows, but its dynamical equation is well known; the swirling strength is a more precise criterion for the identification of vortical flows, but its evolution equation is not known yet. Therefore, we explore the possibility of deriving a dynamical equation for the swirling strength. We then apply the two equations to analyze radiative MHD simulations of the solar atmosphere produced with the CO5BOLD code. We present a detailed review of the swirling strength criterion and the mathematical derivation of its evolution equation. This equation did not exist in the literature before and it constitutes a novel tool that is suitable for the analysis of a wide range of problems in (magneto-)hydrodynamics. By applying this equation to numerical models, we find that hydrodynamical and magnetic baroclinicities are the driving physical processes responsible for vortex generation in the convection zone and the photosphere. Higher up in the chromosphere, the magnetic terms alone dominate. Moreover, we find that the swirling strength is produced at small scales in a chaotic fashion, especially inside magnetic flux concentrations. The swirling strength represents an appropriate criterion for the identification of vortices in turbulent flows, such as those in the solar atmosphere. Moreover, its evolution equation, which is derived in this paper, is pivotal for obtaining precise information about the dynamics of these vortices and the physical mechanisms responsible for their production and evolution. Since this equation is available, the swirling strength is now the ideal quantity to study the dynamics of vortices in (magneto-)hydrodynamics.


title: End-to-End Object Detection with Transformers

abstract: We present a new method that views object detection as a direct set prediction problem. Our approach streamlines the detection pipeline, effectively removing the need for many hand-designed components like a non-maximum suppression procedure or anchor generation that explicitly encode our prior knowledge about the task. The main ingredients of the new framework, called DEtection TRansformer or DETR, are a set-based global loss that forces unique predictions via bipartite matching, and a transformer encoder-decoder architecture. Given a fixed small set of learned object queries, DETR reasons about the relations of the objects and the global image context to directly output the final set of predictions in parallel. The new model is conceptually simple and does not require a specialized library, unlike many other modern detectors. DETR demonstrates accuracy and run-time performance on par with the well-established and highly-optimized Faster RCNN baseline on the challenging COCO object detection dataset. Moreover, DETR can be easily generalized to produce panoptic segmentation in a unified manner. We show that it significantly outperforms competitive baselines. Training code and pretrained models are available at https://github.com/facebookresearch/detr.


title: Benchmarking Graph Data Management and Processing Systems: A Survey

abstract: The development of scalable, representative, and widely adopted benchmarks for graph data systems have been a question for which answers has been sought for decades. We conduct an in-depth study of the existing literature on benchmarks for graph data management and processing, covering 20 different benchmarks developed during the last 15 years. We categorize the benchmarks into three areas focusing on benchmarks for graph processing systems, graph database benchmarks, and bigdata benchmarks with graph processing workloads. This systematic approach allows us to identify multiple issues existing in this area, including i) few benchmarks exist which can produce high workload scenarios, ii) no significant work done on benchmarking graph stream processing as well as graph based machine learning, iii) benchmarks tend to use conventional metrics despite new meaningful metrics have been around for years, iv) increasing number of big data benchmarks appear with graph processing workloads. Following these observations, we conclude the survey by describing key challenges for future research on graph data systems benchmarking.


title: Quantum Divide and Compute: Hardware Demonstrations and Noisy Simulations

abstract: Noisy, intermediate-scale quantum computers come with intrinsic limitations in terms of the number of qubits (circuit "width") and decoherence time (circuit "depth") they can have. Here, for the first time, we demonstrate a recently introduced method that breaks a circuit into smaller subcircuits or fragments, and thus makes it possible to run circuits that are either too wide or too deep for a given quantum processor. We investigate the behavior of the method on one of IBM's 20-qubit superconducting quantum processors with various numbers of qubits and fragments. We build noise models that capture decoherence, readout error, and gate imperfections for this particular processor. We then carry out noisy simulations of the method in order to account for the observed experimental results. We find an agreement within 20% between the experimental and the simulated success probabilities, and we observe that recombining noisy fragments yields overall results that can outperform the results without fragmentation.


title: Comprehensive analysis of the tidal effect in gravitational waves and implication for cosmology

abstract: Detection of gravitational waves (GWs) produced by coalescence of compact binaries provides a novel way to measure the luminosity distance of GW events. Combining their redshift, they can act as standard sirens to constrain cosmological parameters. For various GW detector networks in 2nd-generation (2G), 2.5G and 3G, we comprehensively analyze the method to constrain the equation-of-state (EOS) of binary neutron-stars (BNSs) and extract their redshifts through the imprints of tidal effects in GW waveforms. We find for these events, the observations of electromagnetic counterparts in low-redshift range $z < 0.1$ are important for constraining the tidal effects. Considering 17 different EOSs of NSs or quark-stars, we find GW observations have strong capability to determine the EOS. Applying the events as standard sirens, and considering the constraints of NS's EOS derived from low-redshift observations as prior, we can constrain the dark-energy EOS parameters $w_0$ and $w_a$. In 3G era, the potential constraints are $\Delta w_0\in (0.0006,0.004)$ and $\Delta w_a\in(0.004,0.02)$, which are 1-3 orders smaller than those from traditional methods, including Type Ia supernovas and baryon acoustic oscillations. The constraints are also 1 order smaller than the method of GW standard siren by fixing the redshifts through short-hard $\gamma$-ray bursts, due to more available GW events in this method. Therefore, GW standard sirens, based on the tidal effect measurement, provide a realizable and much more powerful tool in cosmology.


title: A Survey of Languages for Formalizing Mathematics

abstract: In order to work with mathematical content in computer systems, it is necessary to represent it in formal languages. Ideally, these are supported by tools that verify the correctness of the content, allow computing with it, and produce human-readable documents. These goals are challenging to combine and state-of-the-art tools typically have to make difficult compromises. In this paper we discuss languages that have been created for this purpose, including logical languages of proof assistants and other formal systems, semi-formal languages, intermediate languages for exchanging mathematical knowledge, and language frameworks that allow building customized languages. We evaluate their advantages based on our experience in designing and applying languages and tools for formalizing mathematics. We reach the conclusion that no existing language is truly good enough yet and derive ideas for possible future improvements.


title: Progress of Quantum Molecular Dynamics model and its applications in Heavy Ion Collisions

abstract: In this review article, we first briefly introduce the transport theory and quantum molecular dynamics model applied in the study of the heavy ion collisions from low to intermediate energies. The developments of improved quantum molecular dynamics model (ImQMD) and ultra-relativistic quantum molecular dynamics model (UrQMD), are reviewed. The reaction mechanism and phenomena related to the fusion, multinucleon transfer, fragmentation, collective flow and particle production are reviewed and discussed within the framework of the two models. The constraints on the isospin asymmetric nuclear equation of state and in-medium nucleon-nucleon cross sections by comparing the heavy ion collision data with transport models calculations in last decades are also discussed, and the uncertainties of these constraints are analyzed as well. Finally, we discuss the future direction of the development of the transport models for improving the understanding of the reaction mechanism, the descriptions of various observables, the constraint on the nuclear equation of state, as well as for the constraint on in-medium nucleon-nucleon cross sections.


title: The optical signature of few-layer ReSe$_2$

abstract: Optical properties of thin layers of rhenium diselenide (ReSe$_2$) with thickness ranging from mono- (1 ML) to nona-layer (9 MLs) are demonstrated. The photoluminescence (PL) and Raman scattering were measured at low ($T$=5 K) and room ($T$=300 K) temperature, respectively. The PL spectra of ReSe$_2$ layers display two well-resolved emission lines, which blueshift by about 120 meV when the layer thickness decreases from 9 MLs to a monolayer. A rich structure of the observed low-energy Raman scattering modes can be explained within a linear chain model. The two phonon modes of intralayer vibrations, observed in Raman scattering spectra at about 120 cm$^{-1}$, exhibit very sensitive and opposite evolution as a function of layer thickness. It is shown that their energy difference can serve as a convenient and reliable tool to determine the thickness of ReSe$_2$ flakes in the few-layer limit.


title: Global sensitivity analysis informed model reduction and selection applied to a Valsalva maneuver model

abstract: In this study, we develop a methodology for model reduction and selection informed by global sensitivity analysis (GSA) methods. We apply these techniques to a control model that takes systolic blood pressure and thoracic tissue pressure data as inputs and predicts heart rate in response to the Valsalva maneuver (VM). The study compares four GSA methods based on Sobol' indices (SIs) quantifying the parameter influence on the difference between the model output and the heart rate data. The GSA methods include standard scalar SIs determining the average parameter influence over the time interval studied and three time-varying methods analyzing how parameter influence changes over time. The time-varying methods include a new technique, termed limited-memory SIs, predicting parameter influence using a moving window approach. Using the limited-memory SIs, we perform model reduction and selection to analyze the necessity of modeling both the aortic and carotid baroreceptor regions in response to the VM. We compare the original model to three systematically reduced models including (i) the aortic and carotid regions, (ii) the aortic region only, and (iii) the carotid region only. Model selection is done quantitatively using the Akaike and Bayesian Information Criteria and qualitatively by comparing the neurological predictions. Results show that it is necessary to incorporate both the aortic and carotid regions to model the VM.


title: Learning Local Features with Context Aggregation for Visual Localization

abstract: Keypoint detection and description is fundamental yet important in many vision applications. Most existing methods use detect-then-describe or detect-and-describe strategy to learn local features without considering their context information. Consequently, it is challenging for these methods to learn robust local features. In this paper, we focus on the fusion of low-level textual information and high-level semantic context information to improve the discrimitiveness of local features. Specifically, we first estimate a score map to represent the distribution of potential keypoints according to the quality of descriptors of all pixels. Then, we extract and aggregate multi-scale high-level semantic features based by the guidance of the score map. Finally, the low-level local features and high-level semantic features are fused and refined using a residual module. Experiments on the challenging local feature benchmark dataset demonstrate that our method achieves the state-of-the-art performance in the local feature challenge of the visual localization benchmark.


title: Evaluating epidemic forecasts in an interval format

abstract: For practical reasons, many forecasts of case, hospitalization and death counts in the context of the current COVID-19 pandemic are issued in the form of central predictive intervals at various levels. This is also the case for the forecasts collected in the COVID-19 Forecast Hub (https://covid19forecasthub.org/). Forecast evaluation metrics like the logarithmic score, which has been applied in several infectious disease forecasting challenges, are then not available as they require full predictive distributions. This article provides an overview of how established methods for the evaluation of quantile and interval forecasts can be applied to epidemic forecasts in this format. Specifically, we discuss the computation and interpretation of the weighted interval score, which is a proper score that approximates the continuous ranked probability score. It can be interpreted as a generalization of the absolute error to probabilistic forecasts and allows for a decomposition into a measure of sharpness and penalties for over- and underprediction.


title: Factorizations of tropical and sign polynomials

abstract: In this text, we study factorizations of polynomials over the tropical hyperfield and the sign hyperfield, which we call `tropical polynomials' and `sign polynomials', respectively. We classify all irreducible polynomials in either case. We show that tropical polynomials factor uniquely into irreducible factors, but that unique factorization fails for sign polyomials. We describe division algorithms for tropical and sign polynomials by linear terms that correspond to roots of the polynomials.


title: Global Description of Beta Decay with the Axially-Deformed Skyrme Finite Amplitude Method: Extension to Odd-Mass and Odd-Odd Nuclei

abstract: We use the finite amplitude method (FAM), an efficient implementation of the quasiparticle random phase approximation, to compute beta-decay rates with Skyrme energy-density functionals for 3983 nuclei, essentially all the medium-mass and heavy isotopes on the neutron rich side of stability. We employ an extension of the FAM that treats odd-mass and odd-odd nuclear ground states in the equal filling approximation. Our rates are in reasonable agreement both with experimental data where available and with rates from other global calculations.


title: Cross Hashing: Anonymizing encounters in Decentralised Contact Tracing Protocols

abstract: During the COVID-19 (SARS-CoV-2) epidemic, Contact Tracing emerged as an essential tool for managing the epidemic. App-based solutions have emerged for Contact Tracing, including a protocol designed by Apple and Google (influenced by an open-source protocol known as DP3T). This protocol contains two well-documented de-anonymisation attacks. Firstly that when someone is marked as having tested positive and their keys are made public, they can be tracked over a large geographic area for 24 hours at a time. Secondly, whilst the app requires a minimum exposure duration to register a contact, there is no cryptographic guarantee for this property. This means an adversary can scan Bluetooth networks and retrospectively find who is infected. We propose a novel "cross hashing" approach to cryptographically guarantee minimum exposure durations. We further mitigate the 24-hour data exposure of infected individuals and reduce computational time for identifying if a user has been exposed using $k$-Anonymous buckets of hashes and Private Set Intersection. We empirically demonstrate that this modified protocol can offer like-for-like efficacy to the existing protocol.


title: Higgs chameleon

abstract: The existing constraints from particle colliders reveal a suspicious but nonlethal metastability for our current electroweak vacuum of Higgs potential in the standard model of particle physics, which is, however, disfavored in the early Universe if the inflationary Hubble scale is larger than the instability scale when Higgs quartic self-coupling runs into negative value. Alternative to previous trials of acquiring a positive effective mass-squared from Higgs quadratic couplings to Ricci scalar or inflaton field, we propose a third approach to stabilize the Higgs potential in the early Universe by regarding Higgs as chameleon coupled to inflaton alone without conflicting to the present constraints on either Higgs or chameleon.


title: Constraining the weights of Stokes Polytopes using BCFW recursions for Phi^4

abstract: The relationship between certain geometric objects called polytopes and scattering amplitudes has revealed deep structures in QFTs. It has been developed in great depth at the tree- and loop-level amplitudes in $\mathcal{N}=4~\text{SYM}$ theory and has been extended to the scalar $\phi^3$ and $\phi^4$ theories at tree-level. In this paper, we use the generalized BCFW recursion relations for massless planar $\phi^4$ theory to constrain the weights of a class of geometric objects called Stokes polytopes, which manifest in the geometric formulation of $\phi^4$ amplitudes. We see that the weights of the Stokes polytopes are intricately tied to the boundary terms in $\phi^4$ theories. We compute the weights of $N=1,2$, and $3$ dimensional Stokes polytopes corresponding to six-, eight- and ten-point amplitudes respectively. We generalize our results to higher-point amplitudes and show that the generalized BCFW recursions uniquely fix the weights for an $n$-point amplitude.


title: Dynamical History of the Uranian System

abstract: We numerically simulate the past tidal evolution of the five large moons of Uranus (Miranda, Ariel, Umbriel, Titania, and Oberon). We find that the most recent major mean-motion resonance (MMR) between any two moons, the Ariel-Umbriel 5:3 MMR, had a large effect on the whole system. Our results suggest that this resonance is responsible for the current 4.3$^{\circ}$ inclination of Miranda (instead of previously proposed 3:1 Miranda-Umbriel MMR), and that all five moons had their inclinations excited during this resonance. Miranda experienced significant tidal heating during the Ariel-Umbriel 5:3 MMR due to its eccentricity being excited by Ariel's secular perturbations. This tidal heating draws energy from shrinking of Miranda's orbit, rather than Ariel's outward evolution, and can generate heat flows in excess of 100 mW m$^{-2}$, sufficient to produce young coronae on Miranda. We find that this MMR was followed by a sequence of secular resonances, which reshuffled the moons' eccentricities and inclinations. We also find that the precession of Oberon's spin axis is close to a resonance with the precession of Umbriel's orbital plane, and that this spin-orbit resonance was likely excited during the Ariel-Umbriel 5:3 MMR. After the exit from the MMR, subsequent Ariel-Umbriel secular resonance and Oberon-Umbriel spin-orbit resonance may be able to explain the current low inclinations of Ariel and Umbriel. The age of Miranda's surface features tentatively suggests Uranian tidal $Q=15,000-20,000$, which can be further refined in future work.


title: Interactions among intermediate redshift galaxies. The case of SDSSJ134420.86+663717.8

abstract: We present the properties of the central supermassive black holes and the host galaxies of the interacting object SDSSJ134420.86+663717.8. We obtained optical long slit spectroscopy data from the Large Binocular Telescope (LBT) using the Multi Object Double Spectrograph (MODS). Analysing the spectra revealed several strong broad and narrow emission lines of ionised gas in the nuclear region of one galaxy, whereas only narrow emission lines were visible for the second galaxy. The optical spectra were used to plot diagnostic diagrams, deduce rotation curves of the two galaxies, and calculate the masses of the central supermassive black holes. We find that the galaxy with broad emission line features has Seyfert~1 properties, while the galaxy with only narrow emission line features seems to be star-forming in nature. Furthermore, we find that the masses of the central supermassive black holes are almost equal at a few times 10^7 solar mass. Additionally, we present a simple N-body simulation to shed some light on the initial conditions of the progenitor galaxies. We find that for an almost orthogonal approach of the two interacting galaxies, the model resembles the optical image of the system.


title: Refining Implicit Argument Annotation for UCCA

abstract: Predicate-argument structure analysis is a central component in meaning representations of text. The fact that some arguments are not explicitly mentioned in a sentence gives rise to ambiguity in language understanding, and renders it difficult for machines to interpret text correctly. However, only few resources represent implicit roles for NLU, and existing studies in NLP only make coarse distinctions between categories of arguments omitted from linguistic form. This paper proposes a typology for fine-grained implicit argument annotation on top of Universal Conceptual Cognitive Annotation's foundational layer. The proposed implicit argument categorisation is driven by theories of implicit role interpretation and consists of six types: Deictic, Generic, Genre-based, Type-identifiable, Non-specific, and Iterated-set. We exemplify our design by revisiting part of the UCCA EWT corpus, providing a new dataset annotated with the refinement layer, and making a comparative analysis with other schemes.


title: A reinforcement learning approach to rare trajectory sampling

abstract: Very often when studying non-equilibrium systems one is interested in analysing dynamical behaviour that occurs with very low probability, so called rare events. In practice, since rare events are by definition atypical, they are often difficult to access in a statistically significant way. What are required are strategies to "make rare events typical" so that they can be generated on demand. Here we present such a general approach to adaptively construct a dynamics that efficiently samples atypical events. We do so by exploiting the methods of reinforcement learning (RL), which refers to the set of machine learning techniques aimed at finding the optimal behaviour to maximise a reward associated with the dynamics. We consider the general perspective of dynamical trajectory ensembles, whereby rare events are described in terms of ensemble reweighting. By minimising the distance between a reweighted ensemble and that of a suitably parametrised controlled dynamics we arrive at a set of methods similar to those of RL to numerically approximate the optimal dynamics that realises the rare behaviour of interest. As simple illustrations we consider in detail the problem of excursions of a random walker, for the case of rare events with a finite time horizon; and the problem of a studying current statistics of a particle hopping in a ring geometry, for the case of an infinite time horizon. We discuss natural extensions of the ideas presented here, including to continuous-time Markov systems, first passage time problems and non-Markovian dynamics.


title: Dispersion measure components within host galaxies of Fast Radio Bursts: observational constraints from statistical properties of FRBs

abstract: Dispersion measure (DM) of Fast Radio Bursts (FRBs) are commonly used as a indicator of distance assuming that DM in excess of the expected amount within the Milky Way in the direction of each FRB arise mostly from the inter-galactic medium. However, the assumption might not be true if, for example, most FRB progenitors are embedded in ionized circumstellar material (CSM, e.g. supernova remnant). In this study, we jointly analyze distributions of DM, flux density, and fluence of the FRB samples observed by the Parkes telescope and the Australian Square Kilometre Array Pathfinder (ASKAP) using analytical models of FRBs, to constrain fractions of various DM components that shape the overall DM distribution and emission properties of FRBs. Comparing the model predictions with the observations we find that the typical amount of DM in each FRB host galaxy is $\sim 120$ cm$^{-3}$pc which is naturally explained as a combination of interstellar medium (ISM) and halo of an ordinary galaxy, without additional contribution from ionized CSM that is directly associated with an FRB progenitor. Furthermore, we also find that observed flux densities of FRBs do not statistically suffer strong $K$-correction, i.e. the typical luminosity density of FRBs does not significantly change within the range of emitting frequency $\nu_{\rm rest} \sim$ 1--4 GHz.


title: Minimizing Supervision in Multi-label Categorization

abstract: Multiple categories of objects are present in most images. Treating this as a multi-class classification is not justified. We treat this as a multi-label classification problem. In this paper, we further aim to minimize the supervision required for providing supervision in multi-label classification. Specifically, we investigate an effective class of approaches that associate a weak localization with each category either in terms of the bounding box or segmentation mask. Doing so improves the accuracy of multi-label categorization. The approach we adopt is one of active learning, i.e., incrementally selecting a set of samples that need supervision based on the current model, obtaining supervision for these samples, retraining the model with the additional set of supervised samples and proceeding again to select the next set of samples. A crucial concern is the choice of the set of samples. In doing so, we provide a novel insight, and no specific measure succeeds in obtaining a consistently improved selection criterion. We, therefore, provide a selection criterion that consistently improves the overall baseline criterion by choosing the top k set of samples for a varied set of criteria. Using this criterion, we are able to show that we can retain more than 98% of the fully supervised performance with just 20% of samples (and more than 96% using 10%) of the dataset on PASCAL VOC 2007 and 2012. Also, our proposed approach consistently outperforms all other baseline metrics for all benchmark datasets and model combinations.


title: Compositions of pseudo-symmetric integrators with complex coefficients for the numerical integration of differential equations

abstract: In this paper, we are concerned with the construction and analysis of a new class of methods obtained as double jump compositions with complex coefficients and projection on the real axis. It is shown in particular that the new integrators are symmetric and symplectic up to high orders if one uses a symmetric and symplectic basic method. In terms of efficiency, the aforementioned technique requires fewer stages than standard compositions of the same orders and is thus expected to lead to faster methods.


title: Dimension Reduction-based Signal Compression for Uplink Distributed MIMO C-RAN with Limited Fronthaul Capacity

abstract: This paper proposes a dimension reduction-based signal compression scheme for uplink distributed MIMO cloud radio access networks (C-RAN) with an overall excess of receive antennas, in which users are jointly served by distributed multi-antenna receivers connected to a central processor via individual finite-capacity fronthaul links. We first show that, under quantization noise-limited operation, applying linear dimension reduction at each receiver before compressing locally with a uniform quantization noise level results in a sum capacity that scales approximately linearly with fronthaul capacity, and can come within a fixed gap of the cut-set bound. The dimension reduction filters that maximize joint mutual information are then shown to be truncated forms of the conditional Karhunen-Loeve transform, with a block coordinate ascent algorithm for finding a stationary point given. Analysis and numerical results indicate that the signal dimension can be reduced without significant loss of information, particularly at high signal-to-noise ratio, preserving the benefits of using excess antennas. The method is then adapted for the case of imperfect channel state information at the receivers. The scheme significantly outperforms conventional local signal compression at all fronthaul rates, and with complexity linear in network size represents a scalable solution for distributed MIMO C-RAN systems.


title: Communicating over the Torn-Paper Channel

abstract: We consider the problem of communicating over a channel that randomly "tears" the message block into small pieces of different sizes and shuffles them. For the binary torn-paper channel with block length $n$ and pieces of length ${\rm Geometric}(p_n)$, we characterize the capacity as $C = e^{-\alpha}$, where $\alpha = \lim_{n\to\infty} p_n \log n$. Our results show that the case of ${\rm Geometric}(p_n)$-length fragments and the case of deterministic length-$(1/p_n)$ fragments are qualitatively different and, surprisingly, the capacity of the former is larger. Intuitively, this is due to the fact that, in the random fragments case, large fragments are sometimes observed, which boosts the capacity.


title: On reducible non-Weierstrass semigroups

abstract: Weierstrass semigroups are well-known along the literature. We present a new family of non-Weierstrass semigroups which can be written as an intersection of Weierstrass semigroups. In addition, we provide methods for calculating non-Weierstrass semigroups with genus as large as desired.


title: Non-Markovian dynamics of a two-level system in a bosonic bath and a Gaussian fluctuating environment with finite correlation time

abstract: The character of evolution of an open quantum system is often encoded in the correlation function of the environment or, equivalently, in the spectral density function of the interaction. When the environment is heterogeneous, e.g. consists of several independent subenvironments with different spectral functions, one of the subenvironments can be considered auxiliary and used to control decoherence of the open quantum system in the remaining part of the environment. The control can be realized, for example, by adjusting the character of interaction with the subenvironment via suitable parameters of its spectral density. We investigate non-Markovian evolution of a two-level system (qubit) under influence of three independent decoherence channels, two of them have classical nature and originate from interaction with a stochastic field, and the third is a quantum channel formed by interaction with a bosonic bath. By modifying spectral densities of the channels, we study their impact on steady states of the two-level system, evolution of its density matrix and the equilibrium emission spectrums. Additionally, we investigate the impact of the rotation-wave approximation applied to the bath channel on accuracy of the results.


title: Analysis of the Penn Korean Universal Dependency Treebank (PKT-UD): Manual Revision to Build Robust Parsing Model in Korean

abstract: In this paper, we first open on important issues regarding the Penn Korean Universal Treebank (PKT-UD) and address these issues by revising the entire corpus manually with the aim of producing cleaner UD annotations that are more faithful to Korean grammar. For compatibility to the rest of UD corpora, we follow the UDv2 guidelines, and extensively revise the part-of-speech tags and the dependency relations to reflect morphological features and flexible word-order aspects in Korean. The original and the revised versions of PKT-UD are experimented with transformer-based parsing models using biaffine attention. The parsing model trained on the revised corpus shows a significant improvement of 3.0% in labeled attachment score over the model trained on the previous corpus. Our error analysis demonstrates that this revision allows the parsing model to learn relations more robustly, reducing several critical errors that used to be made by the previous model.


title: Effective convergence of coranks of random R\'edei matrices

abstract: We give effective estimates for the $l^1$-distance between the corank distribution of $r \times r$ R\'edei matrices and the measure predicted by the Cohen--Lenstra heuristics. To this end we pinpoint a class of stochastic processes, which we call $c$-transitioning. These stochastic processes are well approximated by Markov processes, and we give an effective ergodic theorem for such processes. With this tool we make effective a theorem of Gerth \cite{Gerth} that initiated the study of the Cohen--Lenstra heuristics for $p = 2$. Gerth's work triggered a series of developments that has recently culminated in the breakthrough of Smith \cite{Smith}. The present work will be used in upcoming work of the authors on further applications of Smith's ideas to the arithmetic of quadratic fields. To this end we extend our main result to several other families of matrix spaces that occur in the study of integral points on the equation $x^2 - dy^2 = l$ as $d$ varies.


title: Breaking the Sample Size Barrier in Model-Based Reinforcement Learning with a Generative Model

abstract: This paper is concerned with the sample efficiency of reinforcement learning, assuming access to a generative model (or simulator). We first consider $\gamma$-discounted infinite-horizon Markov decision processes (MDPs) with state space $\mathcal{S}$ and action space $\mathcal{A}$. Despite a number of prior works tackling this problem, a complete picture of the trade-offs between sample complexity and statistical accuracy is yet to be determined. In particular, all prior results suffer from a severe sample size barrier, in the sense that their claimed statistical guarantees hold only when the sample size exceeds at least $\frac{|\mathcal{S}||\mathcal{A}|}{(1-\gamma)^2}$. The current paper overcomes this barrier by certifying the minimax optimality of two algorithms -- a perturbed model-based algorithm and a conservative model-based algorithm -- as soon as the sample size exceeds the order of $\frac{|\mathcal{S}||\mathcal{A}|}{1-\gamma}$ (modulo some log factor). Moving beyond infinite-horizon MDPs, we further study time-inhomogeneous finite-horizon MDPs, and prove that a plain model-based planning algorithm suffices to achieve minimax-optimal sample complexity given any target accuracy level. To the best of our knowledge, this work delivers the first minimax-optimal guarantees that accommodate the entire range of sample sizes (beyond which finding a meaningful policy is information theoretically infeasible).

categories: cs.LG cs.IT math.IT math.OC math.ST stat.ML stat.TH


title: A Framework for Behavioral Biometric Authentication using Deep Metric Learning on Mobile Devices

abstract: Mobile authentication using behavioral biometrics has been an active area of research. Existing research relies on building machine learning classifiers to recognize an individual's unique patterns. However, these classifiers are not powerful enough to learn the discriminative features. When implemented on the mobile devices, they face new challenges from the behavioral dynamics, data privacy and side-channel leaks. To address these challenges, we present a new framework to incorporate training on battery-powered mobile devices, so private data never leaves the device and training can be flexibly scheduled to adapt the behavioral patterns at runtime. We re-formulate the classification problem into deep metric learning to improve the discriminative power and design an effective countermeasure to thwart side-channel leaks by embedding a noise signature in the sensing signals without sacrificing too much usability. The experiments demonstrate authentication accuracy over 95% on three public datasets, a sheer 15% gain from multi-class classification with less data and robustness against brute-force and side-channel attacks with 99% and 90% success, respectively. We show the feasibility of training with mobile CPUs, where training 100 epochs takes less than 10 mins and can be boosted 3-5 times with feature transfer. Finally, we profile memory, energy and computational overhead. Our results indicate that training consumes lower energy than watching videos and slightly higher energy than playing games.


title: Multiple abrupt phase transitions in urban transport congestion

abstract: During the last decades, the study of cities has been transformed by new approaches combining engineering and complexity sciences. Network theory is playing a central role, facilitating the quantitative analysis of crucial urban dynamics, such as mobility, city growth or urban planning. In this work, we focus on the spatial aspects of congestion. Analyzing a large amount of real city networks, we show that the location of the onset of congestion changes according to the considered urban area, defining, in turn, a set of congestion regimes separated by abrupt transitions. To help unveiling these spatial dependencies of congestion (in terms of network betweenness analysis), we introduce a family of planar road network models composed of a dense urban center connected to an arboreal periphery. These models, coined as GT and DT-MST models, allow us to analytically, numerically and experimentally describe how and why congestion emerges in particular geographical areas of monocentric cities and, subsequently, to describe the congestion regimes and the factors that promote the appearance of their abrupt transitions. We show that the fundamental ingredient behind the observed abrupt transitions is the spatial separation between the urban center and the periphery, and the number of separate areas that form the periphery. Elaborating on the implications of our results, we show that they may have an influence on the design and optimization of road networks regarding urban growth and the management of daily traffic dynamics.


title: On the origin of the weak equivalence principle in a theory of emergent quantum mechanics

abstract: We argue that in a framework for emergent quantum mechanics, the weak equivalence principle is a consequence of concentration of measure in large dimensional spaces of $1$-Lipshitz functions. Furthermore, as a consequence of the emergent framework and the properties that we assume for the fundamental dynamics, it is argued that gravity must be a classical, emergent interaction.


title: Post-Quantum Multi-Party Computation

abstract: We initiate the study of multi-party computation for classical functionalities (in the plain model) with security against malicious polynomial-time quantum adversaries. We observe that existing techniques readily give a polynomial-round protocol, but our main result is a construction of *constant-round* post-quantum multi-party computation. We assume mildly super-polynomial quantum hardness of learning with errors (LWE), and polynomial quantum hardness of an LWE-based circular security assumption. Along the way, we develop the following cryptographic primitives that may be of independent interest: 1. A spooky encryption scheme for relations computable by quantum circuits, from the quantum hardness of an LWE-based circular security assumption. This yields the first quantum multi-key fully-homomorphic encryption scheme with classical keys. 2. Constant-round zero-knowledge secure against multiple parallel quantum verifiers from spooky encryption for relations computable by quantum circuits. To enable this, we develop a new straight-line non-black-box simulation technique against *parallel* verifiers that does not clone the adversary's state. This forms the heart of our technical contribution and may also be relevant to the classical setting. 3. A constant-round post-quantum non-malleable commitment scheme, from the mildly super-polynomial quantum hardness of LWE.


title: Lyra -- containerized microservices for browsing shared biomedical data

abstract: Research papers in the biomedical field come with large and complex data sets that are shared with the scientific community as unstructured data files via public data repositories. Examples are sequencing, microarray, and mass spectroscopy data. The papers discuss and visualize only a small part of the data, the part that is in its research focus. For labs with similar but not identical research interests different parts of the data might be important. They can thus download the full data, preprocess it, integrate it with data from other publications and browse those parts that they are most interested in. This requires substantial work as well as programming and analysis expertise that only few biological labs have on board. In contrast, providing access to published data over web browsers makes all data visible, allows for easy interaction with it, and lowers the barrier to working with data from others.


title: SIRIUS project. I. Star formation models for star-by-star simulations of star clusters and galaxy formation

abstract: Most stars are formed as star clusters in galaxies, which then disperse into galactic disks. Upcoming exascale supercomputational facilities will enable performing simulations of galaxies and their formation by resolving individual stars (star-by-star simulations). This will substantially advance our understanding of star formation in galaxies, star cluster formation, and assembly histories of galaxies. In previous galaxy simulations, a simple stellar population approximation was used. It is, however, difficult to improve the mass resolution with this approximation. Therefore, a model for forming individual stars that can be used in simulations of galaxies must be established. In this first paper of a series of the SIRIUS (SImulations Resolving IndividUal Stars) project, we demonstrate a stochastic star formation model for star-by-star simulations. An assumed stellar initial mass function (IMF) is randomly assigned to newly formed stars. We introduce a maximum search radius to assemble the mass from surrounding gas particles to form star particles. In this study, we perform a series of N-body/smoothed particle hydrodynamics simulations of star cluster formations from turbulent molecular clouds and ultra-faint dwarf galaxies as test cases. The IMF can be correctly sampled if a maximum search radius that is larger than the value estimated from the threshold density for star formation is adopted. In small clouds, the formation of massive stars is highly stochastic because of the small number of stars. We confirm that the star formation efficiency and threshold density do not strongly affect the results. We find that our model can naturally reproduce the relationship between the most massive stars and the total stellar mass of star clusters. Herein, we demonstrate that our models can be applied to simulations varying from star clusters to galaxies for a wide range of resolutions.


title: Coordinated Multicell Beamforming and Power Allocation for Massive MIMO with Low-Resolution ADC/DAC

abstract: In this work, we present a solution for coordinated beamforming and power allocation when base stations employ a massive number of antennas equipped with low-resolution analog-to-digital and digital-to-analog converters. We address total power minimization problems of the coarsely quantized uplink (UL) and downlink (DL) communication systems with target signal-to-interference-plus-noise ratio (SINR) constraints. By combining the UL problem with minimum mean square error combiners and deriving the Lagrangian dual of the DL problem, we prove UL-DL duality and show there is no duality gap even with coarse data converters. Inspired by strong duality, we devise an iterative algorithm to determine the optimal UL transmit powers, and then linearly amplify the UL combiners with proper weights to acquire the optimal DL precoder. Simulation results validate strong duality and evaluate the proposed method in terms of total power consumption and achieved SINR.


title: Learning-Based Synthetic Dual Energy CT Imaging from Single Energy CT for Stopping Power Ratio Calculation in Proton Radiation Therapy

abstract: Purpose: Dual-energy CT (DECT) has been shown to derive stopping power ratio (SPR) map with higher accuracy than conventional single energy CT (SECT) by obtaining the energy dependence of photon interactions. However, DECT is not as widely implemented as SECT in proton radiation therapy simulation. This work presents a learning-based method to synthetize DECT images from SECT for proton radiation therapy. Methods: The proposed method uses a residual attention generative adversarial network. Residual blocks with attention gates were used to force the model focus on the difference between DECT maps and SECT images. To evaluate the accuracy of the method, we retrospectively investigated 20 head-and-neck cancer patients with both DECT and SECT scans available. The high and low energy CT images acquired from DECT acted as learning targets in the training process for SECT datasets and were evaluated against results from the proposed method using a leave-one-out cross-validation strategy. To evaluate our method in the context of a practical application, we generated SPR maps from sDECT using physics-based dual-energy stoichiometric method and compared the maps to those generated from DECT. Results: The synthesized DECT images showed an average mean absolute error around 30 Hounsfield Unit (HU) across the whole-body volume. The corresponding SPR maps generated from synthetic DECT showed an average normalized mean square error of about 1% with reduced noise level and artifacts than those from original DECT. Conclusions: The accuracy of the synthesized DECT image by our machine-learning-based method was evaluated on head and neck patient, and potential feasibility for proton treatment planning and dose calculation was shown by generating SPR map using the synthesized DECT.


abstract: Let $G$ be a simple graph with maximum degree $\Delta$. A classic result of Vizing shows that $\chi'(G)$, the chromatic index of $G$, is either $\Delta$ or $\Delta+1$. We say $G$ is of \emph{Class 1} if $\chi'(G)=\Delta$, and is of \emph{Class 2} otherwise. A graph $G$ is \emph{$\Delta$-critical} if $\chi'(G)=\Delta+1$ and $\chi'(H)<\Delta+1$ for every proper subgraph $H$ of $G$, and is \emph{overfull} if $|E(G)|>\Delta \lfloor (|V(G)|-1)/2 \rfloor$. Clearly, overfull graphs are Class 2. Hilton and Zhao in 1997 conjectured that if $G$ is obtained from an $n$-vertex $\Delta$-regular Class 1 graph with maximum degree greater than $n/3$ by splitting a vertex, then being overfull is the only reason for $G$ to be Class 2. This conjecture was only confirmed when $\Delta\ge \frac{n}{2}(\sqrt{7}-1)\approx 0.82n$. In this paper, we improve the bound on $\Delta$ from $\frac{n}{2}(\sqrt{7}-1)$ to $0.75n$. Considering the structure of $\Delta$-critical graphs with a vertex of degree 2, we also show that for an $n$-vertex $\Delta$-critical graph with $\Delta\ge \frac{3n}{4}$, if it contains a vertex of degree 2, then it is overfull. We actually obtain a more general form of this result, which partially supports the overfull conjecture of Chetwynd and Hilton from 1986, which states that if $G$ is an $n$-vertex $\Delta$-critical graph with $\Delta>n/3$, then $G$ contains an overfull subgraph $H$ with $\Delta(H)=\Delta$. Our proof techniques are new and might shed some light on attacking both of the conjectures when $\Delta$ is large.


title: Representations of the fixed point subalgebra of the vertex algebra associated to a non-degenerate even lattice by an automorphism of order $2$

abstract: Let $V_{L}$ be the vertex algebra associated to a non-degenerate even lattice $L$, $\theta$ the automorphism of $V_{L}$ induced from the $-1$-isometry of $L$, and $V_{L}^{+}$ the fixed point subalgebra of $V_{L}$ under the action of $\theta$. We show that every weak $V_{L}^{+}$-module is completely reducible.


title: Symbolic and Structural Model-Checking

abstract: Brute-force model-checking consists in exhaustive exploration of the state-space of a Petri net, and meets the dreaded state-space explosion problem. In contrast, this paper shows how to solve model-checking problems using a combination of techniques that stay in complexity proportional to the size of the net structure rather than to the state-space size. We combine an SMT based over-approximation to prove that some behaviors are unfeasible, an under-approximation using memory-less sampling of runs to find witness traces or counter-examples, and a set of structural reduction rules that can simplify both the system and the property. This approach was able to win by a clear margin the model-checking contest 2020 for reachability queries as well as deadlock detection, thus demonstrating the practical effectiveness and general applicability of the system of rules presented in this paper.


title: The generalized Lindemann melting coefficient

abstract: Lindemann developed the melting temperature theory over 100 years ago, known as the Lindemann criterion. Its main assumption is that melting occurs when the root-mean-square vibration amplitude of ions and atoms in crystals exceeds a critical fraction, h of the inter-atomic spacing in crystals. The Lindemann coefficient h is undefined and scientific papers report different h values for different elements. Here we present previously unobserved data trends pointing to the fact that the Lindemann coefficient could be linked to the periodic groups of the periodic table, having an exact value for each element belonging to a given periodic group. We report 12 distinctive Lindemann coefficient values corresponding to 12 groups of the periodic table containing solid elements with identifiable melting temperature. Using these vales, the recalculation of the melting temperatures indicates a good match to the experimental values for 39 elements, corresponding to 12 out of 15 periodic groups. This newly observed result opens up the possibility of further refining the Lindemann melting criterion by stimulating analytical studies of the Lindemann coefficient in the light of this newly discovered result.


title: The PVLAS experiment: a 25 year effort to measure vacuum magnetic birefringence

abstract: This paper describes the 25 year effort to measure vacuum magnetic birefringence and dichroism with the PVLAS experiment. The experiment went through two main phases: the first using a rotating superconducting magnet and the second using two rotating permanent magnets. The experiment was not able to reach the predicted value from QED. Nonetheless the experiment set the current best limits on vacuum magnetic birefringence and dichroism for a field of $B_{\rm ext} = 2.5$ T, namely, $\Delta n^{\rm (PVLAS)} = (12\pm17)\times10^{-23}$ and $|\Delta\kappa|^{\rm (PVLAS)} = (10\pm28)\times10^{-23}$. The uncertainty on $\Delta n^{\rm (PVLAS)}$ is about a factor 7 above the predicted value of $\Delta n^{\rm (QED)} = 2.5\times10^{-23}$ @ 2.5 T.


title: Class-Weighted Classification: Trade-offs and Robust Approaches

abstract: We address imbalanced classification, the problem in which a label may have low marginal probability relative to other labels, by weighting losses according to the correct class. First, we examine the convergence rates of the expected excess weighted risk of plug-in classifiers where the weighting for the plug-in classifier and the risk may be different. This leads to irreducible errors that do not converge to the weighted Bayes risk, which motivates our consideration of robust risks. We define a robust risk that minimizes risk over a set of weightings and show excess risk bounds for this problem. Finally, we show that particular choices of the weighting set leads to a special instance of conditional value at risk (CVaR) from stochastic programming, which we call label conditional value at risk (LCVaR). Additionally, we generalize this weighting to derive a new robust risk problem that we call label heterogeneous conditional value at risk (LHCVaR). Finally, we empirically demonstrate the efficacy of LCVaR and LHCVaR on improving class conditional risks.


title: Proportional Choosability of Complete Bipartite Graphs

abstract: Proportional choosability is a list analogue of equitable coloring that was introduced in 2019. The smallest $k$ for which a graph $G$ is proportionally $k$-choosable is the proportional choice number of $G$, and it is denoted $\chi_{pc}(G)$. In the first ever paper on proportional choosability, it was shown that when $2 \leq n \leq m$, $ \max\{ n + 1, 1 + \lceil m / 2 \rceil\} \leq \chi_{pc}(K_{n,m}) \leq n + m - 1$. In this note we improve on this result by showing that $ \max\{ n + 1, \lceil n / 2 \rceil + \lceil m / 2 \rceil\} \leq \chi_{pc}(K_{n,m}) \leq n + m -1- \lfloor m/3 \rfloor$. In the process, we prove some new lower bounds on the proportional choice number of complete multipartite graphs. We also present several interesting open questions.


title: Pressure balance in the multiphase ISM of cosmologically simulated disk galaxies

abstract: Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disk galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyze how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the midplane. Bulk flows (e.g., inflows and fountains) are important at a few disk scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total midplane pressure is well-predicted by the weight of the disk gas, and we show that it also scales linearly with the star formation rate surface density (Sigma_SFR). These results support the notion that the Kennicutt-Schmidt relation arises because Sigma_SFR and the gas surface density (Sigma_g) are connected via the ISM midplane pressure.


title: Structure of two-loop SMEFT anomalous dimensions via on-shell methods

abstract: We describe on-shell methods for computing one- and two-loop anomalous dimensions in the context of effective field theories containing higher-dimension operators. We also summarize methods for computing one-loop amplitudes, which are used as inputs to the computation of two-loop anomalous dimensions, and we explain how the structure of rational terms and judicious renormalization scheme choices can lead to additional vanishing terms in the anomalous dimension matrix at two loops. We describe the two-loop implications for the Standard Model Effective Field Theory (SMEFT). As a by-product of this analysis we verify a variety of one-loop SMEFT anomalous dimensions computed by Alonso, Jenkins, Manohar and Trott.


title: 128 Identical Quantum Sources Integrated on a Single Silica Chip

abstract: Quantum technology is playing an increasingly important role due to the intrinsic parallel processing capabilities endorsed by quantum superposition, exceeding upper limits of classical performances in diverse fields. Integrated photonic chip offers an elegant way to construct large-scale quantum systems in a physically scalable fashion, however, nonuniformity of quantum sources prevents all the elements from being connected coherently for exponentially increasing Hilbert space. Here, we experimentally demonstrate 128 identical quantum sources integrated on a single silica chip. By actively controlling the light-matter interaction in femtosecond laser direct writing, we are able to unify the properties of waveguides comprehensively and therefore the spontaneous four-wave mixing process for quantum sources. We verify the indistinguishability of the on-chip sources by a series of heralded two-source Hong-Ou-Mandel interference, with all the dip visibilities above 90%. In addition, the brightness of the sources is found easily reaching MHz and being applicable to both discrete-variable and continuous-variable platform, showing either clear anti-bunching feature or large squeezing parameter under different pumping regimes. The demonstrated scalability and uniformity of quantum sources, together with integrated photonic network and detection, will enable large-scale all-on-chip quantum processors for real-life applications.


title: Probing the nature of dark matter with accreted globular cluster streams

abstract: The steepness of the central density profiles of dark matter (DM) in low-mass galaxy halos (e.g. dwarf galaxies) is a powerful probe of the nature of DM. We propose a novel scheme to probe the inner profiles of galaxy subhalos using stellar streams. We show that the present day morphological and dynamical properties of accreted globular cluster (GC) streams - those produced from tidal stripping of GCs that initially evolved within satellite galaxies and later merged with the Milky Way (MW) - are sensitive to the central DM density profile and mass of their parent satellites. GCs that accrete within cuspy CDM subhalos produce streams that are physically wider and dynamically hotter than streams that accrete inside cored subhalos. A first comparison of MW streams "GD-1" and "Jhelum" (likely of accreted GC origin) with our simulations indicates a preference for cored subhalos. If these results hold up in future data, the implication is that either the DM cusps were erased by baryonic feedback, or their subhalos naturally possessed cored density profiles implying DM models beyond CDM. Moreover, accreted GC streams are highly structured and exhibit complex morphological features (e.g., parallel structures and "spurs"). This implies that the accretion scenario can naturally explain the recently observed peculiarities in some of the MW streams. We also propose a novel mechanism for forming "gaps" in streams when the remnant of the parent subhalo later passes through the stream. This encounter can last a longer time (and have more of an impact) than the random encounters with DM subhalos previously considered, because the GC stream and its parent subhalo are on similar orbits with small relative velocities. Current and future surveys of the MW halo will uncover numerous faint stellar streams and provide the data needed to substantiate our preliminary tests with this new probe of DM.


title: Illuminating a tadpole's metamorphosis II: observing the on-going transformation with ALMA

abstract: We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the tadpole, a small globule in the Carina Nebula that hosts the HH 900 jet+outflow system. Our data include $^{12}$CO, $^{13}$CO, C$^{18}$O J=2-1, $^{13}$CO, C$^{18}$O J=3-2, and serendipitous detections of DCN J=3-2 and CS J=7-6. With angular resolution comparable to the Hubble Space Telescope (HST), our data reveal for the first time the bipolar molecular outflow in CO, seen only inside the globule, that is launched from the previously unseen jet-driving protostar (the HH 900 YSO). The biconical morphology joins smoothly with the externally irradiated outflow seen in ionized gas tracers outside the globule, tracing the overall morphology of a jet-driven molecular outflow. Continuum emission at the location of the HH 900 YSO appears to be slightly flattened perpendicular to outflow axis. Model fits to the continuum have a best-fit spectral index of $\sim 2$, suggesting cold dust and the onset of grain growth. In position-velocity space, $^{13}$CO and C$^{18}$O gas kinematics trace a C-shaped morphology, similar to infall profiles seen in other sources, although the global dynamical behaviour of the gas remains unclear. Line profiles of the CO isotopologues display features consistent with externally heated gas. We estimate a globule mass of $\sim 1.9$ M$_{\odot}$, indicating a remaining lifetime of $\sim 4$ Myr, assuming a constant photoevaporation rate. This long globule lifetime will shield the disk from external irradiation perhaps prolonging its life and enabling planet formation in regions where disks are typically rapidly destroyed.


title: Outflows, cores and magnetic field orientations in W43-MM1 as seen by ALMA

abstract: It has been proposed that the magnetic field, pervasive in the ISM, plays an important role in the process of massive star formation. To better understand its impact at the pre and protostellar stages, high-angular resolution observations of polarized dust emission toward a large sample of massive dense cores are needed. To this end, we used the Atacama Large Millimeter Array in Band 6 (1.3 mm) in full polarization mode to map the polarized emission from dust grains at a physical scale of $\sim$2700 au in the massive protocluster W43-MM1. We used these data to measure the orientation of the magnetic field at the core scale. Then, we examined the relative orientations of the core-scale magnetic field, of the protostellar outflows determined from CO molecular line emission, and of the major axis of the dense cores determined from 2D Gaussian fit in the continuum emission. We found that the orientation of the dense cores is not random with respect to the magnetic field. Instead, the dense cores are compatible with being oriented 20-50$^\deg$ with respect to the magnetic field. The outflows could be oriented 50-70$^\deg$ with respect to the magnetic field, or randomly oriented with respect to the magnetic field, similar to current results in low-mass star-forming regions. In conclusion, the observed alignment of the position angle of the cores with respect to the magnetic field lines shows that the magnetic field is well coupled with the dense material; however, the 20-50$^\deg$ preferential orientation contradicts the predictions of the magnetically-controlled core-collapse models. The potential correlation of the outflow directions with respect to the magnetic field suggests that, in some cases, the magnetic field is strong enough to control the angular momentum distribution from the core scale down to the inner part of the circumstellar disks where outflows are triggered.


title: SN2019dge: a Helium-rich Ultra-Stripped Envelope Supernova

abstract: We present observations of ZTF18abfcmjw (SN2019dge), a helium-rich supernova with a fast-evolving light curve indicating an extremely low ejecta mass ($\approx 0.3\,M_\odot$) and low kinetic energy ($\approx 1.2\times 10^{50}\,{\rm erg}$). Early-time (<4 d after explosion) photometry reveal evidence of shock cooling from an extended helium-rich envelope of $\sim0.1\,M_\odot$ located at $\sim 3\times 10^{12}\,{\rm cm}$ from the progenitor. Early-time He II line emission and subsequent spectra show signatures of interaction with helium-rich circumstellar material, which extends from $\gtrsim 5\times 10^{13}\,{\rm cm}$ to $\gtrsim 2\times 10^{16}\,{\rm cm}$. We interpret SN2019dge as a helium-rich supernova from an ultra-stripped progenitor, which originates from a close binary system consisting of a mass-losing helium star and a low-mass main sequence star or a compact object (i.e., a white dwarf, a neutron star, or a black hole). We infer that the local volumetric birth rate of 19dge-like ultra-stripped SNe is in the range of 1400--8200$\,{\rm Gpc^{-3}\, yr^{-1}}$ (i.e., 2--12% of core-collapse supernova rate). This can be compared to the observed coalescence rate of compact neutron star binaries that are not formed by dynamical capture.


title: General relativistic non-ideal fluid equations for dark matter from a truncated cumulant expansion

abstract: A new truncation scheme based on the cumulant expansion of the one-particle phase-space distribution function for dark matter particles is developed. Extending the method of moments in relativistic kinetic theory, we derive evolution equations which supplement the covariant conservation of the energy-momentum tensor and particle number current. Truncating the cumulant expansion we obtain a closed, covariant and hyperbolic system of equations which can be used to model the evolution of a general relativistic non-ideal fluid. As a working example we consider a Friedmann-Lema\^itre-Robertson-Walker cosmology with dynamic pressure and solve for the time evolution of the effective equation of state parameter.


title: Gapped Goldstones at the cut-off scale: a non-relativistic EFT

abstract: At finite density, the spontaneous breakdown of an internal non-Abelian symmetry dictates, along with gapless modes, modes whose gap is fixed by the algebra and proportional to the chemical potential: the gapped Goldstones. Generically the gap of these states is comparable to that of other non-universal excitations or to the energy scale where the dynamics is strongly coupled. This makes it non-straightforward to derive a universal ef effective field theory (EFT) description realizing all the symmetries. Focusing on the illustrative example of a fully broken SU (2) group, we demonstrate that such an EFT can be constructed by carving out around the Goldstones, gapless and gapped, at small 3-momentum. The rules governing the EFT, where the gapless Goldstones are soft while the gapped ones are slow, are those of standard nonrelativistic EFTs, like for instance nonrelativistic QED. In particular, the EFT Lagrangian formally preserves gapped Goldstone number, and processes where such number is not conserved are described inclusively by allowing for imaginary parts in the Wilson coefficients. Thus, while the symmetry is manifestly realized in the EFT, unitarity is not. We comment on the application of our construction to the study of the large charge sector of conformal field theories with non-Abelian symmetries.


title: Electrical detection of spin liquids in double moir\'e layers

abstract: Although spin is a fundamental quantum number, measuring spin transport in traditional solid state systems is extremely challenging. This poses a major obstacle to detecting interesting quantum states including certain spin liquids. In this paper we propose a platform that not only allows for the electrical measurement of spin transport, but in which a variety of exotic quantum phases may be stabilized. Our proposal involves two moir\'e superlattices, built from transition metal dichalcogenides (TMD) or graphene, separated from one another by a thin insulating layer. The two Coulomb coupled moir\'e layers, when suitably aligned, give rise to a layer pseudospin degree of freedom. The transport of pseudospin can be accessed from purely electrical measurements of counter-flow or Coulomb drag conductivity. Furthermore, these platforms naturally realize Hubbard models on the triangular lattice with $N = 4\, {\rm or}\, 8$ flavors. The flavor degeneracy motivates a large-N approximation from which we obtain the phase diagram of Mott insulators at different electron fillings and correlation strengths. In addition to conventional phases such as psuedospin superfluids and crystallized insulators, exotic phases including chiral spin liquids and a $U(1)$ spinon Fermi surface spin liquid are also found, all of which will show smoking gun electrical signatures in this setup.


title: The origin of X-ray coronae around simulated disc galaxies

abstract: The existence of hot, accreted gaseous coronae around massive galaxies is a long-standing central prediction of galaxy formation models in the $\Lambda$CDM cosmology. While observations now confirm that extraplanar hot gas is present around late-type galaxies, the origin of the gas is uncertain with suggestions that galactic feedback could be the dominant source of energy powering the emission. We investigate the origin and X-ray properties of the hot gas that surrounds galaxies of halo mass, $(10^{11}-10^{14}) \mathrm{M}_\odot$, in the cosmological hydrodynamical EAGLE simulations. We find that the central X-ray emission, $\leq 0.10 R_{\mathrm{vir}}$, of halos of mass $\leq 10^{13} \mathrm{M}_\odot$ originates from gas heated by supernovae (SNe). However, beyond this region, a quasi-hydrostatic, accreted atmosphere dominates the X-ray emission in halos of mass $\geq 10^{12} \mathrm{M}_\odot$. We predict that a dependence on halo mass of the hot gas to dark matter mass fraction can significantly change the slope of the $L_{\mathrm{X}}-M_{\mathrm{vir}}$ relation (which is typically assumed to be $4/3$ for clusters) and we derive the scaling law appropriate to this case. As the gas fraction in halos increases with halo mass, we find a steeper slope for the $L_{\mathrm{X}}-M_{\mathrm{vir}}$ in lower mass halos, $\leq 10^{14} \mathrm{M}_\odot$. This varying gas fraction is driven by active galactic nuclei (AGN) feedback. We also identify the physical origin of the so-called "missing feedback" problem, the apparently low X-ray luminosities observed from high star-forming, low-mass galaxies. This is explained by the ejection of SNe-heated gas from the central regions of the halo.


title: Application of X-Ray Clumpy Torus Model (XCLUMPY) to 10 Obscured Active Galactic Nuclei Observed with Suzaku and NuSTAR

abstract: We apply XCLUMPY, an X-ray spectral model from a clumpy torus in an active galactic nucleus (AGN), to the broadband X-ray spectra of 10 obscured AGNs observed with both Suzaku and NuSTAR. The infrared spectra of these AGNs were analyzed with the CLUMPY code. Since XCLUMPY adopts the same clump distribution as that in the CLUMPY, we can directly compare the torus parameters obtained from the X-ray spectra and those from the infrared ones. The torus angular widths determined from the infrared spectra ($\sigma_{\mathrm{IR}}$) are systematically larger than those from the X-ray data ($\sigma_{\mathrm{X}}$); the difference ($\sigma_{\mathrm{IR}}-\sigma_{\mathrm{X}}$) correlates with the inclination angle determined from the X-ray spectrum. These results can be explained by the contribution from dusty polar outflows to the observed infrared flux, which becomes more significant at higher inclinations (more edge-on views). The ratio of the hydrogen column density and V-band extinction in the line of sight absorber shows large scatter ($\simeq$1 dex) around the Galactic value, suggesting that a significant fraction of AGNs have dust-rich circumnuclear environments.


title: A hierarchical field-level inference approach to reconstruction from sparse Lyman-$\alpha$ forest data

abstract: We address the problem of inferring the three-dimensional matter distribution from a sparse set of one-dimensional quasar absorption spectra of the Lyman-$\alpha$ forest. Using a Bayesian forward modelling approach, we focus on extending the dynamical model to a fully self-consistent hierarchical field-level prediction of redshift-space quasar absorption sightlines. Our field-level approach rests on a recently developed semiclassical analogue to Lagrangian perturbation theory (LPT), which improves over noise problems and interpolation requirements of LPT. It furthermore allows for a manifestly conservative mapping of the optical depth to redshift space. In addition, this new dynamical model naturally introduces a coarse-graining scale, which we exploited to accelerate the Markov chain Monte-Carlo (MCMC) sampler using simulated annealing. By gradually reducing the effective temperature of the forward model, we were able to allow it to first converge on large spatial scales before the sampler became sensitive to the increasingly larger space of smaller scales. We demonstrate the advantages, in terms of speed and noise properties, of this field-level approach over using LPT as a forward model, and, using mock data, we validated its performance to reconstruct three-dimensional primordial perturbations and matter distribution from sparse quasar sightlines.


title: Non-Simply-Connected Symmetries in 6D SCFTs

abstract: Six-dimensional N=(1,0) superconformal field theories can be engineered geometrically via F-theory on elliptically-fibered Calabi-Yau 3-folds. We include torsional sections in the geometry, which lead to a finite Mordell-Weil group. This allows us to identify the full non-Abelian group structure rather than just the algebra. The presence of torsion also modifies the center of the symmetry groups and the matter representations that can appear. This in turn affects the tensor branch of these theories. We analyze this change for a large class of superconformal theories with torsion and explicitly construct their tensor branches. Finally, we elaborate on the connection to the dual heterotic and M-theory description, in which our configurations are interpreted as generalizations of discrete holonomy instantons.


title: New de Sitter solutions of 10d type IIB supergravity

abstract: We find and study 17 new de Sitter solutions of ten-dimensional (10d) type IIB supergravity with intersecting D5-branes and orientifold O5-planes, as well as a new Minkowski one. These solutions are obtained numerically on 6d group manifolds, the compactness of which is established for 4 of them. We show that all our de Sitter solutions are perturbatively unstable, using a restricted 4d effective theory of four scalar fields. We finally analyse whether our solutions can be promoted to classical string backgrounds. Several of them appear as good candidates, as they satisfy all requirements imposed so far.


title: Determining the systemic redshift of Lyman-alpha emitters with neural networks and improving the measured large-scale clustering

abstract: We explore how to mitigate the clustering distortions in Lyman-$\alpha$ emitters (LAEs) samples caused by the miss-identification of the Lyman-$\alpha$ (Ly$\alpha$) wavelength in their Ly$\alpha$ line profiles. We use the Ly$\alpha$ line profiles from our previous LAE theoretical model that includes radiative transfer in the interstellar and intergalactic mediums. We introduce a novel approach to measure the systemic redshift of LAEs from their Ly$\alpha$ line using neural networks. In detail, we assume that, for a fraction of the whole LAE population their systemic redshift is determined precisely through other spectral features. We then use this subset to train a neural network that predicts the Ly$\alpha$ wavelength given a Ly$\alpha$ line profile. We test two different training sets: i) the LAEs are selected homogeneously and ii) only the brightest LAEs are selected. In comparison with previous approaches in the literature, our methodology improves significantly both accuracy and precision in determining the Ly$\alpha$ wavelength. In fact, after applying our algorithm in ideal Ly$\alpha$ line profiles, we recover the clustering unperturbed down to 1cMpc/h. Then, we test the performance of our methodology in realistic Ly$\alpha$ line profiles by downgrading their quality. The machine learning techniques work well even if the Ly$\alpha$ line profile quality is decreased considerably. We conclude that LAE surveys such as HETDEX would benefit from determining with high accuracy the systemic redshift of a subpopulation and applying our methodology to estimate the systemic redshift of the rest of the galaxy sample.


title: Evidence that 1I/2017 U1 (`Oumuamua) was composed of molecular hydrogen ice

abstract: `Oumuamua (I1 2017) was the first macroscopic ($l\sim100\,{\rm m}$) body observed to traverse the inner solar system on an unbound hyperbolic orbit. Its light curve displayed strong periodic variation, and it showed no hint of a coma or emission from molecular outgassing. Astrometric measurements indicate that 'Oumuamua experienced non-gravitational acceleration on its outbound trajectory, but energy balance arguments indicate this acceleration is inconsistent with a water ice sublimation-driven jet of the type exhibited by solar system comets. We show that all of `Oumaumua's observed properties can be explained if it contained a significant fraction of molecular hydrogen (H$_{2}$) ice. H$_{2}$ sublimation at a rate proportional to the incident solar flux generates a surface-covering jet that reproduces the observed acceleration. Mass wasting from sublimation leads to monotonic increase in the body axis ratio, explaining `Oumuamua's shape. Back-tracing `Oumuamua's trajectory through the Solar System permits calculation of its mass and aspect ratio prior to encountering the Sun. We show that H$_{2}$-rich bodies plausibly form in the coldest dense cores of Giant Molecular Clouds, where number densities are of order $n\sim10^5$, and temperatures approach the $T=3\,{\rm K}$ background. Post-formation exposure to galactic cosmic rays implies a $\tau \sim 100$ Myr age, explaining the kinematics of `Oumuamua's inbound trajectory.


title: Connecting the structure of dark matter haloes to the primordial power spectrum

abstract: A large body of work based on collisionless cosmological N-body simulations going back over two decades has advanced the idea that collapsed dark matter haloes have simple and approximately universal forms for their mass density and pseudo-phase space density (PPSD) distributions. However, a general consensus on the physical origin of these results has not yet been reached. In the present study, we explore to what extent the apparent universality of these forms holds when we vary the initial conditions (i.e., the primordial power spectrum of density fluctuations) away from the standard CMB-normalised case, but still within the context of LCDM with a fixed expansion history. Using simulations that vary the initial amplitude and shape, we show that the structure of dark matter haloes retains a clear memory of the initial conditions. Specifically, increasing (lowering) the amplitude of fluctuations increases (decreases) the concentration of haloes and, if pushed far enough, the density profiles deviate strongly from the NFW form that is a good approximation for the CMB-normalised case. Although, an Einasto form works well. Rather than being universal, the slope of the PPSD (or pseudo-entropy) profile steepens (flattens) with increasing (decreasing) power spectrum amplitude and can exhibit a strong halo mass dependence. Our results therefore indicate that the previously identified universality of the structure of dark matter haloes is mostly a consequence of adopting a narrow range of (CMB-normalised) initial conditions for the simulations. Our new suite provides a useful test-bench against which physical models for the origin of halo structure can be validated.