title: Almost positive kernels on compact Riemannian manifolds

abstract: We show how to build a kernel \[ K_X(x,y)=\sum_{m=0}^Xh(\lambda_m/{\lambda_X})\varphi_m(x)\overline{\varphi_m(y)} \] on a compact Riemannian manifold $M$, which is positive up to a negligible error and such that $K_X(x,x)\approx X$. Here $0=\lambda_0^2\le\lambda_1^2\le\ldots$ are the eigenvalues of the Laplace-Beltrami operator on $M$, listed with repetitions, and $\varphi_0,\,\varphi_1,\ldots$ an associated system of eigenfunctions, forming an orthonormal basis of $L^2(M)$. The function $h$ is smooth up to a certain minimal degree, even, compactly supported in $[-1,1]$ with $h(0)=1$, and $K_X(x,y)$ turns out to be an approximation to the identity.


title: Ab initio simulation of amorphous graphite

abstract: A disorder-order transition from amorphous carbon (aC) to layered amorphous graphene (LAG) has been predicted using ab-initio methods. Amorphous carbon at densities close to the graphitic density show a strong proclivity to layer in NVT simulations near 3000K . Our calculations have shown that graphitization is strongly dependent on the simulation temperature and is not observed under 2500K. The origin of such a disorder-order transition has been studied using various structural analyses. Each layer of LAG is a layer of amorphous graphene (aG) including pentagons and heptagons in addition to hexagons and the planes are separated by about 3.2A. LAG obtained from the NVT simulations were highly stable. The electronic charge density, especially for the pi and pi$^*$ states, was computed with the Hyed-Scuseria-Ernzerhof (HSE) hybrid functional and compared with crystalline graphite (cG). The impact of structural disorder has been studied in detail, especially the consequences of disorder to electronic transport.


title: Tangent spaces and T-invariant curves of Schubert varieties

abstract: The set of T-invariant curves in a Schubert variety through a T-fixed point is relatively easy to characterize in terms of its weights, but the tangent space is more difficult. We prove that the weights of the tangent space are contained in the rational cone generated by the weights of the T-invariant curves. In simply laced types, this remains true if "rational" is replaced by "integral". We also obtain conditions under which every weight of the tangent space is the weight of a T-invariant curve, as well as a smoothness criterion. The results rely on equivariant K-theory, as well as the study of different notions of decomposability of roots.


title: Diffeological Principal Bundles and Principal Infinity Bundles

abstract: In this paper, we study diffeological spaces as certain kinds of discrete simplicial presheaves on the site of cartesian spaces with the coverage of good open covers. The \v{C}ech model structure on simplicial presheaves provides us with a notion of $\infty$-stack cohomology of a diffeological space with values in a diffeological abelian group $A$. We compare $\infty$-stack cohomology of diffeological spaces with two existing notions of \v{C}ech cohomology for diffeological spaces in the literature. Finally, we prove that for a diffeological group $G$, that the nerve of the category of diffeological principal $G$-bundles is weak homotopy equivalent to the nerve of the category of $G$-principal $\infty$-bundles on $X$, bridging the bundle theory of diffeology and higher topos theory.


title: Broad search for gravitational waves from subsolar-mass binaries through LIGO and Virgo's third observing run

abstract: We present a search for gravitational waves from the coalescence of binaries which contain at least one subsolar mass component using data from the LIGO and Virgo observatories through the completion of their third observing run. The observation of a merger with a component below $1\,M_{\odot}$ would be a clear sign of either new physics or the existence of a primordial black hole population; these black holes could also contribute to the dark matter distribution. Our search targets binaries where the primary has mass $M_1$ between 0.1-100$\,M_{\odot}$ and the secondary has mass $M_2$ from 0.1-1$\,M_{\odot}$ for $M_1 < 20\,M_\odot$ and 0.01-1$\,M_{\odot}$ for $M_1 \ge 20\,M_\odot$. Sources with $M_1 < 7\,M_\odot, M_2 > 0.5\,M_\odot$ are also allowed to have orbital eccentricity up to $e_{10} \sim 0.3$. This search region covers from comparable to extreme mass ratio sources up to $10^4:1$. We find no statistically convincing candidates and so place new upper limits on the rate of mergers; our analysis sets the first limits for most subsolar sources with $7\,M_{\odot}< M_1 < 20\,M_{\odot}$ and tightens limits by $\sim 8\times$ $(1.6\times)$ where $M_1 > 20\,M_{\odot}$ ($M_1 < 7\,M_{\odot}$). Using these limits, we constrain the dark matter fraction to below 0.3 (0.7)$\%$ for 1 (0.5)$\,M_{\odot}$ black holes assuming a monochromatic mass function. Due to the high merger rate of primordial black holes beyond the individual source horizon distance, we also use the lack of an observed stochastic background as a complementary probe to limit the dark matter fraction. We find that although the limits are in general weaker than those from the direct search they become comparable at $0.1 \,M_{\odot}$.


title: 3D ToF LiDAR in Mobile Robotics: A Review

abstract: In the past ten years, the use of 3D Time-of-Flight (ToF) LiDARs in mobile robotics has grown rapidly. Based on our accumulation of relevant research, this article systematically reviews and analyzes the use 3D ToF LiDARs in research and industrial applications. The former includes object detection, robot localization, long-term autonomy, LiDAR data processing under adverse weather conditions, and sensor fusion. The latter encompasses service robots, assisted and autonomous driving, and recent applications performed in response to public health crises. We hope that our efforts can effectively provide readers with relevant references and promote the deployment of existing mature technologies in real-world systems.


title: Anharmonic quantum thermal transport across a van der Waals interface

abstract: We investigate the anharmonic phonon scattering across a weakly interacting interface by developing a quantum mechanics-based theory. We find that the contribution from anharmonic three-phonon scatterings to interfacial thermal conductance can be cast into Landauer formula with transmission function being temperature-dependent. Surprisingly, in the weak coupling limit, the transmission due to anharmonic phonon scattering is unbounded with increasing temperature, which is physically impossible for two-phonon processes. We further reveal that the anharmonic contribution in a real heterogeneous interface (e.g., between graphene and monolayer molybdenum disulfide) can dominate over the harmonic process even at room temperature, highlighting the important role of anharmonicity in weakly interacting heterogeneous systems.


title: An Overview on Machine Translation Evaluation

abstract: Since the 1950s, machine translation (MT) has become one of the important tasks of AI and development, and has experienced several different periods and stages of development, including rule-based methods, statistical methods, and recently proposed neural network-based learning methods. Accompanying these staged leaps is the evaluation research and development of MT, especially the important role of evaluation methods in statistical translation and neural translation research. The evaluation task of MT is not only to evaluate the quality of machine translation, but also to give timely feedback to machine translation researchers on the problems existing in machine translation itself, how to improve and how to optimise. In some practical application fields, such as in the absence of reference translations, the quality estimation of machine translation plays an important role as an indicator to reveal the credibility of automatically translated target languages. This report mainly includes the following contents: a brief history of machine translation evaluation (MTE), the classification of research methods on MTE, and the the cutting-edge progress, including human evaluation, automatic evaluation, and evaluation of evaluation methods (meta-evaluation). Manual evaluation and automatic evaluation include reference-translation based and reference-translation independent participation; automatic evaluation methods include traditional n-gram string matching, models applying syntax and semantics, and deep learning models; evaluation of evaluation methods includes estimating the credibility of human evaluations, the reliability of the automatic evaluation, the reliability of the test set, etc. Advances in cutting-edge evaluation methods include task-based evaluation, using pre-trained language models based on big data, and lightweight optimisation models using distillation techniques.


title: Generating Synthetic Mobility Networks with Generative Adversarial Networks

abstract: The increasingly crucial role of human displacements in complex societal phenomena, such as traffic congestion, segregation, and the diffusion of epidemics, is attracting the interest of scientists from several disciplines. In this article, we address mobility network generation, i.e., generating a city's entire mobility network, a weighted directed graph in which nodes are geographic locations and weighted edges represent people's movements between those locations, thus describing the entire mobility set flows within a city. Our solution is MoGAN, a model based on Generative Adversarial Networks (GANs) to generate realistic mobility networks. We conduct extensive experiments on public datasets of bike and taxi rides to show that MoGAN outperforms the classical Gravity and Radiation models regarding the realism of the generated networks. Our model can be used for data augmentation and performing simulations and what-if analysis.


title: In-situ Measurement of the Energy Fraction in Supra-thermal and Energetic Particles at ACE, Wind, and PSP Interplanetary Shocks

abstract: The acceleration of charged particles by interplanetary shocks (IPs) can drain a non-negligible fraction of the plasma pressure. In this study, we have selected 17 IPs observed in-situ at $1\,\text{au}$ by the Advanced Composition Explorer (ACE) and the Wind spacecraft, and 1 shock at $0.8\,\text{au}$ observed by Parker Solar Probe (PSP). We have calculated the time-dependent partial pressure of supra-thermal and energetic particles (smaller and greater than $50\,\text{keV}$ for protons and $30\,\text{keV}$ for electrons, respectively) in both the upstream and downstream regions. The particle fluxes were averaged for 1 hour before and 1 hour after the shock time to remove short time scale effects. Using the MHD Rankine-Hugoniot jump conditions, we find that the fraction of the total upstream energy flux transferred to supra-thermal and energetic downstream particles is typically $\lesssim\!16\%$, in agreement with previous observations and simulations. Notably, by accounting for errors on all measured shock parameters, we have found that for any given fast magnetosonic Mach number, $M_{f}\!<7$, the angle between the shock normal and average upstream magnetic field, $\theta_{Bn}$, is not correlated with the energetic particle pressure; in particular, the partial pressure of energized particles does not decrease for $\theta_{Bn} \gtrsim 45^\circ$. The downstream electron-to-proton energy ratio in the range $\gtrsim\!140\,\text{eV}$ for electrons and $\gtrsim\!70\,\text{keV}$ for protons exceeds the expected $\sim\!1\%$ and nears equipartition ($>\!0.1$) for the Wind events.

