title: Molecules with ALMA at Planet-forming Scales (MAPS) III: Characteristics of Radial Chemical Substructures

abstract: The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a detailed, high resolution (${\sim}$10-20 au) view of molecular line emission in five protoplanetary disks at spatial scales relevant for planet formation. Here, we present a systematic analysis of chemical substructures in 18 molecular lines toward the MAPS sources: IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. We identify more than 200 chemical substructures, which are found at nearly all radii where line emission is detected. A wide diversity of radial morphologies - including rings, gaps, and plateaus - is observed both within each disk and across the MAPS sample. This diversity in line emission profiles is also present in the innermost 50 au. Overall, this suggests that planets form in varied chemical environments both across disks and at different radii within the same disk. Interior to 150 au, the majority of chemical substructures across the MAPS disks are spatially coincident with substructures in the millimeter continuum, indicative of physical and chemical links between the disk midplane and warm, elevated molecular emission layers. Some chemical substructures in the inner disk and most chemical substructures exterior to 150 au cannot be directly linked to dust substructure, however, which indicates that there are also other causes of chemical substructures, such as snowlines, gradients in UV photon fluxes, ionization, and radially-varying elemental ratios. This implies that chemical substructures could be developed into powerful probes of different disk characteristics, in addition to influencing the environments within which planets assemble. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.


title: Deep Hubble Space Telescope Observations of GW170817: Complete Light Curves and the Properties of the Galaxy Merger of NGC 4993

abstract: We present the complete set of {\it Hubble Space Telescope} imaging of the binary neutron star merger GW170817 and its optical counterpart AT 2017gfo. Including deep template imaging in F814W, F110W, F140W, and F160W at 3.4 years post-merger, we re-analyze the full light curve of AT 2017gfo across 12 bands from 5--1273 rest-frame days after merger. We obtain four new detections of the short $\gamma$-ray burst (GRB) 170817A afterglow from 109--170 rest-frame days post-merger. These detections are consistent with the previously observed $\beta=-0.6$ spectral index in the afterglow light curve with no evidence for spectral evolution. We also analyze our limits in the case of novel late-time optical and IR emission signatures, such as a kilonova afterglow or infrared dust echo, but find our limits are not constraining in these contexts. We use the new data to construct deep optical and infrared stacks, reaching limits of $M=-6.3$ to $-4.6$ mag, to analyze the local environment around AT 2017gfo and low surface brightness features in its host galaxy NGC 4993. We rule out the presence of any globular cluster at the position of AT 2017gfo to $2.3 \times 10^{4} L_{\odot}$, including those with the reddest $V-H$ colors. Finally, we analyze the substructure of NGC 4993 in deep residual imaging, and find shell features which extend up to 71.8\arcsec\ (14.2 kpc) from the center of the galaxy. We find that the shells have a cumulative stellar mass of $6.3\times10^{8} M_{\odot}$, roughly 2% the total stellar mass of NGC 4993, and mass-weighted ages of $>$3 Gyr. We conclude that it was unlikely the GW170817 progenitor system formed in the galaxy merger, which based on dynamical signatures and the stellar population in the shells mostly likely occurred 220--685 Myr ago.


title: Fine corrections in the effective string describing SU(2) Yang-Mills theory in three dimensions

abstract: We present a study of the effective string that describes the infrared dynamics of SU(2) Yang-Mills theory in three dimensions. By combining high-precision lattice simulation results for Polyakov-loop correlators at finite temperatures close to (and less than) the deconfinement one with the analytical constraints from renormalization-group arguments, from the exact integrability of the two-dimensional Ising model that describes the universality class of the critical point of the theory, from conformal perturbation theory, and from Lorentz invariance, we derive tight quantitative bounds on the corrections to the effective string action beyond the Nambu-Goto approximation. We show that these corrections are compatible with the predictions derived from a bootstrap analysis of the effective string theory, and have a value which does not allow to prove the Axionic String Ansatz for this model


title: Measuring sloshing, merging and feedback velocities in the Virgo cluster

abstract: We present a detailed analysis of the velocity structure of the Virgo cluster using {\it XMM-Newton} observations. Using a novel technique which uses uses the Cu K$\alpha$ instrumental line to calibrate the EPIC-pn energy scale, we are able to obtain velocity measurements with uncertainties down to $\Delta v \sim 100$ km/s. We created 2D projected maps for the velocity, temperature, metallicity, density, pressure and entropy with an spatial resolution of 0.25$'$. We have found that in the innermost gas there is a high velocity structure, most likely indicating the presence of an outflow from the AGN while our analysis of the cluster cool core using RGS data indicates that the velocity of the gas agrees with the M87 optical redshift. An overall gradient in the velocity is seen, with larger values as we move away from the cluster core. The hot gas located within the western radio flow is redshifted, moving with a velocity $\sim 331$ km/s while the hot gas located within the eastern radio flow is blueshifted, with a velocity $\sim 258$ km/s, suggesting the presence of backflows. Our results reveal the effects of both AGN outflows and gas sloshing, in the complex velocity field of the Virgo cluster.


title: Semiclassical equations of motion for disordered conductors: extrinsic interband velocity, corrected collision integral and spin-orbit torques

abstract: The semiclassical equations of motion are widely used to describe carrier transport in conducting materials. Nevertheless, the substantial challenge of incorporating disorder systematically into the semiclassical model persists, leading to quantitative inaccuracies and occasionally erroneous predictions for the expectation values of physical observables. In the present work we provide a general prescription for reformulating the semiclassical equations of motion for carriers in disordered conductors by taking the quantum mechanical density matrix as the starting point. We focus on external electric fields, without magnetic fields, and spin-independent disorder. The density matrix approach allows averaging over impurity configurations, and the trace of the velocity operator with the disorder-averaged density matrix can be reinterpreted as the semiclassical velocity weighted by the Boltzmann distribution function. Through this rationale the well-known intrinsic group and anomalous velocities are trivially recovered, while we demonstrate the existence of an extrinsic interband velocity, namely a disorder correction to the semiclassical velocity of Bloch electrons, mediated by the interband matrix elements of the Berry connection. A similar correction is present in the non-equilibrium expectation value of the spin operator, contributing to spin-orbit torques. To obtain agreement with diagrammatic approaches the scattering term in the Boltzmann equation is corrected to first order in the electric field, and the Boltzmann equation is solved up to sub-leading order in the disorder potential. Our prescription ensures all vertex corrections present in diagrammatic treatments are taken into account, and to illustrate this we discuss model cases in topological insulators, including the anomalous Hall effect as well as spin-orbit torques.


title: Quenching of satellite galaxies of Milky Way analogues: reconciling theory and observations

abstract: The vast majority of low-mass satellite galaxies around the Milky Way and M31 appear virtually devoid of cool gas and show no signs of recent or ongoing star formation. Cosmological simulations demonstrate that such quenching is expected and is due to the harsh environmental conditions that satellites face when joining the Local Group (LG). However, recent observations of Milky Way analogues in the SAGA survey present a very different picture, showing the majority of observed satellites to be actively forming stars, calling into question the realism of current simulations and the typicality of the LG. Here we use the ARTEMIS suite of high-resolution cosmological hydrodynamical simulations to carry out a careful comparison with observations of dwarf satellites in the LG, SAGA, and the Local Volume (LV) survey. We show that differences between SAGA and the LG and LV surveys, as well as between SAGA and the ARTEMIS simulations, can be strongly reduced by considering differences in the host mass distributions and (more importantly) observational selection effects, specifically that low-mass satellites which have only recently been accreted are more likely to be star-forming, have a higher optical surface brightness, and are therefore more likely to be included in the SAGA survey. This picture is confirmed using data from the deeper LV survey, which shows pronounced quenching at low masses, in accordance with the predictions of LCDM-based simulations.


title: Sulfur abundances in the Galactic bulge and disk

abstract: Context. The measurement of $\alpha$-elements abundances provides a powerful tool to put constraints on chemical evolution and star formation history of galaxies. The majority of studies on the $\alpha$-element Sulfur (S) are focused on local stars, making S behavior in other environments an astronomical topic yet to be analyzed. Aims. The investigation of S in the Galactic bulge has only recently been considered for the first time. This work aims to improve our knowledge on S behavior in this component of the Milky Way. Methods. We present S abundances of 74 dwarf and sub-giant stars in the Galactic bulge, 21 and 30 F and G thick and thin disk stars. We performed local thermodynamic equilibrium analysis and applied corrections for non-LTE on high resolution and high signal-to-noise UVES spectra. S abundances were derived from multiplets 1, 6 and 8 in the metallicity range $-2<$[Fe/H]$<$0.6, by spectrosynthesis or line equivalent widths. Results. We confirm that S behaves like an $\alpha$-element within the Galactic bulge. In the [S/Fe] versus [Fe/H] diagram, S presents a plateau at low metallicity followed by a decreasing of [S/Fe] with the increasing of [Fe/H], until reaching [S/Fe]$\sim0$ at super-solar metallicity. We found that the Galactic bulge is S-rich with respect to both the thick and thin disks at $-1<$[Fe/H]$<0.3$, supporting a more rapid formation and chemical evolution of the Galactic bulge than the disk.


title: Molecules with ALMA at Planet-forming Scales (MAPS) IV: Emission Surfaces and Vertical Distribution of Molecules

abstract: The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a unique opportunity to study the vertical distribution of gas, chemistry, and temperature in the protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. By using the asymmetry of molecular line emission relative to the disk major axis, we infer the emission height ($z$) above the midplane as a function of radius ($r$). Using this method, we measure emitting surfaces for a suite of CO isotopologues, HCN, and C$_2$H. We find that $^{12}$CO emission traces the most elevated regions with $z/r > 0.3$, while emission from the less abundant $^{13}$CO and C$^{18}$O probes deeper into the disk at altitudes of $z/r \lesssim 0.2$. C$_2$H and HCN have lower opacities and SNRs, making surface fitting more difficult, and could only be reliably constrained in AS 209, HD 163296, and MWC 480, with $z/r \lesssim 0.1$, i.e., relatively close to the planet-forming midplanes. We determine peak brightness temperatures of the optically thick CO isotopologues and use these to trace 2D disk temperature structures. Several CO temperature profiles and emission surfaces show dips in temperature or vertical height, some of which are associated with gaps and rings in line and/or continuum emission. These substructures may be due to local changes in CO column density, gas surface density, or gas temperatures, and detailed thermo-chemical models are necessary to better constrain their origins and relate the chemical compositions of elevated disk layers with those of planet-forming material in disk midplanes. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.


