title: Introduction to Monte Carlo for Matrix Models

abstract: We consider a wide range of matrix models and study them using the Monte Carlo technique in the large $N$ limit. The results we obtain agree with exact analytic expressions and recent numerical bootstrap methods for models with one and two matrices. We then present new results for several unsolved multi-matrix models where no other tool is yet available. In order to encourage an exchange of ideas between different numerical approaches to matrix models, we provide programs in Python that can be easily modified to study potentials other than the ones discussed. These programs were tested on a laptop and took between a few minutes to several hours to finish depending on the model, $N$, and the required precision.


title: The THESAN project: predictions for multi-tracer line intensity mapping in the epoch of reionization

abstract: Line intensity mapping (LIM) is rapidly emerging as a powerful technique to study galaxy formation and cosmology in the high-redshift Universe. We present LIM estimates of select spectral lines originating from the interstellar medium (ISM) of galaxies and 21 cm emission from neutral hydrogen gas in the Universe using the large volume, high resolution THESAN reionization simulations. A combination of sub-resolution photo-ionization modelling for HII regions and Monte Carlo radiative transfer calculations is employed to estimate the dust-attenuated spectral energy distributions (SEDs) of high-redshift galaxies ($z\gtrsim5.5$). We show that the derived photometric properties such as the ultraviolet (UV) luminosity function and the UV continuum slopes match observationally inferred values, demonstrating the accuracy of the SED modelling. We provide fits to the luminosity--star formation rate relation (L-SFR) for the brightest emission lines and find that important differences exist between the derived scaling relations and the widely used low-$z$ ones because the interstellar medium of reionization era galaxies is generally less metal-enriched than in their low redshift counterparts. We use these relations to construct line intensity maps of nebular emission lines and cross correlate with the 21 cm emission. Interestingly, the wavenumber at which the correlation switches sign ($k_\mathrm{transition}$) depends heavily on the reionization model and to a lesser extent on the targeted emission line, which is consistent with the picture that $k_\mathrm{transition}$ probes the typical sizes of ionized regions. The derived scaling relations and intensity maps represent a timely state-of-the-art framework for forecasting and interpreting results from current and upcoming LIM experiments.


title: Quantum theory of feedback cooling of an anelastic macro-mechanical oscillator

abstract: Conventional techniques for laser cooling, by coherent scattering off of internal states or through an optical cavity mode, have so far proved inefficient on mechanical oscillators heavier than a few nanograms. That is because larger oscillators vibrate at frequencies much too small compared to the scattering rates achievable by their coupling to auxiliary modes. Decoherence mechanisms typically observed in heavy low frequency elastically suspended oscillators also differ markedly from what is assumed in conventional treatments of laser cooling. We show that for a low-frequency anelastic oscillator forming the mechanically compliant end-mirror of a cavity, detuned optical readout, together with measurement-based feedback to stiffen and dampen it, can harness ponderomotively generated quantum correlations, to realize efficient cooling to the motional ground state. This will pave the way for experiments that call for milligram-scale mechanical oscillators prepared in pure motional states, for example, for tests of gravity's effect on massive quantum systems.


title: Higher-order non-global logarithms from jet calculus

abstract: Non-global QCD observables are characterised by a sensitivity to the full angular distribution of soft radiation emitted coherently in hard scattering processes. This complexity poses a challenge to their all-order resummation, that was formulated at the leading-logarithmic order about two decades ago. In this article we present a solution to the long-standing problem of their resummation beyond this order, and carry out the first complete next-to-leading logarithmic calculation for non-global observables. This is achieved by solving numerically the recently derived set of non-linear differential equations which describe the evolution of soft radiation in the planar, large-$N_c$ limit. As a case study we address the calculation of the transverse energy distribution in the interjet rapidity region in $e^+e^-\to$ dijet production. The calculation is performed by means of an algorithm that we formulate in the language of jet-calculus generating functionals, which also makes the resummation technique applicable to more general non-global problems, such as those that arise in hadronic collisions. We find that NLL corrections are substantial and their inclusion leads to a significant reduction of the perturbative scale uncertainties for these observables. The computer code used in the calculations is made publicly available.


title: Premature Black Hole Death of Population III Stars by Dark Matter

abstract: Population III stars were the first generation of stars, formed in minihalos of roughly primordial element abundances, and therefore metal-free. They are thought to have formed at the cores of dense dark matter clouds. Interactions between baryons and dark matter can therefore have had an important impact on their evolution. In this paper we consider the capture of non- or weakly-annihilating dark matter by these early massive stars. In a wide region of parameter space, interactions of dark matter with baryons lead to premature death of the star as a black hole. We sketch how this modification of the standard evolutionary history of Population III stars might impact the epoch of reionisation, by modifying the amount of UV emission, the transition to Population II star formation, and the X-ray and radio emission from accretion onto the black hole remnants. Signals of massive black holes originating from Population III stars could be observed through gravitational waves from their mergers. Finally, the observation of pair-instability supernovae could effectively preclude premature black hole death across a wide range of parameter space, ranging in mass from $m_{DM} \sim 0.1\text{ GeV}$ to $m_{DM} \sim m_{\rm Pl}$.

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


title: Citizen ASAS-SN Data Release I: Variable Star Classification Using Citizen Science

abstract: We present the first results from Citizen ASAS-SN, a citizen science project for the All-Sky Automated Survey for Supernovae (ASAS-SN) hosted on the Zooniverse platform. Citizen ASAS-SN utilizes the newer, deeper, higher cadence ASAS-SN $g$-band data and tasks volunteers to classify periodic variable star candidates based on their phased light curves. We started from 40,640 new variable candidates from an input list of ${\sim} 7.4$ million stars with $\delta < -60^\circ$ and the volunteers identified 10,420 new discoveries which they classified as 4,234 pulsating variables, 3,132 rotational variables, 2,923 eclipsing binaries, and 131 variables flagged as Unknown. They classified known variable stars with an accuracy of 89% for pulsating variables, 81% for eclipsing binaries, and 49% for rotational variables. We examine user performance, agreement between users, and compare the citizen science classifications with our machine learning classifier updated for the $g$-band light curves. In general, user activity correlates with higher classification accuracy and higher user agreement. We used the user's "Junk" classifications to develop an effective machine learning classifier to separate real from false variables, and there is a clear path for using this "Junk" training set to significantly improve our primary machine learning classifier. We also illustrate the value of Citizen ASAS-SN for identifying unusual variables with several examples.


title: TOI-712: a system of adolescent mini-Neptunes extending to the habitable zone

abstract: As an all-sky survey, NASA's $TESS$ mission is able to detect the brightest and rarest types of transiting planetary systems, including young planets that enable study of the evolutionary processes that occur within the first billion years. Here, we report the discovery of a young, multi-planet system orbiting the bright K4.5V star, TOI-712 ($V = 10.838$, $M_\star = 0.733_{-0.025}^{+0.026} M_\odot$, $R_\star = 0.674\pm0.016 R_\odot$, $T_{\rm eff} = 4622_{-60}^{+61}$ K). From the $TESS$ light curve, we measure a rotation period of 12.48 days, and derive an age between about $500$ Myr and 1.1 Gyr. The photometric observations reveal three transiting mini-Neptunes ($R_b = 2.049^{+0.12}_{-0.080} R_\oplus$, $R_c = 2.701^{+0.092}_{-0.082} R_\oplus$, $R_d = 2.474^{+0.090}_{-0.082} R_\oplus $), with orbital periods of $P_b = 9.531$ days, $P_c = 51.699$ days, and $P_d = 84.839$ days. After modeling the three-planet system, an additional Earth-sized candidate is identified, TOI-712.05 ($P = 4.32$ days, $R_P = 0.81 \pm 0.11 R_\oplus$). We calculate that the habitable zone falls between 0.339 and 0.844 au (82.7 and 325.3 days), placing TOI-712 d near its inner edge. Among planetary systems harboring temperate planets, TOI-712 ($T = 9.9$) stands out as a relatively young star bright enough to motivate further characterization.


title: Patterns of $C$- and $CP$-violation in hadronic $\eta$ and $\eta'$ three-body decays

abstract: We construct hadronic amplitudes for the three-body decays $\eta^{(\prime)}\to\pi^+\pi^-\pi^0$ and $\eta'\to\eta\pi^+\pi^-$ in a non-perturbative fashion, allowing for $C$- and $CP$-violating asymmetries in the $\pi^+\pi^-$ distributions. These amplitudes are consistent with the constraints of analyticity and unitarity. We find that the currently most accurate Dalitz-plot distributions taken by the KLOE-2 and BESIII collaborations confine the patterns of these asymmetries to a relative per mille level. Our dispersive representation allows us to extract the individual coupling strengths of the $C$- and $CP$-violating contributions arising from effective isoscalar and isotensor operators in $\eta^{(\prime)}\to\pi^+\pi^-\pi^0$ and an effective isovector operator in $\eta'\to\eta\pi^+\pi^-$, while the strongly different sensitivities to these operators can be understood from chiral power counting arguments.


title: Proof of the universal density of charged states in QFT

abstract: We prove a recent conjecture by Harlow and Ooguri concerning a universal formula for the charged density of states in QFT at high energies for global symmetries associated with finite groups. An equivalent statement, based on the entropic order parameter associated with charged operators in the thermofield double state, was proven in a previous article by Casini, Huerta, Pontello, and the present author. Here we describe how the statement about the entropic order parameter arises, and how it gets transformed into the universal density of states. The use of the certainty principle, relating the entropic order and disorder parameters, is crucial for the proof. We remark that although the immediate application of this result concerns charged states, the origin and physics of such density can be understood by looking at the vacuum sector only. We also describe how these arguments lie at the origin of the so-called entropy equipartition in these type of systems, and how they generalize to QFT's on non-compact manifolds.


title: Cosmological inference from the emulator based halo model II: Joint analysis of galaxy-galaxy weak lensing and galaxy clustering from HSC-Y1 and SDSS