categories: astro-ph.SR astro-ph.HE physics.plasm-ph physics.space-ph


title: Birational classification of pointless del Pezzo surfaces of degree 8

abstract: Let k be a perfect field. Recently J.-L. Colliot-Th\'el\`ene showed that two pointless quadric surfaces over k are birationally equivalent if and only if they are isomorphic. We show that this result holds for arbitrary del Pezzo surfaces of degree $8$ with the Picard number $1$, and describe minimal surfaces birationally equivalent to a given pointless del Pezzo surface of degree $8$.


title: A Unified Nonparametric Test of Transformations on Distribution Functions with Nuisance Parameters

abstract: This paper proposes a simple unified approach to testing transformations on cumulative distribution functions (CDFs) in the presence of nuisance parameters. The proposed test is constructed based on a new characterization that avoids the estimation of nuisance parameters. The critical values are obtained through a numerical bootstrap method which can easily be implemented in practice. Under suitable conditions, the proposed test is shown to be asymptotically size controlled and consistent. The local power property of the test is established. Finally, Monte Carlo simulations and an empirical study show that the test performs well on finite samples.


title: Classification of Planar Monomials Over Finite Fields of Small Order

abstract: For all finite fields of order up to $2^{30}$, we computationally prove that there are no planar monomials besides the ones already known.


title: Characterizing generalized axisymmetric quantum states in $d\times d$ systems

abstract: We introduce a family of highly symmetric bipartite quantum states in arbitrary dimensions. It consists of all states that are invariant under local phase rotations and local cyclic permutations of the basis. We solve the separability problem for a subspace of these states and show that a sizable part of the family is bound entangled. We also calculate some of the Schmidt numbers for the family in $d = 3$, thereby characterizing the dimensionality of entanglement. Our results allow us to estimate entanglement properties of arbitrary states, as general states can be symmetrized to the considered family by local operations.


title: Some minimal bimolecular mass-action systems with limit cycles

abstract: We discuss three examples of bimolecular mass-action systems with three species, due to Feinberg, Berner, Heinrich, and Wilhelm. Each system has a unique positive equilibrium which is unstable for certain rate constants and then exhibits stable limit cycles, but no chaotic behaviour. For some rate constants in the Feinberg--Berner system, a stable equilibrium, an unstabe limit cycle, and a stable limit cycle coexist. All three networks are minimal in some sense. By way of homogenising the above three examples, we construct bimolecular mass-conserving mass-action systems with four species that admit a stable limit cycle. The homogenised Feinberg--Berner system and the homogenised Wilhelm--Heinrich system admit the coexistence of a stable equilibrium, an unstable limit cycle, and a stable limit cycle.


title: Fast geostatistical inference under positional uncertainty: Analysing DHS household survey data

abstract: Household survey data from the Demographic and Health Surveys (DHS) Program is published with GPS coordinates. However, almost all geostatistical analyses of such data ignore that the published GPS coordinates are randomly displaced (jittered). In this short report, we develop a geostatistical model that accounts for the positional uncertainty when analysing DHS surveys, and provide a fast implementation using Template Model Builder. The key focus is inference with Gaussian random fields under positional uncertainty, and our approach works for both Gaussian and non-Gaussian likelihoods. A simulation study with a binomial observation model shows that the new approach performs equally or better than the common approach of ignoring jittering, both in terms of more accurate parameter estimates and improved predictive measures. We demonstrate that the improvement would be larger under stronger jittering. An analysis of contraceptive use in Kenya shows that the approach is fast and easy to use in practice.


title: Gradient-type estimates for the dynamic $\varphi^4_2$-model

abstract: We prove gradient bounds for the Markov semigroup of the dynamic $\varphi^4_2$-model on a torus of fixed size $L>0$. For sufficiently large mass $m>0$ these estimates imply exponential contraction of the Markov semigroup. Our method is based on pathwise estimates of the linearized equation. To compensate the lack of exponential integrability of the stochastic drivers we use a stopping time argument and the strong Markov property in the spirit of Cass--Litterer--Lyons. Following the classical approach of Bakry-\'Emery, as a corollary we prove a Poincar\'e/spectral gap inequality for the $\varphi^4_2$-measure of sufficiently large mass $m>0$ with almost optimal carr\'e du champ.


title: Effects of high order interatomic potential on elastic phonon scatterings

abstract: Interatomic potentials beyond quadratic order provide scattering sources for phonon transport in lattice. By using a weakly-interacting interface model, we investigated the relation between the order of interatomic potential and the multiple-phonon scattering process. We find that high order interatomic potential not only causes multiple-phonon scattering processes, but also has significant impacts on elastic phonon scattering processes. Using fourth order potential as an example, we show that it can significantly affects elastic phonon scatterings, through the formation of localized phonons. Such impact is closely related to the correlations of interfacial atoms and it becomes more significant with increasing temperature. Our work suggests that it is insufficient to consider only quadratic potential to investigate elastic phonon transport.


title: Banding vs. Quality: Perceptual Impact and Objective Assessment

abstract: Staircase-like contours introduced to a video by quantization in flat areas, commonly known as banding, have been a long-standing problem in both video processing and quality assessment communities. The fact that even a relatively small change of the original pixel values can result in a strong impact on perceived quality makes banding especially difficult to be detected by objective quality metrics. In this paper, we study how banding annoyance compares to more commonly studied scaling and compression artifacts with respect to the overall perceptual quality. We further propose a simple combination of VMAF and the recently developed banding index, CAMBI, into a banding-aware video quality metric showing improved correlation with overall perceived quality.


title: Galaxy Rotation Curve Measurements with Low Cost 21 cm Radio Telescope

abstract: Probing the Universe with atomic hydrogen 21 cm emission is a fascinating and challenging work in astronomy. Radio telescopes play a vital role in detecting and imaging these faint signals. Powerful radio telescopes are complex to construct and operate. We have built a simple, low-cost 21 cm radio telescope primarily for educational training purposes. The design uses a custom horn antenna, ready-to-use radio-frequency components, and a software-defined radio module. The telescope operates efficiently from a rooftop in a city environment. Using this telescope, we have conducted observations and successfully detected the 21 cm line emissions from the different directions of our galactic plane. Based on the Doppler-shift observed in these measurements, we have successfully derived the Galactic rotation velocity (rotation curve) in those directions. The paper presents the details of the telescope construction, 21 cm observation, and the Galactic rotation curve derivation.


title: Hubble Space Telescope Imaging of Luminous Extragalactic Infrared Transients and Variables from the SPIRITS Survey

abstract: SPIRITS--the SPitzer InfraRed Intensive Transients Survey--searched for luminous infrared (IR) transients and variables in nearly 200 nearby galaxies from 2014 to 2019, using the warm Spitzer telescope at 3.6 and 4.5 microns. Among the SPIRITS variables are IR-bright objects that are undetected in ground-based optical surveys. We classify them as (1) transients, (2) periodic variables, and (3) irregular variables. The transients include "SPRITE"s (eSPecially Red Intermediate-luminosity Transient Events), having maximum luminosities fainter than supernovae, red IR colors, and a wide range of outburst durations (days to years). Here we report deep optical and near-IR imaging with the Hubble Space Telescope (HST) of 21 SPIRITS variables. They were initially considered SPRITE transients, but many eventually proved instead to be periodic or irregular variables as more data were collected. HST images show most of these cool and dusty variables are associated with star-forming regions in late-type galaxies, implying an origin in massive stars. Two SPRITEs lacked optical progenitors in deep pre-outburst HST images; however, one was detected during eruption at J and H, indicating a dusty object with an effective temperature of ~1050 K. One faint SPRITE turned out to be a dusty classical nova. About half the HST targets proved to be periodic variables, with pulsation periods of 670-2160 days; they are likely dusty asymptotic-giant-branch (AGB) stars with masses of ~5-10 Msun. A few of them were warm enough to be detected in deep HST frames, but most are too cool. Out of six irregular variables, two were red supergiants with optical counterparts in HST images; four were too enshrouded for HST detection.


title: Validating First-Principles Phonon Lifetimes via Inelastic Neutron Scattering

abstract: Phonon lifetimes are a key component of quasiparticle theories of transport, yet first-principles lifetimes are rarely directly compared to inelastic neutron scattering (INS) results. Existing comparisons show discrepancies even at temperatures where perturbation theory is expected to be reliable. In this work, we demonstrate that the reciprocal space voxel ($q$-voxel), which is the finite region in reciprocal space required in INS data analysis, must be explicitly accounted for within theory in order to draw a meaningful comparison. We demonstrate accurate predictions of peak widths of the scattering function when accounting for the $q$-voxel in CaF$_2$ and ThO$_2$. Passing this test implies high fidelity of the phonon interactions and the approximations used to compute the Green's function, serving as critical benchmark of theory, and indicating that other material properties should be accurately predicted; which we demonstrate for thermal conductivity.


title: FASURA: A Scheme for Quasi-Static Massive MIMO Unsourced Random Access Channels

abstract: This article considers the massive MIMO unsourced random access problem on a quasi-static Rayleigh fading channel. Given a fixed message length and a prescribed number of channel uses, the objective is to construct a coding scheme that minimizes the energy-per-bit subject to a fixed probability of error. The proposed scheme differs from other state-of-the-art schemes in that it blends activity detection, single-user coding, pilot-aided and temporary decisions-aided iterative channel estimation and decoding, minimum-mean squared error (MMSE) estimation, and successive interference cancellation (SIC). We show that an appropriate combination of these ideas can substantially outperform state-of-the-art coding schemes when the number of active users is more than 100, making this the best performing scheme known for this regime.


title: Differentially Private Estimation of Heterogeneous Causal Effects