title: Molecules with ALMA at Planet-forming Scales (MAPS XVIII): Kinematic Substructures in the Disks of HD 163296 and MWC 480

abstract: We explore the dynamical structure of the protoplanetary disks surrounding HD 163296 and MWC 480 as part of the Molecules with ALMA at Planet Forming Scales (MAPS) large program. Using the $J = 2-1$ transitions of $^{12}$CO, $^{13}$CO and C$^{18}$O imaged at spatial resolutions of $\sim 0.^{\prime \prime}15$ and with a channel spacing of $200$ ${\rm m\,s^{-1}}$, we find perturbations from Keplerian rotation in the projected velocity fields of both disks ($\lesssim\!5\%$ of the local Keplerian velocity), suggestive of large-scale (10s of au in size), coherent flows. By accounting for the azimuthal dependence on the projection of the velocity field, the velocity fields were decomposed into azimuthally averaged orthogonal components, $v_{\phi}$, $v_r$ and $v_z$. Using the optically thick $^{12}$CO emission as a probe of the gas temperature, local variations of $\approx\! 3$ K ($\approx\! 5 \%$ relative changes) were observed and found to be associated with the kinematic substructures. The MWC 480 disk hosts a suite of tightly wound spiral arms. The spirals arms, in conjunction with the highly localized perturbations in the gas velocity structure (kinematic planetary signatures), indicate a giant planet, $\sim\! 1$ $M_{\rm Jup}$, at a radius of $\approx 245$ au. In the disk of HD 163296, the kinematic substructures were consistent with previous studies of Pinte et al. (2018a) and Teague et al. (2018a) advocating for multiple $\sim\! 1$ $M_{\rm Jup}$ planets embedded in the disk. These results demonstrate that molecular line observations that characterize the dynamical structure of disks can be used to search for the signatures of embedded planets. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.


title: The ZTF-BTS Type Ia supernovae luminosity function is consistent with a single progenitor channel for the explosions

abstract: We construct the Type Ia supernovae (SNe Ia) luminosity function (LF) using the Zwicky Transient Facility Bright Transient Survey (BTS) catalogue. While this magnitude-limited survey has an unprecedented number of objects, it suffers from large distance uncertainties and lacks an estimation of host extinction. We bypass these issues by calculating the intrinsic luminosities from the shape parameters of the light curve's $ g $ and $ r $ bands, with the luminosities calibrated from the well observed SNe Ia sample of the Carnegie Supernova Project, allowing us to construct, for the first time, the intrinsic LF of SNe Ia. We then use a novel tight relation between the color stretch and the synthesized $^{56}$Ni mass, $M_\mathrm{Ni56}$, to determine the $M_\mathrm{Ni56}$ distribution of SNe Ia. We find that the LFs are unimodal, with their peaks in line with previous results, but have a much lower rate of dim events and luminous events. We show that the features on top of the unimodal LF-derived distributions are all compatible with statistical noise, consistent with a single progenitor channel for the explosions. We further derive, for the first time, the SNe Ia distribution of host galaxy extinction, and find a mean selective extinction of $E(B-V)\approx0.1$ and a non-negligible fraction with large, $ >1\,\text{mag} $, extinction in the optical bands. The high extinction is typical for luminous SNe, supporting their young population origin.


title: Scar states in a system of interacting chiral fermions

abstract: We study the nature of many-body eigenstates of a system of interacting chiral spinless fermions on a ring. We find a coexistence of fermionic and bosonic types of eigenstates in parts of the many-body spectrum. Some bosonic eigenstates, native to the strong interaction limit, persist at intermediate and weak couplings, enabling persistent density oscillations in the system, despite it being far from integrability.


title: Molecules with ALMA at Planet-forming Scales (MAPS). VII. Sub-stellar O/H and C/H and super-stellar C/O in planet feeding gas

abstract: The elemental composition of the gas and dust in a protoplanetary disk influences the compositions of the planets that form in it. We use the Molecules with ALMA at Planet-forming Scales (MAPS) data to constrain the elemental composition of the gas at the locations of potentially forming planets. The elemental abundances are inferred by comparing source-specific gas-grain thermochemical models, with variable C/O ratios and small-grain abundances, from the DALI code with CO and C2H column densities derived from the high-resolution observations of the disks of AS 209, HD 163296, and MWC 480. Elevated C/O ratios (~2.0), even within the CO ice line, are necessary to match the inferred C2H column densities, over most of the pebble disk. Combined with constraints on the CO abundances in these systems, this implies that both the O/H and C/H ratios in the gas are substellar by a factor of 4-10, with the O/H depleted by a factor of 20-50, resulting in the high C/O ratios. This necessitates that even within the CO ice line, most of the volatile carbon and oxygen is still trapped on grains in the midplane. Planets accreting gas in the gaps of the AS 209, HD 163296, and MWC 480 disks will thus acquire very little carbon and oxygen after reaching the pebble isolation mass. In the absence of atmosphere-enriching events, these planets would thus have a strongly substellar O/H and C/H and superstellar C/O atmospheric composition.


title: The Cosmic Evolution of Binary Black Holes in Young, Globular and Nuclear Star Clusters: Rates, Masses, Spins and Mixing Fractions

abstract: The growing population of binary black holes (BBHs) observed by gravitational wave detectors is a potential Rosetta stone for understanding their formation channels. Here, we use an upgraded version of our semi-analytic codes {\sc fastcluster} and {\sc cosmo$\mathcal{R}$ate} to investigate the cosmic evolution of four different BBH populations: isolated BBHs and dynamically formed BBHs in nuclear star clusters (NSCs), globular clusters (GCs), and young star clusters (YSCs). With our approach, we can study different channels assuming the same stellar and binary input physics. We find that the merger rate density of BBHs in GCs and NSCs is barely affected by stellar metallicity ($Z$), while the rate of isolated BBHs changes wildly with $Z$. BBHs in YSCs behave in an intermediate way between isolated and GC/NSC BBHs. The local merger rate density of Nth-generation black holes (BHs), obtained by summing up hierarchical mergers in GCs, NSCs and YSCs, ranges from $\sim{1}$ to $\sim{4}$ Gpc$^{-3}$ yr$^{-1}$ and is mostly sensitive to the spin parameter. We find that the mass function of primary BHs evolves with redshift in GCs and NSCs, becoming more top-heavy at higher $z$. In contrast, the primary BH mass function almost does not change with redshift in YSCs and in the field. This signature of the BH mass function has relevant implications for Einstein Telescope and Cosmic Explorer. Finally, our analysis suggests that multiple channels contribute to the BBH population of the second gravitational-wave transient catalog.

categories: astro-ph.HE astro-ph.GA astro-ph.SR gr-qc


title: Molecules with ALMA at Planet-forming Scales (MAPS). XV. Tracing protoplanetary disk structure within 20 au

abstract: Constraining the distribution of gas and dust in the inner 20 au of protoplanetary disks is difficult. At the same time, this region is thought to be responsible for most planet formation, especially around the water ice line at 3-10 au. Under the assumption that the gas is in a Keplerian disk, we use the exquisite sensitivity of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA large program to construct radial surface brightness profiles with a ~3 au effective resolution for the CO isotopologue J=2-1 lines using the line velocity profile. IM Lup reveals a central depression in 13CO and C18O that is ascribed to a pileup of ~500 $M_\oplus$ of dust in the inner 20 au, leading to a gas-to-dust ratio of around <10. This pileup is consistent with efficient drift of grains ($\gtrsim$ 100 $M_\oplus$ Myr$^{-1}$) and a local gas-to-dust ratio that suggests that the streaming instability could be active. The CO isotopologue emission in the GM Aur disk is consistent with a small (~15 au), strongly depleted gas cavity within the ~40 au dust cavity. The radial surface brightness profiles for both the AS 209 and HD 163296 disks show a local minimum and maximum in the C18O emission at the location of a known dust ring (~14 au) and gap (~10 au), respectively. This indicates that the dust ring has a low gas-to-dust ratio ($>$ 10) and that the dust gap is gas-rich enough to have optically thick C18O.


title: Molecules with ALMA at Planet-forming Scales (MAPS) XIX. Spiral Arms, a Tail, and Diffuse Structures Traced by CO around the GM Aur Disk

abstract: The concentric gaps and rings commonly observed in protoplanetary disks in millimeter continuum emission have lent the impression that planet formation generally proceeds within orderly, isolated systems. While deep observations of spatially resolved molecular emission have been comparatively limited, they are increasingly suggesting that some disks interact with their surroundings while planet formation is underway. We present an analysis of complex features identified around GM Aur in $^{12}$CO $J=2-1$ images at a spatial resolution of $\sim40$ au. In addition to a Keplerian disk extending to a radius of $\sim550$ au, the CO emission traces flocculent spiral arms out to radii of $\sim$1200 au, a tail extending $\sim1800$ au southwest of GM Aur, and diffuse structures extending from the north side of the disk up to radii of $\sim1900$ au. The diffuse structures coincide with a "dust ribbon" previously identified in scattered light. The large-scale asymmetric gas features present a striking contrast with the mostly axisymmetric, multi-ringed millimeter continuum tracing the pebble disk. We hypothesize that GM Aur's complex gas structures result from late infall of remnant envelope or cloud material onto the disk. The morphological similarities to the SU Aur and AB Aur systems, which are also located in the L1517 cloud, provide additional support to a scenario in which interactions with the environment are playing a role in regulating the distribution and transport of material in all three of these Class II disk systems. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.


title: Detection of the Crab Nebula by the prototype Schwarzschild-Couder Telescope