abstract: We present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing ($\Delta\!\Sigma$) and projected galaxy clustering ($w_{\rm p}$) measured from the Hyper Suprime-Cam Year-1 (HSC-Y1) data and spectroscopic Sloan Digital Sky Survey (SDSS) galaxy catalogs in the redshift range $0.15<z<0.7$. We define luminosity-limited samples of SDSS galaxies to serve as the tracers of $w_{\rm p}$ in three spectroscopic redshift bins, and as the lens samples for $\Delta\!\Sigma$. For the $\Delta\!\Sigma$ measurements, we select a single sample of 4 million source galaxies over 140 deg$^2$ from HSC-Y1 with photometric redshifts (photo-$z$) greater than 0.75, enabling a better handle of photo-$z$ errors by comparing the $\Delta\!\Sigma$ amplitudes for the three lens redshift bins. For cosmological parameter inference, we use an input galaxy-halo connection model built on the {\tt Dark Emulator} package with a halo occupation distribution that includes nuisance parameters to marginalize over modeling uncertainties. We model the $\Delta\!\Sigma$ and $w_{\rm p}$ measurements on scales from $R\simeq 3$ and $2\,h^{-1}{\rm Mpc}$, respectively, up to $30\,h^{-1}{\rm Mpc}$ assuming a flat $\Lambda$CDM cosmology. With various tests using mock catalogs described in Miyatake et al. (2021), we show that any bias in the clustering amplitude $S_8\equiv \sigma_8(\Omega_{\rm m}/0.3)^{0.5}$ due to uncertainties in the galaxy-halo connection is less than $\sim50$\% of the statistical uncertainty of $S_8$, {\it unless} the assembly bias effect is unexpectedly large. Our best-fit models have $S_8=0.795^{+0.049}_{-0.042}$ (mode and 68\% credible interval) for the flat $\Lambda$CDM model; we find tighter constraints on the quantity $S_8(\alpha=0.17)\equiv\sigma_8(\Omega_{\rm m}/0.3)^{0.17} =0.745^{+0.039}_{-0.031}$. (abriged)


title: Quasar cosmology: dark energy evolution and spatial curvature

abstract: We analyse some open debates in cosmology in light of the most updated quasar (QSO) sample, covering a wide redshift range up to $z\sim7.5$, combined with type Ia supernovae (SNe) and baryon acoustic oscillations (BAO). Indeed, extending the cosmological analyses with high-redshift data is key to distinguishing between different cosmological models that are degenerate at low redshifts, and allowing better constraints on a possible dark energy (DE) evolution. Also, we discuss combinations of BAO, SNe, and QSO data to understand their compatibility and implications for extensions of the standard cosmological model. Specifically, we consider a flat and non-flat $\Lambda \mathrm{CDM}$ cosmology, a flat and non-flat DE model with a constant DE equation of state parameter ($w$), and four flat DE models with variable $w$, namely the Chevallier-Polarski-Linder and Jassal-Bagla-Padmanabhan models, and an "exponential" and Barboza-Alcaniz parameterisations. We find that a joint analysis of QSO+SNe with BAO is only possible in the context of a flat Universe. Indeed BAO confirms the flatness condition assuming a curved geometry, whilst SNe+QSO show evidence of a closed space. We also find $\Omega_{M,0}=0.3$ in all data sets assuming a flat $\Lambda \mathrm{CDM}$ model. Yet, all the other models show a statistically significant deviation at 2-3$\sigma$ with the combined SNe+QSO+BAO data set. In the models where DE density evolves with time, SNe+QSO+BAO data always prefer $\Omega_{M,0}>0.3$, $w_{0}<-1$ and $w_{a}>0$. This DE phantom behaviour is mainly driven by SNe+QSO, while BAO are closer to the flat $\Lambda \mathrm{CDM}$ model.


title: $pp$ Solar Neutrinos at DARWIN

abstract: The DARWIN collaboration recently argued that DARWIN (DARk matter WImp search with liquid xenoN) can collect, via neutrino--electron scattering, a large, useful sample of solar $pp$-neutrinos, and measure their survival probability with sub-percent precision. We explore the physics potential of such a sample in more detail. We estimate that, with 300 ton-years of data, DARWIN can also measure, with the help of current solar neutrino data, the value of $\sin^2\theta_{13}$, with the potential to exclude $\sin^2\theta_{13}=0$ close to the three-sigma level. We explore in some detail how well DARWIN can constrain the existence of a new neutrino mass-eigenstate $\nu_4$ that is quasi-mass-degenerate with $\nu_1$ and find that DARWIN's sensitivity supersedes that of all current and near-future searches for new, very light neutrinos. In particular, DARWIN can test the hypothesis that $\nu_1$ is a pseudo-Dirac fermion as long as the induced mass-squared difference is larger than $10^{-13}$ eV$^2$, one order of magnitude more sensitive than existing constraints. Throughout, we allowed for the hypotheses that DARWIN is filled with natural xenon or $^{136}$Xe-depleted xenon.


title: Modeling the galaxy-halo connection with machine learning

abstract: To extract information from the clustering of galaxies on non-linear scales, we need to model the connection between galaxies and halos accurately and in a flexible manner. Standard halo occupation distribution (HOD) models make the assumption that the galaxy occupation in a halo is a function of only its mass, however, in reality, the occupation can depend on various other parameters including halo concentration, assembly history, environment, spin, etc. Using the IllustrisTNG hydrodynamic simulation as our target, we show that machine learning tools can be used to capture this high-dimensional dependence and provide more accurate galaxy occupation models. Specifically, we use a random forest regressor to identify which secondary halo parameters best model the galaxy-halo connection and symbolic regression to augment the standard HOD model with simple equations capturing the dependence on those parameters, namely the local environmental overdensity and shear, at the location of a halo. This not only provides insights into the galaxy-formation relationship but, more importantly, improves the clustering statistics of the modeled galaxies significantly. Our approach demonstrates that machine learning tools can help us better understand and model the galaxy-halo connection, and are therefore useful for galaxy formation and cosmology studies from upcoming galaxy surveys.


title: Closed form expression for the 5-point Veneziano amplitude in terms of 4-point amplitudes

abstract: We propose a number system covariance principle between $p$-adic and Archimedean frameworks. We use it to derive several closed-form expressions for the five-point open string tachyon scattering amplitude.


title: NRPyElliptic: A Fast Hyperbolic Relaxation Elliptic Solver for Numerical Relativity, I: Conformally Flat, Binary Puncture Initial Data

abstract: We introduce NRPyElliptic, an elliptic solver for numerical relativity (NR) built within the NRPy+ framework. As its first application, NRPyElliptic sets up conformally flat, binary black hole (BBH) puncture initial data (ID) on a single numerical domain, similar to the widely used TwoPunctures code. Unlike TwoPunctures, NRPyElliptic employs a hyperbolic relaxation scheme, whereby arbitrary elliptic PDEs are trivially transformed into a hyperbolic system of PDEs. As consumers of NR ID generally already possess expertise in solving hyperbolic PDEs, they will generally find NRPyElliptic easier to tweak and extend than other NR elliptic solvers. When evolved forward in (pseudo)time, the hyperbolic system exponentially reaches a steady state that solves the elliptic PDEs. Notably NRPyElliptic accelerates the relaxation waves, which makes it many orders of magnitude faster than the usual constant-wavespeed approach. While it is still ${\sim}12$x slower than TwoPunctures at setting up full-3D BBH ID, NRPyElliptic requires only ${\approx}0.3\%$ of the runtime for a full BBH simulation in the Einstein Toolkit. Future work will focus on improving performance and generating other types of ID, such as binary neutron star.


title: The Simons Observatory: Galactic Science Goals and Forecasts

abstract: Observing in six frequency bands from 27 to 280 GHz over a large sky area, the Simons Observatory (SO) is poised to address many questions in Galactic astrophysics in addition to its principal cosmological goals. In this work, we provide quantitative forecasts on astrophysical parameters of interest for a range of Galactic science cases. We find that SO can: constrain the frequency spectrum of polarized dust emission at a level of $\Delta\beta_d \lesssim 0.01$ and thus test models of dust composition that predict that $\beta_d$ in polarization differs from that measured in total intensity; measure the correlation coefficient between polarized dust and synchrotron emission with a factor of two greater precision than current constraints; exclude the non-existence of exo-Oort clouds at roughly 2.9$\sigma$ if the true fraction is similar to the detection rate of giant planets; map more than 850 molecular clouds with at least 50 independent polarization measurements at 1 pc resolution; detect or place upper limits on the polarization fractions of CO(2-1) emission and anomalous microwave emission at the 0.1% level in select regions; and measure the correlation coefficient between optical starlight polarization and microwave polarized dust emission in $1^\circ$ patches for all lines of sight with $N_{\rm H} \gtrsim 2\times10^{20}$ cm$^{-2}$. The goals and forecasts outlined here provide a roadmap for other microwave polarization experiments to expand their scientific scope via Milky Way astrophysics.


title: Weighted Quantum Channel Compiling through Proximal Policy Optimization

abstract: We propose a general and systematic strategy to compile arbitrary quantum channels without using ancillary qubits, based on proximal policy optimization -- a powerful deep reinforcement learning algorithm. We rigorously prove that, in sharp contrast to the case of compiling unitary gates, it is impossible to compile an arbitrary channel to arbitrary precision with any given finite elementary channel set, regardless of the length of the decomposition sequence. However, for a fixed accuracy $\epsilon$ one can construct a universal set with constant number of $\epsilon$-dependent elementary channels, such that an arbitrary quantum channel can be decomposed into a sequence of these elementary channels followed by a unitary gate, with the sequence length bounded by $O(\frac{1}{\epsilon}\log\frac{1}{\epsilon})$. Through a concrete example concerning topological compiling of Majorana fermions, we show that our proposed algorithm can conveniently and effectively reduce the use of expensive elementary gates through adding the weighted cost into the reward function of the proximal policy optimization.


title: Stabilizing Disorder-Free Localization

abstract: Disorder-free localization is a paradigm of nonergodicity in translation-invariant quantum many-body systems hosting gauge symmetries. The quench dynamics starting from simple initial states, which correspond to extensive superpositions of gauge superselection sectors, exhibits many-body localization with the system dynamically inducing its own disorder. An open question concerns the stability of disorder-free localization in the presence of gauge-breaking errors, and whether processes due to the latter can be controllably suppressed. Here, we show that translation-invariant \textit{single-body gauge terms} induce a quantum Zeno effect that reliably protects disorder-free localization against errors up to times at least polynomial in the protection strength. Our experimentally feasible scheme not only shows that disorder-free localization can be reliably stabilized, but also opens promising prospects for its observation in quantum simulators.