abstract: Estimating heterogeneous treatment effects in domains such as healthcare or social science often involves sensitive data where protecting privacy is important. We introduce a general meta-algorithm for estimating conditional average treatment effects (CATE) with differential privacy (DP) guarantees. Our meta-algorithm can work with simple, single-stage CATE estimators such as S-learner and more complex multi-stage estimators such as DR and R-learner. We perform a tight privacy analysis by taking advantage of sample splitting in our meta-algorithm and the parallel composition property of differential privacy. In this paper, we implement our approach using DP-EBMs as the base learner. DP-EBMs are interpretable, high-accuracy models with privacy guarantees, which allow us to directly observe the impact of DP noise on the learned causal model. Our experiments show that multi-stage CATE estimators incur larger accuracy loss than single-stage CATE or ATE estimators and that most of the accuracy loss from differential privacy is due to an increase in variance, not biased estimates of treatment effects.


title: Nested spheroidal figures of equilibrium I. Approximate solutions for rigid rotations

abstract: We discuss the equilibrium conditions for a body made of two homogeneous components separated by oblate spheroidal surfaces and in relative motion. While exact solutions are not permitted for rigid rotation (unless a specific ambient pressure), approximations can be obtained for configurations involving a small confocal parameter. The problem then admits two families of solutions, depending on the pressure along the common interface (constant or quadratic with the cylindrical radius). We give in both cases the pressure and the rotation rates as a function of the fractional radius, ellipticities and mass-density jump. Various degrees of flattening are allowed but there are severe limitations for global rotation, as already known from classical theory (e.g. impossibility of confocal and coelliptical solutions, gradient of ellipticity outward). States of relative rotation are much less constrained, but these require a mass-density jump. This analytical approach compares successfully with the numerical solutions obtained from the Self-Consistent-Field method. Practical formula are derived in the limit of small ellipticities appropriate for slowly-rotating star/planet interiors.


title: SupermarQ: A Scalable Quantum Benchmark Suite

abstract: The emergence of quantum computers as a new computational paradigm has been accompanied by speculation concerning the scope and timeline of their anticipated revolutionary changes. While quantum computing is still in its infancy, the variety of different architectures used to implement quantum computations make it difficult to reliably measure and compare performance. This problem motivates our introduction of SupermarQ, a scalable, hardware-agnostic quantum benchmark suite which uses application-level metrics to measure performance. SupermarQ is the first attempt to systematically apply techniques from classical benchmarking methodology to the quantum domain. We define a set of feature vectors to quantify coverage, select applications from a variety of domains to ensure the suite is representative of real workloads, and collect benchmark results from the IBM, IonQ, and AQT@LBNL platforms. Looking forward, we envision that quantum benchmarking will encompass a large cross-community effort built on open source, constantly evolving benchmark suites. We introduce SupermarQ as an important step in this direction.


title: A policy gradient approach for optimization of smooth risk measures

abstract: We propose policy gradient algorithms for solving a risk-sensitive reinforcement learning (RL) problem in on-policy as well as off-policy settings. We consider episodic Markov decision processes, and model the risk using the broad class of smooth risk measures of the cumulative discounted reward. We propose two template policy gradient algorithms that optimize a smooth risk measure in on-policy and off-policy RL settings, respectively. We derive non-asymptotic bounds that quantify the rate of convergence of our proposed algorithms to a stationary point of the smooth risk measure. As special cases, we establish that our algorithms apply to optimization of mean-variance and distortion risk measures, respectively.


title: Szeg\"{o}-Weinberger type inequalities for symmetric domains in simply connected space forms

abstract: We consider the Neumann eigenvalue problem for Laplacian on a bounded multi-connected domain contained in simply connected space forms. Under certain symmetry assumptions on the domain, we prove Szeg\"{o}-Weinberger type inequalities for the first $n$ positive Neumann eigenvalues.


title: Listening to the Universe with Next Generation Ground-Based Gravitational-Wave Detectors

abstract: In this study, we use simple performance metrics to assess the science capabilities of future ground-based gravitational-wave detector networks -- composed of A+ or Voyager upgrades to the LIGO, Virgo, and KAGRA observatories and proposed next generation observatories such as Cosmic Explorer and Einstein Telescope. These metrics refer to coalescences of binary neutron stars (BNSs) and binary black holes (BBHs) and include: (i) network detection efficiency and detection rate of cosmological sources as a function of redshift, (ii) signal-to-noise ratios and the accuracy with which intrinsic and extrinsic parameters would be measured, and (iii) enabling multimessenger astronomy with gravitational waves by accurate 3D localization and early warning alerts. We further discuss the science enabled by the small population of rare and extremely loud events. While imminent upgrades will provide impressive advances in all these metrics, next generation observatories will deliver an improvement of an order-of-magnitude or more in most metrics. In fact, a network containing two or three such facilities will detect half of all the BNS and BBH mergers up to a redshift of $z=1$ and $z=20$, respectively, give access to hundreds of BNSs and ten thousand BBHs with signal-to-noise ratios exceeding 100, readily localize hundreds to thousands of mergers to within $1\,{\rm deg^2}$ on the sky and better than 10% in luminosity distance, respectively, and consequently, enable mutlimessenger astronomy through follow-up surveys in the electromagnetic spectrum several times a week. Such networks will further shed light on potential cosmological merger populations and detect an abundance of high-fidelity BNS and BBH signals which will allow investigations of the high-density regime of matter at an unprecedented level and enable precision tests of general relativity in the strong-filed regime, respectively.


title: Wastewater Pipe Condition Rating Model Using K- Nearest Neighbors

abstract: Risk-based assessment in pipe condition mainly focuses on prioritizing the most critical assets by evaluating the risk of pipe failure. This paper's goal is to classify a comprehensive pipe rating model which is obtained based on a series of pipe physical, external, and hydraulic characteristics that are identified for the proposed methodology. The traditional manual method of assessing sewage structural conditions takes a long time. By building an automated process using K-Nearest Neighbors (K-NN), this study presents an effective technique to automate the identification of the pipe defect rating using the pipe repair data. First, we performed the Shapiro Wilks Test for 1240 data from the Dept. of Engineering & Environmental Services, Shreveport, Louisiana Phase 3 with 12 variables to determine if factors could be incorporated in the final rating. We then developed a K-Nearest Neighbors model to classify the final rating from the statistically significant factors identified in Shapiro Wilks Test. This classification process allows recognizing the worst condition of wastewater pipes that need to be replaced immediately. This comprehensive model is built according to the industry-accepted and used guidelines to estimate the overall condition. Finally, for validation purposes, the proposed model is applied to a small portion of a US wastewater collection system in Shreveport, Louisiana. Keywords: Pipe rating, Shapiro Wilks Test, K-Nearest Neighbors (KNN), Failure, Risk analysis


title: Stacked Periodograms as a Probe of Exoplanetary Populations

abstract: Ongoing, extreme-precision Doppler radial velocity surveys seek planets with masses less than several M$_\oplus$; population-level studies to determine the distribution of planetary masses, however, remain difficult due to the required observational time investment, as well as challenges associated with robustly detecting the lowest mass planets. We outline a novel approach that leverages extensive, existing RV datasets to constrain masses of exoplanet populations: stacking periodograms of RV timeseries across many targets. We show that an exoplanet population may be statistically identifiable in the stacked periodogram, even when individual planets do not pass the threshold of detection. We discuss analytical, statistical properties of the stacked periodogram, perform simulations to demonstrate the efficacy of the method, and investigate the influence of semi-structured window functions and stellar activity. Analysis of the Lick-Carnegie Exoplanet Survey data set reveals a marginally significant ($1.6\sigma$) signal consistent with a population of exoplanets occupying $3-7$ day periods with typical $K$ between $1.6-5.1$ m s$^{-1}$. More detailed investigation of signals associated with stellar activity and yearly systematics may be necessary to confirm this result or detect other underlying Keplerian contributions.


title: The mass profile of NGC 3377 from a Bayesian approach

abstract: The mass profile for the moderately bright elliptical NGC 3377 is studied through an spherical Jeans analysis, combined with a Bayesian approach. The prior distributions are generated from dark matter simulations and observational constraints. The observational data set consist of Gemini/GMOS long-slit observations aligned with the major and minor axes of the galaxy, and are supplemented with data from the literature for the diffuse stellar component, globular clusters and planetary nebulae. Although the galaxy is assumed to alternatively reside in central and satellite haloes, the comparison with literature results prefer the latter option. Several options of constant anisotropy are considered, as well as both NFW and Einasto mass profiles. The analysis points to an intermediate mass halo, presenting a virial mass around $(3.6\pm0.6)\times 10^{11}\,{\rm M_{\odot}}$.


title: Retraction based Direct Search Methods for Derivative Free Riemannian Optimization

abstract: Direct search methods represent a robust and reliable class of algorithms for solving black-box optimization problems. In this paper, we explore the application of those strategies to Riemannian optimization, wherein minimization is to be performed with respect to variables restricted to lie on a manifold. More specifically, we consider classic and line search extrapolated variants of direct search, and, by making use of retractions, we devise tailored strategies for the minimization of both smooth and nonsmooth functions. As such we analyze, for the first time in the literature, a class of retraction based algorithms for minimizing nonsmooth objectives on a Riemannian manifold without having access to (sub)derivatives. Along with convergence guarantees we provide a set of numerical performance illustrations on a standard set of problems.


title: Adsorption, self-assembly and self-metalation of tetra-cyanophenyl porphyrins on semiconducting CoO(100) films

abstract: The adsorption properties of free base 5,10,15,20-tetrakis(p-cyanophenyl)porphyrin (2H-TCNPP) on thin films of rock salt (rs) CoO(100) on Au(111) was studied in ultra-high vacuum by a combination of low-temperature scanning tunneling microscopy and spectroscopy (STM/STS) and density functional theory (DFT). Films of rs-CoO(100) on Au(111) are prepared with excellent quality in a suitable thickness range. Particularly, we found that films of only 1 nm thickness show a semiconducting energy gap of $E_\mathrm{g}=(2.5\pm 0.2)\,\textrm{eV}$. Upon deposition at 300 K, 2H-TCNPP adsorbs flat-lying and self-assembles in a long-range ordered superstructure that is stable at 80 K. The adsorption geometry of the molecules on the surface and within the self-assembly is analyzed by DFT. We find that the self-assemblies are stabilized by hydrogen bridge bonding via the functional cyano groups. Our STS data shows molecular states within the fundamental gap of the CoO. By comparison with the calculated DOS we determine the energetic positions of the frontier orbitals and find that the first three LUMO states 2H-TCNPP are located within the band gap, whereas the HOMO is shifted 1 eV below the CoO conduction band edge. Upon annealing to 420 K the molecules change their appearance in STM images and a new prominent electronic state located at the center of the molecule is formed. We interpret this changed configuration as Co-TCNPP created by self-metalation on the oxide surface.