abstract: The Schwarzschild-Couder Telescope (SCT) is a medium-sized telescope technology proposed for the Cherenkov Telescope Array. It uses a novel dual-mirror optical design that removes comatic aberrations across its entire field of view. The SCT camera employs high-resolution silicon photomultiplier (SiPM) sensors with a pixel size of 4 arcminutes. A prototype SCT (pSCT) has been constructed at the Fred Lawrence Whipple Observatory in Arizona, USA. An observing campaign in 2020, with a partial camera of 1600 pixels (2.7 degrees by 2.7 degrees field of view) resulted in detection of the Crab Nebula at 8.6 sigma statistical significance. Work on the pSCT camera and optical system is ongoing to improve performance and prepare for an upcoming camera upgrade. The pSCT camera upgrade will replace the current camera modules with improved SiPMs and readout electronics and will expand the camera to its full design field of view of 8 degrees in diameter (11,328 pixels). The fully upgraded pSCT will enable next-generation very-high-energy gamma-ray astrophysics through excellent background rejection and angular resolution. In this presentation we describe first results from the successful operation of the pSCT and future plans.


title: Density-matrix renormalization group study of optical conductivity of the Mott insulator for two-dimensional clusters

abstract: The real part of optical conductivity $\text{Re}\sigma(\omega)$ of the Mott insulators has a large amount of information on how spin and charge degrees of freedom interact with each other. By using the time-dependent density-matrix renormalization group, we study $\text{Re}\sigma(\omega)$ of the two-dimensional Hubbard model on a square lattice at half filling. We find an excitonic peak at the Mott-gap edge of $\text{Re}\sigma(\omega)$ not only for the two-dimensional square lattice but also for two- and four-leg ladders. For the square lattice, however, we do not clearly find a gap between an excitonic peak and continuum band, which indicates that a bound state is not well defined. The emergence of an excitonic peak in $\text{Re}\sigma(\omega)$ implies the formation of a spin polaron. Examining the dependence of $\text{Re}\sigma(\omega)$ on the on-site Coulomb interaction and next-nearest neighbor hoppings, we confirm that an excitonic peak is generated from a magnetic effect. Electron scattering due to an electron-phonon interaction is expected to easily suppress an excitonic peak since spectral width of an excitonic peak is very narrow. Introducing a large broadening in $\text{Re}\sigma(\omega)$ by modeling the electron-phonon coupling present in La$_{2}$CuO$_{4}$ and Nd$_{2}$CuO$_{4}$, we obtain $\text{Re}\sigma(\omega)$ comparable with experiments.


title: Phase transitions of the Dicke model: a unified perspective

abstract: The Dicke model exhibits a variety of phase transitions. The quantum phase transition from the normal phase to the super-radiant phase is marked by a dramatic change in the scaling of the participation ratio. We find that the ground state in the super-radiant phase exhibits multifractality manifest in the participation ratio scaling as the square root of the full Hilbert space dimension. The thermal phase transition temperature, for which we obtain an exact analytical expression, is strikingly captured by the mutual information between two spins. In the excited state quantum phase transition within the super-radiant phase, we discover a new upper cut-off energy; the central energy band between the lower and upper cut-off energies shows distinctly different behaviour. This finding is corroborated with the aid of several eigenvector properties: von Neumann entanglement entropy between spins and bosons, the mean photon number, concurrence between two spins, and participation ratio. Thus we obtain a unified picture for the three different kinds of phase transitions.


title: Molecules with ALMA at Planet-forming Scales. XX. The Massive Disk Around GM Aurigae

abstract: Gas mass remains one of the most difficult protoplanetary disk properties to constrain. With much of the protoplanetary disk too cold for the main gas constituent, H2, to emit, alternative tracers such as dust, CO, or the H2 isotopolog HD are used. However, relying on disk mass measurements from any single tracer requires assumptions about the tracer's abundance relative to \hh\ and the disk temperature structure. Using new Atacama Large Millimeter/submillimeter Array (ALMA) observations from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program as well as archival ALMA observations, we construct a disk physical/chemical model of the protoplanetary disk GM Aur. Our model is in good agreement with the spatially resolved CO isotopolog emission from eleven rotational transitions with spatial resolution ranging from 0.15'' to 0.46'' (24-73 au at 159 pc) and the spatially unresolved HD J=1-0 detection from Herschel. Our best-fit model favors a cold protoplanetary disk with a total gas mass of approximately 0.2 solar masses, a factor of 10 reduction in CO gas inside roughly 100 au and a factor of 100 reduction outside of 100 au. Despite its large mass, the disk appears to be on the whole gravitationally stable based on the derived Toomre Q parameter. However, the region between 70 and 100 au, corresponding to one of the millimeter dust rings, is close to being unstable based on the calculated Toomre Q of <1.7. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.


title: Dark Energy at early times and ACT: a larger Hubble constant without late-time priors

abstract: In this paper we fit two models of Early Dark Energy (EDE) (an increase in the expansion rate before recombination) to the combination of Atacama Cosmology Telescope (ACT) measurements of the Cosmic Microwave Background (CMB) with data from either the WMAP or the Planck satellite, along with measurements of the baryon acoustic oscillations and uncalibrated supernovae luminosity distance. We study a phenomenological axion-like potential ('axEDE') and a scalar field experiencing a first-order phase-transition ('NEDE'). We find that for both models the 'Planck-free' analysis yields non-zero EDE at > 2 sigma and an increased value for $H_0 \sim 70-74$ km/s/Mpc, compatible with local measurements, without the inclusion of any prior on $H_0$. On the other hand, the inclusion of Planck data restricts the EDE contribution to an upper-limit only at 95% C.L. For axEDE, the combination of Planck and ACT leads to constraints 30% weaker than with Planck alone, and there is no residual Hubble tension. On the other hand, NEDE is more strongly constrained in a Planck+ACT analysis, and the Hubble tension remains at $\sim 3\sigma$, illustrating the ability for CMB data to distinguish between EDE models. We explore the apparent inconsistency between the Planck and ACT data and find that it comes (mostly) from a slight tension between the temperature power spectrum at multipoles around $\sim 1000$ and $\sim 1500$. Finally, through a mock analysis of ACT data, we demonstrate that the preference for EDE is not driven by a lack of information at high-$\ell$ when removing Planck data, and that a LCDM fit to the fiducial EDE cosmology results in a significant bias on $\{H_0,\omega_{\rm cdm}\}$. More accurate measurements of the TT power spectra above $\ell\sim 2500$ and EE between $\ell \sim 300-500$ will play a crucial role in differentiating EDE models.


title: Potential of the Sterile Insect Technique for Control of Deer Ticks, $\textit{Ixodes scapularis}$

abstract: The deer tick, $\textit{Ixodes scapularis}$, is a vector for numerous human diseases, including Lyme disease, anaplasmosis, and babesiosis. Concern is rising in the US and abroad as the population and range of this species grow and new diseases emerge. Herein I consider the potential for control of $\textit{I. scapularis}$ using the Sterile Insect Technique (SIT), which acts by reducing net fertility through release of sterile males. I construct a population model with density-dependent and -independent growth, migration, and an Allee effect (decline of the population when it is small), and use this model to simulate sterile tick release in both single- and multi-patch frameworks. I test two key concerns with implementing $\textit{I. scapularis}$ SIT: that the ticks' lengthy life course could make control take too long and that low migration might mean sterile males need thorough manual dispersal to all parts of the control area. Results suggest that typical $\textit{I. scapularis}$ SIT programs will take about eight years, a prediction near the normal range for the technique, but that thorough distribution of sterile ticks over the control area is indeed critical, increasing expense substantially by necessitating aerial release. With particularly high rearing costs also expected for $\textit{I. scapularis}$, the latter finding suggests that cost-effectiveness improvements to aerial release may be a prerequisite to $\textit{I. scapularis}$ SIT.


title: Categories of graphs for operadic structures

abstract: We recall several categories of graphs which are useful for describing homotopy-coherent versions of generalized operads (e.g. cyclic operads, modular operads, properads, and so on), and give new, uniform definitions for their morphisms. This allows for straightforward comparisons, and we use this to show that certain free-forgetful adjunctions between categories of generalized operads can be realized at the level of presheaves. This includes adjunctions between operads and cyclic operads, between dioperads and augmented cyclic operads, and betweeen wheeled properads and modular operads.


title: The Emergence of the Shape Bias Results from Communicative Efficiency

abstract: By the age of two, children tend to assume that new word categories are based on objects' shape, rather than their color or texture; this assumption is called the shape bias. They are thought to learn this bias by observing that their caregiver's language is biased towards shape based categories. This presents a chicken and egg problem: if the shape bias must be present in the language in order for children to learn it, how did it arise in language in the first place? In this paper, we propose that communicative efficiency explains both how the shape bias emerged and why it persists across generations. We model this process with neural emergent language agents that learn to communicate about raw pixelated images. First, we show that the shape bias emerges as a result of efficient communication strategies employed by agents. Second, we show that pressure brought on by communicative need is also necessary for it to persist across generations; simply having a shape bias in an agent's input language is insufficient. These results suggest that, over and above the operation of other learning strategies, the shape bias in human learners may emerge and be sustained by communicative pressures.


title: Molecules with ALMA at Planet-forming Scales (MAPS) V: CO gas distributions

abstract: Here we present high resolution (15-24 au) observations of CO isotopologue lines from the Molecules with ALMA on Planet-forming Scales (MAPS) ALMA Large Program. Our analysis employs $^{13}$CO and C$^{18}$O ($J$=2-1), (1-0), and C$^{17}$O (1-0) line observations of five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermo-chemical modeling framework based on the CO data, the broadband spectral energy distribution, and the mm-continuum emission; (2) an empirical temperature distribution based on optically thick CO lines; and (3) a direct fit to the C$^{17}$O hyperfine lines. Results from these methods generally show excellent agreement. The CO gas column density profiles of the five disks show significant variations in the absolute value and the radial shape. Assuming a gas-to-dust mass ratio of 100, all five disks have a global CO-to-H$_2$ abundance of 10-100 times lower than the ISM ratio. The CO gas distributions between 150-400 au match well with models of viscous disks, supporting the long-standing theory. CO gas gaps appear to be correlated with continuum gap locations, but some deep continuum gaps do not have corresponding CO gaps. The relative depths of CO and dust gaps are generally consistent with predictions of planet-disk interactions, but some CO gaps are 5-10 times shallower than predictions based on dust gaps. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.