categories: cond-mat.dis-nn cond-mat.quant-gas cond-mat.stat-mech cond-mat.str-el quant-ph


title: Characterization of Two 2mm-detected Optically-Obscured Dusty Star-Forming Galaxies

abstract: The 2mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high redshift ($z\gtrsim4$), massive, dusty star-forming galaxies (DSFGs). Here we present two, likely high redshift sources, identified in the survey whose physical characteristics are consistent with a class of optical/near-infrared (OIR) invisible DSFGs found elsewhere in the literature. We first perform a rigorous analysis of all available photometric data to fit spectral energy distributions and estimate redshifts before deriving physical properties based on our findings. Our results suggest the two galaxies, called MORA-5 and MORA-9, represent two extremes of the "OIR-dark" class of DSFGs. MORA-5 ($z_{\rm phot}=4.3^{+1.5}_{-1.3}$) is a significantly more active starburst with a star-formation rate of 830$^{+340}_{-190}$M$_\odot$yr$^{-1}$ compared to MORA-9 ($z_{\rm phot}=4.3^{+1.3}_{-1.0}$) whose star-formation rate is a modest 200$^{+250}_{-60}$M$_\odot$yr$^{-1}$. Based on the stellar masses (M$_{\star}\approx10^{10-11}$M$_\odot$), space density ($n\sim(5\pm2)\times10^{-6}$Mpc$^{-3}$, which incorporates two other spectroscopically confirmed OIR-dark DSFGs in the MORA sample at $z=4.6$ and $z=5.9$), and gas depletion timescales ($<1$Gyr) of these sources, we find evidence supporting the theory that OIR-dark DSFGs are the progenitors of recently discovered $3<z<4$ massive quiescent galaxies.


title: Does Complexity Equal Anything?

abstract: We present a new infinite class of gravitational observables in asymptotically Anti-de Sitter space living on codimension-one slices of the geometry, the most famous of which is the volume of the maximal slice. We show that these observables display universal features for the thermofield-double state: they grow linearly in time at late times and reproduce the switch-back effect in shock wave geometries. We argue that any member of this class of observables is an equally viable candidate as the extremal volume for a gravitational dual of complexity.


title: VLA Radio Study of a Sample of Nearby X-ray and Optically Bright Early-Type Galaxies

abstract: Many massive early-type galaxies host central radio sources and hot X-ray atmospheres indicating the presence of radio-mechanical active galactic nucleus (AGN) feedback. The duty cycle and detailed physics of the radio-mode AGN feedback is still a matter of debate. To address these questions, we present 1-2 GHz Karl G. Jansky Very Large Array (VLA) radio observations of a sample of the 42 nearest optically and X-ray brightest early-type galaxies. We detect radio emission in 41/42 galaxies. However, the galaxy without a radio source, NGC 499, has recently been detected at lower frequencies by the Low-Frequency Array (LOFAR). Furthermore, 27/42 galaxies in our sample host extended radio structures and 34/42 sources show environmental interactions in the form of X-ray cavities. We find a significant correlation between the radio flux density and the largest linear size of the radio emission and between the radio power and the luminosity of the central X-ray point-source. The central radio spectral indices of the galaxies span a wide range of values, with the majority of the systems having steep spectra and the rest flat spectra. These results are consistent with AGN activity, where the central radio sources are mostly switched on, thus the duty cycle is very high. 7/14 galaxies with point-like radio emission (Fanaroff-Riley Class 0; FR 0) also show X-ray cavities indicating that, despite the lack of extended radio structures at 1-2 GHz, these AGN do launch jets capable of inflating lobes and cavities.


title: Creep effects on the Campbell response in type II superconductors

abstract: Applying the strong pinning formalism to the mixed state of a type II superconductor, we study the effect of thermal fluctuations (or creep) on the penetration of an ac magnetic field as quantified by the so-called Campbell length $\lambda_\textrm{C}$. Within strong pinning theory, vortices get pinned by individual defects, with the jumps in the pinning energy ($\Delta e_\mathrm{pin}$) and force ($\Delta f_\mathrm{pin}$) between bistable pinned and free states quantifying the pinning process. We find that the evolution of the Campbell length $\lambda_{\rm C}(t)$ as a function of time $t$ is the result of two competing effects, the change in the force jumps $\Delta f_\mathrm{pin}(t)$ and a change in the trapping area $S_\mathrm{trap}(t)$ of vortices; the latter describes the area around the defect where a nearby vortex gets and remains trapped. Contrary to naive expectation, we find that during the decay of the critical state in a zero-field cooled (ZFC) experiment, the Campbell length $\lambda_{\rm C}(t)$ is usually nonmonotonic, first decreasing with time $t$ and then increasing for long waiting times. Field cooled (FC) experiments exhibit hysteretic effects in $\lambda_\textrm{C}$; relaxation then turns out to be predominantly monotonic, but its magnitude and direction depends on the specific phase of the cooling--heating cycle. Furthermore, when approaching equilibrium, the Campbell length relaxes to a finite value, different from the persistent current which vanishes at long waiting times $t$, e.g., above the irreversibility line. Finally, measuring the Campbell length $\lambda_\textrm{C}(t)$ for different states, zero-field cooled, field cooled, and relaxed, as a function of different waiting times $t$ and temperatures $T,$ allows to "spectroscopyse" the pinning potential of the defects.


title: Consistency conditions and primordial black holes in single field inflation

abstract: We discuss new consistency relations for single field models of inflation capable of generating primordial black holes (PBH), and their observational implications for CMB $\mu$-space distortions. These inflationary models include a short period of non-attractor evolution: the scale-dependent profile of curvature perturbation is characterized by a pronounced dip, followed by a rapid growth leading to a peak responsible for PBH formation. We investigate the squeezed and the collapsed limits of three and four point functions of curvature perturbation around the dip, showing that they satisfy consistency relations connecting their values to the total amplification of the curvature spectrum, and to the duration of the non-attractor era. Moreover, the corresponding non-Gaussian parameters are scale-dependent in proximity of the dip, with features that again depend on the amplification of the spectrum. For typical PBH scenarios requiring an order ${\cal O}(10^7)$ enhancement of the spectrum from large towards small scales, we generally find values $f_{\rm NL}^{\rm sq}\,=\,{\cal O}(10)$ and $\tau_{\rm NL}^{\rm col}\,=\,{\cal O}(10^3)$ in a range of scales that can be probed by CMB $\mu$-space distortions. Using these consistency relations, we carefully analyze how the scale-dependence of non-Gaussian parameters leads to characteristic features in $\langle \mu T \rangle$ and $\langle \mu \mu \rangle$ correlators, providing distinctive probes of inflationary PBH scenarios that can be tested using well-understood CMB physics.


title: Three-Dimensional Real Space Invariants, Obstructed Atomic Insulators and A New Principle for Active Catalytic Sites

abstract: Topologically trivial insulators come in two kinds: atomic, where the Wannier charge centers (WCCs) are localized on the atoms, and obstructed atomic, where the WCCs are located away from the atoms. The latter, which can exhibit interesting surface states and possibly have much larger band gaps than the topological insulators, have so far not been classified in three-dimensional (3D) crystalline materials. In this paper, we developed the 3D real space invariants (RSIs) for the 1651 Shubnikov space groups with the spin-orbit coupling and provide the full classification of 3D obstructed atomic insulators (OAIs) by the RSIs. We then apply the theory to the entire database of materials on the Topological Quantum Chemistry website and Topological Magnetic Materials website, obtaining all the OAIs so far existing in nature. We find that, out of the 34013 paramagnetic and 296 magnetic topologically trivial insulators, there are 3383 paramagnetic and 30 magnetic OAIs. All of them present a filling anomaly under certain open-boundary conditions and exhibit obstructed surface states (OSSs). We further refine the atomic insulator concept to obtain the orbital-selected OAIs (OOAIs), where the WCC of the system is located at a Wyckoff position occupied by an atom but forms a symmetric representation that does not belong to the outer-shell electrons of the given atom. In such a way, we obtain a further 121 OOAIs. Furthermore, we analyze the catalytic properties of one of the OAIs in a "proof of principle" experiment. By using the single crystal of 2H-MoS2 as a hydrogen evolution catalyst, we directly prove that the catalytic activities arise from the surfaces with OSSs. Additional potential applications of the 3D RSIs and OAIs in, for example, electrochemistry, asymmetric catalysis, superconductivity, and Josephson diode will be discussed.


title: Hamiltonian Dynamics with Non-Newtonian Momentum for Rapid Sampling

abstract: Sampling from an unnormalized probability distribution is a fundamental problem in machine learning with applications including Bayesian modeling, latent factor inference, and energy-based model training. After decades of research, variations of MCMC remain the default approach to sampling despite slow convergence. Auxiliary neural models can learn to speed up MCMC, but the overhead for training the extra model can be prohibitive. We propose a fundamentally different approach to this problem via a new Hamiltonian dynamics with a non-Newtonian momentum. In contrast to MCMC approaches like Hamiltonian Monte Carlo, no stochastic step is required. Instead, the proposed deterministic dynamics in an extended state space exactly sample the target distribution, specified by an energy function, under an assumption of ergodicity. Alternatively, the dynamics can be interpreted as a normalizing flow that samples a specified energy model without training. The proposed Energy Sampling Hamiltonian (ESH) dynamics have a simple form that can be solved with existing ODE solvers, but we derive a specialized solver that exhibits much better performance. ESH dynamics converge faster than their MCMC competitors enabling faster, more stable training of neural network energy models.


title: Obstructed surface states as the origin of catalytic activity in inorganic heterogeneous catalysts