title: A Parkes "Murriyang" Search for Pulsars and Transients in the Large Magellanic Cloud

abstract: We report on the latest results of a Parkes multibeam survey for pulsars and dispersed radio bursts in the Large Magellanic Cloud (LMC). We conducted both periodicity and single-pulse searches at a much larger range of trial dispersion measures (DMs) than previously searched. We detected 229 single pulses with signal-to-noise ratio ($\rm{S/N}) > 7$ that were classified by the deep learning network FETCH as being real (with $> 90$% likelihood), of which 9 were from the known giant-pulse-emitting pulsar PSR B0540$-$69. Two possibly repeating sources were detected with DMs suggesting that they lie within the LMC, but these require confirmation. Only 3 of the 220 unknown pulses had S/N greater than 8, and the DM distribution for these detected pulses follows an exponential falloff with increasing DM and does not show any excess of signals at DM values expected for the LMC. These features suggest that the detected pulses are not likely to be real, although they are visually compelling. We also report the discovery of a new pulsar (PSR J0556$-$67) in our periodicity search. This pulsar has a spin period of 791 ms, a DM of 71 cm$^{-3}$ pc, an estimated 1400 MHz flux density of $\sim 0.12$ mJy, and shows no evidence of binary motion. Future observations may be able to confirm whether any of the weak but promising candidates in our single pulse and periodicity searches, including our two possible repeaters, are real or not.


title: RMF-Owl: A Collision-Tolerant Flying Robot for Autonomous Subterranean Exploration

abstract: This work presents the design, hardware realization, autonomous exploration and object detection capabilities of RMF-Owl, a new collision-tolerant aerial robot tailored for resilient autonomous subterranean exploration. The system is custom built for underground exploration with focus on collision tolerance, resilient autonomy with robust localization and mapping, alongside high-performance exploration path planning in confined, obstacle-filled and topologically complex underground environments. Moreover, RMF-Owl offers the ability to search, detect and locate objects of interest which can be particularly useful in search and rescue missions. A series of results from field experiments are presented in order to demonstrate the system's ability to autonomously explore challenging unknown underground environments.


title: Quantum regression in dephasing phenomena

abstract: We investigate the validity of quantum regression for a family of quantum Hamiltonians on a multipartite system leading to phase-damping reduced dynamics. After finding necessary and sufficient conditions for the CP-divisibility of the corresponding channel, we evaluate a hierarchy of equations equivalent to the validity of quantum regression under arbitrary interventions; in particular, we find necessary conditions for a nontrivial dephasing to be compatible with quantum regression. In this framework, we study a class of dephasing-type generalized spin-boson (GSB) models, investigating the existence of qubit-environment coupling functions that ensure the exact validity of quantum regression.


title: Survey of Orion Disks with ALMA (SODA) I: Cloud-level demographics of 873 protoplanetary disks

abstract: (Abridged) Surveys of protoplanetary disks in nearby star-forming regions (SFRs) have provided important information on their demographics, but due to their sample sizes, they cannot be used to study how disk properties vary with the environment. We conduct a survey of the unresolved millimeter continuum emission of 873 protoplanetary disks identified by Spitzer in the L1641 and L1647 regions of the Orion A cloud. This is the largest such survey yet, allowing us to identify even weak trends in the median disk mass as a function of position in the cloud and cluster membership. The sample detection rates and median masses are also compared to those of nearby (<300 pc) SFRs. The sample was observed with ALMA at 225 GHz, with a median rms of 0.08 mJy/beam (1.5 M$_{\oplus}$). The data were reduced and imaged using an innovative parallel data processing approach. We detect 58% (502/873) of the observed disks. This includes 20 disks with dust masses >100 M$_{\oplus}$, and two objects associated with extended dust emission. We infer a median disk dust mass in the full sample of $2.2^{+0.2}_{-0.2}$ M$_{\oplus}$. In L1641 and L1647 median dust masses are $2.1^{+0.2}_{-0.2}$ M$_{\oplus}$ and $2.6^{+0.4}_{-0.5}$ M$_{\oplus}$, respectively. The disk mass distribution of the full sample is similar to that of nearby low-mass SFRs at similar ages of 1-3 Myr. We find only weak trends in disk (dust) masses with galactic longitude and between the YSO clusters identified in the sample, with median masses varying by $\lesssim$ 50%. Age differences may explain these median disk mass variations. Apart from this, disk masses are essentially constant at scales of ~100 pc. This also suggests that the majority of disks, even in different SFRs, are formed with similar initial masses and evolve at similar rates, assuming no external irradiation, with disk mass loss rates of $\sim 10^{-8}$ M$_{\odot}$/yr.


title: Transfer matrices for AdS3/CFT2

abstract: We work out the algebraic Bethe ansatz for the worldsheet theory of the $AdS_3\times S^3\times T^4$ superstring, and use it to derive the transfer matrices for fundamental particles and bound states of the string and mirror model. We also show how the Bethe equations and transfer matrices are modified in the presence of an Abelian twist. These will be an important ingredient in the exploration of the mirror thermodynamic Bethe ansatz equations recently proposed by Frolov and Sfondrini, and their generalisation to twisted and deformed models.


title: Fundamental invariants of tensors, Latin hypercubes, and rectangular Kronecker coefficients

abstract: We study polynomial SL-invariants of tensors, mainly focusing on fundamental invariants which are of smallest degrees. In particular, we prove that certain 3-dimensional analogue of the Alon--Tarsi conjecture on Latin cubes considered previously by B\"urgisser and Ikenmeyer, implies positivity of (generalized) Kronecker coefficients at rectangular partitions and as a result provides values for degree sequences of fundamental invariants.


title: Universal approximation of credit portfolio losses using Restricted Boltzmann Machines

abstract: We introduce a new portfolio credit risk model based on Restricted Boltzmann Machines (RBMs), which are stochastic neural networks capable of universal approximation of loss distributions. We test the model on an empirical dataset of default probabilities of 1'012 US companies and we show that it outperforms commonly used parametric factor copula models -- such as the Gaussian or the t factor copula models -- across several credit risk management tasks. In particular, the model leads to better fits for the empirical loss distribution and more accurate risk measure estimations. We introduce an importance sampling procedure which allows risk measures to be estimated at high confidence levels in a computationally efficient way and which is a substantial improvement over the Monte Carlo techniques currently available for copula models. Furthermore, the statistical factors extracted by the model admit an interpretation in terms of the underlying portfolio sector structure and provide practitioners with quantitative tools for the management of concentration risk. Finally, we show how to use the model for stress testing by estimating stressed risk measures (e.g. stressed VaR) under various macroeconomic stress test scenarios, such as those specified by the FRB's Dodd-Frank Act stress test.


title: In This Apportionment Lottery, the House Always Wins

abstract: Apportionment is the problem of distributing $h$ indivisible seats across states in proportion to the states' populations. In the context of the US House of Representatives, this problem has a rich history and is a prime example of interactions between mathematical analysis and political practice. Grimmett (2004) suggested to apportion seats in a randomized way such that each state receives exactly their proportional share $q_i$ of seats in expectation (ex ante proportionality) and receives either $\lfloor q_i \rfloor$ or $\lceil q_i \rceil$ many seats ex post (quota). However, there is a vast space of randomized apportionment methods satisfying these two axioms, and so we additionally consider prominent axioms from the apportionment literature. Our main result is a randomized method satisfying quota, ex ante proportionality and house monotonicity - a property that prevents paradoxes when the number of seats changes and which we require to hold ex post. This result is based on a generalization of dependent rounding on bipartite graphs, which we call cumulative rounding and which might be of independent interest, as we demonstrate via applications beyond apportionment.


title: The relative heat content for submanifolds in sub-Riemannian geometry

abstract: We study the small-time asymptotics of the relative heat content for submanifolds in sub-Riemannian geometry. First, we prove the existence of a smooth tubular neighborhood for submanifolds of any codimension, assuming they do not have characteristic points. Next, we propose a definition of relative heat content for submanifolds of codimension $k\geq 1$ and we build an approximation of this quantity, via smooth tubular neighborhoods. Finally, we show that this approximation fails to recover the asymptotic expansion of the relative heat content of the submanifold, by studying an explicit example.


title: Nested spheroidal figures of equilibrium II. Generalization to L layers

abstract: We present a vectorial formalism to determine the approximate solutions to the problem of a composite body made of $L$ homogeneous, rigidly rotating layers bounded by spheroidal surfaces. The method is based on the 1st-order expansion of the gravitational potential over confocal parameters, thereby generalizing the method described in Paper I for $L=2$. For a given relative geometry of the ellipses and a given set of mass-density jumps at the interfaces, the sequence of rotation rates and interface pressures is obtained analytically by recursion. A wide range of equilibria result when layers rotate in an asynchronous manner, although configurations with a negative oblateness gradient are more favorable. In contrast, states of global rotation (all layers move at the same rate), found by solving a linear system of $L-1$ equations, are much more constrained. In this case, we mathematically demonstrate that confocal and coelliptical configurations are not permitted. Approximate formula for small ellipticities are derived. These results reinforce and prolongate known results and classical theorems restricted to small elliptiticities. Comparisons with the numerical solutions computed from the Self-Consistent-Field method are successful.


title: Occupation similarity through bipartite graphs