title: Machine Learning for Online Algorithm Selection under Censored Feedback

abstract: In online algorithm selection (OAS), instances of an algorithmic problem class are presented to an agent one after another, and the agent has to quickly select a presumably best algorithm from a fixed set of candidate algorithms. For decision problems such as satisfiability (SAT), quality typically refers to the algorithm's runtime. As the latter is known to exhibit a heavy-tail distribution, an algorithm is normally stopped when exceeding a predefined upper time limit. As a consequence, machine learning methods used to optimize an algorithm selection strategy in a data-driven manner need to deal with right-censored samples, a problem that has received little attention in the literature so far. In this work, we revisit multi-armed bandit algorithms for OAS and discuss their capability of dealing with the problem. Moreover, we adapt them towards runtime-oriented losses, allowing for partially censored data while keeping a space- and time-complexity independent of the time horizon. In an extensive experimental evaluation on an adapted version of the ASlib benchmark, we demonstrate that theoretically well-founded methods based on Thompson sampling perform specifically strong and improve in comparison to existing methods.


title: Hierarchical Robust Adaptive Control for Wind Turbines with Actuator Fault

abstract: This paper solves the problem of regulating the rotor speed tracking error for wind turbines in the full-load region by an effective robust-adaptive control strategy. The developed controller compensates for the uncertainty in the control input effectiveness caused by a pitch actuator fault, unmeasurable wind disturbance, and nonlinearity in the model. Wind turbines have multi-layer structures such that the high-level structure is nonlinearly coupled through an aggregation of the low-level control authorities. Hence, the control design is divided into two stages. First, an $\mathcal{L}_2$ controller is designed to attenuate the influence of wind disturbance fluctuations on the rotor speed. Then, in the low-level layer, a controller is designed using a proposed adaptation mechanism to compensate for actuator faults. The theoretical results show that the closed-loop equilibrium point of the regulated rotor speed tracking error dynamics in the high level is finite-gain $\mathcal{L}_2$ stable, and the closed-loop error dynamics in the low level is globally asymptotically stable. Simulation results show that the developed controller significantly reduces the root mean square of the rotor speed error compared to some well-known works, despite the largely fluctuating wind disturbance, and the time-varying uncertainty in the control input effectiveness.


title: Chaos in the Bose-Hubbard model and random two-body Hamiltonians

abstract: We investigate the chaotic phase of the Bose-Hubbard model [L. Pausch et al, Phys. Rev. Lett. 126, 150601 (2021)] in relation to the bosonic embedded random matrix ensemble, which mirrors the dominant few-body nature of many-particle interactions, and hence the Fock space sparsity of quantum many-body systems. The energy dependence of the chaotic regime is well described by the bosonic embedded ensemble, which also reproduces the Bose-Hubbard chaotic eigenvector features, quantified by the expectation value and eigenstate-to-eigenstate fluctuations of fractal dimensions. Despite this agreement, in terms of the fractal dimension distribution, these two models depart from each other and from the Gaussian orthogonal ensemble as Hilbert space grows. These results provide further evidence of a way to discriminate among different many-body Hamiltonians in the chaotic regime.


title: Stationary solution and $H$ theorem for a generalized Fokker-Planck equation

abstract: We investigate a family of generalized Fokker-Planck equations that contains Richardson and porous media equations as members. Considering a confining drift term that is related to an effective potential, we show that each equation of this family has a stationary solution that depends on this potential. This stationary solution encompasses several well-known probability distributions. Moreover, we verify an $H$ theorem for the generalized Fokker-Planck equations using free-energy-like functionals. We show that the energy-like part of each functional is based on the effective potential and the entropy-like part is a generalized Tsallis entropic form, which has an unusual dependence on the position and can be related to a generalization of the Kullback-Leibler divergence. We also verify that the optimization of this entropic-like form subjected to convenient constraints recovers the stationary solution. The analysis presented here includes several studies about $H$ theorems for other generalized Fokker-Planck equations as particular cases.


title: Extended Version of GTGraffiti: Spray Painting Graffiti Art from Human Painting Motions with a Cable Driven Parallel Robot

abstract: We present GTGraffiti, a graffiti painting system from Georgia Tech that tackles challenges in art, hardware, and human-robot collaboration. The problem of painting graffiti in a human style is particularly challenging and requires a system-level approach because the robotics and art must be designed around each other. The robot must be highly dynamic over a large workspace while the artist must work within the robot's limitations. Our approach consists of three stages: artwork capture, robot hardware, and planning & control. We use motion capture to capture collaborator painting motions which are then composed and processed into a time-varying linear feedback controller for a cable-driven parallel robot (CDPR) to execute. In this work, we will describe the capturing process, the design and construction of a purpose-built CDPR, and the software for turning an artist's vision into control commands. Our work represents an important step towards faithfully recreating human graffiti artwork by demonstrating that we can reproduce artist motions up to 2m/s and 20m/s$^2$ within 9.3mm RMSE to paint artworks. Changes to the submitted manuscript are colored in blue.


title: Radiation GRMHD Simulations of the Hard State of Black Hole X-ray Binaries and the Collapse of a Hot Accretion Flow

abstract: We present global radiation GRMHD simulations of strongly magnetized accretion onto a spinning, stellar mass black hole at sub-Eddington rates. Using a frequency-dependent Monte Carlo procedure for Compton scattering, we self-consistently evolve a two-temperature description of the ion-electron fluid and its radiation field. For an Eddington ratio $L/L_{\rm Edd} \gtrsim 10^{-3}$, the emergent spectrum forms an apparent power law shape from thermal Comptonization up to a cutoff at $\simeq 100$ keV, characteristic of that seen in the hard spectral states of black hole X-ray binary systems. At these luminosities, the radiative efficiency is high ($\approx 24\%$) and results in a denser midplane region where magnetic fields are dynamically important. For $L/L_{\rm Edd} \sim 10^{-2}$, our hot accretion flow appears to undergo thermal runaway and collapse. Our simulations demonstrate that hot accretion flows can be radiatively efficient and provide an estimate of their maximum luminosity.


title: Singularities of Ricci flow and diffeomorphisms

abstract: Comparing and recognizing metrics can be extraordinarily difficult because of the group of diffeomorphisms. Two metrics, that could even be the same, could look completely different in different coordinates. This is the gauge problem. The general gauge problem is extremely subtle for non-compact spaces. Often it can be avoided if one uses some additional structure of the particular situation. However, in many problems there is no additional structure. Instead we solve the gauge problem directly in great generality. The techniques and ideas apply to many problems. We use them to solve a well-known open problem in Ricci flow: Strong rigidity of cylinders. Strong rigidity is an illustration of a {\it shrinker principle} that uniqueness radiates out from a compact set. It implies that if one tangent flow at a future singular point is a cylinder, then all tangent flows are. We solve the gauge problem by solving a nonlinear system of PDEs. The PDE produces a diffeomorphism that fixes an appropriate gauge in the spirit of the slice theorem for group actions. We then show optimal bounds for the displacement function of the diffeomorphism. Strong rigidity relies on gauge fixing and several other new ideas. One of these is "propagation of almost splitting", another is quadratic rigidity in the right gauge, and a third is an optimal polynomial growth bound for PDEs that holds in great generality.


title: Incremental Abstraction in Distributed Probabilistic SLAM Graphs

abstract: Scene graphs represent the key components of a scene in a compact and semantically rich way, but are difficult to build during incremental SLAM operation because of the challenges of robustly identifying abstract scene elements and optimising continually changing, complex graphs. We present a distributed, graph-based SLAM framework for incrementally building scene graphs based on two novel components. First, we propose an incremental abstraction framework in which a neural network proposes abstract scene elements that are incorporated into the factor graph of a feature-based monocular SLAM system. Scene elements are confirmed or rejected through optimisation and incrementally replace the points yielding a more dense, semantic and compact representation. Second, enabled by our novel routing procedure, we use Gaussian Belief Propagation (GBP) for distributed inference on a graph processor. The time per iteration of GBP is structure-agnostic and we demonstrate the speed advantages over direct methods for inference of heterogeneous factor graphs. We run our system on real indoor datasets using planar abstractions and recover the major planes with significant compression.


title: Upper tail of the spectral radius of sparse Erd\H{o}s-R\'{e}nyi graphs

abstract: We consider an Erd\H{o}s-R\'{e}nyi graph $\mathbb{G}(n,p)$ on $n$ vertices with edge probability $p$ such that \[ \sqrt{\frac{\log n}{\log \log n}} \ll np \le n^{1/2-o(1)}, \label{eq:abs} \tag{$\dagger$} \] and derive the upper tail large deviations of $\lambda(\mathbb{G}(n,p))$, the largest eigenvalue of its adjacency matrix. Within this regime we show that, for $p \gg n^{-2/3}$ the $\log$-probability of the upper tail event of $\lambda(\mathbb{G}(n,p))$ equals to that of planting a clique of an appropriate size (upon ignoring smaller order terms), while for $p \ll n^{-2/3}$ the same is given by that of the existence of a high degree vertex. This, in particular, shows an emergence of {\em non-planted localized structure} in the latter regime. We also confirm that in the entire regime \eqref{eq:abs} the large deviation probability is asymptotically approximated by the solution of the mean-field variational problem, and further identify the typical structure of $\mathbb{G}(n,p)$ conditioned on the upper tail event of $\lambda(\mathbb{G}(n,p))$ in a certain sub-regime of $p$. For $p$ such that $\log(np) \gtrsim \log n$ the large deviations of $\lambda(\mathbb{G}(n,p))$ is deduced from those of ${\rm Hom}(C_{2t}, \mathbb{G}(n,p))$, the homomorphism counts of cycle of length $2t$, for $t \ge 3$ and $p$ such that $n^{1/2-o(1)} \ge np \gg n^{1/t}$. In this latter regime the typical structure of $\mathbb{G}(n,p)$ conditioned on the upper tail of ${\rm Hom}(C_{2t}, \mathbb{G}(n,p))$ is identified and asymptotic tightness of the mean-field approximation is also established.


title: KroneckerBERT: Learning Kronecker Decomposition for Pre-trained Language Models via Knowledge Distillation