abstract: The discovery of new catalysts that are efficient, sustainable, and low-cost is a major research endeavor for many industrial chemical processes. This requires an understanding and determination of the catalytic origins for the given catalysts, which still remains a challenge. Here we describe a novel method to identify new catalysts based on searching for crystalline symmetry-protected obstructed atomic insulators (OAIs) that have metallic surface states on otherwise semiconducting or insulating compounds. The Wannier charge centers in OAIs are pinned by symmetries at some empty Wyckoff positions so that surfaces that accommodate these sites are guaranteed to have metallic obstructed surface states (OSSs). Beyond the well-studied 2H-MoS2, we further verified our theory on the catalysts, 2H-MoTe2, and 1T'-MoTe2, whose catalytic active sites are consistent with our calculations of obstructed Wannier charge centers (OWCCs) and OSSs. In addition, we have predicted the location of catalytic active sites and confirmed these predictions by exploring the hydrogen evolution reaction on NiPS3 bulk single crystals, which we find to be one of the most promising new catalysts with high activity and, moreover, of low cost. Most importantly, we successfully identified several high-efficient catalysts just by considering the number of OWCCs and the crystal symmetry of the OAIs. Using the real space invariant theory and high-throughput computational methods applied to a database of 34013 topologically trivial insulators, we have identified 1788 unique OAIs (3383 ICSDs), of which 465 are potential high-quality catalysts for heterogeneous reactions. The Miller indices of the active surfaces are also obtained. Our new methodology will facilitate and accelerate the discovery of new catalysts for a wide range of heterogeneous redox reactions, where sustainability, toxicity, and cost must be considered.


title: The Unanticipated Phenomenology of the Blazar PKS~2131$-$021: A Unique Super-Massive Black Hole Binary Candidate

abstract: Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS~2131$-$021, which strongly suggest an SMBHB with an orbital separation of $\sim 0.001-0.01$ pc. Our 45.1-year radio light curve shows two epochs of strong sinusoidal variation with the same period and phase to within $<2\%$ and $\sim 10\%$, respectively, straddling a 20-year period when this variation was absent. Our simulated light curves accurately reproduce the ``red noise'' of this object, and Lomb-Scargle, weighted wavelet Z-transform, and least-squares sine wave analyses demonstrate conclusively, at the $4.6\sigma$ significance level, that the periodicity in this object is not due to random fluctuations in flux density. The observed period translates to $2.082\pm 0.003$ years in the rest frame at the $z=1.285$ redshift of PKS~2131$-$021. The periodic variation in PKS~2131$-$021 is remarkably sinusoidal. We present a model in which orbital motion, combined with the strong Doppler boosting of the approaching relativistic jet, produces a sine-wave modulation in the flux density which easily fits the observations. Given the rapidly-developing field of gravitational wave experiments with pulsar timing arrays, closer counterparts to PKS~2131$-$021 and searches using the techniques we have developed are strongly motivated. These results constitute a compelling demonstration that the phenomenology, not the theory, must provide the lead in this field.


title: Long-Lived Light Mediators from Higgs boson Decay at HL-LHC, FCC-hh and a Proposal of Dedicated LLP Detectors for FCC-hh

abstract: We study the pair production of the long-lived mediator particles from the decay of the SM Higgs boson and their subsequent decay into standard model particles. We compute the projected sensitivity, both model-independently and with a minimal model, of using the muon spectrometer of the CMS detector at the HL-LHC experiment for ggF, VBF, and Vh production modes of the Higgs boson and various decay modes of the mediator particle, along with dedicated detectors for LLP searches like CODEX-b and MATHUSLA. Subsequently, we study the improvement with the FCC-hh detector at the 100\,TeV collider experiment for such long-lived mediators, again focusing on the muon spectrometer. We propose dedicated LLP detector designs for the 100\,TeV collider experiment, DELIGHT (\textbf{De}tector for \textbf{l}ong-l\textbf{i}ved particles at hi\textbf{gh} energy of 100\,\textbf{T}eV), and study their sensitivities.


title: No second law of entanglement manipulation after all

abstract: Many fruitful analogies have emerged between the theories of quantum entanglement and thermodynamics, motivating the pursuit of an axiomatic description of entanglement akin to the laws of thermodynamics. A long-standing open problem has been to establish a true second law of entanglement, and in particular a unique function which governs all transformations between entangled systems, mirroring the role of entropy in thermodynamics. Contrary to previous promising evidence, here we show that this is impossible, and no direct counterpart to the second law of thermodynamics can be established. This is accomplished by demonstrating the irreversibility of entanglement theory from first principles -- assuming only the most general microscopic physical constraints of entanglement manipulation, we show that entanglement theory is irreversible under all non-entangling transformations. We furthermore rule out reversibility without significant entanglement expenditure, showing that reversible entanglement transformations require the generation of macroscopically large amounts of entanglement according to certain measures. Our results not only reveal fundamental differences between quantum entanglement transformations and thermodynamic processes, but also showcase a unique property of entanglement which distinguishes it from other known quantum resources.

categories: quant-ph cond-mat.stat-mech math-ph math.MP


title: What really makes an accretion disc MAD

abstract: Magnetically arrested accretion discs (MADs) around black holes (BH) have the potential to stimulate the production of powerful jets and account for recent ultra-high-resolution observations of BH environments. Their main properties are usually attributed to the accumulation of dynamically significant net magnetic (vertical) flux throughout the arrested region, which is then regulated by interchange instabilities. Here we propose instead that it is mainly a dynamically important {\it toroidal} field -- the result of dynamo action triggered by the significant but still relatively weak vertical field -- that defines and regulates the properties of MADs. We suggest that rapid convection-like instabilities, involving interchange of toroidal flux tubes and operating concurrently with the magnetorotatonal instability (MRI), can regulate the structure of the disc and the escape of net flux. We generalize the convective stability criteria and disc structure equations to include the effects of a strong toroidal field and show that convective flows could be driven towards two distinct marginally stable states, one of which we associate with MADs. We confirm the plausibility of our theoretical model by comparing its quantitative predictions to simulations of both MAD and SANE (strongly magnetized but not "arrested") discs, and suggest a set of criteria that could help to distinguish MADs from other accretion states. Contrary to previous claims in the literature, we argue that MRI is not suppressed in MADs and is probably responsible for the existence of the strong toroidal field.


title: Mirror Symmetry for Five-Parameter Hulek-Verrill Manifolds

abstract: We study the mirrors of five-parameter Calabi-Yau threefolds first studied by Hulek and Verrill in the context of observed modular behaviour of the zeta functions for Calabi-Yau manifolds. Toric geometry allows for a simple explicit construction of these mirrors, which turn out to be familiar manifolds. These are elliptically fibred in multiple ways. By studying the singular fibres, we are able to identify the rational curves of low degree on the mirror manifolds. This verifies the mirror symmetry prediction obtained by studying the mirror map near large complex structure points. We undertake also an extensive study of the periods of the Hulek-Verrill manifolds and their monodromies. On the mirror, we compute the genus-zero and -one instanton numbers, which are labelled by 5 indices, as $h^{1,1}=5$. There is an obvious permutation symmetry on these indices, but in addition there is a surprising repetition of values. We trace this back to an $S_{6}$ symmetry made manifest by certain constructions of the complex structure moduli space of the Hulek-Verrill manifold. Among other consequences, we see in this way that the moduli space has six large complex structure limits. It is the freedom to expand the prepotential about any one of these points that leads to this symmetry in the instanton numbers. An intriguing fact is that the group that acts on the instanton numbers is larger than $S_6$ and is in fact an infinite hyperbolic Coxeter group, that we study. The group orbits have a 'web' structure, and with certain qualifications the instanton numbers are only nonzero if they belong to what we term 'positive webs'. This structure has consequences for instanton numbers at all genera.


title: NECOLA: Towards a Universal Field-level Cosmological Emulator

abstract: We train convolutional neural networks to correct the output of fast and approximate N-body simulations at the field level. Our model, Neural Enhanced COLA --NECOLA--, takes as input a snapshot generated by the computationally efficient COLA code and corrects the positions of the cold dark matter particles to match the results of full N-body Quijote simulations. We quantify the accuracy of the network using several summary statistics, and find that NECOLA can reproduce the results of the full N-body simulations with sub-percent accuracy down to $k\simeq1~h{\rm Mpc}^{-1}$. Furthermore, the model, that was trained on simulations with a fixed value of the cosmological parameters, is also able to correct the output of COLA simulations with different values of $\Omega_{\rm m}$, $\Omega_{\rm b}$, $h$, $n_s$, $\sigma_8$, $w$, and $M_\nu$ with very high accuracy: the power spectrum and the cross-correlation coefficients are within $\simeq1\%$ down to $k=1~h{\rm Mpc}^{-1}$. Our results indicate that the correction to the power spectrum from fast/approximate simulations or field-level perturbation theory is rather universal. Our model represents a first step towards the development of a fast field-level emulator to sample not only primordial mode amplitudes and phases, but also the parameter space defined by the values of the cosmological parameters.


title: The Effective Vector Boson Approximation in High-Energy Muon Collisions

abstract: Due to the inclination for forward gauge radiation, lepton colliders beyond a few TeV are effectively electroweak (EW) boson colliders, suggesting the treatment of EW bosons as constituents of high-energy leptons. In the context of a muon collider, we revisit the validity of $W$ and $Z$ parton distribution functions (PDFs) at leading order in $2\to n$ process. We systematically investigate universal and quasi-universal power-law and logarithmic corrections that arise when deriving (polarized) weak boson PDFs in the collinear limit. We go on to survey a multitude of $2\to n$ processes at $\sqrt{s}=2-30$ TeV via polarized and unpolarized EW boson fusion/scattering. To conduct this study, we report a public implementation of the Effective $W/Z$ and Weizs\"acker-Williams Approximations, which we collectively call the Effective Vector Boson Approximation, into the Monte Carlo event generator \texttt{MadGraph5\_aMC@NLO}. This implementation lays the groundwork for developing matrix-element matching prescriptions involving EW parton showers and renormalized EW PDFs. To further with this agenda, we give recommendations on using $W/Z$ PDFs.


title: Time-Resolved Photometry of the High-Energy Radiation of M Dwarfs with the Star-Planet Activity Research CubeSat (SPARCS)

abstract: Know thy star, know thy planet,... especially in the ultraviolet (UV). Over the past decade, that motto has grown from mere wish to necessity in the M dwarf regime, given that the intense and highly variable UV radiation from these stars is suspected of strongly impacting their planets' habitability and atmospheric loss. This has led to the development of the Star-Planet Activity Research CubeSat (SPARCS), a NASA-funded 6U CubeSat observatory fully devoted to the photometric monitoring of the UV flaring of M dwarfs hosting potentially habitable planets. The SPARCS science imaging system uses a 9-cm telescope that feeds two delta-doped UV-optimized CCDs through a dichroic beam splitter, enabling simultaneous monitoring of a target field in the near-UV and far-UV. A dedicated onboard payload processor manages science observations and performs near-real time image processing to sustain an autonomous dynamic exposure control algorithm needed to mitigate pixel saturation during flaring events. The mission is currently half-way into its development phase. We present an overview of the mission's science drivers and its expected contribution to our understanding of star-planet interactions. We also present the expected performance of the autonomous dynamic exposure control algorithm, a first-of-its-kind on board a space-based stellar astrophysics observatory.