abstract: Similarity between occupations is a crucial piece of information when making career decisions. However, the notion of a single and unified occupation similarity measure is more of a limitation than an asset. The goal of the study is to assess multiple explainable occupation similarity measures that can provide different insights into inter-occupation relations. Several such measures are derived using the framework of bipartite graphs. Their viability is assessed on more than 450,000 job transitions occurring in Slovenia in the period between 2012 and 2021. The results support the hypothesis that several similarity measures are plausible and that they present different feasible career paths. The complete implementation and part of the datasets are available at https://repo.ijs.si/pboskoski/bipartite_job_similarity_code.


title: Radio spectra of narrow-line Seyfert 1 galaxies observed with ATCA and VLASS

abstract: We present radio spectral analyses for a sample of 29 radio-quiet (RQ) and three radio-loud (RL) narrow-line Seyfert 1 galaxies (NLS1s) detected with the Australia Telescope Compact Array at both 5.5 GHz and 9.0 GHz. The sample is characterized by L_bol/L_Edd > 0.15. The radio slopes in 25 of the 29 RQ NLS1s are steep (alpha_5.5-9.0 < -0.5), as found in earlier studies of RQ high L_bol/L_Edd AGN. This steep radio emission may be related to AGN-driven outflows, which are likely more prevalent in high L_bol/L_Edd AGN. In two of the three RL NLS1s, the radio slopes are flat or inverted (alpha_5.5-9.0 > -0.5), indicating a compact optically-thick source, likely a relativistic jet. Archival data at 3.0 GHz, 1.4 GHz, and 0.843 GHz are also compiled, yielding a sample of 17 NLS1s detected in three bands or more. In nine objects, the radio spectra flatten at lower frequencies, with median slopes of alpha_5.5-9.0 = -1.21 +/- 0.17, flattening to alpha_3.0-5.5 = -0.97 +/- 0.27, and to alpha_1.4-3.0 = -0.63 +/- 0.16. A parabolic fit suggests a median spectral turnover of ~ 1 GHz, which implies synchrotron self-absorption in a source with a size of only a fraction of 1 pc, possibly a compact wind or a weak jet. Two objects show significant spectral steepening to alpha < -2 above 3 GHz or 5 GHz, which may suggest relic emission from past ejection of radio emitting plasma, of the order of a few years to a few decades ago. Finally, two objects present a single spectral slope consistent with star-forming activity.


title: Outing Power Outages: Real-time and Predictive Socio-demographic Analytics for New York City

abstract: Electrical outages continue to occur despite technological innovations and improvements to electric power distribution infrastructure. In this paper, we describe a tool that was designed to acquire and collect data on electric power outages in New York City since July 2020. The electrical outages are then displayed on a front-end application, which is publicly available. We use the collected outage data to analyze these outages and their socio-economic impacts on electricity vulnerable population groups. We determined that there was a slightly negative linear relationship between income and number of outages. Finally, a Markov Influence Graph was created to better understand the spatial and temporal relationships between outages.


title: Information-driven transitions in projections of underdamped dynamics

abstract: Low-dimensional representations of underdamped systems often provide insightful grasps and analytical tractability. Here, we build such representations via information projections, obtaining an optimal model that captures the most information on observed spatial trajectories. We show that, in paradigmatic systems, the minimization of the information loss drives the appearance of a discontinuous transition in the optimal model parameters. Our results raise serious warnings for general inference approaches and unravel fundamental properties of effective dynamical representations, impacting several fields, from biophysics to dimensionality reduction.


title: Development of a real-time tritium-in-water monitor

abstract: In this paper, we report the development and performance of a detector module envisaging a tritium-in-water real-time activity monitor. The monitor is based on modular detection units whose number can be chosen according to the required sensitivity. The full system is being designed to achieve a \emph{Minimum Detectable Activity} (\emph{MDA}) of 100\,Bq/L of tritium-in-water activity which is the limit established by the E.U. Council Directive 2013/51/Euratom for water intended for human consumption. The same system can be used as a real-time pre-alert system for nuclear power plant regarding tritium-in water environmental surveillance. The first detector module was characterized, commissioned and installed immediately after the discharge channel of the Arrocampo dam (Almaraz nuclear power plant, Spain) on the Tagus river. Due to the high sensitivity of the single detection modules, the system requires radioactive background mitigation techniques through the use of active and passive shielding. We have extrapolated a \emph{MDA} of 3.6\,kBq/L for a single module being this value limited by the cosmic background. The obtained value for a single module is already compatible with a real-time environmental surveillance and pre-alert system. Further optimization of the single-module sensitivity will imply the reduction of the number of modules and the cost of the detector system.


title: Generalised bayesian sample copula of order $m$

abstract: In this work we propose a semiparametric bivariate copula whose density is defined by a piecewise constant function on disjoint squares. We obtain the maximum likelihood estimators of model parameters and prove that they reduce to the sample copula under specific conditions. We further propose to carry out a full Bayesian analysis of the model and introduce a spatial dependent prior distribution for the model parameters. This prior allows the parameters to borrow strength across neighbouring regions to produce smooth posterior estimates. To characterise the posterior distribution, via the full conditional distributions, we propose a data augmentation technique. A Metropolis-Hastings step is required and we propose a novel adaptation scheme for the random walk proposal distribution. We implement a simulation study and an analysis of a real dataset to illustrate the performance of our model and inference algorithms.


title: Surface-tension-driven evolution of a viscoplastic liquid coating the interior of a cylindrical tube

abstract: One mechanism for airway closure in the lung is the surface-tension-driven instability of the mucus layer which lines the airway wall. We study the instability of an axisymmetric layer of viscoplastic Bingham liquid coating the interior of a rigid tube, which is a simple model for an airway that takes into account the yield stress of mucus. An evolution equation for the thickness of the liquid layer is derived using long-wave theory, from which we also derive a simpler thin-film evolution equation. In the thin-film case, we show that two branches of marginally-yielded static solutions of the evolution equation can be used to both predict the size of initial perturbation required to trigger instability and quantify how increasing the capillary Bingham number (a parameter measuring yield stress relative to surface tension) reduces the final deformation of the layer. Using numerical solutions of the long-wave evolution equation, we quantify how the critical layer thickness required to form a liquid plug in the tube increases as the capillary Bingham number is increased. We discuss the significance of these findings for modelling airway closure in obstructive conditions such as cystic fibrosis, where the mucus layer is often thicker and has a higher yield stress.


title: Few-body Bose gases in low dimensions -- a laboratory for quantum dynamics

abstract: Cold atomic gases have become a paradigmatic system for exploring fundamental physics, which at the same time allows for applications in quantum technologies. The accelerating developments in the field have led to a highly advanced set of engineering techniques that, for example, can tune interactions, shape the external geometry, select among a large set of atomic species with different properties, or control the number of atoms. In particular, it is possible to operate in lower dimensions and drive atomic systems into the strongly correlated regime. In this review, we discuss recent advances in few-body cold atom systems confined in low dimensions from a theoretical viewpoint. We mainly focus on bosonic systems in one dimension and provide an introduction to the static properties before we review the state-of-the-art research into quantum dynamical processes stimulated by the presence of correlations. Besides discussing the fundamental physical phenomena arising in these systems, we also provide an overview of the calculational and numerical tools and methods that are commonly used, thus delivering a balanced and comprehensive overview of the field. We conclude by giving an outlook on possible future directions that are interesting to explore in these correlated systems.

categories: cond-mat.quant-gas cond-mat.other physics.atom-ph quant-ph


title: Kolmogorov equations on the space of probability measures associated to the nonlinear filtering equation: the viscosity approach

abstract: We study the backward Kolmogorov equation on the space of probability measures associated to the Kushner-Stratonovich equation of nonlinear filtering. We prove existence and uniqueness in the viscosity sense and, in particular, we provide a comparison theorem. In the context of stochastic filtering it is natural to consider measure-valued processes that satisfy stochastic differential equations. In the literature, a classical way to address this problem is by assuming that these measure-valued processes admit a density. Our approach is different and we work directly with measures. Thus, the backward Kolmogorov equation we study is a second-order partial differential equation of parabolic type on the space of probability measures with compact support. In the literature only few results are available on viscosity solutions for this kind of problems and in particular the uniqueness is a very challenging issue. Here we find a viscosity solution and then we prove that it is unique via comparison theorem.


title: Magnetic properties of equiatomic CrMnFeCoNi

abstract: Magnetic, specific heat, and structural properties of the equiatomic Cantor alloy system are reported for temperatures between 5 kelvin and 300 kelvin, and up to fields of 70 kilo-oersted. Magnetization measurements performed on as-cast, annealed, and cold-worked samples reveal a strong processing history dependence and that high-temperature annealing after cold-working does not restore the alloy to a pristine state. Measurements on known precipitates show that the two transitions, detected at 43 kelvin and 85 kelvin, are intrinsic to the Cantor alloy and not the result of an impurity phase. Experimental and ab initio density functional theory (DFT) computational results suggest that these transitions are a weak ferrimagnetic transition and a spin-glass-like transition, respectively, and magnetic and specific heat measurements provide evidence of significant Stoner enhancement and electron-electron interactions within the material.


title: Stochastic trust-region and direct-search methods: A weak tail bound condition and reduced sample sizing

abstract: Using tail bounds, we introduce a new probabilistic condition for function estimation in stochastic derivative-free optimization which leads to a reduction in the number of samples and eases algorithmic analyses. Moreover, we develop simple stochastic direct-search and trust-region methods for the optimization of a potentially non-smooth function whose values can only be estimated via stochastic observations. For trial points to be accepted, these algorithms require the estimated function values to yield a sufficient decrease measured in terms of a power larger than 1 of the algoritmic stepsize. Our new tail bound condition is precisely imposed on the reduction estimate used to achieve such a sufficient decrease. This condition allows us to select the stepsize power used for sufficient decrease in such a way to reduce the number of samples needed per iteration. In previous works, the number of samples necessary for global convergence at every iteration $k$ of this type of algorithms was $O(\Delta_{k}^{-4})$, where $\Delta_k$ is the stepsize or trust-region radius. However, using the new tail bound condition, and under mild assumptions on the noise, one can prove that such a number of samples is only $O(\Delta_k^{-2 - \varepsilon})$, where $\varepsilon > 0$ can be made arbitrarily small by selecting the power of the stepsize in the sufficient decrease test arbitrarily close to $1$. The global convergence properties of the stochastic direct-search and trust-region algorithms are established under the new tail bound condition.