abstract: The development of over-parameterized pre-trained language models has made a significant contribution toward the success of natural language processing. While over-parameterization of these models is the key to their generalization power, it makes them unsuitable for deployment on low-capacity devices. We push the limits of state-of-the-art Transformer-based pre-trained language model compression using Kronecker decomposition. We use this decomposition for compression of the embedding layer, all linear mappings in the multi-head attention, and the feed-forward network modules in the Transformer layer. We perform intermediate-layer knowledge distillation using the uncompressed model as the teacher to improve the performance of the compressed model. We present our KroneckerBERT, a compressed version of the BERT_BASE model obtained using this framework. We evaluate the performance of KroneckerBERT on well-known NLP benchmarks and show that for a high compression factor of 19 (5% of the size of the BERT_BASE model), our KroneckerBERT outperforms state-of-the-art compression methods on the GLUE. Our experiments indicate that the proposed model has promising out-of-distribution robustness and is superior to the state-of-the-art compression methods on SQuAD.


title: Neutrino interactions in liquid scintillators including active-sterile neutrino mixing

abstract: Neutrinos play an important role in core-collapse supernova events since they are a key piece to understand the explosion mechanisms. The analysis of the neutrino fluxes can bring answers to neutrino's related problems e.g.: mass hierarchy, spectral splitting, sterile neutrinos, etc. In this work we study the impact of neutrino oscillations and the possible existence of eV sterile neutrinos upon the supernova neutrino flux $F_\nu(E)$. We have calculated the energy distribution of the neutrino flux from a supernova and the total number of events which would be detected in a liquid scintillator. We also present an analysis for the conversion probabilities as a function of the active-sterile neutrino mixing parameters. Finally, we have carried out a statistical analysis to extract values for the mixing parameters of the model.


title: Student networks on online teaching due to COVID-19: Academic effects of strong friendship ties and perceived academic prestige in physics and mathematics courses

abstract: Collaboration among students is fundamental for knowledge building and competency development. Nonetheless, the effectiveness of student collaboration depends on the extent that these interactions take place under conditions that favor commitment, trust, and decision-making among those who interact. The sanitary situation and the transition to remote teaching has added new challenges for collaboration given that students' interactions are mediated by Information and Communication Technologies (ICTs). In this study we explore the effectiveness of different collaborative relationships on physics and mathematics, from a sample of secondary students from two schools located in rural and urban areas in southern Chile. We used Social Network Analysis (SNA) to map students' friendships relations, academic prestige, and collaboration on both courses. Later we combined the collaboration network with friendship and academic prestige on the course, to separate strong from weak friendship working ties, and those among students who enjoy or not academic prestige. Multiple linear regression models showed, on average, positive effects of collaboration on grades. Yet, when isolating the effects of the types of collaboration, the positive effects are observed only between those who display more strong friendship ties. Also, we found differences on the social networks and their effects over grades between both courses, presumably due to their pedagogical nature. With these results we contribute to the literature of collaboration and its effectiveness based on the nature of students' relationships, and advocate for the importance of instructional design in fostering appropriate motivations and guidelines for constructive collaboration in the classroom.


title: Resolving dichotomy in compact objects through continuous gravitational waves observation

abstract: More than two dozen soft gamma-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs) have been detected so far. These are isolated compact objects. Many of them are either found to be associated with supernova remnants or their surface magnetic fields are directly measured, confirming that they are neutron stars (NSs). However, it has been argued that some SGRs and AXPs are highly magnetized white dwarfs (WDs). Meanwhile, the existence of super-Chandrasekhar WDs has remained to be a puzzle. However, not even a single such massive WD has been observed directly. Moreover, some WD pulsars are detected in electromagnetic surveys and some of their masses are still not confirmed. Here we calculate the signal-to-noise ratio for all these objects, considering different magnetic field configurations and thereby estimate the required time for their detection by various gravitational wave (GW) detectors. For SGRs and AXPs, we show that, if these are NSs, they can hardly be detected by any of the GW detectors, while if they are WDs, Big Bang Observer (BBO), DECi-hertz Interferometer Gravitational wave Observatory (DECIGO) and Advanced Laser Interferometer Antenna (ALIA) would be able to detect them within a few days to a year of integration, depending on the magnetic field strength and its configuration. Similarly, if a super-Chandrasekhar WD has a dominant toroidal field, we show that even Laser Interferometer Space Antenna (LISA) and TianQin would be able to detect it within one year of integration. We also discuss how GWs can confirm the masses of the WD pulsars.


title: DIISC-I: The Discovery of Kinematically Anomalous HI Clouds in M 100

abstract: We report the discovery of two kinematically anomalous atomic hydrogen (HI) clouds in M 100 (NGC 4321), which was observed as part of the Deciphering the Interplay between the Interstellar medium, Stars, and the Circumgalactic medium (DIISC) survey in HI 21 cm at 3.3 km s$^{-1}$ spectroscopic and 44 arcsec$\times$30 arcsec spatial resolution using the Karl G. Jansky Very Large Array. These clouds were identified as structures that show significant kinematic offsets from the rotating disk of M100. The velocity offsets of 40 km s$^{-1}$ observed in these clouds are comparable to the offsets seen in intermediate-velocity clouds (IVCs) in the circumgalactic medium (CGM) of the Milky Way and nearby galaxies. We find that one anomalous cloud in M 100 is associated with star-forming regions detected in H$\alpha$ and far-ultraviolet imaging. Our investigation shows that anomalous clouds in M 100 may originate from multiple mechanisms, such as star formation feedback-driven outflows, ram-pressure stripping, and tidal interactions with satellite galaxies. Moreover, we do not detect any cool CGM at 38.8 kpc from the center of M 100, giving an upper limit of N(HI) $\le$ $1.7\times10^{13}$ cm$^{-2}$ (3$\sigma$). Since M 100 is in the Virgo cluster, the non-existence of neutral/cool CGM is a likely pathway for turning it into a red galaxy.


title: Distilling GHZ States using Stabilizer Codes

abstract: Entanglement distillation is a well-studied problem in quantum information, where one typically starts with $n$ noisy Bell pairs and distills $k$ Bell pairs of higher fidelity. While distilling Bell pairs is the canonical setting, it is important to study the distillation of multipartite entangled states because these can be useful for realizing distributed algorithms on quantum networks. In this paper, we study the distillation of GHZ states using quantum error correcting codes (QECCs). Using the stabilizer formalism, we begin by explaining the QECC-based Bell pair distillation protocol in arXiv:0708.3699, which relies particularly on the transpose symmetry between Alice's and Bob's qubits in Bell states. Extending this idea, we show that, given $n$ GHZ states, performing a matrix on Alice's qubits is equivalent to performing a "stretched" version of the transpose of the matrix on the qubits of Bob and Charlie. We call this mapping to the stretched version of the matrix the GHZ-map, and show that it is an algebra homomorphism. Using this property, we show that Alice projecting her qubits onto an $[[n,k]]$ stabilizer code implies the simultaneous projection of Bob's and Charlie's qubits onto an induced $[[2n,k]]$ stabilizer code. Guided by this insight, we develop a GHZ distillation protocol based on local operations and classical communication that uses any stabilizer code. Inspired by stabilizer measurements on GHZ states, we also develop a new algorithm to generate logical Pauli operators of any stabilizer code and use it in the protocol. Since quantum codes with finite rate and almost linear minimum distance have recently been discovered, this paper paves the way for high-rate high-output-fidelity GHZ distillation. We provide simulation results on the $5$-qubit perfect code to emphasize the importance of the placement of a certain local Clifford operation in the protocol.


title: Ricochets on Asteroids II: Sensitivity of laboratory experiments of low velocity grazing impacts on substrate grain size

abstract: We compare low velocity impacts that ricochet with the same impact velocity and impact angle into granular media with similar bulk density, porosity and friction coefficient but different mean grain size. The ratio of projectile diameter to mean grain length ranges from 4 in our coarsest medium to 50 in our finest sand. Using high speed video and fluorescent markers, we measure the ratio of pre- to post-impact horizontal and vertical velocity components, which we refer to as coefficients of restitution, and the angle of deflection caused by the impact in the horizontal plane. Coefficients of restitution are sensitive to mean grain size with the ratio associated with the horizontal velocity component about twice as large for our coarsest gravel as that for our finest sand. This implies that coefficients for hydro-static-like, drag-like and lift-like forces, used in empirical force laws, are sensitive to mean grain size. The coefficient that is most strongly sensitive to grain size is the lift coefficient which decreases by a factor of 3 between our coarsest and finest media. The deflection angles are largest in the coarser media and their size approximately depends on grain size to the 3/2 power. This scaling is matched with a model where momentum transfer takes place via collisions with individual grains. The dependence of impact mechanics on substrate size distribution should be considered in future models for populations of objects that impact granular asteroid surfaces.


title: TNS: Terrain Traversability Mapping and Navigation System for Autonomous Excavators

abstract: We present a terrain traversability mapping and navigation system (TNS) for autonomous excavator applications in an unstructured environment. We use an efficient approach to extract terrain features from RGB images and 3D point clouds and incorporate them into a global map for planning and navigation. Our system can adapt to changing environments and update the terrain information in real-time. Moreover, we present a novel dataset, the Complex Worksite Terrain (CWT) dataset, which consists of RGB images from construction sites with seven categories based on navigability. Our novel algorithms improve the mapping accuracy over previous SOTA methods by 4.17-30.48% and reduce MSE on the traversability map by 13.8-71.4%. We have combined our mapping approach with planning and control modules in an autonomous excavator navigation system and observe 49.3% improvement in the overall success rate. Based on TNS, we demonstrate the first autonomous excavator that can navigate through unstructured environments consisting of deep pits, steep hills, rock piles, and other complex terrain features.


title: Bootstrapping More QM Systems

abstract: We test the bootstrap approach for determining the spectrum of one dimensional Hamiltonians. In this paper we focus on problems that have a two parameter search space in the bootstrap approach: the double well and a periodic potential associated to the Mathieu equation. For the double well, we compare the energies with contributions from perturbative and non-perturbative results, finding good agreement. For the periodic potentials, we notice that the bootstrap approach gives the band structure of the periodic potential, but it has trouble finding the quasi-momentum of the system. To make further progress on the dispersion relation of the bands, new techniques are needed.