categories: astro-ph.IM astro-ph.EP astro-ph.HE astro-ph.SR


title: Panoptic 3D Scene Reconstruction From a Single RGB Image

abstract: Understanding 3D scenes from a single image is fundamental to a wide variety of tasks, such as for robotics, motion planning, or augmented reality. Existing works in 3D perception from a single RGB image tend to focus on geometric reconstruction only, or geometric reconstruction with semantic segmentation or instance segmentation. Inspired by 2D panoptic segmentation, we propose to unify the tasks of geometric reconstruction, 3D semantic segmentation, and 3D instance segmentation into the task of panoptic 3D scene reconstruction - from a single RGB image, predicting the complete geometric reconstruction of the scene in the camera frustum of the image, along with semantic and instance segmentations. We thus propose a new approach for holistic 3D scene understanding from a single RGB image which learns to lift and propagate 2D features from an input image to a 3D volumetric scene representation. We demonstrate that this holistic view of joint scene reconstruction, semantic, and instance segmentation is beneficial over treating the tasks independently, thus outperforming alternative approaches.


title: Autonomous Attack Mitigation for Industrial Control Systems

abstract: Defending computer networks from cyber attack requires timely responses to alerts and threat intelligence. Decisions about how to respond involve coordinating actions across multiple nodes based on imperfect indicators of compromise while minimizing disruptions to network operations. Currently, playbooks are used to automate portions of a response process, but often leave complex decision-making to a human analyst. In this work, we present a deep reinforcement learning approach to autonomous response and recovery in large industrial control networks. We propose an attention-based neural architecture that is flexible to the size of the network under protection. To train and evaluate the autonomous defender agent, we present an industrial control network simulation environment suitable for reinforcement learning. Experiments show that the learned agent can effectively mitigate advanced attacks that progress with few observable signals over several months before execution. The proposed deep reinforcement learning approach outperforms a fully automated playbook method in simulation, taking less disruptive actions while also defending more nodes on the network. The learned policy is also more robust to changes in attacker behavior than playbook approaches.


title: Optical demultiplexing of fractal-structured beams in turbulent atmospheric environments

abstract: When information is spatially repeated in self-similar fractal beam patterns, only a portion of the diffracted beam is needed to reconstruct the kernel data. What is unique to a fractal-encoding scheme is that the image demultiplexing process can be, to a first approximation, easily performed optically. In prior work, we experimentally and numerically study fractal-encoded optical beams and their mid- and far-field propagation without added turbulence. Here, we present preliminary simulations of fractal-encoded beams with high turbulence ($C_n^2 \geq 10^{-14}$ m$^{-2/3}$) where we achieve respectable bit error rates of $10^{-3}$. These results are impressive given that: data with low fractal orders is shown, simple threshold-algorithms are used (i.e., no machine learning), and only a third of the beam, off-axis, is needed. More robust channel encoding is associated with increased fractal orders, larger collection areas, and higher kernel singular value decomposition entropy.


title: On the Frequency Bias of Generative Models

abstract: The key objective of Generative Adversarial Networks (GANs) is to generate new data with the same statistics as the provided training data. However, multiple recent works show that state-of-the-art architectures yet struggle to achieve this goal. In particular, they report an elevated amount of high frequencies in the spectral statistics which makes it straightforward to distinguish real and generated images. Explanations for this phenomenon are controversial: While most works attribute the artifacts to the generator, other works point to the discriminator. We take a sober look at those explanations and provide insights on what makes proposed measures against high-frequency artifacts effective. To achieve this, we first independently assess the architectures of both the generator and discriminator and investigate if they exhibit a frequency bias that makes learning the distribution of high-frequency content particularly problematic. Based on these experiments, we make the following four observations: 1) Different upsampling operations bias the generator towards different spectral properties. 2) Checkerboard artifacts introduced by upsampling cannot explain the spectral discrepancies alone as the generator is able to compensate for these artifacts. 3) The discriminator does not struggle with detecting high frequencies per se but rather struggles with frequencies of low magnitude. 4) The downsampling operations in the discriminator can impair the quality of the training signal it provides. In light of these findings, we analyze proposed measures against high-frequency artifacts in state-of-the-art GAN training but find that none of the existing approaches can fully resolve spectral artifacts yet. Our results suggest that there is great potential in improving the discriminator and that this could be key to match the distribution of the training data more closely.


title: A prescreening method for variational quantum state eigensolver

abstract: The development of Fault-Tolerant Quantum Computer (FTQC) gradually raises a possibility to implement the Quantum Phase Estimation (QPE) algorithm. However, QPE works only for normalized systems. This requires the minimum and maximum of eigenvalues of a Hamiltonian. Variational Quantum Eigensolver (VQE) well developed in the Noisy Intermediate Scale Quantum (NISQ) era is necessary for preparing the initial eigenvectors that close to the exact states. In this paper, we propose a method to derive all of the states with high accuracy by using the Variational Quantum State Eigensolver (VQSE) and Subspace-Search VQE (SSVQE) methods. We show that by using the VQSE and the SSVQE prescreening methods, we can derive all of the hydrogen molecules states correctly.


title: Effective Resistance for Pandemics: Mobility Network Sparsification for High-Fidelity Epidemic Simulation

abstract: Network science has increasingly become central to the field of epidemiology and our ability to respond to infectious disease threats. However, many networks derived from modern datasets are not just large, but dense, with a high ratio of edges to nodes. This includes human mobility networks where most locations have a large number of links to many other locations. Simulating large-scale epidemics requires substantial computational resources and in many cases is practically infeasible. One way to reduce the computational cost of simulating epidemics on these networks is sparsification, where a representative subset of edges is selected based on some measure of their importance. We test several sparsification strategies, ranging from naive thresholding to random sampling of edges, on mobility data from the U.S. Following recent work in computer science, we find that the most accurate approach uses the effective resistances of edges, which prioritizes edges that are the only efficient way to travel between their endpoints. The resulting sparse network preserves many aspects of the behavior of an SIR model, including both global quantities, like the epidemic size, and local details of stochastic events, including the probability each node becomes infected and its distribution of arrival times. This holds even when the sparse network preserves fewer than $10\%$ of the edges of the original network. In addition to its practical utility, this method helps illuminate which links of a weighted, undirected network are most important to disease spread.


title: Unified 3D Mesh Recovery of Humans and Animals by Learning Animal Exercise

abstract: We propose an end-to-end unified 3D mesh recovery of humans and quadruped animals trained in a weakly-supervised way. Unlike recent work focusing on a single target class only, we aim to recover 3D mesh of broader classes with a single multi-task model. However, there exists no dataset that can directly enable multi-task learning due to the absence of both human and animal annotations for a single object, e.g., a human image does not have animal pose annotations; thus, we have to devise a new way to exploit heterogeneous datasets. To make the unstable disjoint multi-task learning jointly trainable, we propose to exploit the morphological similarity between humans and animals, motivated by animal exercise where humans imitate animal poses. We realize the morphological similarity by semantic correspondences, called sub-keypoint, which enables joint training of human and animal mesh regression branches. Besides, we propose class-sensitive regularization methods to avoid a mean-shape bias and to improve the distinctiveness across multi-classes. Our method performs favorably against recent uni-modal models on various human and animal datasets while being far more compact.


title: Simulation Study of the Space Charge Limit in Heavy-ion Synchrotrons