title: Constrained Visual-Inertial Localization With Application And Benchmark in Laparoscopic Surgery

abstract: We propose a novel method to tackle the visual-inertial localization problem for constrained camera movements. We use residuals from the different modalities to jointly optimize a global cost function. The residuals emerge from IMU measurements, stereoscopic feature points, and constraints on possible solutions in SE(3). In settings where dynamic disturbances are frequent, the residuals reduce the complexity of the problem and make localization feasible. We verify the advantages of our method in a suitable medical use case and produce a dataset capturing a minimally invasive surgery in the abdomen. Our novel clinical dataset MITI is comparable to state-of-the-art evaluation datasets, contains calibration and synchronization and is available at https://mediatum.ub.tum.de/1621941.


title: Topological Universality of the Art Gallery Problem

abstract: We prove that any compact semi-algebraic set is homeomorphic to the solution space of some art gallery problem. Previous works have established similar universality theorems, but holding only up to homotopy equivalence, rather than homeomorphism, and prior to this work, the existence of art galleries even for simple spaces such as the M\"obius strip or the three-holed torus were unknown. Our construction relies on an elegant and versatile gadget to copy guard positions with minimal overhead. It is simpler than previous constructions, consisting of a single rectangular room with convex slits cut out from the edges. We show that both the orientable and non-orientable surfaces of genus $n$ admit galleries with only $O(n)$ vertices.


title: Kinetic Field Theory for Cosmic Structure Formation

abstract: We apply kinetic field theory to non-linear cosmic structure formation. Kinetic field theory decomposes the cosmic density field into particles and follows their trajectories through phase space. We assume that initial particle momenta are drawn from a Gaussian random field. We place particular emphasis on the late-time, asymptotic behaviour on small spatial scales of low-order statistical measures for the distribution of particles in configuration and velocity space. Our main result is that the power spectra for density and velocity fluctuations in ensembles of particles freely streaming along Zel'dovich trajectories asymptotically fall off with wave number $k$ like $k^{-3}$ for $k\to\infty$, irrespective of the cosmological model and the type of dark matter assumed, with the exponent set only by the number of spatial dimensions. This conclusion remains valid for density-fluctuation power spectra if particle interactions are taken into account in a mean-field approximation. We also show that the bispectrum of freely-streaming particles falls off asymptotically like $k^{-11/2}$ under the same general conditions.


title: An interpolating particle method for the Vlasov-Poisson equation

abstract: In this paper we present a novel particle method for the Vlasov--Poisson equation. Unlike in conventional particle methods, the particles are not interpreted as point charges, but as point values of the distribution function. In between the particles, the distribution function is reconstructed using mesh-free interpolation. Our numerical experiments confirm that this approach results in significantly increased accuracy and noise reduction. At the same time, many benefits of the conventional schemes are preserved.


title: Reward-Free Policy Space Compression for Reinforcement Learning

abstract: In reinforcement learning, we encode the potential behaviors of an agent interacting with an environment into an infinite set of policies, the policy space, typically represented by a family of parametric functions. Dealing with such a policy space is a hefty challenge, which often causes sample and computation inefficiencies. However, we argue that a limited number of policies are actually relevant when we also account for the structure of the environment and of the policy parameterization, as many of them would induce very similar interactions, i.e., state-action distributions. In this paper, we seek for a reward-free compression of the policy space into a finite set of representative policies, such that, given any policy $\pi$, the minimum R\'enyi divergence between the state-action distributions of the representative policies and the state-action distribution of $\pi$ is bounded. We show that this compression of the policy space can be formulated as a set cover problem, and it is inherently NP-hard. Nonetheless, we propose a game-theoretic reformulation for which a locally optimal solution can be efficiently found by iteratively stretching the compressed space to cover an adversarial policy. Finally, we provide an empirical evaluation to illustrate the compression procedure in simple domains, and its ripple effects in reinforcement learning.


title: The R136 star cluster dissected with Hubble Space Telescope/STIS. III. The most massive stars and their clumped winds

abstract: Context: The star cluster R136 inside the LMC hosts a rich population of massive stars, including the most massive stars known. The strong stellar winds of these very luminous stars impact their evolution and the surrounding environment. We currently lack detailed knowledge of the wind structure that is needed to quantify this impact. Aims: To observationally constrain the stellar and wind properties of the massive stars in R136, in particular the parameters related to wind clumping. Methods: We simultaneously analyse optical and UV spectroscopy of 53 O-type and 3 WNh-stars using the FASTWIND model atmosphere code and a genetic algorithm. The models account for optically thick clumps and effects related to porosity and velocity-porosity, as well as a non-void interclump medium. Results: We obtain stellar parameters, surface abundances, mass-loss rates, terminal velocities and clumping characteristics and compare these to theoretical predictions and evolutionary models. The clumping properties include the density of the interclump medium and the velocity-porosity of the wind. For the first time, these characteristics are systematically measured for a wide range of effective temperatures and luminosities. Conclusions: We confirm a cluster age of 1.0-2.5 Myr and derive an initial stellar mass of $\geq 250 {\rm M}_\odot$ for the most massive star in our sample, R136a1. The winds of our sample stars are highly clumped, with an average clumping factor of $f_{\rm cl}=29\pm15$. We find tentative trends in the wind-structure parameters as a function of mass-loss rate, suggesting that the winds of stars with higher mass-loss rates are less clumped. We compare several theoretical predictions to the observed mass-loss rates and terminal velocities and find that none satisfactorily reproduces both quantities. The prescription of Krti\v{c}ka & Kub\'at (2018) matches best the observed mass-loss rates.


title: Constraints on the mass and self-coupling of Ultra-Light Scalar Field Dark Matter using observational limits on galactic central mass

abstract: It is well known that Ultra-Light Dark Matter (ULDM), usually scalar fields of mass $m \sim 10^{-22}~{\rm eV}$, can solve some of the outstanding problems of the Cold Dark Matter (CDM) paradigm. Such a scalar field could have non-negligible self-coupling $\lambda$. In this work, using the known observational upper limit on the amount of centrally concentrated dark matter in a galaxy, we arrive at the observational constraints in the $\lambda - m$ (self coupling $-$ mass) parameter space. It is found that the observational limit on the mass $m$ of the ULDM depends upon the sign and strength of the self-interactions. We demonstrate that, for $m \sim 10^{-22}~{\rm eV}$, self-coupling values of ${\cal O}(10^{-96})$ (corresponding to a scattering length of $a_s \sim 10^{-82}~{\rm m}$) can be probed using limits on the dark matter mass within 10 pc of the centre of M87 galaxy. Our analysis suggests that if Ultra Light Axions (ULAs) form all of dark matter, its mass has to be less than $\sim 6 \times 10^{-23}$ eV.


title: A Survey on Scalable LoRaWAN for Massive IoT: Recent Advances, Potentials, and Challenges

abstract: Long-range (LoRa) technology is most widely used for enabling low-power wide area networks (WANs) on unlicensed frequency bands. Despite its modest data rates, it provides extensive coverage for low-power devices, making it an ideal communication system for many internet of things (IoT) applications. In general, LoRa is considered as the physical layer, whereas LoRaWAN is the medium access control (MAC) layer of the LoRa stack that adopts a star topology to enable communication between multiple end devices (EDs) and the network gateway. The chirp spread spectrum modulation deals with LoRa signal interference and ensures long-range communication. At the same time, the adaptive data rate mechanism allows EDs to dynamically alter some LoRa features, such as the spreading factor (SF), code rate, and carrier frequency to address the time variance of communication conditions in dense networks. Despite the high LoRa connectivity demand, LoRa signals interference and concurrent transmission collisions are major limitations. Therefore, to enhance LoRaWAN capacity, the LoRa Alliance released many LoRaWAN versions, and the research community has provided numerous solutions to develop scalable LoRaWAN technology. Hence, we thoroughly examine LoRaWAN scalability challenges and state-of-the-art solutions in both the physical and MAC layers. These solutions primarily rely on SF, logical, and frequency channel assignment, whereas others propose new network topologies or implement signal processing schemes to cancel the interference and allow LoRaWAN to connect more EDs efficiently. A summary of the existing solutions in the literature is provided at the end of the paper, describing the advantages and disadvantages of each solution and suggesting possible enhancements as future research directions.


title: Higher regularity for singular K\"ahler-Einstein metrics

abstract: We study singular K\"ahler-Einstein metrics that are obtained as non-collapsed limits of polarized K\"ahler-Einstein manifolds. Our main result is that if the metric tangent cone at a point is locally isomorphic to the germ of the singularity, then the metric converges to the metric on its tangent cone at a polynomial rate on the level of K\"ahler potentials. When the tangent cone at the point has a smooth cross section, then the result implies polynomial convergence of the metric in the usual sense, generalizing a result due to Hein-Sun. We show that a similar result holds even in certain cases where the tangent cone is not locally isomorphic to the germ of the singularity. Finally we prove a rigidity result for complete $\partial\bar\partial$-exact Calabi-Yau metrics with maximal volume growth. This generalizes a result of Conlon-Hein, which applies to the case of asymptotically conical manifolds.


title: Artificial Viscosity -- Then and Now

abstract: In this paper, we recount the history of artificial viscosity, beginning with its origin in previously unpublished and unavailable documents, continuing on to current research and ending with recent work describing its physical basis that suggests new directions for improvement. We focus on the underlying ideas that recognize the finiteness of scale and of measurement.

categories: physics.comp-ph astro-ph.HE astro-ph.IM math-ph math.MP


title: Dynamical scaling as a signature of multiple phase competition in Yb$_2$Ti$_2$O$_7$