title: Affine mobi spaces

abstract: The category of mobi algebras has been introduced as a model to the unit interval of real numbers. The notion of mobi space over a mobi algebra has been proposed as a model for spaces with geodesic paths. In this paper we analyse the particular case of affine mobi spaces and show that there is an isomorphism of categories between R-modules and pointed affine mobi spaces over a mobi algebra R as soon as R is a unitary ring in which 2 is an invertible element.


title: Multi-Sentence Resampling: A Simple Approach to Alleviate Dataset Length Bias and Beam-Search Degradation

abstract: Neural Machine Translation (NMT) is known to suffer from a beam-search problem: after a certain point, increasing beam size causes an overall drop in translation quality. This effect is especially pronounced for long sentences. While much work was done analyzing this phenomenon, primarily for autoregressive NMT models, there is still no consensus on its underlying cause. In this work, we analyze errors that cause major quality degradation with large beams in NMT and Automatic Speech Recognition (ASR). We show that a factor that strongly contributes to the quality degradation with large beams is \textit{dataset length-bias} - \textit{NMT datasets are strongly biased towards short sentences}. To mitigate this issue, we propose a new data augmentation technique -- \textit{Multi-Sentence Resampling (MSR)}. This technique extends the training examples by concatenating several sentences from the original dataset to make a long training example. We demonstrate that MSR significantly reduces degradation with growing beam size and improves final translation quality on the IWSTL$15$ En-Vi, IWSTL$17$ En-Fr, and WMT$14$ En-De datasets.


title: Parametric Modeling Approach to COVID-19 Pandemic Data

abstract: The problem of skewness is common among clinical trials and survival data which has being the research focus derivation and proposition of different flexible distributions. Thus, a new distribution called Extended Rayleigh Lomax distribution is constructed from Rayleigh Lomax distribution to capture the excessiveness of some survival data. We derive the new distribution by using beta logit function proposed by Jones (2004). Some statistical properties of the distribution such as probability density function, cumulative density function, reliability rate, hazard rate, reverse hazard rate, moment generating functions, likelihood functions, skewness, kurtosis and coefficient of variation are obtained. We also performed the expected estimation of model parameters by maximum likelihood; goodness of fit and model selection criteria including Anderson Darling (AD), CramerVon Misses (CVM), Kolmogorov Smirnov (KS), Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC) and Consistent Akaike Information Criterion (CAIC) are employed to select the better distribution from those models considered in the work. The results from the statistics criteria show that the proposed distribution performs better with better representation of the States in Nigeria COVID-19 death cases data than other competing models.


title: Implicit Regularization Effects of the Sobolev Norms in Image Processing

abstract: In this paper, we propose to use the general $L^2$-based Sobolev norms, i.e., $H^s$ norms where $s\in \mathbb{R}$, to measure the data discrepancy due to noise in image processing tasks that are formulated as optimization problems. As opposed to a popular trend of developing regularization methods, we emphasize that an implicit regularization effect can be achieved through the class of Sobolev norms as the data-fitting term. Specifically, we analyze that the implicit regularization comes from the weights that the $H^s$ norm imposes on different frequency contents of an underlying image. We further analyze the underlying noise assumption of using the Sobolev norm as the data-fitting term from a Bayesian perspective, build the connections with the Sobolev gradient-based methods and discuss the preconditioning effects on the convergence rate of the gradient descent algorithm, leading to a better understanding of functional spaces/metrics and the optimization process involved in image processing. Numerical results in full waveform inversion, image denoising and deblurring demonstrate the implicit regularization effects.


title: The Hubble Tension, The $M$ Crisis of Late Time $H(z)$ Deformation Models and the Reconstruction of Quintessence Lagrangians

abstract: We present a detailed and pedagogical analysis of recent cosmological data, including CMB, BAO, SnIa and the recent local measurement of $H_0$. We thus obtain constraints on the parameters of standard dark energy parameterizations, including $\Lambda CDM$, and $H(z)$ deformation models such as $wCDM$ (constant equation of state $w$ of dark energy), and the CPL model (corresponding to the evolving dark energy equation-of-state parameter $w(z) = w_0 + w_a \frac{z}{1+z}$). The fitted parameters include the dark matter density $\Omega_{0m}$, the SnIa absolute magnitude $M$, the Hubble constant $H_0$ and the dark energy parameters (e.g., $w$ for $wCDM$). All models considered lead to a best-fit value of $M$ that is inconsistent with the locally determined value obtained by Cepheid calibrators ($M$ tension). We then use the best-fit dark energy parameters to reconstruct the quintessence Lagrangian that would be able to reproduce these best-fit parameterizations. Due to the derived late phantom behavior of the best-fit dark energy equation-of-state parameter $w(z)$, the reconstructed quintessence models have a negative kinetic term and are therefore plagued with instabilities.


title: Hadrons as QCD Bound States

abstract: Bound state perturbation theory is well established for QED atoms. Today the hyperfine splitting of Positronium is known to $O(\alpha^7\log\alpha)$. Whereas standard expansions of scattering amplitudes start from free states, bound states are expanded around eigenstates of the Hamiltonian including a binding potential. The eigenstate wave functions have all powers of $\alpha$, requiring a choice in the ordering of the perturbative expansion. Temporal $(A^0=0)$ gauge permits an expansion starting from valence Fock states, bound by their instantaneous gauge field. This formulation is applicable in any frame and seems promising even for hadrons in QCD. The $O(\alpha_s^0)$ confining potential is determined (up to a universal scale) by a homogeneous solution of Gauss' law.


title: Unprovability in Mathematics: A First Course on Ordinal Analysis

abstract: These are the lecture notes of an introductory course on ordinal analysis. Our selection of topics is guided by the aim to give a complete and direct proof of a mathematical independence result: Kruskal's theorem for binary trees is unprovable in conservative extensions of Peano arithmetic (note that much stronger results of this type are due to Harvey Friedman). Concerning prerequisites, we assume a solid introduction to mathematical logic but no specialized knowledge of proof theory. The material in these notes is intended for 12 lectures and 6 exercise sessions of 90 minutes each.


title: A new look at Ternary Boolean Algebras

abstract: We present a new approach to ternary Boolean algebras in which negation is derived from the ternary operation. The key aspect is the replacement of complete commutativity by other axioms that do not require the ternary operation to be symmetric.


title: Quantum anonymous veto: A set of new protocols

abstract: We propose a set of protocols for quantum anonymous veto (QAV) broadly categorized under the probabilistic, iterative, and deterministic schemes. The schemes are based upon different types of quantum resources. Specifically, they may be viewed as single photon-based, bipartite and multipartite entangled states-based, orthogonal state-based and conjugate coding-based. The set of the proposed schemes is analyzed for all the requirements of a valid QAV scheme (e.g., privacy, verifiability, robustness, binding, eligibility and correctness). The proposed schemes are observed to be more efficient in comparison to the existing QAV schemes and robust up to the moderate decoherence rate. In addition, a trade-off between correctness and robustness of the probabilistic QAV schemes is observed. Further, the multipartite dense coding based determinsitic QAV scheme is most efficient scheme among the set of schemes proposed here. A bipartite entanglement based iterative scheme employing dense coding is yet another efficient and practical scheme. The intrinsic connections between dining cryptographer-net with anonymous veto-net is also explored in the process of designing new protocols.


title: Asymptotic nonlocality in gauge theories

abstract: Asymptotically nonlocal field theories represent a sequence of higher-derivative theories whose limit point is a ghost-free, infinite-derivative theory. Here we extend this framework, developed previously in a theory of real scalar fields, to gauge theories. We focus primarily on asymptotically nonlocal scalar electrodynamics, first identifying equivalent gauge-invariant formulations of the Lagrangian, one with higher-derivative terms and the other with auxiliary fields instead. We then study mass renormalization of the complex scalar field in each formulation, showing that an emergent nonlocal scale (i.e., one that does not appear as a fundamental parameter in the Lagrangian of the finite-derivative theories) regulates loop integrals as the limiting theory is approached, so that quadratic divergences can be hierarchically smaller than the lightest Lee-Wick partner. We conclude by making preliminary remarks on the generalization of our approach to non-Abelian theories, including an asymptotically nonlocal standard model.


title: Evaluating Multiway Multilingual NMT in the Turkic Languages

abstract: Despite the increasing number of large and comprehensive machine translation (MT) systems, evaluation of these methods in various languages has been restrained by the lack of high-quality parallel corpora as well as engagement with the people that speak these languages. In this study, we present an evaluation of state-of-the-art approaches to training and evaluating MT systems in 22 languages from the Turkic language family, most of which being extremely under-explored. First, we adopt the TIL Corpus with a few key improvements to the training and the evaluation sets. Then, we train 26 bilingual baselines as well as a multi-way neural MT (MNMT) model using the corpus and perform an extensive analysis using automatic metrics as well as human evaluations. We find that the MNMT model outperforms almost all bilingual baselines in the out-of-domain test sets and finetuning the model on a downstream task of a single pair also results in a huge performance boost in both low- and high-resource scenarios. Our attentive analysis of evaluation criteria for MT models in Turkic languages also points to the necessity for further research in this direction. We release the corpus splits, test sets as well as models to the public.


title: Molecules with ALMA at Planet-forming Scales (MAPS) VIII: CO Gap in AS 209--Gas Depletion or Chemical Processing?

abstract: Emission substructures in gas and dust are common in protoplanetary disks. Such substructures can be linked to planet formation or planets themselves. We explore the observed gas substructures in AS 209 using thermochemical modeling with RAC2D and high-spatial resolution data from the Molecules with ALMA at Planet-forming Scales(MAPS) program. The observations of C$^{18}$O J=2-1 emission exhibit a strong depression at 88 au overlapping with the positions of multiple gaps in millimeter dust continuum emission. We find that the observed CO column density is consistent with either gas surface-density perturbations or chemical processing, while C$_2$H column density traces changes in the C/O ratio rather than the H$_2$ gas surface density. However, the presence of a massive planet (> 0.2 M$_{Jup}$) would be required to account for this level of gas depression, which conflicts with constraints set by the dust emission and the pressure profile measured by gas kinematics. Based on our models, we infer that a local decrease of CO abundance is required to explain the observed structure in CO, dominating over a possible gap-carving planet present and its effect on the H$_2$ surface density. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement Series.