abstract: The SIS100 synchrotron as a part of the new FAIR accelerator facility at GSI should be operated at the "space charge limit" for light and heavy ion beams. Beam losses due to space charge induced resonance crossing should not exceed a few percent during a full cycle. The recent advances in the performance of particle tracking tools with self-consistent solvers for the 3D space charge forces now allow us to reliably identify low-loss areas in tune space, considering the full 1s (160'000 turns) accumulation plateau in SIS100. A realistic magnet error model, extracted from precise bench measurements of the SIS100 main dipole and quadrupole magnets, is included in the simulations. Previously such beam dynamics simulations required non-self-consistent space charge models. By comparing to the self-consistent simulations results we are now able to demonstrate that the predictions from such faster space charge models can be used to identify low-loss regions with sufficient accuracy. The findings are applied by identifying a low-loss working point region in SIS100 for the design FAIR beam parameters. The bunch intensity at the space charge limit is determined. Several counter-measures to space charge are proposed to enlarge the low-loss area and to further increase the space charge limit.


title: Slapping Cats, Bopping Heads, and Oreo Shakes: Understanding Indicators of Virality in TikTok Short Videos

abstract: Short videos have become one of the leading media used by younger generations to express themselves online and thus a driving force in shaping online culture. In this context, TikTok has emerged as a platform where viral videos are often posted first. In this paper, we study what elements of short videos posted on TikTok contribute to their virality. We apply a mixed-method approach to develop a codebook and identify important virality features. We do so vis-\`a-vis three research hypotheses; namely, that: 1) the video content, 2) TikTok's recommendation algorithm, and 3) the popularity of the video creator contribute to virality. We collect and label a dataset of 400 TikTok videos and train classifiers to help us identify the features that influence virality the most. While the number of followers is the most powerful predictor, close-up and medium-shot scales also play an essential role. So does the lifespan of the video, the presence of text, and the point of view. Our research highlights the characteristics that distinguish viral from non-viral TikTok videos, laying the groundwork for developing additional approaches to create more engaging online content and proactively identify possibly risky content that is likely to reach a large audience.


title: Bootstrapping (D, D) Conformal Matter

abstract: We use the numerical conformal bootstrap to study six-dimensional $\mathcal{N}=(1,0)$ superconformal field theories with flavor symmetry $\mathfrak{so}_{4k}$. We present evidence that minimal $(D_k, D_k)$ conformal matter saturates the unitarity bounds for arbitrary $k$. Furthermore, using the extremal-functional method, we check that the chiral-ring relations are correctly reproduced, extract the anomalous dimensions of low-lying long superconformal multiplets, and find hints for novel OPE selection rules involving type-$\mathcal{B}$ multiplets.


title: Testing the Potential for Radio Variability in Disks around T Tauri Stars with Observations and Chemical Modeling

abstract: A multiwavelength observing campaign of the T Tauri star (TTS) GM Aur was undertaken in 2019 December that obtained Swift X-ray and NUV fluxes, HST NUV spectra, LCOGT ugri and TESS photometry, CHIRON Halpha spectra, ALMA 13CO and C18O line fluxes, and VLA 3 cm continuum fluxes taken contemporaneously over one month. The X-ray to optical observations were presented previously. Here we present the ALMA and VLA data and make comparisons to GM Aur's accretion and X-ray properties. We report no variability in the observed millimeter CO emission. Using disk chemistry models, we show that the magnitude of the changes seen in the FUV luminosity of GM Aur could lead to variation of up to ~6% in CO line emission and changes in the X-ray luminosity could lead to larger changes of ~25%. However, the FUV and X-ray luminosity increases must last at least 100 years in order to induce changes, which seems implausible in the TTS stage; also, these changes would be too small to be detectable by ALMA. We report no variability in the 3 cm emission observed by the VLA, showing that changes of less than a factor of ~3 in the accretion rates of TTSs do not lead to detectable changes in the mass-loss rate traced by the jet at centimeter wavelengths. We conclude that typically seen changes in the FUV and X-ray luminosities of TTSs do not lead to observable changes in millimeter CO line emission or jet centimeter continuum emission.


title: Understanding the Use of e-Prints on Reddit and 4chan's Politically Incorrect Board

abstract: The dissemination and reach of scientific knowledge have increased at a blistering pace. In this context, e-Print servers have played a central role by providing scientists with a rapid and open mechanism for disseminating research without waiting for the (lengthy) peer review process. While helping the scientific community in several ways, e-Print servers also provide scientific communicators and the general public with access to a wealth of knowledge without paying hefty subscription fees. This motivates us to study how e-Prints are positioned within Web community discussions. In this paper, we analyze data from two Web communities: 14 years of Reddit data and over 4 from 4chan's Politically Incorrect board. Our findings highlight the presence of e-Prints in both science-enthusiast and general-audience communities. Real-world events and distinct factors influence the e-Prints people's discussions; e.g., there was a surge of COVID-19-related research publications during the early months of the outbreak and increased references to e-Prints in online discussions. Text in e-Prints and in online discussions referencing them has a low similarity, suggesting that the latter are not exclusively talking about the findings in the former. Further, our analysis of a sample of threads highlights: 1) misinterpretation and generalization of research findings, 2) early research findings being amplified as a source for future predictions, and 3) questioning findings from a pseudoscientific e-Print. Overall, our work emphasizes the need to quickly and effectively validate non-peer-reviewed e-Prints that get substantial press/social media coverage to help mitigate wrongful interpretations of scientific outputs.


title: On Johnson's "sufficientness" postulates for features-sampling models

abstract: In the 1920's, the English philosopher W.E. Johnson introduced a characterization of the symmetric Dirichlet prior distribution in terms of its predictive distribution. This is typically referred to as Johnson's "sufficientness" postulate, and it has been the subject of many contributions in Bayesian statistics, leading to predictive characterization for infinite-dimensional generalizations of the Dirichlet distribution, i.e. species-sampling models. In this paper, we review "sufficientness" postulates for species-sampling models, and then investigate analogous predictive characterizations for the more general features-sampling models. In particular, we present a "sufficientness" postulate for a class of features-sampling models referred to as Scaled Processes (SPs), and then discuss analogous characterizations in the general setup of features-sampling models.


title: Symmetries and local transformations of translationally invariant Matrix Product States

abstract: We determine the local symmetries and local transformation properties of translationally invariant matrix product states (MPS). We focus on physical dimension $d=2$ and bond dimension $D=3$ and use the procedure introduced in D. Sauerwein et al., Phys. Rev. Lett. 123, 170504 (2019) to determine all (including non--global) symmetries of those states. We identify and classify the stochastic local transformations (SLOCC) that are allowed among MPS. We scrutinize two very distinct sets of MPS and show the big diversity (also compared to the case $D=2$) occurring in both, their symmetries and the possible SLOCC transformations. These results reflect the variety of local properties of MPS, even if restricted to translationally invariant states with low bond dimension. Finally, we show that states with non-trivial local symmetries are of measure zero for $d = 2$ and $D > 3$.


title: Perturb-and-max-product: Sampling and learning in discrete energy-based models

abstract: Perturb-and-MAP offers an elegant approach to approximately sample from a energy-based model (EBM) by computing the maximum-a-posteriori (MAP) configuration of a perturbed version of the model. Sampling in turn enables learning. However, this line of research has been hindered by the general intractability of the MAP computation. Very few works venture outside tractable models, and when they do, they use linear programming approaches, which as we will show, have several limitations. In this work we present perturb-and-max-product (PMP), a parallel and scalable mechanism for sampling and learning in discrete EBMs. Models can be arbitrary as long as they are built using tractable factors. We show that (a) for Ising models, PMP is orders of magnitude faster than Gibbs and Gibbs-with-Gradients (GWG) at learning and generating samples of similar or better quality; (b) PMP is able to learn and sample from RBMs; (c) in a large, entangled graphical model in which Gibbs and GWG fail to mix, PMP succeeds.


title: Smart sensors network for accurate indirect heat accounting in apartment buildings

abstract: A new method for accurate indirect heat accounting in apartment buildings has been recently developed by the Centre Suisse d'Electronique et de Microtechnique (CSEM). It is based on a data driven approach aimed to the smart networking of any type of indirect heat allocation devices, which can provide, for each heat delivery point of an apartment building, measurements or estimations of the temperature difference between the heat transfer fluid and the indoor environment. The analysis of the data gathered from the devices installed on the heating bodies, together with the measurements of the overall building heat consumption provided by direct heat metering, allows the evaluation of the characteristic thermal model parameters of heating bodies at actual installation and working conditions. Thus overcoming the negative impact on accuracy of conventional indirect heat accounting due to off-design operation, in which these measurement systems normally operate. The method has been tested on conventional heat cost allocators (HCA), and on innovative smart radiator thermostatic valves developed by CSEM. The evaluations were carried out at the centralized heating system mock-up of the Istituto Nazionale di Ricerca Metrologica (INRIM), and also in a real building in Neuchatel, Switzerland. The method has proven to be an effective tool to improve the accuracy of indirect heat metering systems; compared to conventional HCA systems, the error on the individual heating bill is reduced by 20% to 50%.


title: Joint Species Distribution Modeling of Percentage Cover Data with Exclusive Competition for Space

abstract: Joint species distribution models (JSDM) are among the most important statistical tools in community ecology. They are routinely used for inference and various prediction tasks, such as to build species distribution maps or biomass estimation over spatial areas. Existing JSDM's cannot, however, model mutual exclusion between species, which may happen in some species groups, such as mosses in the bottom layer of a peatland site. We tackle this deficiency in the context of modeling plant percentage cover data, where mutual exclusion arises from limited growing space and competition for light. We propose a hierarchical JSDM where multivariate latent Gaussian variable model describes species' niche preferences and Dirichlet-Multinomial distribution models the observation process and exclusive competition for space between species. We use both stationary and non-stationary multivariate Gaussian processes to model residual phenomena. We also propose a decision theoretic model comparison and validation approach to assess the goodness of JSDMs in four different types of predictive tasks. We apply our models and methods to a case study on modeling vegetation cover in a boreal peatland. Our results show that ignoring the interspecific interactions and competition for space significantly reduces models' predictive performance and leads to biased estimates for total percentage cover both for individual species and over all species combined. A model's relative predictive performance also depends on the model comparison methods highlighting that model comparison and assessment should resemble the true predictive task. Our results also demonstrate that the proposed joint species distribution model can be used to simultaneously infer interspecific correlations in niche preference as well as mutual exclusive competition for space and through that provide novel insight into ecological research.


title: Automatic ultrasound vessel segmentation with deep spatiotemporal context learning

abstract: Accurate, real-time segmentation of vessel structures in ultrasound image sequences can aid in the measurement of lumen diameters and assessment of vascular diseases. This, however, remains a challenging task, particularly for extremely small vessels that are difficult to visualize. We propose to leverage the rich spatiotemporal context available in ultrasound to improve segmentation of small-scale lower-extremity arterial vasculature. We describe efficient deep learning methods that incorporate temporal, spatial, and feature-aware contextual embeddings at multiple resolution scales while jointly utilizing information from B-mode and Color Doppler signals. Evaluating on femoral and tibial artery scans performed on healthy subjects by an expert ultrasonographer, and comparing to consensus expert ground-truth annotations of inner lumen boundaries, we demonstrate real-time segmentation using the context-aware models and show that they significantly outperform comparable baseline approaches.


title: Gravitational Decoherence: A Thematic Overview

abstract: Gravitational decoherence (GD) refers to the effects of gravity in actuating the classical appearance of a quantum system. Because the underlying processes involve issues in general relativity (GR), quantum field theory (QFT) and quantum information, GD has fundamental theoretical significance. There is a great variety of GD models, many of them involving physics that diverge from GR and/or QFT. This overview has two specific goals along one central theme: (i) present theories of GD based on GR and QFT and explore their experimental predictions; (ii) place other theories of GD under the scrutiny of GR and QFT, and point out their theoretical differences. We also describe how GD experiments in space in the coming decades can provide evidences at two levels: a) discriminate alternative quantum theories and non-GR theories; b) discern whether gravity is a fundamental or an effective theory.


title: V1405 Cas Nova 2021: About the nature of a multi-maxima nova

abstract: CONTEXT: The Nova V1405 Cas, which erupted on 18 March 2021, showed a pre-maximum phase of almost two months. After the main maximum, during several month the nova only weakened slightly and showed clear fluctuations in brightness. AIMS: The aim of this work is to investigate the relationship between brightness and half-width (FWHM) and equivalent width (EW) of the H-alpha emission line, and from this to infer possible causes for the abnormal behavior of the Nova. METHODS: Magnitudes and spectra from online databases were used, the H-alpha line of the spectra was measured and 2-day normal values were formed in order to establish a temporal comparison. RESULTS: After the main maximum, during the period investigated by JD 2459292 to 2459305 at intervals of 2-5 weeks secondary maxima with different characteristics. The maxima of the equivalent width were 6-14 days after the respective brightness maxima. The maxima are interpreted as subsequent outbreaks of the nova. Possibly these are subsequent weaker eruptions of the nova or other mechanisms that lead to an increase in density in the envelope.


title: Cosmological implications of EW vacuum instability: constraints on the Higgs-curvature coupling from inflation

abstract: The current experimentally measured parameters of the Standard Model (SM) suggest that our Universe lies in a metastable electroweak vacuum, where the Higgs field is prone to vacuum decay to a lower state with catastrophic consequences. Our measurements dictate that such an event has not taken place yet, despite the many different mechanisms that could have triggered it in our past light-cone. The focus of our work has been to calculate the probability of the false vacuum to decay during the period of inflation and use it to constrain the last unknown renormalisable SM parameter $\xi$, which couples the Higgs field with space-time curvature. More specifically, we derived lower $\xi$-bounds from vacuum stability in three inflationary models: quadratic and quartic chaotic inflation, and Starobinsky-like power-law inflation. We also took the time-dependence of the Hubble rate into account both in the geometry of our past light-cone and in the Higgs effective potential, which is approximated with three-loop renormalisation group improvement supplemented with one-loop curvature corrections.


title: Concerted Rolling and Membrane Penetration Revealed by Atomistic Simulations of Antimicrobial Peptides

abstract: Short peptides with antimicrobial activity have therapeutic potential for treating bacterial infections. Mechanisms of actions for antimicrobial peptides require binding the biological membrane of their target, which often represents a key mechanistic step. A multitude of data-driven approaches have been developed to predict potential antimicrobial peptide sequences; however, these methods are usually agnostic to the physical interactions between the peptide and the membrane. Towards developing higher throughput screening methodologies, here we use Markov State Modeling and all-atom molecular dynamics simulations to quantify the membrane binding and insertion kinetics of three prototypical and antimicrobial peptides (alpha-helical magainin 2 and PGLa and beta-hairpin tachyplesin 1). By leveraging a set of collective variables that capture the essential physics of the amphiphilic and cationic peptide-membrane interactions we reveal how the slowest kinetic process of membrane insertion is the dynamic rolling of the peptide from a prebound to fully inserted state. These results add critical details to how antimicrobial peptides insert into bacterial membranes.


title: Interpolation of rational scrolls

abstract: We show in many cases that there exist rational scrolls which are balanced, i.e. they contain the expected number of general linear spaces as rulings. For example, there exist balanced scrolls of degree $mk+1$ and fibre dimension $k$ in $\P^{2k+1}$ for all $m\geq 1$.


title: Fundamental limitations on the device-independent quantum conference key agreement

abstract: We provide several general upper bounds on the rate of a key secure against a quantum adversary in the device-independent conference key agreement (DI-CKA) scenario. They include bounds by reduced entanglement measures and those based on multipartite secrecy monotones such as a multipartite squashed entanglement-based measure, which we refer to as reduced c-squashed entanglement. We compare the latter bound with the known lower bound for the protocol of conference key distillation based on the parity Clauser-Horne-Shimony-Holt game. We also show that the gap between the DI-CKA rate and the device-dependent rate is inherited from the bipartite gap between device-independent and device-dependent key rates, giving examples that exhibit the strict gap.


title: Robust Auction Design in the Auto-bidding World

abstract: In classic auction theory, reserve prices are known to be effective for improving revenue for the auctioneer against quasi-linear utility maximizing bidders. The introduction of reserve prices, however, usually do not help improve total welfare of the auctioneer and the bidders. In this paper, we focus on value maximizing bidders with return on spend constraints -- a paradigm that has drawn considerable attention recently as more advertisers adopt auto-bidding algorithms in advertising platforms -- and show that the introduction of reserve prices has a novel impact on the market. Namely, by choosing reserve prices appropriately the auctioneer can improve not only the total revenue but also the total welfare. Our results also demonstrate that reserve prices are robust to bidder types, i.e., reserve prices work well for different bidder types, such as value maximizers and utility maximizers, without using bidder type information. We generalize these results for a variety of auction mechanisms such as VCG, GSP, and first-price auctions. Moreover, we show how to combine these results with additive boosts to improve the welfare of the outcomes of the auction further. Finally, we complement our theoretical observations with an empirical study confirming the effectiveness of these ideas using data from online advertising auctions.


title: Production of 177Lu with deuterons at IFMIF-DONES facility

abstract: The International Fusion Materials Irradiation Facility - Demo Oriented NEutron Source (IFMIF-DONES) is a single-sited novel Research Infrastructure for testing, validation and qualification of the materials to be used in a fusion reactor. The main purpose of IFMIF-DONES is related to fusion technology and neutron irradiation of the materials to be used in the future fusion power plants. However, there is an important effort to take advantage of the outstanding characteristics of the facility in terms of neutrons and deuterons. One of the applications could be radioisotopes production with deuterons and neutrons. We discuss here the possible production of radioisotopes with deuterons at DONES. In this work, we have focused on the production of 177Lu with deuterons. The study has been carried out through the design and simulation of a device as cooling system for the sample producing the radioisotope. The results show the viability of using DONES for such production. In addition, the study suggests that new nuclear energy data above 20 MeV for deuterons is mandatory for an accurate study in this field.


title: A priori estimates for finite-energy sign-changing blowing-up solutions of critical elliptic equations

abstract: We prove sharp pointwise blow-up estimates for finite-energy sign-changing solutions of critical equations of Schr\"odinger-Yamabe type on a closed Riemannian manifold $(M,g)$ of dimension $n \ge 3$. This is a generalisation of the so-called $C^0$-theory for positive solutions of Schr\"odinger-Yamabe type equations. To deal with the sign-changing case we develop a method of proof that combines an \emph{a priori} bubble-tree analysis with a finite-dimensional reduction, and reduces the proof to obtaining sharp \emph{a priori} blow-up estimates for a linear problem.


title: Generalized Integer Splines on Arbitrary Graphs

abstract: Generalized integer splines on a graph $G$ with integer edge weights are integer vertex labelings such that if two vertices share an edge in $G$, the vertex labels are congruent modulo the edge weight. We introduce collapsing operations that reduce any simple graph to a single vertex, carrying with it the edge weight information. This corresponds to a sequence of surjective maps between the associated spline modules, leading to an explicit construction of a module basis in terms of the edge weights.


title: Cross-polarized surface lattice resonances in a rectangular lattice plasmonic metasurface

abstract: Multiresonant metasurfaces could enable many applications in filtering, sensing and nonlinear optics. However, developing a metasurface with more than one high-quality-factor or high-Q resonance at designated resonant wavelengths is challenging. Here, we experimentally demonstrate a plasmonic metasurface exhibiting different, narrow surface lattice resonances by exploiting the polarization degree of freedom where different lattice modes propagate along different dimensions of the lattice. The surface consists of aluminum nanostructures in a rectangular periodic lattice. The resulting surface lattice resonances were measured around 630 nm and 1160 nm with Q-factors of ~50 and ~800, respectively. The latter is a record-high plasmonic Q-factor within the near-infrared type-II window. Such metasurfaces could benefit applications such as frequency conversion and all-optical switching.


title: Yu-Shiba-Rusinov states in 2D superconductors with arbitrary Fermi contours

abstract: Magnetic impurities on a superconductor induce sub-gap Yu-Shiba-Rusinov (YSR) bound states, localized at the impurity site and fading away from it for distances up to several nanometers. In this article, we present a theoretical method to calculate the spatial distribution of the YSR spectrum of a two-dimensional superconductor with arbitrary Fermi contours (FCs) in the presence of magnetic impurities. Based on the Green's Function (GF) formalism, we obtain a general analytical expression by approximating an arbitrary contour shape to a regular polygon. This method allows us to show the connection between the spatial decay (and, hence, the extension) of YSR states and the shape of the FC of the host superconductor. We further apply this formalism to compute the evolution of YSR states in the presence of a nearby impurity atom, and compare the results with Scanning Tunneling Microscopy (STM) measurements on interacting manganese dimers on the $\beta$-Bi2Pd superconductor. The method can be easily extended to any arbitrary number of magnetically coupled impurities, thus, providing a useful tool for simulating the spectral properties of interacting YSR states in artificial atomic nanostructures.


title: A model of persistent breaking of continuous symmetry

abstract: We consider a UV-complete field-theoretic model in general dimensions, including $d=2+1$, that exhibits spontaneous breaking of continuous symmetry, persisting to arbitrarily large temperatures. Our model consists of two copies of the long-range vector models, with $O(m)$ and $O(N-m)$ global symmetry groups, perturbed by double-trace operators. Using conformal perturbation theory we find weakly-coupled IR fixed points for $N\geq 6$ that reveal a spontaneous breaking of global symmetry. Namely, at finite temperature the lower rank group is broken, with the pattern persisting at all temperatures due to scale-invariance. We provide evidence that the models in question are unitary and invariant under full conformal symmetry. Our work generalizes recent results, which considered the particular case of $m=1$ and reported persistent breaking of the discrete $\mathbb{Z}_2=O(1)$. Furthermore, we show that this model exhibits a continuous family of weakly interacting field theories at finite $N$.


title: Does the $J^{PC}=1^{+-}$ counterpart of the $X(3872)$ exist?

abstract: We explore the possible existence of the $J^{PC}=1^{+-}$ counterpart of the $X(3872)$ state in a coupled-channels calculation within a constituent quark model, with the aim of confirming the existence of the so-called $\tilde{X}(3872)$ state observed by the COMPASS Collaboration. Two states are found in the energy region of the $\tilde X(3872)$ signal, both with almost equal mixture of $c\bar c$ $2^1P_1$ state and $D^*\bar D^{(*)}$ channels: One that can be identified as the dressed $c\bar c$ $2^1P_1$ and a bound state below the $D\bar D^*$ threshold. We provide predictions of strong and radiative decays that could help to clarify the existence of such structures.


title: Resistive switching in HfO2-x/La0.67Sr0.33MnO3 heterostructure: An intriguing case of low H-field susceptibility of an E-field controlled active interface

abstract: High-performance non-volatile resistive random access memories (ReRAM) and their small stimuli control are of immense interest for high-speed computation and big-data processing in the emergent Internet of Things (IOT) arena. Here, we examine the resistive switching (RS) behavior in growth controlled HfO2/La0.67Sr0.33MnO3 heterostructures and their tunability under low magnetic field. It is demonstrated that oxygen-deficient HfO2 films show bipolar switching with high on/off ratio, stable retention, as well as good endurance owing to the orthorhombic-rich phase constitution and charge (de)-trapping-enabled Schottky-type conduction. Most importantly, we have demonstrated that, the RS can be tuned by a very low externally applied magnetic field (0-30 mT). Remarkably, application of a magnetic field of 30 mT causes the RS to be fully quenched and frozen in the high resistance state (HRS) even after the removal of magnetic field. However, the quenched state could be resurrected by applying higher bias voltage than the one for initial switching. This is argued to be a consequence of the electronically and ionically active nature of the HfO2-x/LSMO interface on both sides, and its susceptibility to the electric and low magnetic field effects. This result could pave the way for new designs of interface engineered high-performance oxitronic ReRAM devices.


title: Measurement of scintillation response of CsI[Na] to low-energy nuclear recoils by COHERENT

abstract: We present results of several measurements of CsI[Na] scintillation response to 3-60 keV energy nuclear recoils performed by the COHERENT collaboration using tagged neutron elastic scattering experiments and an endpoint technique. Earlier results, used to estimate the coherent elastic neutrino-nucleus scattering (CEvNS) event rate for the first observation of this process achieved by COHERENT at the Spallation Neutron Source (SNS), have been reassessed. We discuss corrections for the identified systematic effects and update the respective uncertainty values. The impact of updated results on future precision tests of CEvNS is estimated. We scrutinize potential systematic effects that could affect each measurement. In particular we confirm the response of the H11934-200 Hamamatsu photomultiplier tube (PMT) used for the measurements presented in this study to be linear in the relevant signal scale region.


title: HOLZ: High-Order Entropy Encoding of Lempel-Ziv Factor Distances

abstract: We propose a new representation of the offsets of the Lempel-Ziv (LZ) factorization based on the co-lexicographic order of the processed prefixes. The selected offsets tend to approach the k-th order empirical entropy. Our evaluations show that this choice of offsets is superior to the rightmost LZ parsing and the bit-optimal LZ parsing on datasets with small high-order entropy.


title: Intrinsic quantum correlations for Gaussian localized Dirac cat states in phase space

abstract: Following the information-based approach to Dirac spinors under a constant magnetic field, the phase-space representation of symmetric and anti-symmetric localized Dirac cat states is obtained. The intrinsic entanglement profile implied by the Dirac Hamiltonian is then investigated so as to shed a light on quantum states as carriers of qubits correlated by phase-space variables. Corresponding to the superposition of Gaussian states, cat states exhibit non-trivial elementary information dynamics which include the interplay between intrinsic entanglement and quantum superposition as reported by the corresponding Dirac archetypes. Despite the involved time-evolution as non-stationary states, the Wigner function constrains the elementary information quantifiers according to a robust framework which can be consistently used for quantifying the time-dependent $SU(2) \otimes SU(2)$ (spin projection and intrinsic parity) correlation profile of phase-space localized Dirac spinor states. Our results show that the Dirac Wigner functions for cat states -- described in terms of generalized Laguerre polynomials -- exhibit an almost maximized timely persistent mutual information profile which is engendered by either classical- or quantum-like spin-parity correlations, depending on the magnetic field intensity.


title: Linear-time Minimization of Wheeler DFAs

abstract: Wheeler DFAs (WDFAs) are a sub-class of finite-state automata which is playing an important role in the emerging field of compressed data structures: as opposed to general automata, WDFAs can be stored in just $\log\sigma + O(1)$ bits per edge, $\sigma$ being the alphabet's size, and support optimal-time pattern matching queries on the substring closure of the language they recognize. An important step to achieve further compression is minimization. When the input $\mathcal A$ is a general deterministic finite-state automaton (DFA), the state-of-the-art is represented by the classic Hopcroft's algorithm, which runs in $O(|\mathcal A|\log |\mathcal A|)$ time. This algorithm stands at the core of the only existing minimization algorithm for Wheeler DFAs, which inherits its complexity. In this work, we show that the minimum WDFA equivalent to a given input WDFA can be computed in linear $O(|\mathcal A|)$ time. When run on de Bruijn WDFAs built from real DNA datasets, an implementation of our algorithm reduces the number of nodes from 14% to 51% at a speed of more than 1 million nodes per second.


title: Secure Namespaced Kernel Audit for Containers

abstract: Despite the wide usage of container-based cloud computing, container auditing for security analysis relies mostly on built-in host audit systems, which often lack the ability to capture high-fidelity container logs. State-of-the-art reference-monitor-based audit techniques greatly improve the quality of audit logs, but their system-wide architecture is too costly to be adapted for individual containers. Moreover, these techniques typically require extensive kernel modifications, making them difficult to deploy in practical settings. In this paper, we present saBPF (secure audit BPF), an extension of the eBPF framework capable of deploying secure system-level audit mechanisms at the container granularity. We demonstrate the practicality of saBPF in Kubernetes by designing an audit framework, an intrusion detection system, and a lightweight access control mechanism. We evaluate saBPF and show that it is comparable in performance and security guarantees to audit systems from the literature that are implemented directly in the kernel.


title: Introducing MAMBO: Materials And Molecules Basic Ontology

abstract: Recent advances in computational and experimental technologies applied to the design and development of novel materials have brought out the need for systematic, rational and efficient methods for the organization of knowledge in the field. In this work, we present the initial steps carried out in the development of MAMBO - an ontology focused on the organization of concepts and knowledge in the field of materials based on molecules and targeted to applications. Our approach is guided by the needs of the communities involved in the development of novel molecular materials with functional properties at the nanoscale. As such, MAMBO aims at bridging the gaps of ongoing efforts in the development of ontologies in the materials science domain. By extending current work in the field, the modular nature of MAMBO also allows straightforward extension of concepts and relations to neighboring domains. Our work is expected to enable the systematic integration of computational and experimental data in specific domains of interest (nanomaterials, molecular materials, organic an polymeric materials, supramolecular and bio-organic systems, etc.). Moreover, MAMBO can be applied to the development of data-driven integrated predictive frameworks for the design of novel materials with tailored functional properties.


title: On the clique behavior of graphs of low degree

abstract: To any simple graph $G$, the clique graph operator $K$ associates the graph $K(G)$ which is the intersection graph of the maximal complete subgraphs of $G$. The iterated clique graphs are defined by $K^{0}(G)=G$ and $K^{n}(G)=K(K^{n-1}(G))$ for $n\geq 1$. If there are $m<n$ such that $K^{m}(G)$ is isomorphic to $K^{n}(G)$ we say that $G$ is convergent, otherwise, $G$ is divergent. The first example of a divergent graph was shown by Neumann-Lara in the 1970s, and is the graph of the octahedron. In this paper, we prove that among the connected graphs with maximum degree 4, the octahedron is the only one that is divergent.


title: Accelerated replica exchange stochastic gradient Langevin diffusion enhanced Bayesian DeepONet for solving noisy parametric PDEs

abstract: The Deep Operator Networks~(DeepONet) is a fundamentally different class of neural networks that we train to approximate nonlinear operators, including the solution operator of parametric partial differential equations (PDE). DeepONets have shown remarkable approximation and generalization capabilities even when trained with relatively small datasets. However, the performance of DeepONets deteriorates when the training data is polluted with noise, a scenario that occurs very often in practice. To enable DeepONets training with noisy data, we propose using the Bayesian framework of replica-exchange Langevin diffusion. Such a framework uses two particles, one for exploring and another for exploiting the loss function landscape of DeepONets. We show that the proposed framework's exploration and exploitation capabilities enable (1) improved training convergence for DeepONets in noisy scenarios and (2) attaching an uncertainty estimate for the predicted solutions of parametric PDEs. In addition, we show that replica-exchange Langeving Diffusion (remarkably) also improves the DeepONet's mean prediction accuracy in noisy scenarios compared with vanilla DeepONets trained with state-of-the-art gradient-based optimization algorithms (e.g. Adam). To reduce the potentially high computational cost of replica, in this work, we propose an accelerated training framework for replica-exchange Langevin diffusion that exploits the neural network architecture of DeepONets to reduce its computational cost up to 25% without compromising the proposed framework's performance. Finally, we illustrate the effectiveness of the proposed Bayesian framework using a series of experiments on four parametric PDE problems.


title: Three-dimensional Skyrme Hartree-Fock-Bogoliubov solver in coordinate-space representation

abstract: The coordinate-space representation of the Hartree-Fock-Bogoliubov theory is the method of choice to study weakly bound nuclei whose properties are affected by the quasiparticle continuum space. To describe such systems, we developed a three-dimensional Skyrme-Hartree-Fock-Bogoliubov solver HFBFFT based on the existing, highly optimized and parallelized Skyrme-Hartree-Fock code Sky3D. The code does not impose any self-consistent spatial symmetries such as mirror inversions or parity. The underlying equations are solved in HFBFFT directly in the canonical basis using the fast Fourier transform. To remedy the problems with pairing collapse, we implemented the soft energy cutoff and pairing annealing. The convergence of HFB solutions was improved by a sub-iteration method. The Hermiticity violation of differential operators brought by Fourier-transform-based differentiation has also been solved. The accuracy and performance of HFBFFT were tested by benchmarking it against other HFB codes, both spherical and deformed, for a set of nuclei, both well-bound and weakly-bound.


title: Convex Chance-Constrained Programs with Wasserstein Ambiguity

abstract: Chance constraints yield non-convex feasible regions in general. In particular, when the uncertain parameters are modeled by a Wasserstein ball, arXiv:1806.07418 and arXiv:1809.00210 showed that the distributionally robust (pessimistic) chance constraint admits a mixed-integer conic representation. This paper identifies sufficient conditions that lead to convex feasible regions of chance constraints with Wasserstein ambiguity. First, when uncertainty arises from the left-hand side of a pessimistic individual chance constraint, we derive a convex and conic representation if the Wasserstein ball is centered around a Gaussian distribution. Second, when uncertainty arises from the right-hand side of a pessimistic joint chance constraint, we show that the ensuing feasible region is convex if the Wasserstein ball is centered around a log-concave distribution (or, more generally, an $\alpha$-concave distribution with $\alpha \geq -1$). In addition, we propose a block coordinate ascent algorithm for this class of chance constraints and prove its convergence to global optimum. Furthermore, we extend the convexity results and conic representation to optimistic chance constraints.