abstract: $\rm Yb_2Ti_2O_7$ is a celebrated example of a pyrochlore magnet with highly-frustrated, anisotropic exchange interactions. To date, attention has largely focused on its unusual, static properties, many of which can be understood as coming from the competition between different types of magnetic order. Here we use inelastic neutron scattering with exceptionally high energy resolution to explore the dynamical properties of $\rm Yb_2Ti_2O_7$. We find that spin correlations exhibit dynamical scaling, analogous to behavior found near to a quantum critical point. We show that the observed scaling collapse can be explained within a phenomenological theory of multiple--phase competition, and confirm that a scaling collapse is also seen in semi--classical simulations of a microscopic model of $\rm Yb_2Ti_2O_7$. These results suggest a general picture for dynamics in systems with competing ground states.


title: StickyLand: Breaking the Linear Presentation of Computational Notebooks

abstract: How can we better organize code in computational notebooks? Notebooks have become a popular tool among data scientists, as they seamlessly weave text and code together, supporting users to rapidly iterate and document code experiments. However, it is often challenging to organize code in notebooks, partially because there is a mismatch between the linear presentation of code and the non-linear process of exploratory data analysis. We present StickyLand, a notebook extension for empowering users to freely organize their code in non-linear ways. With sticky cells that are always shown on the screen, users can quickly access their notes, instantly observe experiment results, and easily build interactive dashboards that support complex visual analytics. Case studies highlight how our tool can enhance notebook users's productivity and identify opportunities for future notebook designs. StickyLand is available at https://github.com/xiaohk/stickyland.


title: Semi-Blind Joint Channel and Symbol Estimation in IRS-Assisted Multi-User MIMO Networks

abstract: Intelligent reflecting surface (IRS) is a promising technology for beyond 5th Generation of the wireless communications. In fully passive IRS-assisted systems, channel estimation is challenging and should be carried out only at the base station or at the terminals since the elements of the IRS are incapable of processing signals. In this letter, we formulate a tensor-based semi-blind receiver that solves the joint channel and symbol estimation problem in an IRS-assisted multi-user multiple-input multiple-output system. The proposed approach relies on a generalized PARATUCK tensor model of the signals reflected by the IRS, based on a two-stage closed-form semi-blind receiver using Khatri-Rao and Kronecker factorizations. Simulation results demonstrate the superior performance of the proposed semi-blind receiver, in terms of the normalized mean squared error and symbol error rate, as well as a lower computational complexity, compared to recently proposed parallel factor analysis-based receivers.


title: A gradient-free subspace-adjusting ensemble sampler for infinite-dimensional Bayesian inverse problems

abstract: Sampling of sharp posteriors in high dimensions is a challenging problem, especially when gradients of the likelihood are unavailable. In low to moderate dimensions, affine-invariant methods, a class of ensemble-based gradient-free methods, have found success in sampling concentrated posteriors. However, the number of ensemble members must exceed the dimension of the unknown state in order for the correct distribution to be targeted. Conversely, the preconditioned Crank-Nicolson (pCN) algorithm succeeds at sampling in high dimensions, but samples become highly correlated when the posterior differs significantly from the prior. In this article we combine the above methods in two different ways as an attempt to find a compromise. The first method involves inflating the proposal covariance in pCN with that of the current ensemble, whilst the second performs approximately affine-invariant steps on a continually adapting low-dimensional subspace, while using pCN on its orthogonal complement.


title: Counterfactual Phenotyping with Censored Time-to-Events

abstract: Estimation of treatment efficacy of real-world clinical interventions involves working with continuous outcomes such as time-to-death, re-hospitalization, or a composite event that may be subject to censoring. Counterfactual reasoning in such scenarios requires decoupling the effects of confounding physiological characteristics that affect baseline survival rates from the effects of the interventions being assessed. In this paper, we present a latent variable approach to model heterogeneous treatment effects by proposing that an individual can belong to one of latent clusters with distinct response characteristics. We show that this latent structure can mediate the base survival rates and helps determine the effects of an intervention. We demonstrate the ability of our approach to discover actionable phenotypes of individuals based on their treatment response on multiple large randomized clinical trials originally conducted to assess appropriate treatments to reduce cardiovascular risk.


title: Total insecurity of communication via strong converse for quantum privacy amplification

abstract: Quantum privacy amplification is a central task in quantum cryptography. Given shared randomness, which is initially correlated with a quantum system held by an eavesdropper, the goal is to extract uniform randomness which is decoupled from the latter. The optimal rate for this task is known to satisfy the strong converse property and we provide a lower bound on the corresponding strong converse exponent. In the strong converse region, the distance of the final state of the protocol from the desired decoupled state converges exponentially fast to its maximal value, in the asymptotic limit. We show that this necessarily leads to totally insecure communication by establishing that the eavesdropper can infer any sent messages with certainty, when given very limited extra information. In fact, we prove that in the strong converse region, the eavesdropper has an exponential advantage in inferring the sent message correctly, compared to the achievability region. Additionally we establish the following technical result, which is central to our proofs, and is of independent interest: the smoothing parameter for the smoothed max-relative entropy satisfies the strong converse property.


title: Efficient and Differentiable Conformal Prediction with General Function Classes

abstract: Quantifying the data uncertainty in learning tasks is often done by learning a prediction interval or prediction set of the label given the input. Two commonly desired properties for learned prediction sets are \emph{valid coverage} and \emph{good efficiency} (such as low length or low cardinality). Conformal prediction is a powerful technique for learning prediction sets with valid coverage, yet by default its conformalization step only learns a single parameter, and does not optimize the efficiency over more expressive function classes. In this paper, we propose a generalization of conformal prediction to multiple learnable parameters, by considering the constrained empirical risk minimization (ERM) problem of finding the most efficient prediction set subject to valid empirical coverage. This meta-algorithm generalizes existing conformal prediction algorithms, and we show that it achieves approximate valid population coverage and near-optimal efficiency within class, whenever the function class in the conformalization step is low-capacity in a certain sense. Next, this ERM problem is challenging to optimize as it involves a non-differentiable coverage constraint. We develop a gradient-based algorithm for it by approximating the original constrained ERM using differentiable surrogate losses and Lagrangians. Experiments show that our algorithm is able to learn valid prediction sets and improve the efficiency significantly over existing approaches in several applications such as prediction intervals with improved length, minimum-volume prediction sets for multi-output regression, and label prediction sets for image classification.

categories: cs.LG cs.AI math.ST stat.ME stat.ML stat.TH


title: ReorientBot: Learning Object Reorientation for Specific-Posed Placement

abstract: Robots need the capability of placing objects in arbitrary, specific poses to rearrange the world and achieve various valuable tasks. Object reorientation plays a crucial role in this as objects may not initially be oriented such that the robot can grasp and then immediately place them in a specific goal pose. In this work, we present a vision-based manipulation system, ReorientBot, which consists of 1) visual scene understanding with pose estimation and volumetric reconstruction using an onboard RGB-D camera; 2) learned waypoint selection for successful and efficient motion generation for reorientation; 3) traditional motion planning to generate a collision-free trajectory from the selected waypoints. We evaluate our method using the YCB objects in both simulation and the real world, achieving 93% overall success, 81% improvement in success rate, and 22% improvement in execution time compared to a heuristic approach. We demonstrate extended multi-object rearrangement showing the general capability of the system.


title: The Surface Group Conjectures for groups with two generators

abstract: The Surface Group Conjectures are statements about recognising surface groups among one-relator groups, using either the structure of their finite-index subgroups, or all subgroups. We resolve these conjectures in the two generator case. More generally, we prove that every two-generator one-relator group with every infinite-index subgroup free is itself either free or a surface group.


title: GroupViT: Semantic Segmentation Emerges from Text Supervision

abstract: Grouping and recognition are important components of visual scene understanding, e.g., for object detection and semantic segmentation. With end-to-end deep learning systems, grouping of image regions usually happens implicitly via top-down supervision from pixel-level recognition labels. Instead, in this paper, we propose to bring back the grouping mechanism into deep networks, which allows semantic segments to emerge automatically with only text supervision. We propose a hierarchical Grouping Vision Transformer (GroupViT), which goes beyond the regular grid structure representation and learns to group image regions into progressively larger arbitrary-shaped segments. We train GroupViT jointly with a text encoder on a large-scale image-text dataset via contrastive losses. With only text supervision and without any pixel-level annotations, GroupViT learns to group together semantic regions and successfully transfers to the task of semantic segmentation in a zero-shot manner, i.e., without any further fine-tuning. It achieves a zero-shot accuracy of 52.3% mIoU on the PASCAL VOC 2012 and 22.4% mIoU on PASCAL Context datasets, and performs competitively to state-of-the-art transfer-learning methods requiring greater levels of supervision. We open-source our code at https://github.com/NVlabs/GroupViT .


title: The Dichotomous Affiliate Stable Matching Problem: Approval-Based Matching with Applicant-Employer Relations

abstract: While the stable marriage problem and its variants model a vast range of matching markets, they fail to capture complex agent relationships, such as the affiliation of applicants and employers in an interview marketplace. To model this problem, the existing literature on matching with externalities permits agents to provide complete and total rankings over matchings based off of both their own and their affiliates' matches. This complete ordering restriction is unrealistic, and further the model may have an empty core. To address this, we introduce the Dichotomous Affiliate Stable Matching (DASM) Problem, where agents' preferences indicate dichotomous acceptance or rejection of another agent in the marketplace, both for themselves and their affiliates. We also assume the agent's preferences over entire matchings are determined by a general weighted valuation function of their (and their affiliates') matches. Our results are threefold: (1) we use a human study to show that real-world matching rankings follow our assumed valuation function; (2) we prove that there always exists a stable solution by providing an efficient, easily-implementable algorithm that finds such a solution; and (3) we experimentally validate the efficiency of our algorithm versus a linear-programming-based approach.


title: Radio Signals from Axion Star-Neutron Star Binaries