title: Post-OCR Document Correction with large Ensembles of Character Sequence-to-Sequence Models

abstract: In this paper, we propose a novel method based on character sequence-to-sequence models to correct documents already processed with Optical Character Recognition (OCR) systems. The main contribution of this paper is a set of strategies to accurately process strings much longer than the ones used to train the sequence model while being sample- and resource-efficient, supported by thorough experimentation. The strategy with the best performance involves splitting the input document in character n-grams and combining their individual corrections into the final output using a voting scheme that is equivalent to an ensemble of a large number of sequence models. We further investigate how to weigh the contributions from each one of the members of this ensemble. We test our method on nine languages of the ICDAR 2019 competition on post-OCR text correction and achieve a new state-of-the-art performance in five of them. Our code for post-OCR correction is shared at https://github.com/jarobyte91/post_ocr_correction.


title: Program-to-Circuit: Exploiting GNNs for Program Representation and Circuit Translation

abstract: Circuit design is complicated and requires extensive domain-specific expertise. One major obstacle stuck on the way to hardware agile development is the considerably time-consuming process of accurate circuit quality evaluation. To significantly expedite the circuit evaluation during the translation from behavioral languages to circuit designs, we formulate it as a Program-to-Circuit problem, aiming to exploit the representation power of graph neural networks (GNNs) by representing C/C++ programs as graphs. The goal of this work is four-fold. First, we build a standard benchmark containing 40k C/C++ programs, each of which is translated to a circuit design with actual hardware quality metrics, aiming to facilitate the development of effective GNNs targeting this high-demand circuit design area. Second, 14 state-of-the-art GNN models are analyzed on the Program-to-Circuit problem. We identify key design challenges of this problem, which should be carefully handled but not yet solved by existing GNNs. The goal is to provide domain-specific knowledge for designing GNNs with suitable inductive biases. Third, we discuss three sets of real-world benchmarks for GNN generalization evaluation, and analyze the performance gap between standard programs and the real-case ones. The goal is to enable transfer learning from limited training data to real-world large-scale circuit design problems. Fourth, the Program-to-Circuit problem is a representative within the Program-to-X framework, a set of program-based analysis problems with various downstream tasks. The in-depth understanding of strength and weaknesses in applying GNNs on Program-to-Circuit could largely benefit the entire family of Program-to-X. Pioneering in this direction, we expect more GNN endeavors to revolutionize this high-demand Program-to-Circuit problem and to enrich the expressiveness of GNNs on programs.


title: Automatic Tuning of Tensorflow's CPU Backend using Gradient-Free Optimization Algorithms

abstract: Modern deep learning (DL) applications are built using DL libraries and frameworks such as TensorFlow and PyTorch. These frameworks have complex parameters and tuning them to obtain good training and inference performance is challenging for typical users, such as DL developers and data scientists. Manual tuning requires deep knowledge of the user-controllable parameters of DL frameworks as well as the underlying hardware. It is a slow and tedious process, and it typically delivers sub-optimal solutions. In this paper, we treat the problem of tuning parameters of DL frameworks to improve training and inference performance as a black-box optimization problem. We then investigate applicability and effectiveness of Bayesian optimization (BO), genetic algorithm (GA), and Nelder-Mead simplex (NMS) to tune the parameters of TensorFlow's CPU backend. While prior work has already investigated the use of Nelder-Mead simplex for a similar problem, it does not provide insights into the applicability of other more popular algorithms. Towards that end, we provide a systematic comparative analysis of all three algorithms in tuning TensorFlow's CPU backend on a variety of DL models. Our findings reveal that Bayesian optimization performs the best on the majority of models. There are, however, cases where it does not deliver the best results.


title: Are Group Acknowledgements Worth Anything in IEEE 802.15.4 DSME: A Comparative Analysis

abstract: For data collection scenarios in the Industrial Internet of Things, wireless communication provides a cost-effective and easy-to-deploy alternative to wired networks. The main focus lies on energy efficiency and reliability, as many devices are battery operated. IEEE 802.15.4 DSME enhances reliability by acknowledging each packet individually, imposing an overhead for each transmitted packet, and increasing energy consumption. In networks with little interference, it may be beneficial to aggregate the acknowledgments for multiple nodes and broadcast them in a compressed format to all nodes in the neighborhood. The IEEE 802.15.4 2012 standard describes such a group acknowledgment scheme which, however, disappears in later iterations of the standard. This paper compares different group acknowledgment schemes and proposes a novel group acknowledgment scheme with the goal to examine whether group acknowledgments constitute a viable alternative to regular acknowledgments in reliable data-collection scenarios. Our analysis suggests that apart from a few cases, GACKs do not constitute a valid alternative to the direct acknowledgement of data packets.


title: Molecules with ALMA at Planet-forming Scales (MAPS) I: Program Overview and Highlights

abstract: Planets form and obtain their compositions in dust and gas-rich disks around young stars, and the outcome of this process is intimately linked to the disk chemical properties. The distributions of molecules across disks regulate the elemental compositions of planets, including C/N/O/S ratios and metallicity (O/H and C/H), as well as access to water and prebiotically relevant organics. Emission from molecules also encodes information on disk ionization levels, temperature structures, kinematics, and gas surface densities, which are all key ingredients of disk evolution and planet formation models. The Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program was designed to expand our understanding of the chemistry of planet formation by exploring disk chemical structures down to 10 au scales. The MAPS program focuses on five disks - around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480 - in which dust substructures are detected and planet formation appears to be ongoing. We observed these disks in 4 spectral setups, which together cover ~50 lines from over 20 different species. This paper introduces the ApJS MAPS Special Issue by presenting an overview of the program motivation, disk sample, observational details, and calibration strategy. We also highlight key results, including discoveries of links between dust, gas, and chemical sub-structures, large reservoirs of nitriles and other organics in the inner disk regions, and elevated C/O ratios across most disks. We discuss how this collection of results is reshaping our view of the chemistry of planet formation.

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


title: A bound for the $p$-domination number of a graph in terms of its eigenvalue multiplicities

abstract: Let $G$ be a connected graph of order $n$ with domination number $\gamma(G)$. Wang, Yan, Fang, Geng and Tian [Linear Algebra Appl. 607 (2020), 307-318] showed that for any Laplacian eigenvalue $\lambda$ of $G$ with multiplicity $m_G(\lambda)$, it holds that $\gamma(G)\leq n-m_G(\lambda)$. Using techniques from the theory of star sets, in this work we prove that the same bound holds when $\lambda$ is an arbitrary adjacency eigenvalue of a non-regular graph, and we characterize the cases of equality. Moreover, we show a result that gives a relationship between start sets and the $p$-domination number, and we apply it to extend the aforementioned spectral bound to the $p$-domination number using the adjacency and Laplacian eigenvalue multiplicities.


title: STraTA: Self-Training with Task Augmentation for Better Few-shot Learning

abstract: Despite their recent successes in tackling many NLP tasks, large-scale pre-trained language models do not perform as well in few-shot settings where only a handful of training examples are available. To address this shortcoming, we propose STraTA, which stands for Self-Training with Task Augmentation, an approach that builds on two key ideas for effective leverage of unlabeled data. First, STraTA uses task augmentation, a novel technique that synthesizes a large amount of data for auxiliary-task fine-tuning from target-task unlabeled texts. Second, STraTA performs self-training by further fine-tuning the strong base model created by task augmentation on a broad distribution of pseudo-labeled data. Our experiments demonstrate that STraTA can substantially improve sample efficiency across 12 few-shot benchmarks. Remarkably, on the SST-2 sentiment dataset, STraTA, with only 8 training examples per class, achieves comparable results to standard fine-tuning with 67K training examples. Our analyses reveal that task augmentation and self-training are both complementary and independently effective.


title: Comparison of $pp$ and $p\bar{p}$ differential elastic cross sections and observation of colourless C-odd gluonic exchanges

abstract: The TOTEM 2.76, 7, 8, and 13 TeV pp elastic differential cross sections are extrapolated using a data-driven approach to obtain the 1.96 TeV pp elastic cross section. A difference with a 3.4$\sigma$ significance is observed between the extrapolated pp elastic cross section and the D0 p$\bar{\rm p}$ elastic cross section at 1.96 TeV in the region of the diffractive minimum and the second maximum of the pp cross section, providing evidence for colourless C-odd gluonic exchanges, also denoted as odderon. These results are combined with a TOTEM analysis of the same $t$-channel C-odd exchanges based on forward pp elastic scattering at TeV scale. The combined significance is larger than 5$\sigma$ and is interpreted as the first observation of odderon exchange.


title: Bipartite dimer model: perfect t-embeddings and Lorentz-minimal surfaces

abstract: This is the second paper in the series devoted to the study of the dimer model on t-embeddings of planar bipartite graphs. We introduce the notion of perfect t-embeddings and assume that the graphs of the associated origami maps converge to a Lorentz-minimal surface $\mathrm{S}_\xi$ as $\delta\to 0$. In this setup we prove (under very mild technical assumptions) that the gradients of the height correlation functions converge to those of the Gaussian Free Field defined in the intrinsic metric of the surface $\mathrm{S}_\xi$. We also formulate several open questions motivated by our work.


title: The role of neutral hydrogen in setting the abundances of molecular species in the Milky Way's diffuse interstellar medium. I. Observational constraints from ALMA and NOEMA

abstract: We have complemented existing observations of HI absorption with new observations of HCO$^+$, C$_2$H, HCN, and HNC absorption from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA) in the direction of 20 background radio continuum sources with $4^\circ \leq |b| \leq 81^\circ$ to constrain the atomic gas conditions that are suitable for the formation of diffuse molecular gas. We find that these molecular species form along sightlines where $A_V \gtrsim 0.25$, consistent with the threshold for the HI-to-H$_2$ transition at solar metallicity. Moreover, we find that molecular gas is associated only with structures that have an HI optical depth $> 0.1$, a spin temperature $< 80$ K, and a turbulent Mach number $\gtrsim 2$. We also identify a broad, faint component to the HCO$^+$ absorption in a majority of sightlines. Compared to the velocities where strong, narrow HCO$^+$ absorption is observed, the HI at these velocities has a lower cold neutral medium (CNM) fraction and negligible CO emission. The relative column densities and linewidths of the different molecular species observed here are similar to those observed in previous experiments over a range of Galactic latitudes, suggesting that gas in the solar neighborhood and gas in the Galactic plane are chemically similar. For a select sample of previously-observed sightlines, we show that the absorption line profiles of HCO$^+$, HCN, HNC, and C$_2$H are stable over periods of $\sim 3$ years and $\sim 25$ years, likely indicating that molecular gas structures in these directions are at least $\gtrsim 100$ AU in size


title: Cross-Modality Domain Adaptation for Vestibular Schwannoma and Cochlea Segmentation

abstract: Automatic methods to segment the vestibular schwannoma (VS) tumors and the cochlea from magnetic resonance imaging (MRI) are critical to VS treatment planning. Although supervised methods have achieved satisfactory performance in VS segmentation, they require full annotations by experts, which is laborious and time-consuming. In this work, we aim to tackle the VS and cochlea segmentation problem in an unsupervised domain adaptation setting. Our proposed method leverages both the image-level domain alignment to minimize the domain divergence and semi-supervised training to further boost the performance. Furthermore, we propose to fuse the labels predicted from multiple models via noisy label correction. Our results on the challenge validation leaderboard showed that our unsupervised method has achieved promising VS and cochlea segmentation performance with mean dice score of 0.8261 $\pm$ 0.0416; The mean dice value for the tumor is 0.8302 $\pm$ 0.0772. This is comparable to the weakly-supervised based method.


title: MindCraft: Theory of Mind Modeling for Situated Dialogue in Collaborative Tasks

abstract: An ideal integration of autonomous agents in a human world implies that they are able to collaborate on human terms. In particular, theory of mind plays an important role in maintaining common ground during human collaboration and communication. To enable theory of mind modeling in situated interactions, we introduce a fine-grained dataset of collaborative tasks performed by pairs of human subjects in the 3D virtual blocks world of Minecraft. It provides information that captures partners' beliefs of the world and of each other as an interaction unfolds, bringing abundant opportunities to study human collaborative behaviors in situated language communication. As a first step towards our goal of developing embodied AI agents able to infer belief states of collaborative partners in situ, we build and present results on computational models for several theory of mind tasks.


title: The generalized Ermakov conservative system: A discussion

abstract: Using older and recent results on the integrability of two-dimensional (2d) dynamical systems, we prove that the results obtained in a recent publication concerning the 2d generalized Ermakov system can be obtained as special cases of a more general approach. This approach is geometric and can be used to study efficiently similar dynamical systems.


title: Machine-learning-based identification for initial clustering structure in relativistic heavy-ion collisions

abstract: $\alpha$-clustering structure is a significant topic in light nuclei. A Bayesian convolutional neural network (BCNN) is applied to classify initial non-clustered and clustered configurations, namely Woods-Saxon distribution and three-$\alpha$ triangular (four-$\alpha$ tetrahedral) structure for $^{12}$C ($^{16}$O), from heavy-ion collision events generated within a multi-phase transport (AMPT) model. Azimuthal angle and transverse momentum distributions of charged pions are taken as inputs to train the classifier. On multiple-event basis, the overall classification accuracy can reach $95\%$ for $^{12}$C/$^{16}$O + $^{197}$Au events at $\sqrt{S_{NN}} =$ 200 GeV. With proper constructions of samples, the predicted deviations on mixed samples with different proportions of both configurations could be within $5\%$. In addition, setting a simple confidence threshold can further improve the predictions on the mixed dataset. Our results indicate promising and extensive possibilities of application of machine-learning-based techniques to real data and some other problems in physics of heavy-ion collisions.


title: Semi-biproducts of monoids

abstract: It is shown that the category of \emph{semi-biproducts} of monoids is equivalent to the category of \emph{pseudo-actions}. A semi-biproduct of monoids is a new notion, obtained through generalizing a biproduct of commutative monoids. By dropping commutativity and requiring some of the homomorphisms in the biproduct diagram to be merely identity-preserving maps, we obtain a semi-biproduct. A pseudo-action is a new notion as well. It consists of three ingredients: a pre-action, a factor system and a correction system. In the category of groups all correction systems are trivial. This is perhaps the reason why this notion, to the author's best knowledge, has never been considered before.


title: Interactive All-Hex Meshing via Cuboid Decomposition

abstract: Standard PolyCube-based hexahedral (hex) meshing methods aim to deform the input domain into an axis-aligned PolyCube volume with integer corners; if this deformation is bijective, then applying the inverse map to the voxelized PolyCube yields a valid hex mesh. A key challenge in these methods is to maintain the bijectivity of the PolyCube deformation, thus reducing the robustness of these algorithms. In this work, we present an interactive pipeline for hex meshing that sidesteps this challenge by using a new representation of PolyCubes as unions of cuboids. We begin by deforming the input tetrahedral mesh into a near-PolyCube domain whose faces are close but not perfectly aligned to the major axis directions. We then build a PolyCube by optimizing the layout of a set of cuboids with user guidance to closely fit the deformed domain. Finally, we construct an inversion-free pullback map from the voxelized PolyCube to the input domain while optimizing for mesh quality metrics. We allow extensive user control over each stage, such as editing the voxelized PolyCube, positioning surface vertices, and exploring the trade-off among competing quality metrics, while also providing automatic alternatives. We validate our method on over one hundred shapes, including models that are challenging for past PolyCube-based and frame-field-based methods. Our pipeline reliably produces hex meshes with quality on par with or better than state-of-the-art. We additionally conduct a user study with 20 participants in which the majority prefer hex meshes they make using our tool to the ones from automatic state-of-the-art methods. This demonstrates the need for intuitive interactive hex meshing tools where the user can dictate the priorities of their mesh.


title: Solar differential rotation reproduced with high-resolution simulation

abstract: The Sun rotates differentially with a fast equator and slow pole. Convection in the solar interior is thought to maintain the differential rotation. However, although many numerical simulations have been conducted to reproduce the solar differential rotation, previous high-resolution calculations with solar parameters fall into the anti-solar (fast pole) differential rotation regime. Consequently, we still do not know the true reason why the Sun has a fast-rotating equator. While the construction of the fast equator requires a strong rotational influence on the convection, the previous calculations have not been able to achieve the situation without any manipulations. The problem is called convective conundrum. The convection and the differential rotation in numerical simulations were different from the observations. Here, we show that a high-resolution calculation succeeds in reproducing the solar-like differential rotation. Our calculations indicate that the strong magnetic field generated by a small-scale dynamo has a significant impact on thermal convection. The successful reproduction of the differential rotation, convection, and magnetic field achieved in our calculation is an essential step to understanding the cause of the most basic nature of solar activity, specifically, the 11-year cycle of sunspot activity.


title: Trions in MoS$_2$ are quantum superpositions of intra- and intervalley spin states

abstract: We report magneto-photoluminescence spectroscopy of gated MoS$_2$ monolayers in high magnetic fields to 28 T. At B = 0T and electron density $n_s\sim 10^{12}cm^-2$, we observe three trion resonances that cannot be explained within a single-particle picture. Employing ab initio calculations that take into account three-particle correlation effects as well as local and non-local electron-hole exchange interaction, we identify those features as quantum superpositions of inter- and intravalley spin states. We experimentally investigate the mixed character of the trion wave function via the filling factor dependent valley Zeeman shift in positive and negative magnetic fields. Our results highlight the importance of exchange interactions for exciton physics in monolayer MoS$_2$ and provide new insights into the microscopic understanding of trion physics in 2D multi-valley semiconductors for low excess carrier densities.


title: Phase Transitions of Zirconia: Machine-Learned Force Fields Beyond Density Functional Theory

abstract: We present an approach to generate machine-learned force fields (MLFF) with beyond density functional theory (DFT) accuracy. Our approach combines on-the-fly active learning and $\Delta$-machine learning in order to generate an MLFF for zirconia based on the random phase approximation (RPA). Specifically, an MLFF trained on-the-fly during DFT based molecular dynamics simulations is corrected by another MLFF that is trained on the differences between RPA and DFT calculated energies, forces and stress tensors. Thanks to the relatively smooth nature of the differences, the expensive RPA calculations are performed only on a small number of representative structures of small unit cells. These structures are determined by a singular value decomposition rank compression of the kernel matrix with low spatial resolution. This dramatically reduces the computational cost and allows us to generate an MLFF fully capable of reproducing high-level quantum-mechanical calculations beyond DFT. We carefully validate our approach and demonstrate its success in studying the phase transitions of zirconia.


title: Graph Algorithms for Multiparallel Word Alignment

abstract: With the advent of end-to-end deep learning approaches in machine translation, interest in word alignments initially decreased; however, they have again become a focus of research more recently. Alignments are useful for typological research, transferring formatting like markup to translated texts, and can be used in the decoding of machine translation systems. At the same time, massively multilingual processing is becoming an important NLP scenario, and pretrained language and machine translation models that are truly multilingual are proposed. However, most alignment algorithms rely on bitexts only and do not leverage the fact that many parallel corpora are multiparallel. In this work, we exploit the multiparallelity of corpora by representing an initial set of bilingual alignments as a graph and then predicting additional edges in the graph. We present two graph algorithms for edge prediction: one inspired by recommender systems and one based on network link prediction. Our experimental results show absolute improvements in $F_1$ of up to 28% over the baseline bilingual word aligner in different datasets.


title: Unruh-DeWitt detector's response to a non-relativistic particle

abstract: Unruh-DeWitt (UDW) detector that couples locally to a quantum field is an important tool for operationally studying field properties. Here, we study the response of an inertial UDW detector to a non-relativistic particle state of a massive scalar field and we find that the transition probability of the detector splits into the vacuum contribution and the matter contribution. We show that the matter part oscillates with the interaction time duration and the oscillation period depends on the difference between the mass of the particle and the energy gap of the detector; a strong resonance pattern for the transition probability is found when such difference equals to zero. We compare the matter part contribution with the non-relativistic probability density of the particle and find that they provide qualitatively similar description of the particle.