abstract: Axion stars could form binaries with neutron stars. Given the extremely strong external magnetic field exhibited by individual neutron stars, there can be a substantial conversion of axions to photons in these binaries. The photon emission is doubly modulated due to the neutron star spinning and the axion star orbiting, yielding a unique discovery signal. Similar features are also generated in binaries between a neutron star and an axion-clouded black hole. Encouragingly, such binaries are found to be within the reach of ongoing and upcoming experiments (e.g., the Five hundred meter Aperture Spherical Telescope and the future Square Kilometer Array) for certain parameter regions. They thus provide a promising astronomical laboratory for detecting axions and axion dark matter.


title: Message passing all the way up

abstract: The message passing framework is the foundation of the immense success enjoyed by graph neural networks (GNNs) in recent years. In spite of its elegance, there exist many problems it provably cannot solve over given input graphs. This has led to a surge of research on going "beyond message passing", building GNNs which do not suffer from those limitations -- a term which has become ubiquitous in regular discourse. However, have those methods truly moved beyond message passing? In this position paper, I argue about the dangers of using this term -- especially when teaching graph representation learning to newcomers. I show that any function of interest we want to compute over graphs can, in all likelihood, be expressed using pairwise message passing -- just over a potentially modified graph, and argue how most practical implementations subtly do this kind of trick anyway. Hoping to initiate a productive discussion, I propose replacing "beyond message passing" with a more tame term, "augmented message passing".


title: Hybrid Learning for Orchestrating Deep Learning Inference in Multi-user Edge-cloud Networks

abstract: Deep-learning-based intelligent services have become prevalent in cyber-physical applications including smart cities and health-care. Collaborative end-edge-cloud computing for deep learning provides a range of performance and efficiency that can address application requirements through computation offloading. The decision to offload computation is a communication-computation co-optimization problem that varies with both system parameters (e.g., network condition) and workload characteristics (e.g., inputs). Identifying optimal orchestration considering the cross-layer opportunities and requirements in the face of varying system dynamics is a challenging multi-dimensional problem. While Reinforcement Learning (RL) approaches have been proposed earlier, they suffer from a large number of trial-and-errors during the learning process resulting in excessive time and resource consumption. We present a Hybrid Learning orchestration framework that reduces the number of interactions with the system environment by combining model-based and model-free reinforcement learning. Our Deep Learning inference orchestration strategy employs reinforcement learning to find the optimal orchestration policy. Furthermore, we deploy Hybrid Learning (HL) to accelerate the RL learning process and reduce the number of direct samplings. We demonstrate efficacy of our HL strategy through experimental comparison with state-of-the-art RL-based inference orchestration, demonstrating that our HL strategy accelerates the learning process by up to 166.6x.


title: Rotationally Equivariant Super-Resolution of Velocity Fields in Two-Dimensional Fluids Using Convolutional Neural Networks

abstract: This paper investigates the super-resolution (SR) of velocity fields in two-dimensional fluids from the viewpoint of rotational equivariance. SR refers to techniques that estimate high-resolution images from those in low resolution and has lately been applied in fluid mechanics. The rotational equivariance of SR models is defined as the property in which the super-resolved velocity field is rotated according to a rotation of the input, which leads to the inference covariant to the orientation of fluid systems. Generally, the covariance in physics is related to symmetries. To clarify a relationship to symmetries, the rotational consistency of datasets for SR is newly introduced as the invariance of pairs of low- and high-resolution velocity fields with respect to rotation. This consistency is sufficient and necessary for SR models to acquire rotational equivariance from large datasets with supervised learning. Such a large dataset is not required when rotational equivariance is imposed on SR models through weight sharing of convolution kernels as prior knowledge. Even if a fluid system has rotational symmetry, this symmetry may not carry over to a velocity dataset, which is not rotationally consistent. This inconsistency can occur when the rotation does not commute with the generation of low-resolution velocity fields. These theoretical suggestions are supported by the results from numerical experiments, where two existing convolutional neural networks (CNNs) are converted into rotationally equivariant CNNs and the inferences of the four CNNs are compared after the supervised training.


title: Triangular lattice quantum dimer model with variable dimer density

abstract: Quantum dimer models are known to host topological quantum spin liquid phases, and it has recently become possible to simulate such models with Rydberg atoms trapped in arrays of optical tweezers. Here, we present large-scale quantum Monte Carlo simulation results on an extension of the triangular lattice quantum dimer model with terms in the Hamiltonian annihilating and creating single dimers. We find distinct odd and even $\mathbb{Z}_2$ spin liquids, along with several phases with no topological order: a staggered crystal, a nematic phase, and a trivial symmetric phase with no obvious broken symmetry. We also present dynamic spectra of the phases, and note implications for experiments on Rydberg atoms.

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


title: TDCOSMO. IX. Systematic comparison between lens modelling software programs: time delay prediction for WGD 2038-4008

abstract: The importance of alternative methods for measuring the Hubble constant, such as time-delay cosmography, is highlighted by the recent Hubble tension. It is paramount to thoroughly investigate and rule out systematic biases in all measurement methods before we can accept new physics as the source of this tension. In this study, we perform a check for systematic biases in the lens modelling procedure of time-delay cosmography by comparing independent and blind time-delay predictions of the system WGD 2038-4008 from two teams using two different software programs: GLEE and LENSTRONOMY. The predicted time delays from the two teams incorporate the stellar kinematics of the deflector and the external convergence from line-of-sight structures. The un-blinded time-delay predictions from the two teams agree within $1.2\sigma$, implying that once the time delay is measured the inferred Hubble constant will also be mutually consistent. However, there is a $\sim$4$\sigma$ discrepancy between the power-law model slope and external shear, which is a significant discrepancy at the level of lens models before the stellar kinematics and the external convergence are incorporated. We identify the difference in the reconstructed point spread function (PSF) to be the source of this discrepancy. When the same reconstructed PSF was used by both teams, we achieved excellent agreement, within $\sim$0.6$\sigma$, indicating that potential systematics stemming from source reconstruction algorithms and investigator choices are well under control. We recommend that future studies supersample the PSF as needed and marginalize over multiple algorithms or realizations for the PSF reconstruction to mitigate the systematics associated with the PSF. A future study will measure the time delays of the system WGD 2038-4008 and infer the Hubble constant based on our mass models.


title: The Boring Monopole

abstract: We study false vacuum decays catalysed by magnetic monopoles which act as tunnelling sites with exponentially enhanced decay rates. The field profiles describing the decay do not have the typically assumed $O(3)/O(4)$ symmetry, thus requiring an extension of the usual decay rate calculation. To numerically determine the saddle point solutions which describe the tunnelling process we use a new algorithm based on the mountain pass theorem. This method can be applied more widely to phase transitions with reduced symmetry, such as decays away from the zero and infinite temperature limits. Our results show that monopole-catalysed tunnelling can dominate over the homogeneous false vacuum decay for a wide range of parameters, significantly modify the gravitational wave signal or trigger phase transitions which would otherwise never complete. A single boring monopole in our Hubble patch may determine the lifetime of the Standard Model vacuum.


title: Nonlinear optical responses in nodal line semimetals

abstract: Shift current and second harmonic generation (SHG) are nonlinear optical responses that are often used for photovoltaic effect and the detection of subtle symmetry breaking patterns in materials. It has been recently shown that topological semimetals such as Weyl semimetals with broken inversion symmetry can generate large nonlinear optical responses. In this paper, we investigate both shift current and SHG in a class of topological nodal line semimetal (NLSM) systems with broken inversion symmetry, which may arise intrinsically or via an applied electric field. It is shown that certain classes of NLSM systems offer more singular shift current and SHG in comparison to the well-known Weyl semimetal systems.


title: A duality connecting neural network and cosmological dynamics

abstract: We demonstrate that the dynamics of neural networks trained with gradient descent and the dynamics of scalar fields in a flat, vacuum energy dominated Universe are structurally profoundly related. This duality provides the framework for synergies between these systems, to understand and explain neural network dynamics and new ways of simulating and describing early Universe models. Working in the continuous-time limit of neural networks, we analytically match the dynamics of the mean background and the dynamics of small perturbations around the mean field, highlighting potential differences in separate limits. We perform empirical tests of this analytic description and quantitatively show the dependence of the effective field theory parameters on hyperparameters of the neural network. As a result of this duality, the cosmological constant is matched inversely to the learning rate in the gradient descent update.


title: Extremely Buried Nucleus of IRAS 17208$-$0014 Observed at Sub-Millimeter and Near-Infrared Wavelengths

abstract: The ultraluminous infrared galaxy IRAS 17208$-$0014 is a late-stage merger that hosts a buried active galactic nucleus (AGN). To investigate its nuclear structure, we performed high spatial resolution ($\sim0.\!\!^{\prime\prime}04\sim32\,\mathrm{pc}$) Atacama Large Millimeter/submillimeter Array (ALMA) observations in Band 9 ($\sim$450\,\micron\ or $\sim$660\,GHz), along with near-infrared AKARI spectroscopy in 2.5--5.0\,\micron. The Band 9 dust continuum peaks at the AGN location, and toward this position CO($J$=6--5) and CS($J$=14--13) are detected in absorption. Comparison with non-local thermal equilibrium calculations indicates that, within the central beam ($r\sim20\,\mathrm{pc}$), there exists a concentrated component that is dense ($10^7\,\mathrm{cm}^{-2}$) and warm ($>$200\,K) and has a large column density ($N_\mathrm{H_2}>10^{23}\,\mathrm{cm}^{-2}$). The AKARI spectrum shows deep and broad CO ro-vibrational absorption at 4.67\,\micron. Its band profile is well reproduced with a similarly dense and large column but hotter ($\sim$1000\,K) gas. The region observed through absorption in the near-infrared is highly likely in the nuclear direction, as in the sub-millimeter, but with a narrower beam including a region closer to the nucleus. The central component is considered to possess a hot structure where vibrationally excited HCN emission originates. The most plausible heating source for the gas is X-rays from the AGN. The AKARI spectrum does not show other AGN signs in 2.5--4\,\micron, but this absence may be usual for AGNs buried in a hot mid-infrared core. Besides, based on our ALMA observations, we relate various nuclear structures of IRAS 17208$-$0014 that have been proposed in the literature.