title: How and When Random Feedback Works: A Case Study of Low-Rank Matrix Factorization

abstract: The success of gradient descent in ML and especially for learning neural networks is remarkable and robust. In the context of how the brain learns, one aspect of gradient descent that appears biologically difficult to realize (if not implausible) is that its updates rely on feedback from later layers to earlier layers through the same connections. Such bidirected links are relatively few in brain networks, and even when reciprocal connections exist, they may not be equi-weighted. Random Feedback Alignment (Lillicrap et al., 2016), where the backward weights are random and fixed, has been proposed as a bio-plausible alternative and found to be effective empirically. We investigate how and when feedback alignment (FA) works, focusing on one of the most basic problems with layered structure -- low-rank matrix factorization. In this problem, given a matrix $Y_{n\times m}$, the goal is to find a low rank factorization $Z_{n \times r}W_{r \times m}$ that minimizes the error $\|ZW-Y\|_F$. Gradient descent solves this problem optimally. We show that FA converges to the optimal solution when $r\ge \mbox{rank}(Y)$. We also shed light on how FA works. It is observed empirically that the forward weight matrices and (random) feedback matrices come closer during FA updates. Our analysis rigorously derives this phenomenon and shows how it facilitates convergence of FA*, a closely related variant of FA. We also show that FA can be far from optimal when $r < \mbox{rank}(Y)$. This is the first provable separation result between gradient descent and FA. Moreover, the representations found by gradient descent and FA can be almost orthogonal even when their error $\|ZW-Y\|_F$ is approximately equal. As a corollary, these results also hold for training two-layer linear neural networks when the training input is isotropic, and the output is a linear function of the input.


title: A Hierarchical Multi-Task Approach to Gastrointestinal Image Analysis

abstract: A large number of different lesions and pathologies can affect the human digestive system, resulting in life-threatening situations. Early detection plays a relevant role in the successful treatment and the increase of current survival rates to, e.g., colorectal cancer. The standard procedure enabling detection, endoscopic video analysis, generates large quantities of visual data that need to be carefully analyzed by an specialist. Due to the wide range of color, shape, and general visual appearance of pathologies, as well as highly varying image quality, such process is greatly dependent on the human operator experience and skill. In this work, we detail our solution to the task of multi-category classification of images from the gastrointestinal (GI) human tract within the 2020 Endotect Challenge. Our approach is based on a Convolutional Neural Network minimizing a hierarchical error function that takes into account not only the finding category, but also its location within the GI tract (lower/upper tract), and the type of finding (pathological finding/therapeutic intervention/anatomical landmark/mucosal views' quality). We also describe in this paper our solution for the challenge task of polyp segmentation in colonoscopies, which was addressed with a pretrained double encoder-decoder network. Our internal cross-validation results show an average performance of 91.25 Mathews Correlation Coefficient (MCC) and 91.82 Micro-F1 score for the classification task, and a 92.30 F1 score for the polyp segmentation task. The organization provided feedback on the performance in a hidden test set for both tasks, which resulted in 85.61 MCC and 86.96 F1 score for classification, and 91.97 F1 score for polyp segmentation. At the time of writing no public ranking for this challenge had been released.


title: Automated skin lesion segmentation using multi-scale feature extraction scheme and dual-attention mechanism

abstract: Segmenting skin lesions from dermoscopic images is essential for diagnosing skin cancer. But the automatic segmentation of these lesions is complicated due to the poor contrast between the background and the lesion, image artifacts, and unclear lesion boundaries. In this work, we present a deep learning model for the segmentation of skin lesions from dermoscopic images. To deal with the challenges of skin lesion characteristics, we designed a multi-scale feature extraction module for extracting the discriminative features. Further in this work, two attention mechanisms are developed to refine the post-upsampled features and the features extracted by the encoder. This model is evaluated using the ISIC2018 and ISBI2017 datasets. The proposed model outperformed all the existing works and the top-ranked models in two competitions.


title: Faint objects in motion: the new frontier of high precision astrometry

abstract: Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of Earth-mass habitable worlds around the nearest stars, to distant Milky Way objects, and out to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local Universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry missions: NEAT proposed for the ESA M3 opportunity, micro-NEAT for the S1 opportunity, and Theia for the M4 and M5 opportunities. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this White Paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review instrumentation and mission profiles.

categories: astro-ph.IM astro-ph.CO astro-ph.EP astro-ph.GA astro-ph.HE


title: CNN Filter Learning from Drawn Markers for the Detection of Suggestive Signs of COVID-19 in CT Images

abstract: Early detection of COVID-19 is vital to control its spread. Deep learning methods have been presented to detect suggestive signs of COVID-19 from chest CT images. However, due to the novelty of the disease, annotated volumetric data are scarce. Here we propose a method that does not require either large annotated datasets or backpropagation to estimate the filters of a convolutional neural network (CNN). For a few CT images, the user draws markers at representative normal and abnormal regions. The method generates a feature extractor composed of a sequence of convolutional layers, whose kernels are specialized in enhancing regions similar to the marked ones, and the decision layer of our CNN is a support vector machine. As we have no control over the CT image acquisition, we also propose an intensity standardization approach. Our method can achieve mean accuracy and kappa values of $0.97$ and $0.93$, respectively, on a dataset with 117 CT images extracted from different sites, surpassing its counterpart in all scenarios.


title: Two-step adversarial debiasing with partial learning -- medical image case-studies

abstract: The use of artificial intelligence (AI) in healthcare has become a very active research area in the last few years. While significant progress has been made in image classification tasks, only a few AI methods are actually being deployed in hospitals. A major hurdle in actively using clinical AI models currently is the trustworthiness of these models. More often than not, these complex models are black boxes in which promising results are generated. However, when scrutinized, these models begin to reveal implicit biases during the decision making, such as detecting race and having bias towards ethnic groups and subpopulations. In our ongoing study, we develop a two-step adversarial debiasing approach with partial learning that can reduce the racial disparity while preserving the performance of the targeted task. The methodology has been evaluated on two independent medical image case-studies - chest X-ray and mammograms, and showed promises in bias reduction while preserving the targeted performance.


title: Automatic Semantic Segmentation of the Lumbar Spine: Clinical Applicability in a Multi-parametric and Multi-centre Study on Magnetic Resonance Images

abstract: One of the major difficulties in medical image segmentation is the high variability of these images, which is caused by their origin (multi-centre), the acquisition protocols (multi-parametric), as well as the variability of human anatomy, the severity of the illness, the effect of age and gender, among others. The problem addressed in this work is the automatic semantic segmentation of lumbar spine Magnetic Resonance images using convolutional neural networks. The purpose is to assign a class label to each pixel of an image. Classes were defined by radiologists and correspond to different structural elements like vertebrae, intervertebral discs, nerves, blood vessels, and other tissues. The proposed network topologies are variants of the U-Net architecture. Several complementary blocks were used to define the variants: Three types of convolutional blocks, spatial attention models, deep supervision and multilevel feature extractor. This document describes the topologies and analyses the results of the neural network designs that obtained the most accurate segmentations. Several of the proposed designs outperform the standard U-Net used as baseline, especially when used in ensembles where the output of multiple neural networks is combined according to different strategies.


title: Exploration of a Singular Fluid Spacetime

abstract: We investigate the properties of a special class of singular solutions for a self-gravitating perfect fluid in general relativity: the singular isothermal sphere. For arbitrary constant equation-of-state parameter $w=p/\rho$, there exist static, spherically-symmetric solutions with density profile $\propto 1/r^2$, with the constant of proportionality fixed to be a special function of $w$. Like black holes, singular isothermal spheres possess a fixed mass-to-radius ratio independent of size, but no horizon cloaking the curvature singularity at $r=0$. For $w=1$, these solutions can be constructed from a homogeneous dilaton background, where the metric spontaneously breaks spatial homogeneity. We study the perturbative structure of these solutions, finding the radial modes and tidal Love numbers, and also find interesting properties in the geodesic structure of this geometry. Finally, connections are discussed between these geometries and dark matter profiles, the double copy, and holographic entropy, as well as how the swampland distance conjecture can obscure the naked singularity.


title: Early-type galaxy density profiles from IllustrisTNG: III. Effects on outer kinematic structure

abstract: Early-type galaxies (ETGs) possess total density profiles close to isothermal, which can lead to non-Gaussian line-of-sight velocity dispersion (LOSVD) under anisotropic stellar orbits. However, recent observations of local ETGs in the MASSIVE Survey reveal outer kinematic structures at $1.5 R_{\mathrm{eff}}$ (effective radius) that are inconsistent with fixed isothermal density profiles; the authors proposed varying density profiles as an explanation. We aim to verify this conjecture and understand the influence of stellar assembly on these kinematic features through mock ETGs in IllustrisTNG. We create mock Integral-Field-Unit observations to extract projected stellar kinematic features for 207 ETGs with stellar mass $M_{\ast}\geqslant 10^{11} \mathrm{M_{\odot}}$ in TNG100-1. The mock observations reproduce the key outer ($1.5R_{\mathrm{eff}}$) kinematic structures in the MASSIVE ETGs, including the puzzling positive correlation between velocity dispersion profile outer slope $\gamma_{\mathrm{outer}}$ and the kurtosis $h_{4}$'s gradient. We find that $h_{4}$ is uncorrelated with stellar orbital anisotropy beyond $R_{\mathrm{eff}}$; instead we find that the variations in $\gamma_{\mathrm{outer}}$ and outer $h_{4}$ (a good proxy for $h_{4}$ gradient) are both driven by variations of the density profile at the outskirts across different ETGs. These findings corroborate the proposed conjecture and rule out velocity anisotropy as the origin of non-Gaussian outer kinematic structure in ETGs. We also find that the outer kurtosis and anisotropy correlate with different stellar assembly components, with the former related to minor mergers or flyby interactions while the latter is mainly driven by major mergers, suggesting distinct stellar assembly origins that decorrelates the two quantities.


title: Enhancing disorder-free localization through dynamically emergent local symmetries

abstract: Disorder-free localization is a recently discovered phenomenon of nonergodicity that can emerge in quantum many-body systems hosting gauge symmetries when the initial state is prepared in a superposition of gauge superselection sectors. Thermalization is then prevented up to all accessible evolution times despite the model being nonintegrable and translation-invariant. In a recent work [Halimeh, Zhao, Hauke, and Knolle, arXiv:2111.02427], it has been shown that terms linear in the gauge-symmetry generator stabilize disorder-free localization in $\mathrm{U}(1)$ gauge theories against gauge errors that couple different superselection sectors. Here, we show in the case of $\mathbb{Z}_2$ gauge theories that disorder-free localization can not only be stabilized, but also \textit{enhanced} by the addition of translation-invariant terms linear in a local $\mathbb{Z}_2$ \textit{pseudogenerator} that acts identically to the full generator in a single superselection sector, but not necessarily outside of it. We show analytically and numerically how this leads through the quantum Zeno effect to the dynamical emergence of a renormalized gauge theory with an enhanced local symmetry, which contains the $\mathbb{Z}_2$ gauge symmetry of the ideal model, associated with the $\mathbb{Z}_2$ pseudogenerator. The resulting proliferation of superselection sectors due to this dynamically emergent gauge theory creates an effective disorder greater than that in the original model, thereby enhancing disorder-free localization. We demonstrate the experimental feasibility of the $\mathbb{Z}_2$ pseudogenerator by providing a detailed readily implementable experimental proposal for the observation of disorder-free localization in a Rydberg setup.

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


title: Secular chaos in white-dwarf planetary systems: origins of metal pollution and short-period planetary companions

abstract: Secular oscillations in multi-planet systems can drive chaotic evolution of a small inner body through non-linear resonant perturbations. This "secular chaos" readily pushes the inner body to an extreme eccentricity, triggering tidal interactions or collision with the central star. We present a numerical study of secular chaos in systems with two planets and test particles using the ring-averaging method, with emphasis on the relationship between the planets' properties and the time-scale and efficiency of chaotic diffusion. We find that secular chaos can excite extreme eccentricities on time-scales spanning several orders of magnitude in a given system. We apply our results to the evolution of planetary systems around white dwarfs (WDs), specifically the tidal disruption and high-eccentricity migration of planetesimals and planets. We find that secular chaos in a planetesimal belt driven by large ($\gtrsim 10 M_{\oplus}$), distant ($\gtrsim 10 \, {\rm au}$) planets can sustain metal accretion onto a WD over Gyr time-scales. We constrain the total mass of planetesimals initially present within the chaotic zone by requiring that the predicted mass delivery rate to the Roche limit be consistent with the observed metal accretion rates of WDs with atmospheric pollution throughout the cooling sequence. Based on the occurrence of long-period exoplanets and exo-asteroid belts, we conclude that secular chaos can be a significant (perhaps dominant) channel for polluting solitary WDs. Secular chaos can also produce short-period planets and planetesimals around WDs in concert with various circularization mechanisms. We discuss prospects for detecting exoplanets driving secular chaos around WDs using direct imaging and microlensing.


title: LEGWORK: A python package for computing the evolution and detectability of stellar-origin gravitational-wave sources with space-based detectors

abstract: We present LEGWORK (LISA Evolution and Gravitational Wave Orbit Kit), an open-source Python package for making predictions about stellar-origin gravitational wave sources and their detectability in LISA or other space-based gravitational wave detectors. LEGWORK can be used to evolve the orbits of sources due to gravitational wave emission, calculate gravitational wave strains (using post-Newtonian approximations), compute signal-to-noise ratios and visualise the results. It can be applied to a variety of potential sources, including binaries consisting of white dwarfs, neutron stars and black holes. Although we focus on double compact objects, in principle LEGWORK can be used for any system with a user-specified orbital evolution, such as those affected by a third object or gas drag. We optimised the package to make it efficient for use in population studies which can contain tens-of-millions of sources. This paper describes the package and presents several potential use cases. We explain in detail the derivations of the expressions behind the package as well as identify and clarify some discrepancies currently present in the literature. We hope that LEGWORK will enable and accelerate future studies triggered by the rapidly growing interest in gravitational wave sources.


title: Detecting low-mass haloes with strong gravitational lensing I: the effect of data quality and lensing configuration

abstract: This paper aims to quantify how the lowest halo mass that can be detected with galaxy-galaxy strong gravitational lensing depends on the quality of the observations and the characteristics of the observed lens systems. Using simulated data, we measure the lowest detectable NFW mass at each location of the lens plane, in the form of detailed \emph{sensitivity maps}. In summary, we find that: (i) the lowest detectable mass $M_{\rm low}$ decreases linearly as the signal-to-noise ratio (SNR) increases and the sensitive area is larger when we decrease the noise; (ii) a moderate increase in angular resolution (0.07" vs 0.09") and pixel scale (0.01" vs 0.04") improves the sensitivity by on average 0.25 dex in halo mass, with more significant improvement around the most sensitive regions; (iii) the sensitivity to low-mass objects is largest for bright and complex lensed galaxies located inside the caustic curves and lensed into larger Einstein rings (i.e $r_{E}\geq1.0"$). We find that for the sensitive mock images considered in this work, the minimum mass that we can detect at the redshift of the lens lies between $1.5\times10^{8}$ and $3\times10^{9}M_{\odot}$. We derive analytic relations between $M_{\rm low}$, the SNR and resolution and discuss the impact of the lensing configuration and source structure. Our results start to fill the gap between approximate predictions and real data and demonstrate the challenging nature of calculating precise forecasts for gravitational imaging. In light of our findings, we discuss possible strategies for designing strong lensing surveys and the prospects for HST, Keck, ALMA, Euclid and other future observations.


title: MASCOT -- An ESO-ARO legacy survey of molecular gas in nearby SDSS-MaNGA galaxies: I. first data release, and global and resolved relations between H_2 and stellar content

abstract: We present the first data release of the MaNGA-ARO Survey of CO Targets (MASCOT), an ESO Public Spectroscopic Survey conducted at the Arizona Radio Observatory (ARO). We measure the CO(1-0) line emission in a sample of 187 nearby galaxies selected from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey that has obtained integral field unit (IFU) spectroscopy for a sample of ~ 10,000 galaxies at low redshift. The main goal of MASCOT is to probe the molecular gas content of star-forming galaxies with stellar masses > 10^9.5 M_solar and with associated MaNGA IFU observations and well-constrained quantities like stellar masses, star formation rates and metallicities. In this paper we present the first results of the MASCOT survey, providing integrated CO(1-0) measurements that cover several effective radii of the galaxy and present CO luminosities, CO kinematics, and estimated H2 gas masses. We observe that the decline of galaxy star formation rate with respect to the star formation main sequence (SFMS) increases with the decrease of molecular gas and with a reduced star formation efficiency, in agreement with results of other integrated studies. Relating the molecular gas mass fractions with the slope of the stellar age gradients inferred from the MaNGA observations, we find that galaxies with lower molecular gas mass fractions tend to show older stellar populations close to the galactic center, while the opposite is true for galaxies with higher molecular gas mass fractions, providing tentative evidence for inside-out quenching.


title: Detection of a $\sim$100,000 M$_\odot$ black hole in M31's most massive globular cluster: A tidally stripped nucleus

abstract: We investigate the presence of a central black hole (BH) in B023-G078, M31's most massive globular cluster. We present high-resolution, adaptive-optics assisted, integral-field spectroscopic kinematics from Gemini/NIFS that shows a strong rotation ($\sim$20 km/s) and a velocity dispersion rise towards the center (37 km/s). We combine the kinematic data with a mass model based on a two-component fit to $HST$ ACS/HRC data of the cluster to estimate the mass of a putative BH. Our dynamical modeling suggests a $>$3$\sigma$ detection of a BH component of 9.1$^{+2.6}_{-2.8}\times$10$^4$ M$_\odot$ (1$\sigma$ uncertainties). The inferred stellar mass of the cluster is 6.22$^{+0.03}_{-0.05}\times$10$^6$ M$_\odot$, consistent with previous estimates, thus the BH makes up 1.5% of its mass. We examine whether the observed kinematics are caused by a collection of stellar mass BHs by modeling an extended dark mass as a Plummer profile. The upper limit on the size scale of the extended mass is 0.56 pc (95% confidence), which does not rule out an extended mass. There is compelling evidence that B023-G078 is the tidally stripped nucleus of a galaxy with a stellar mass $>$10$^9$ M$_{\odot}$, including its high mass, two-component luminosity profile, color, metallicity gradient, and spread in metallicity. Given the emerging evidence that the central BH occupation fraction of $>$10$^9$ M$_{\odot}$ galaxies is high, the most plausible interpretation of the kinematic data is that B023-G078 hosts a central BH. This makes it the strongest BH detection in a lower mass ($<$10$^7$ M$_{\odot}$) stripped nucleus, and one of the few dynamically detected intermediate-mass BHs.


title: Galaxy shapes of Light (GaLight): a 2D modeling of galaxy images

abstract: Galight is a Python-based open-source package that can be used to perform two-dimensional model fitting of optical and near-infrared images to characterize the light distribution of galaxies with components including a disk, bulge, bar, and quasar. The decomposition of stellar components has been demonstrated in published studies of inactive galaxies and quasar host galaxies observed by the Hubble Space Telescope and Subaru's Hyper Suprime-Cam. Galight utilizes the image modeling capabilities of lenstronomy while redesigning the user interface for the analysis of large samples of extragalactic sources. The package is user-friendly with some automatic features such as determining the cutout size of the modeling frame, searching for PSF-stars in field-of-view, estimating the noise map of the data, identifying all the objects to set the initial model, and associated parameters to fit them simultaneously. These features minimize the manpower and allow the automatic fitting tasks. The software is distributed under the MIT license. The source code, installation guidelines, and example notebooks code can be found at https://galight.readthedocs.io/en/latest/


title: Cascade Appearance Signatures of Sterile Neutrinos at 1-100 TeV

abstract: Neutrino telescopes provide strong sensitivity to sterile neutrino oscillations through matter-enhanced oscillation, occurring in the few TeV energy range for eV$^{2}$-scale neutrino mass-squared splittings. Prior searches have focused on $\nu_\mu$ disappearance, which has a particularly strong sensitivity to the mixing angle $\theta_{24}$ via $\nu_\mu\rightarrow\nu_s$ transitions. Nowadays, the $\nu_\mu\rightarrow\nu_e$ and $\nu_\mu\rightarrow\nu_\tau$ appearance channels have been considered less promising at neutrino telescopes, due in part to the much smaller target volume for cascades relative to tracks. This work explores the detectability of these signatures at neutrino telescopes given present constraints on sterile neutrino mixing, and as an example, forecasts the sensitivity of the IceCube Neutrino Observatory to the mixing angles $\theta_{14}$, $\theta_{24}$, and $\theta_{34}$ in the 3+1 sterile neutrino model using the cascade channel with ten years of data. We find that $\nu_\tau$ appearance signatures consistent with the existing IceCube $\nu_\mu$ disappearance best-fit point are discoverable for values of $\theta_{34}$ consistent with world constraints, and that the sterile neutrino parameters favored by the BEST and gallium anomalies are expected to be testable at the 95\% confidence level.


title: Who Decides if AI is Fair? The Labels Problem in Algorithmic Auditing

abstract: Labelled "ground truth" datasets are routinely used to evaluate and audit AI algorithms applied in high-stakes settings. However, there do not exist widely accepted benchmarks for the quality of labels in these datasets. We provide empirical evidence that quality of labels can significantly distort the results of algorithmic audits in real-world settings. Using data annotators typically hired by AI firms in India, we show that fidelity of the ground truth data can lead to spurious differences in performance of ASRs between urban and rural populations. After a rigorous, albeit expensive, label cleaning process, these disparities between groups disappear. Our findings highlight how trade-offs between label quality and data annotation costs can complicate algorithmic audits in practice. They also emphasize the need for development of consensus-driven, widely accepted benchmarks for label quality.


title: Hydrodynamical transport of angular momentum in accretion disks in the presence of nonlinear perturbations due to noise

abstract: The origin of hydrodynamical instability and turbulence in the Keplerian accretion disk is a long-standing puzzle. The flow therein is linearly stable. Here we explore the evolution of perturbation in this flow in the presence of an additional force. Such a force, which is expected to be stochastic in nature hence behaving as noise, could result from thermal fluctuations (however small be), grain-fluid interactions, feedback from outflows in astrophysical disks, etc. We essentially establish the evolution of nonlinear perturbation in the presence of Coriolis and external forces, which is the modified Landau equation. We obtain that even in the linear regime, under suitable forcing and Reynolds number, the otherwise least stable perturbation evolves to a very large saturated amplitude, leading to nonlinearity and plausible turbulence. Hence, forcing essentially leads a linear stable mode to unstable. We further show that nonlinear perturbation diverges at a shorter time-scale in the presence of force, leading to a fast transition to turbulence. Interestingly, the emergence of nonlinearity depends only on the force but not on the initial amplitude of perturbation, unlike the original Landau equation-based solution.


title: Determining the Dark Matter distribution in galaxies with Deep Learning

abstract: We present a novel method to infer the Dark Matter (DM) content and spatial distribution within galaxies, based on convolutional neural networks trained within state-of-the-art hydrodynamical simulations (Illustris TNG100). The framework we have developed is capable of inferring the DM mass distribution within galaxies of mass $~10^{11}-10^{13}M_{\odot}$ with very high performance from the gravitationally baryon dominated internal regions to the DM-rich, baryon-depleted outskirts of the galaxies. With respect to traditional methods, the one presented here also possesses the advantages of not relying on a pre-assigned shape for the DM distribution, to be applicable to galaxies not necessarily in isolation, and to perform very well even in the absence of spectroscopic observations


title: The geometry of geometries: matroid theory, old and new

abstract: The theory of matroids or combinatorial geometries originated in linear algebra and graph theory, and has deep connections with many other areas, including field theory, matching theory, submodular optimization, Lie combinatorics, and total positivity. Matroids capture the combinatorial essence that these different settings share. In recent years, the (classical, polyhedral, algebraic, and tropical) geometric roots of the field have grown much deeper, bearing new fruits. We survey some recent successes, stemming from three geometric models of a matroid: the matroid polytope, the Bergman fan, and the conormal fan.


title: Operator Spreading in the Memory Matrix Formalism

abstract: The spread and scrambling of quantum information is a topic of considerable current interest. Numerous studies suggest that quantum information evolves according to hydrodynamical equations of motion, even though it is a starkly different quantity to better-known hydrodynamical variables such as charge and energy. In this work we show that the well-known memory matrix formalism for traditional hydrodynamics can be applied, with relatively little modification, to the question of operator growth in many-body quantum systems. On a conceptual level, this shores up the connection between information scrambling and hydrodynamics. At a practical level, it provides a framework for calculating quantities related to operator growth like the butterfly velocity and front diffusion constant, and for understanding how these quantities are constrained by microscopic symmetries. We apply this formalism to calculate operator-hydrodynamical coefficients perturbatively in a family of Floquet models. Our formalism allows us to identify the processes affecting information transport that arise from the spatiotemporal symmetries of the model.

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


title: Electromagnetic radiation from axion condensates in a time dependent magnetic field

abstract: Time dependent magnetic fields can be sourced by spinning neutron stars, orbiting binaries and merging neutron stars. We consider electromagnetic radiation from axion condensates in the background of an alternating magnetic field. We find that a resonant peak in radiation can occur when the frequency of the alternating magnetic field is comparable with the axion mass scale. More interestingly, in situations where the frequency of the alternating magnetic field itself changes with time, as can be the case in binary mergers due to steady increase in orbital frequency, the resonant peak in radiation may occur for a range of axion mass scales scanned by the time-varying magnetic field frequency.


title: Operator Growth and Symmetry-Resolved Coefficient Entropy in Quantum Maps

abstract: Operator growth, or operator spreading, describes the process where a "simple" operator acquires increasing complexity under the Heisenberg time evolution of a chaotic dynamics, therefore has been a key concept in the study of quantum chaos in both single-particle and many-body systems. An explicit way to quantify the complexity of an operator is the Shannon entropy of its operator coefficients over a chosen set of operator basis, dubbed "coefficient entropy". However, it remains unclear if the basis-dependency of the coefficient entropy may result in a false diagnosis of operator growth, or the lack thereof. In this paper, we examine the validity of coefficient entropy in the presence of hidden symmetries. Using the quantum cat map as an example, we show that under a generic choice of operator basis, the coefficient entropy fails to capture the suppression of operator growth caused by the symmetries. We further propose "symmetry-resolved coefficient entropy" as the proper diagnosis of operator complexity, which takes into account robust unknown symmetries, and demonstrate its effectiveness in the case of quantum cat map.


title: The Type II AGN-host galaxy connection: insights from the VVDS and VIPERS surveys

abstract: We present a study of optically-selected Type II AGN at 0.5 < z < 0.9 from the VIPERS and VVDS surveys, to investigate the connection between AGN activity and physical properties of their host galaxies. The host stellar mass is estimated through spectral energy distribution fitting with the CIGALE code, and star formation rates are derived from the [OII]$\lambda$3727 $\r{A}$ line luminosity. We find that 49% of the AGN host galaxies are on or above the main sequence (MS), 40% lie in the sub-MS locus, and 11% in the quiescent locus. Using the [OIII]$\lambda$5007 $\r{A}$ line luminosity as a proxy of the AGN power, we find that at fixed AGN power Type II AGN host galaxies show a bimodal behaviour: systems with host galaxy stellar mass <10$^{10}$ M$_{\odot}$, reside along the MS or in the starbursts locus (high-SF Type II AGN), while systems residing in massive host-galaxies (>10$^{10}$ M$_{\odot}$) show a lower level of star formation (low-SF Type II AGN). At all stellar masses, the offset from the MS is positively correlated with the AGN power. We interpret this correlation as evidence of co-evolution between the AGN and the host, possibly due to the availability of cold gas. In the most powerful AGN with host galaxies below the MS we find a hint, though weak, of asymmetry in the [OIII] line profile, likely due to outflowing gas, consistent with a scenario in which AGN feedback removes the available gas and halts the star formation in the most massive hosts.


title: Phonon-induced rotation of the electronic nematic director in superconducting Bi$_{2}$Se$_{3}$

abstract: The doped topological insulator $A_{x}\mathrm{Bi_{2}Se_{3}}$, with $A=\{\mathrm{Cu},\mathrm{Sr},\mathrm{Nb}\}$, becomes a nematic superconductor below $T_{c}\sim3-4\,\mathrm{K}$. The associated electronic nematic director is described by an angle $\alpha$ and is experimentally manifested in the elliptical shape of the in-plane critical magnetic field $H_{c2}$. Because of the threefold rotational symmetry of the lattice, $\alpha$ is expected to align with one of three high-symmetry directions corresponding to the in-plane nearest-neighbor bonds, consistent with a $Z_{3}$-Potts nematic transition. Here, we show that the nematic coupling to the acoustic phonons, which makes the nematic correlation length tend to diverge along certain directions only, can fundamentally alter this phenomenology in trigonal lattices. Compared to hexagonal lattices, the former possesses a sixth independent elastic constant $c_{14}$ due to the fact that the in-plane shear strain doublet $(\epsilon_{xx}-\epsilon_{yy},-2\epsilon_{xy})$ and the out-of-plane shear strain doublet $(2\epsilon_{yz},-2\epsilon_{xz})$ transform as the same irreducible representation. We find that, when $c_{14}$ overcomes a threshold value, which is expected to be the case in doped $\mathrm{Bi_{2}Se_{3}}$, the nematic director $\alpha$ unlocks from the high-symmetry directions due to the competition between the quadratic phonon-mediated interaction and the cubic nematic anharmonicity. This implies the breaking of the residual in-plane twofold rotational symmetry ($C_{2x}$), resulting in a triclinic phase. We discuss the implications of these findings to the structure of nematic domains and to the shape of the in-plane $H_{c2}$ in $A_{x}\mathrm{Bi_{2}Se_{3}}$, and to presence of nodes inside the superconducting state.


title: Neutrino meets ultralight dark matter: $\boldsymbol{0\nu\beta\beta}$ decay and cosmology

abstract: We explore the neutrinoless double beta ($0\nu \beta\beta$) decay induced by an ultralight dark matter field coupled to neutrinos. The effect on $0\nu\beta\beta$ decay is significant if the coupling violates the lepton number, for which the $\Delta L=2$ transition is directly driven by the dark matter field without further suppression of small neutrino masses. As the ultralight dark matter can be well described by a classical field, the effect features a periodic modulation pattern in decay events. However, we find that in the early Universe such coupling will be very likely to alter the standard cosmological results. In particular, the requirement of neutrino free-streaming before the matter-radiation equality severely constrains the parameter space, such that the future $0\nu \beta\beta$ decay experiments can hardly see any signal even with a meV sensitivity to the effective neutrino mass.


title: Learning Provably Robust Motion Planners Using Funnel Libraries

abstract: This paper presents an approach for learning motion planners that are accompanied with probabilistic guarantees of success on new environments that hold uniformly for any disturbance to the robot's dynamics within an admissible set. We achieve this by bringing together tools from generalization theory and robust control. First, we curate a library of motion primitives where the robustness of each primitive is characterized by an over-approximation of the forward reachable set, i.e., a "funnel". Then, we optimize probably approximately correct (PAC)-Bayes generalization bounds for training our planner to compose these primitives such that the entire funnels respect the problem specification. We demonstrate the ability of our approach to provide strong guarantees on two simulated examples: (i) navigation of an autonomous vehicle under external disturbances on a five-lane highway with multiple vehicles, and (ii) navigation of a drone across an obstacle field in the presence of wind disturbances.


title: Formal Synthesis of Controllers for Uncertain Linear Systems against $\omega$-Regular Properties: A Set-based Approach

abstract: In this paper, we present how to synthesize controllers to enforce $\omega$-regular properties over linear control systems affected by bounded disturbances. In particular, these controllers are synthesized based on so-called hybrid controlled invariant (HCI) sets. To compute these sets, we first construct a product system between the linear control system and the deterministic Streett automata (DSA) modeling the desired property. Then, we compute the maximal HCI set over the state set of the product system by leveraging a set-based approach. To ensure termination of the computation of the HCI sets within a finite number of iterations, we also propose two iterative schemes to compute approximations of the maximal HCI set. Finally, we show the effectiveness of our results via two case studies.


title: Revealing the ultra-sensitive calorimetric properties of supercon-ducting magic-angle twisted bilayer graphene

abstract: The allegedly unconventional superconducting phase of magic-angle twisted bilayer graphene (MATBG)1 has been predicted to possess extraordinary thermal properties, as it is formed from a highly diluted electron ensemble with both a record-low carrier density n ~ 10^11 cm-2 and electronic heat capacity Ce < 100 kB. While these attributes position MATBG as a ground-breaking material platform for revolutionary calorimetric applications2, these properties have so far not been experimentally shown. Here we reveal the ultra-sensitive calorimetric properties of a superconducting MATBG device, by monitoring its temperature dependent critical current Ic under continuous laser heating with a wavelength of 1550nm. From the bolometric effect, we are able to extract the temperature dependence of the electronic thermal conductance Gth, which remarkably has a non-zero value Gth = 0.19 pW/K at 35mK and in the low temperature limit is consistent with a power law dependence, as expected for nodal superconductors. Photo-voltage measurements on this non-optimized device reveal a peak responsivity of S = 5.8 x 10^7 V/W when the device is biased close to Ic, with a noise-equivalent power of NEP = 5.5 x 10^-16 WHz^-1/2. Analysis of the intrinsic perfor-mance shows that a theoretically achievable limit is defined by thermal fluctuations and can be as low as NEPTEF < 10^-20 WHz-1/2, with operation speeds as fast as ~ 500 ns. This establishes superconducting MATBG as a revolutionizing active material for ultra-sensitive photon-detection applications, which could enable currently unavailable technologies such as THz photon-number-resolving single-photon-detectors.


title: Ocean Mover's Distance: Using Optimal Transport for Analyzing Oceanographic Data

abstract: Remote sensing observations from satellites and global biogeochemical models have combined to revolutionize the study of ocean biogeochemical cycling, but comparing the two data streams to each other and across time remains challenging due to the strong spatial-temporal structuring of the ocean. Here, we show that the Wasserstein distance provides a powerful metric for harnessing these structured datasets for better marine ecosystem and climate predictions. Wasserstein distance complements commonly used point-wise difference methods such as the root mean squared error, by quantifying differences in terms of spatial displacement in addition to magnitude. As a test case we consider Chlorophyll (a key indicator of phytoplankton biomass) in the North-East Pacific Ocean, obtained from model simulations, in situ measurements, and satellite observations. We focus on two main applications: 1) Comparing model predictions with satellite observations, and 2) temporal evolution of Chlorophyll both seasonally and over longer time frames. Wasserstein distance successfully isolates temporal and depth variability and quantifies shifts in biogeochemical province boundaries. It also exposes relevant temporal trends in satellite Chlorophyll consistent with climate change predictions. Our study shows that optimal transport vectors underlying Wasserstein distance provide a novel visualization tool for testing models and better understanding temporal dynamics in the ocean.


title: Stellar proper motions in the outskirts of classical dwarf spheroidal galaxies with Gaia EDR3

abstract: We use Gaia EDR3 data to identify stars associated with six classical dwarf spheroidals (Draco, Ursa Minor, Sextans, Sculptor, Fornax, Carina) at their outermost radii, beyond their nominal King stellar limiting radius. For all of the dSphs examined, we find radial velocity matches with stars residing beyond the King limiting radius and with $> 50\%$ astrometric probability (four in Draco, two in Ursa Minor, eight in Sextans, two in Sculptor, twelve in Fornax, and five in Carina), indicating that these stars are associated with their respective dwarf spheroidals (dSphs) at high probability. We compare the positions of our candidate "extra-tidal" stars with the orbital tracks of the galaxies, and identify stars, both with and without radial velocity matches, that are consistent with lying along the orbital track of the satellites. However, given the small number of candidate stars, we cannot make any conclusive statements about the significance of these spatially correlated stars. Cross matching with publicly available catalogs of RR Lyrae, we find one RR Lyrae candidate with $> 50\%$ astrometric probability outside the limiting radius in each of Sculptor and Fornax, two such candidates in Draco, nine in Ursa Minor, seven in Sextans, and zero in Carina. Follow-up spectra on all of our candidates, including possible metallicity information, will help confirm association with their respective dSphs, and could represent evidence for extended stellar halos or tidal debris around these classical dSphs.


title: Synthesis-Guided Feature Learning for Cross-Spectral Periocular Recognition

abstract: A common yet challenging scenario in periocular biometrics is cross-spectral matching - in particular, the matching of visible wavelength against near-infrared (NIR) periocular images. We propose a novel approach to cross-spectral periocular verification that primarily focuses on learning a mapping from visible and NIR periocular images to a shared latent representational subspace, and supports this effort by simultaneously learning intra-spectral image reconstruction. We show the auxiliary image reconstruction task (and in particular the reconstruction of high-level, semantic features) results in learning a more discriminative, domain-invariant subspace compared to the baseline while incurring no additional computational or memory costs at test-time. The proposed Coupled Conditional Generative Adversarial Network (CoGAN) architecture uses paired generator networks (one operating on visible images and the other on NIR) composed of U-Nets with ResNet-18 encoders trained for feature learning via contrastive loss and for intra-spectral image reconstruction with adversarial, pixel-based, and perceptual reconstruction losses. Moreover, the proposed CoGAN model beats the current state-of-art (SotA) in cross-spectral periocular recognition. On the Hong Kong PolyU benchmark dataset, we achieve 98.65% AUC and 5.14% EER compared to the SotA EER of 8.02%. On the Cross-Eyed dataset, we achieve 99.31% AUC and 3.99% EER versus SotA EER of 4.39%.


title: GAP Enhancing Semantic Interoperability of Genomic Datasets and Provenance Through Nanopublications

abstract: While the publication of datasets in scientific repositories has become broadly recognised, the repositories tend to have increasing semantic-related problems. For instance, they present various data reuse obstacles for machine-actionable processes, especially in biological repositories, hampering the reproducibility of scientific experiments. An example of these shortcomings is the GenBank database. We propose GAP, an innovative data model to enhance the semantic data meaning to address these issues. The model focuses on converging related approaches like data provenance, semantic interoperability, FAIR principles, and nanopublications. Our experiments include a prototype to scrape genomic data and trace them to nanopublications as a proof of concept. For this, (meta)data are stored in a three-level nanopub data model. The first level is related to a target organism, specifying data in terms of biological taxonomy. The second level focuses on the biological strains of the target, the central part of our contribution. The strains express information related to deciphered (meta)data of the genetic variations of the genomic material. The third level stores related scientific papers (meta)data. We expect it will offer higher data storage flexibility and more extensive interoperability with other data sources by incorporating and adopting associated approaches to store genomic data in the proposed model.


title: Interpreting Angle Dependent Magnetoresistance in Layered Materials: Application to Cuprates

abstract: The evolution of the low temperature electronic structure of the cuprate metals from the overdoped to the underdoped side has recently been addressed through Angle-Dependant Magneto-Resistance (ADMR) experiments in La$_{1.6-x}$Nd$_{0.4}$Sr$_x$CuO$_4$. The results show a striking difference between hole dopings $p = 0.24$ and $p = 0.21$ which lie on either side of a putative quantum critical point at intermediate $p$. Motivated by this, we here study the theory of ADMR in correlated layered materials, paying special attention to the role of angle dependent quasiparticle weights $Z_{\mathbf{k}}$. Such a $Z_{\mathbf{k}}$ is expected to characterize a number of popular models of the cuprate materials, particularly when underdoped. Further, in the limit of weak interlayer hopping the quasiparticle weight will affect the $c$-axis transport measured in ADMR experiments. We show that proper inclusion of the quasiparticle weight does not support an interpretation of the data in terms of a $(\pi, \pi)$ spin density wave ordered state, in agreement with the lack of direct evidence for such order. We show that a simple model of Fermi surface reconfiguring across a van Hove point captures many of the striking differences seen between $p = 0.21$ and $p = 0.24$. We comment on why such a model may be appropriate for interpreting the ADMR data, despite having a large Fermi surface at $p = 0.21$, seemingly in contradiction with other evidence for a small Fermi surafce at that doping level.


title: Practical Guidance on Modeling Choices for the Virtual Twins Method

abstract: Individuals can vary drastically in their response to the same treatment, and this heterogeneity has driven the push for more personalized medicine. Accurate and interpretable methods to identify subgroups that respond to the treatment differently from the population average are necessary to achieving this goal. The Virtual Twins (VT) method by Foster et al. \cite{Foster} is a highly cited and implemented method for subgroup identification because of its intuitive framework. However, since its initial publication, many researchers still rely heavily on the authors' initial modeling suggestions without examining newer and more powerful alternatives. This leaves much of the potential of the method untapped. We comprehensively evaluate the performance of VT with different combinations of methods in each of its component steps, under a collection of linear and nonlinear problem settings. Our simulations show that the method choice for step 1 of VT is highly influential in the overall accuracy of the method, and Superlearner is a promising choice. We illustrate our findings by using VT to identify subgroups with heterogeneous treatment effects in a randomized, double-blind nicotine reduction trial.


title: Scaling of population resilience with dispersal length and habitat size

abstract: Environmental fluctuations can create population-depleted areas and even extinct areas for the population. This effect is more severe in the presence of the Allee effect (decreasing growth rate at low population densities). Dispersal inside the habitat provides a rescue effect on population-depleted areas, enhancing the population resilience to environmental fluctuations. Habitat reduction decreases the effectiveness of the dispersal rescue mechanism. We report here how the population resilience to environmental fluctuations decreases when the dispersal length or the habitat size are reduced. The resilience reduction is characterized by a decrease of the extinction threshold for environmental fluctuations. The extinction threshold is shown to scale with the ratio between the dispersal length and the scale of environmental synchrony, i.e., it is the dispersal connection between non-environmentally-correlated regions that provides resilience to environmental fluctuations. Habitat reduction also decreases the resilience to environmental fluctuations, when the habitat size is similar to or smaller than the characteristic dispersal distances. The power laws of these scaling behaviors are characterized here. Alternative scaling functions with spatial scales of population synchrony are found to fit the simulations worse. These results support the dispersal length as the critical scale for extinction induced by habitat reduction.


title: Flexible Bayesian Product Mixture Models for Vector Autoregressions

abstract: Bayesian non-parametric methods based on Dirichlet process mixtures have seen tremendous success in various domains and are appealing in being able to borrow information by clustering samples that share identical parameters. However, such methods can face hurdles in heterogeneous settings where objects are expected to cluster only along a subset of axes or where clusters of samples share only a subset of identical parameters. We overcome such limitations by developing a novel class of product of Dirichlet process location-scale mixtures that enable independent clustering at multiple scales, which result in varying levels of information sharing across samples. First, we develop the approach for independent multivariate data. Subsequently we generalize it to multivariate time-series data under the framework of multi-subject Vector Autoregressive (VAR) models that is our primary focus, which go beyond parametric single-subject VAR models. We establish posterior consistency and develop efficient posterior computation for implementation. Extensive numerical studies involving VAR models show distinct advantages over competing methods, in terms of estimation, clustering, and feature selection accuracy. Our resting state fMRI analysis from the Human Connectome Project reveals biologically interpretable connectivity differences between distinct intelligence groups, while another air pollution application illustrates the superior forecasting accuracy compared to alternate methods.


title: Magrathea-Pathfinder: A 3D adaptive-mesh code for geodesic ray tracing in $N$-body simulations

abstract: We introduce Magrathea-Pathfinder, a relativistic ray-tracing framework that can reconstruct the past light cone of observers in cosmological simulations. The code directly computes the 3D trajectory of light rays through the null geodesic equations, with the weak-field limit as its only approximation. This approach offers high levels of versatility while removing the need for many of the standard ray-tracing approximations such as plane-parallel, Born, or multiple-lens. Moreover, the use of adaptive integration steps and interpolation strategies based on adaptive-mesh refinement (AMR) grids allows Magrathea-Pathfinder to accurately account for the non-linear regime of structure formation and fully take advantage of the small-scale gravitational clustering. To handle very large N-body simulations, the framework has been designed as a high-performance computing post-processing tool relying on a hybrid parallelization that combines MPI tasks with C++11 std::threads. In this paper, we describe how realistic cosmological observables can be computed from numerical simulation using ray-tracing techniques. We discuss in particular the production of simulated catalogues and sky maps that account for all the observational effects considering first-order metric perturbations (such as peculiar velocities, gravitational potential, integrated Sachs-Wolfe, time-delay, and gravitational lensing). We perform convergence tests of our gravitational lensing algorithms and conduct performance benchmarks of the null geodesic integration procedures. Magrathea-Pathfinder introduces sophisticated ray-tracing tools to make the link between the space of N-body simulations and light-cone observables. This should provide new ways of exploring existing cosmological probes and building new ones beyond standard assumptions in order to prepare for the next generation of large-scale structure surveys.


title: The RayGalGroupSims cosmological simulation suite for the study of relativistic effects: an application to lensing-matter clustering statistics

abstract: General Relativistic effects on the clustering of matter in the universe provide a sensitive probe of cosmology and gravity theories that can be tested with the upcoming generation of galaxy surveys. Here, we present a suite of large volume high-resolution N-body simulations specifically designed to generate light-cone data for the study of relativistic effects on lensing-matter observables. RayGalGroupSims (or in short RayGal) consists of two N-body simulations of $(2625\,h^{-1}\,{\rm Mpc})^3$ volume with $4096^3$ particles of a standard flat $\Lambda$CDM model and a non-standard $w$CDM phantom dark energy model. Light-cone data from the simulations have been generated using a parallel ray-tracing algorithm that has accurately solved billion geodesic equations. Catalogues and maps with relativistic weak-lensing which include post-Born effects, magnification bias (MB) and redshift space distortions (RSD) due to gravitational redshift, Doppler, transverse Doppler, Integrated Sachs-Wolfe/Rees-Sciama effects, are publicly released. Using this dataset, we are able to reproduce the linear and quasi-linear predictions from the Class relativistic code for the 10 (cross-)power spectra (3$\times$2 points) of the matter density fluctuation field and the gravitational convergence at $z=0.7$ and $z=1.8$. We find $1-30\%$ level contribution from both MB and RSD to the matter power spectrum, while the Fingers-of-God effect is visible at lower redshift in the non-linear regime. MB contributes at the $10-30\%$ level to the convergence power spectrum leading to a deviation between the shear power-spectrum and the convergence power-spectrum. MB also plays a significant role in the galaxy-galaxy lensing by decreasing the density-convergence spectra by $20\%$, while coupling non-trivial configurations (such as the one with the convergence at the same or even lower redshift than the density field).


title: Adaptive Transmit Waveform Design

abstract: Recent research efforts in the Anti-Submarine Warfare (ASW) community have focused on developing sonar systems that adapt to their acoustic environment, referred to as "cognitive" sonars. Cognitive active sonar systems utilize principles of the perception action cycle of cognition to leverage information gathered from earlier sensing interactions with the underwater acoustic environment. This in turn informs the selection of system parameters to optimize target detection, classification, localization, and tracking performance in that acoustic environment. Of the many system parameters such a cognitive sonar system could potentially adapt, the acoustic signal transmitted into the medium, also known as the transmit waveform, has a profound impact on system performance. Many of the physical characteristics of the acoustic environment are contained in the return echo signal that is composed of amplitude scaled (target strength), time-delayed (target range) and Doppler shifted (target range-rate) echoes of the transmit waveform. This paper briefly describes a spectrally compact adaptive FM waveform model using Multi-Tone Sinusoidal Frequency Modulation (MTSFM). The MTSFM waveform's frequency and phase modulation functions are composed of a finite set of weighted sinusoidal harmonics. The weights for each harmonic are utilized as a discrete set of design coefficients. Adjusting these coefficients results in constant amplitude, spectrally compact FM waveforms with unique characteristics. The adaptability of the MTSFM combined with its transmitter friendly properties make it an attractive waveform type for a variety of active sonar applications and may provide a cognitive sonar system the ability to generate a complementary set of finely tuned waveforms for the novel scenarios and environments that it may encounter.


title: Unsupervised Learning of Symmetry Protected Topological Phase Transitions

abstract: Symmetry-protected topological (SPT) phases are short-range entangled phases of matter with a non-local order parameter which are preserved under a local symmetry group. Here, by using unsupervised learning algorithm, namely the diffusion maps, we demonstrate that can differentiate between symmetry broken phases and topologically ordered phases, and between non-trivial topological phases in different classes. In particular, we show that the phase transitions associated with these phases can be detected in different bosonic and fermionic models in one dimension. This includes the interacting SSH model, the AKLT model and its variants, and weakly interacting fermionic models. Our approach serves as an inexpensive computational method for detecting topological phases transitions associated with SPT systems which can be also applied to experimental data obtained from quantum simulators.


title: Kernel-based diffusion approximated Markov decision processes for autonomous navigation and control on unstructured terrains

abstract: We propose a diffusion approximation method to the continuous-state Markov Decision Processes (MDPs) that can be utilized to address autonomous navigation and control in unstructured off-road environments. In contrast to most decision-theoretic planning frameworks that assume fully known state transition models, we design a method that eliminates such a strong assumption that is often extremely difficult to engineer in reality. We first take the second-order Taylor expansion of the value function. The Bellman optimality equation is then approximated by a partial differential equation, which only relies on the first and second moments of the transition model. By combining the kernel representation of the value function, we design an efficient policy iteration algorithm whose policy evaluation step can be represented as a linear system of equations characterized by a finite set of supporting states. We first validate the proposed method through extensive simulations in 2D obstacle avoidance and 2.5D terrain navigation problems. The results show that the proposed approach leads to a much superior performance over several baselines. We then develop a system that integrates our decision-making framework with onboard perception and conduct real-world experiments in both cluttered indoor and unstructured outdoor environments. The results from the physical systems further demonstrate the applicability of our method in challenging real-world environments.


title: SMACE: A New Method for the Interpretability of Composite Decision Systems

abstract: Interpretability is a pressing issue for decision systems. Many post hoc methods have been proposed to explain the predictions of a single machine learning model. However, business processes and decision systems are rarely centered around a unique model. These systems combine multiple models that produce key predictions, and then apply decision rules to generate the final decision. To explain such decisions, we propose the Semi-Model-Agnostic Contextual Explainer (SMACE), a new interpretability method that combines a geometric approach for decision rules with existing interpretability methods for machine learning models to generate an intuitive feature ranking tailored to the end user. We show that established model-agnostic approaches produce poor results on tabular data in this setting, in particular giving the same importance to several features, whereas SMACE can rank them in a meaningful way.


title: Gravitational Wave Gastronomy

abstract: The symmetry breaking of grand unified gauge groups in the early Universe often leaves behind relic topological defects such as cosmic strings, domain walls, or monopoles. For some symmetry breaking chains, hybrid defects can form where cosmic strings attach to domain walls or monopoles attach to strings. In general, such hybrid defects are unstable, with one defect "eating" the other via the conversion of its rest mass into the other's kinetic energy and subsequently decaying via gravitational waves. In this work, we determine the gravitational wave spectrum from 1) the destruction of a cosmic string network by the nucleation of monopoles which cut up and "eat" the strings, 2) the collapse and decay of a monopole-string network by strings that "eat" the monopoles, 3) the destruction of a domain wall network by the nucleation of string-bounded holes on the wall that expand and "eat" the wall, and 4) the collapse and decay of a string-bounded wall network by walls that "eat" the strings. We call the gravitational wave signals produced from the "eating" of one topological defect by another gravitational wave gastronomy. We find that the four gravitational wave gastronomy signals considered yield unique spectra that can be used to narrow down the SO(10) symmetry breaking chain to the Standard Model and the scales of symmetry breaking associated with the consumed topological defects. Moreover, the systems we consider are unlikely to have a residual monopole or domain wall problem.


title: Lecture Notes on Berry Phases and Topology

abstract: In these notes, we review the role of Berry phases and topology in noninteracting electron systems. Topics including the adiabatic theorem, parallel transport, and Wannier functions are reviewed, with a focus on the connection to topological insulators.


title: The relativistic 5f electronic structure of delocalized $\alpha$-U and localized $\delta$-Pu from the self consistent vertex corrected GW approach and X-ray Emission Spectroscopy

abstract: The recently developed self-consistent vertex corrected GW method is used to calculate the 5f electronic structure in delocalized $\alpha$-U and localized $\delta$-Pu, each of which is confirmed by the historical experimental approaches of direct and inverse photoemission. Tender X-Ray Emission Spectroscopy (XES), in a novel application to 5f electronic structure, is used to experimentally prove the existence of 5f delocalization in $\alpha$-U.


title: Scintillation Timescales of Bright FRBs Detected by CHIME/FRB

abstract: We describe a pipeline to measure scintillation in fast radio bursts (FRBs) detected by CHIME/FRB in the 400-800 MHz band by analyzing the frequency structure of the FRB's spectrum. We use the pipeline to measure the characteristic frequency bandwidths of scintillation between $4-100$ kHz in 12 FRBs corresponding to timescales of $\sim$2-40 $\mu$s for 10 FRBs detected by CHIME/FRB. For the other two FRBs, we did not detect scintillation in the region our analysis is sensitive. We compared the measured scintillation timescales to the NE2001 predictions for the scintillation timescales from the Milky Way. We find a strong correlation to be an indication that in most instances, the observed scintillation of FRBs can be explained by the Milky Way.


title: $\mathrm{GL}_n$-structure and principal $\mathfrak{sl}_2$-triple on the cohomology ring of complex Grassmannian

abstract: In this note we describe the cohomology ring of the Grassmannian of $k$-planes in $n$-dimensional complex vector space as an $\mathrm{GL}_n$-module. We give explicit formulas for the operators of its principal $\mathfrak{sl}_2$-triple. It is proved that one of these operators corresponds to the shifted cohomology degree operator and the second operator coincides with the Lefschetz map on cohomology (as in the hard Lefschetz theorem). We check that the cohomology ring of the complex Grassmannian as a $\mathrm{GL}_n$-representation is isomorphic to the $k$-th exterior power of the standard $n$-dimensional representation.


title: Learning Scene Dynamics from Point Cloud Sequences

abstract: Understanding 3D scenes is a critical prerequisite for autonomous agents. Recently, LiDAR and other sensors have made large amounts of data available in the form of temporal sequences of point cloud frames. In this work, we propose a novel problem -- sequential scene flow estimation (SSFE) -- that aims to predict 3D scene flow for all pairs of point clouds in a given sequence. This is unlike the previously studied problem of scene flow estimation which focuses on two frames. We introduce the SPCM-Net architecture, which solves this problem by computing multi-scale spatiotemporal correlations between neighboring point clouds and then aggregating the correlation across time with an order-invariant recurrent unit. Our experimental evaluation confirms that recurrent processing of point cloud sequences results in significantly better SSFE compared to using only two frames. Additionally, we demonstrate that this approach can be effectively modified for sequential point cloud forecasting (SPF), a related problem that demands forecasting future point cloud frames. Our experimental results are evaluated using a new benchmark for both SSFE and SPF consisting of synthetic and real datasets. Previously, datasets for scene flow estimation have been limited to two frames. We provide non-trivial extensions to these datasets for multi-frame estimation and prediction. Due to the difficulty of obtaining ground truth motion for real-world datasets, we use self-supervised training and evaluation metrics. We believe that this benchmark will be pivotal to future research in this area. All code for benchmark and models will be made accessible.


title: Achieving Short-Blocklength RCU bound via CRC List Decoding of TCM with Probabilistic Shaping

abstract: This paper applies probabilistic amplitude shaping (PAS) to a cyclic redundancy check (CRC) aided trellis coded modulation (TCM) to achieve the short-blocklength random coding union (RCU) bound. In the transmitter, the equally likely message bits are first encoded by distribution matcher to generate amplitude symbols with the desired distribution. The binary representations of the distribution matcher outputs are then encoded by a CRC. Finally, the CRC-encoded bits are encoded and modulated by Ungerboeck's TCM scheme, which consists of a $\frac{k_0}{k_0+1}$ systematic tail-biting convolutional code and a mapping function that maps coded bits to channel signals with capacity-achieving distribution. This paper proves that, for the proposed transmitter, the CRC bits have uniform distribution and that the channel signals have symmetric distribution. In the receiver, the serial list Viterbi decoding (S-LVD) is used to estimate the information bits. Simulation results show that, for the proposed CRC-TCM-PAS system with 87 input bits and 65-67 8-AM coded output symbols, the decoding performance under additive white Gaussian noise channel achieves the RCU bound with properly designed CRC and convolutional codes.


title: Strategies and performance of the CMS silicon tracker alignment during LHC Run 2

abstract: The strategies for and the performance of the CMS silicon tracking system alignment during the 2015-2018 data-taking period of the LHC are described. The alignment procedures during and after data taking are explained. Alignment scenarios are also derived for use in the simulation of the detector response. Systematic effects, related to intrinsic symmetries of the alignment task or to external constraints, are discussed and illustrated for different scenarios.


title: Blandford-Znajek Process in Quadratic Gravity

abstract: The Blandford-Znajek process, which uses a magnetized plasma to extract energy from a rotating black hole, is one of the leading candidates for powering relativistic jets. In this work, we investigate the Blandford-Znajek process in two well-motivated quadratic gravity theories: scalar Gauss-Bonnet and dynamical Chern-Simons gravity. We solve analytically for a split-monopole magnetosphere to first order in the small-coupling approximation and second relative order in the slow-rotation approximation. The extracted power at fixed spin and magnetic flux is enhanced in scalar Gauss-Bonnet and reduced in dynamical Chern-Simons gravity, compared to general relativity. We find that there is a degeneracy between spin and the coupling constants of the theories at leading order in the slow rotation approximation that is broken at higher orders.


title: On the Potential of Execution Traces for Batch Processing Workload Optimization in Public Clouds

abstract: With the growing amount of data, data processing workloads and the management of their resource usage becomes increasingly important. Since managing a dedicated infrastructure is in many situations infeasible or uneconomical, users progressively execute their respective workloads in the cloud. As the configuration of workloads and resources is often challenging, various methods have been proposed that either quickly profile towards a good configuration or determine one based on data from previous runs. Still, performance data to train such methods is often lacking and must be costly collected. In this paper, we propose a collaborative approach for sharing anonymized workload execution traces among users, mining them for general patterns, and exploiting clusters of historical workloads for future optimizations. We evaluate our prototype implementation for mining workload execution graphs on a publicly available trace dataset and demonstrate the predictive value of workload clusters determined through traces only.


title: Field Synthesis with Azimuthally-Varying, Cascaded, Cylindrical Metasurfaces using a Wave Matrix Approach

abstract: In recent years, there has been extensive research on planar metasurfaces capable of arbitrarily controlling scattered fields. However, rigorous studies on conformal metasurfaces, such as those that are cylindrical, have been few in number likely due to their more complex geometry. Here, wave propagation in cascaded cylindrical structures consisting of layers of dielectric spacers and azimuthally-varying metasurfaces (subwavelength patterned metallic claddings) is investigated. A wave matrix approach, which incorporates the advantages of both ABCD matrices and scattering matrices (S matrices), is adopted. Wave matrices are used to model the higher order coupling between metasurface layers, overcoming fabrication difficulties associated with previous works. The proposed framework provides an efficient approach to synthesize the inhomogeneous sheet admittances that realize a desired cylindrical field transformation. Design examples are reported to illustrate the power and potential applications of the proposed method in antenna design and stealth technology.


title: Stronger Generalization Guarantees for Robot Learning by Combining Generative Models and Real-World Data

abstract: We are motivated by the problem of learning policies for robotic systems with rich sensory inputs (e.g., vision) in a manner that allows us to guarantee generalization to environments unseen during training. We provide a framework for providing such generalization guarantees by leveraging a finite dataset of real-world environments in combination with a (potentially inaccurate) generative model of environments. The key idea behind our approach is to utilize the generative model in order to implicitly specify a prior over policies. This prior is updated using the real-world dataset of environments by minimizing an upper bound on the expected cost across novel environments derived via Probably Approximately Correct (PAC)-Bayes generalization theory. We demonstrate our approach on two simulated systems with nonlinear/hybrid dynamics and rich sensing modalities: (i) quadrotor navigation with an onboard vision sensor, and (ii) grasping objects using a depth sensor. Comparisons with prior work demonstrate the ability of our approach to obtain stronger generalization guarantees by utilizing generative models. We also present hardware experiments for validating our bounds for the grasping task.


title: Propagating irreversibility fronts in cyclically-sheared suspensions

abstract: The interface separating a liquid from its vapor phase is diffuse: the composition varies continuously from one phase to the other over a finite length. Recent experiments on dynamic jamming fronts in two dimensions [Waitukaitis et al., Europhysics Letters 102, 44001 (2013)] identified a diffuse interface between jammed and unjammed discs. In both cases, the thickness of the interface diverges as a critical transition is approached. We investigate the generality of this behavior using a third system: a model of cyclically-sheared non-Brownian suspensions. As we sediment the particles towards a boundary, we observe a diffuse traveling front that marks the interface between irreversible and reversible phases. We argue that the front width is linked to a diverging correlation lengthscale in the bulk, which we probe by studying avalanches near criticality. Our results show how diffuse interfaces may arise generally when an incompressible phase is brought to a critical point.


title: Effect of decreasing population growth-rate on deforestation and population sustainability

abstract: We consider the effect of non-constant parameters on the human-forest interaction logistic model coupled with human technological growth introduced in "Deforestation and world population sustainability: a quantitative analysis"[1]. In recent years in fact, a decrease in human population growth rate has emerged which can be measured to about 1.7% drop per year since 1960 value which coincides with latest UN projections for next decades up to year 2100 [2]. We therefore consider here the effect of decreasing human population growth-rate on the aforementioned model and we evaluate its effect on the probability of survival of human civilisation without going through a catastrophic collapse in population. We find that for realistic values of the human population carrying capacity of the earth (measured by parameter beta) this decrease would not affect previous results leading to a low probability of avoiding a catastrophic collapse. For larger more optimistic values of beta instead, a decrease in growth-rate would tilt the probability in favour of a positive outcome, i.e. from 10-20% up to even 95% likelihood of avoiding collapse.


title: The HD 137496 system: A dense, hot super-Mercury and a cold Jupiter

abstract: Most of the currently known planets are small worlds with radii between that of the Earth and that of Neptune. The characterization of planets in this regime shows a large diversity in compositions and system architectures, with distributions hinting at a multitude of formation and evolution scenarios. Using photometry from the K2 satellite and radial velocities measured with the HARPS and CORALIE spectrographs, we searched for planets around the bright and slightly evolved Sun-like star HD 137496. We precisely estimated the stellar parameters, $M_*$ = 1.035 +/- 0.022 $M_\odot$, $R_*$ = 1.587 +/- 0.028 $R_\odot$, $T_\text{eff}$ = 5799 +/- 61 K, together with the chemical composition of the slightly evolved star. We detect two planets orbiting HD 137496. The inner planet, HD 137496 b, is a super-Mercury (an Earth-sized planet with the density of Mercury) with a mass of $M_b$ = 4.04 +/- 0.55 $M_\oplus$, a radius of $R_b = 1.31_{-0.05}^{+0.06} R_\oplus,$ and a density of $\rho_b = 10.49_{-1.82}^{+2.08}$ $\mathrm{g cm^{-3}}$. With an interior modeling analysis, we find that the planet is composed mainly of iron, with the core representing over 70% of the planet's mass ($M_{core}/M_{total} = 0.73^{+0.11}_{-0.12}$). The outer planet, HD 137496 c, is an eccentric ($e$ = 0.477 +/- 0.004), long period ($P$ = $479.9_{-1.1}^{+1.0}$ days) giant planet ($M_c\sin i_c$ = 7.66 +/- 0.11 $M_{Jup}$) for which we do not detect a transit. HD 137496 b is one of the few super-Mercuries detected to date. The accurate characterization reported here enhances its role as a key target to better understand the formation and evolution of planetary systems. The detection of an eccentric long period giant companion also reinforces the link between the presence of small transiting inner planets and long period gas giants.


title: On amphichirality of symmetric unions

abstract: It is still unknown whether there is a nontrivial knot with Jones polynomial equal to that of the unknot. Tanaka shows that if an amphichiral knot is a symmetric union of the unknot with one twist region, then its Jones polynomial is trivial. Hence, he proposes that if any of these knots were nontrivial, a nontrivial knot with trivial Jones polynomial would exist. We first show such a knot is always trivial and hence cannot be used to answer the above question. We then generalize the argument to symmetric unions of any knots and show that if a symmetric union of a knot $J$ with one twist region is amphichiral, then it is the connected sum of $J$ and its mirror image $-J$.


title: New results on the p-adic valuation of Stirling numbers

abstract: We generalize results on the $p$-adic valuations of $S(n,k)$, the Stirling number of the second kind and $s(n,k)$ the Stirling number of the first kind. We have several new estimates for these valuations, along with criteria for when the estimates are sharp. The primary foci are the explicit evaluation of $\nu_2(S(n,k))$ with $n=c2^h$, $k=b2^h+a$, $a, b, c, h, k, n \in Z^+$, and $1\le a \le 2^{h-1}$, and $\nu_p(S(n,k))$ when $n=cp^h$ for an odd prime $p$. We have strong new results, which generalize and strengthen previous results, for all primes. We also have some new results on the $p$-adic valuations $\nu_p(s(n,k))$ for all primes. We generally assume that $p-1|n-k$ for exact values of $\nu_p(S(n,k))$ or $\nu_p(s(n,k))$. In addition, we have proved some new Amdeberhan-type identities for Stirling numbers of both kinds. We also extend some recent results and propose two new conjectures, as well as proofs and extensions of previous ones.


title: Constraining general U(1) interactions from neutrino-electron scattering measurements at DUNE near detector

abstract: The neutrino-electron scattering process is a powerful tool to explore new physics beyond the standard model. Recently the possibility of DUNE Near Detector (ND) to constrain various new physics scenarios using this process have been highlighted in the literature. In this work, we consider the most general U(1) model and probe the constraints on the mass and coupling strength of the additional $Z'$ from $\nu-e$ scattering at DUNE ND. The presence of the $Z^\prime$ gives rise to extra interference effects. In the context of the general U(1) model, the destructive interference can occur in either neutrino or anti-neutrino channel or for both or none. This opens up the possibilities of getting four different type of signal in the neutrino and ant-neutrino runs of DUNE. We perform the analysis using both the total rate and binned events spectrum. Our results show that in a bin by bin analysis the effect of destructive interference is less compared to the analysis using total rate. We present the bounds on the $ m_{Z^\prime} - g_X$ plane from $\nu-e$ scattering measurements at DUNE ND and compare these with those obtained from other $\nu-e$ scattering, COHERENT, and beam dump experiments. We show that the DUNE ND can give the best bound for certain mass ranges of $Z^\prime$.


title: Ultrathin ferrimagnetic GdFeCo films with very low damping

abstract: Ferromagnetic materials dominate as the magnetically active element in spintronic devices, but come with drawbacks such as large stray fields, and low operational frequencies. Compensated ferrimagnets provide an alternative as they combine the ultrafast magnetization dynamics of antiferromagnets with a ferromagnet-like spin-orbit-torque (SOT) behavior. However to use ferrimagnets in spintronic devices their advantageous properties must be retained also in ultrathin films (t < 10 nm). In this study, ferrimagnetic Gdx(Fe87.5Co12.5)1-x thin films in the thickness range t = 2-20 nm were grown on high resistance Si(100) substrates and studied using broadband ferromagnetic resonance measurements at room temperature. By tuning their stoichiometry, a nearly compensated behavior is observed in 2 nm Gdx(Fe87.5Co12.5)1-x ultrathin films for the first time, with an effective magnetization of Meff = 0.02 T and a low effective Gilbert damping constant of {\alpha} = 0.0078, comparable to the lowest values reported so far in 30 nm films. These results show great promise for the development of ultrafast and energy efficient ferrimagnetic spintronic devices.


title: Individual Pulse Monitoring and Dose Control System for Pre-Clinical Implementation of FLASH-RT

abstract: Ultra-high dose rate electron sources require dose rate independent dosimeters and a calibrated dose control system for accurate delivery. In this study, we developed a single-pulse dose monitoring and a real-time dose-based control system for a converted clinical linear accelerator (LINAC). A point scintillator detector was coupled to a gated amplifier and a real-time controller for dose monitoring and feedback control loop. The controller was programmed to integrate dose and measure pulse width of each radiation pulse and gate the LINAC beam when the prescribed dose was delivered. The scintillator was mounted in solid water phantom and placed underneath mice skin for in vivo dose monitoring. Additionally, the scintillator was characterized in terms of its radiation stability, mean dose-rate, and dose per pulse dependence. Dose integration was performed for each radiation pulse and displayed in real-time. The scintillator was shown to be linear with mean dose-rate (40-380 Gy/s) and dose per pulse (0.3-1.3 Gy/Pulse) to within +/- 3%. However, the plastic scintillator was subject to significant radiation damage (16%/kGy) and would need to be calibrated frequently. Pulse-counting control was accurately implemented with direct correspondence between the intended and the actual delivered pulses. The dose-based control was sufficient to gate on any pulse of the LINAC. In-vivo dosimetry monitoring with a 1 cm circular cut-out revealed that a ramp-up of 4-5 pulses was present during which the average dose per pulse was ~0.045 +/- 0.004 Gy/Pulse, whereas after the ramp-up it stabilized at 0.65 +/- 0.01 Gy/Pulse. The tools presented in this study can be used to determine the beam parameter space pertinent to the FLASH effect. Additionally, this study is the first instance of real-time dose-based control for a modified LINAC at ultra-high dose rates.


title: The spheromak tilting and how it affects modelling coronal mass ejections

abstract: Spheromak type flux ropes are increasingly used for modelling coronal mass ejections (CMEs). Many models aim in accurately reconstructing the magnetic field topology of CMEs, considering its importance in assessing their impact on modern technology and human activities in space and on ground. However, so far there is little discussion about how the details of the magnetic structure of a spheromak affect its evolution through the ambient field in the modelling domain, and what impact this has on the accuracy of magnetic field topology predictions. If the spheromak has its axis of symmetry (geometric axis) at an angle with respect to the direction of the ambient field, then the spheromak starts rotating so that its symmetry axis finally aligns with the ambient field. When using the spheromak in space weather forecasting models this tilting can happen already during insertion and significantly affects the results. In this paper we highlight this issue previously not examined in the field of space weather and we estimate the angle by which the spheromak rotates under different conditions. To do this we generated simple purely radial ambient magnetic field topologies (weak/strong positive/negative) and inserted spheromaks with varying initial speed and tilt, and magnetic helicity sign. We employ different physical and geometric criteria to locate the magnetic centre of mass and axis of symmetry of the spheromak. We confirm that spheromaks rotate in all investigated conditions and their direction and angle of rotation depend on the spheromak's initial properties and ambient magnetic field strength and orientation.


title: Variational Adiabatic Gauge Transformation on real quantum hardware for effective low-energy Hamiltonians and accurate diagonalization

abstract: Effective low-energy theories represent powerful theoretical tools to reduce the complexity in modeling interacting quantum many-particle systems. However, common theoretical methods rely on perturbation theory, which limits their applicability to weak interactions. Here we introduce the Variational Adiabatic Gauge Transformation (VAGT), a non-perturbative hybrid quantum algorithm that can use nowadays quantum computers to learn the variational parameters of the unitary circuit that brings the Hamiltonian to either its block-diagonal or full-diagonal form. If a Hamiltonian can be diagonalized via a shallow quantum circuit, then VAGT can learn the optimal parameters using a polynomial number of runs. The accuracy of VAGT is tested trough numerical simulations, as well as simulations on Rigetti and IonQ quantum computers.


title: Computer Vision for Supporting Image Search

abstract: Computer vision and multimedia information processing have made extreme progress within the last decade and many tasks can be done with a level of accuracy as if done by humans, or better. This is because we leverage the benefits of huge amounts of data available for training, we have enormous computer processing available and we have seen the evolution of machine learning as a suite of techniques to process data and deliver accurate vision-based systems. What kind of applications do we use this processing for ? We use this in autonomous vehicle navigation or in security applications, searching CCTV for example, and in medical image analysis for healthcare diagnostics. One application which is not widespread is image or video search directly by users. In this paper we present the need for such image finding or re-finding by examining human memory and when it fails, thus motivating the need for a different approach to image search which is outlined, along with the requirements of computer vision to support it.


title: Disformal Transition of a Black Hole to a Wormhole in Scalar-Tensor Horndeski Theory

abstract: We consider disformal transformations in a subclass of Horndeski theory in which a scalar field is kinetically coupled to the Einstein tensor. We apply a disformal transformation on a seed hairy black hole solution of this theory and we show that there is a transition of a black hole to a wormhole. We also show that the null energy condition is violated in the wormhole configuration and we study the stability of the wormhole solution by calculating the time evolution of scalar perturbations in this geometry.


title: Film Trailer Generation via Task Decomposition

abstract: Movie trailers perform multiple functions: they introduce viewers to the story, convey the mood and artistic style of the film, and encourage audiences to see the movie. These diverse functions make automatic trailer generation a challenging endeavor. We decompose it into two subtasks: narrative structure identification and sentiment prediction. We model movies as graphs, where nodes are shots and edges denote semantic relations between them. We learn these relations using joint contrastive training which leverages privileged textual information (e.g., characters, actions, situations) from screenplays. An unsupervised algorithm then traverses the graph and generates trailers that human judges prefer to ones generated by competitive supervised approaches.


title: On two congruence conjectures of Z.-W. Sun involving Franel numbers

abstract: In this paper, we mainly prove the following conjectures of Z.-W. Sun \cite{S13}: Let $p>2$ be a prime. If $p=x^2+3y^2$ with $x,y\in\mathbb{Z}$ and $x\equiv1\pmod 3$, then $$x\equiv\frac14\sum_{k=0}^{p-1}(3k+4)\frac{f_k} {2^k}\equiv\frac12\sum_{k=0}^{p-1}(3k+2)\frac{f_k}{(-4)^k}\pmod{p^2},$$ and if $p\equiv1\pmod3$, then $$\sum_{k=0}^{p-1}\frac{f_k}{2^k}\equiv\sum_{k=0}^{p-1}\frac{f_k}{(-4)^k}\pmod{p^3},$$ where $f_n=\sum_{k=0}^n\binom{n}k^3$ stands for the $n$th Franel number.


title: Models of the in situ formation of detected extrasolar giant planets

abstract: (Abridged) We present numerical simulations of the formation of the planetary companions to 47 UMa, rho CrB, and 51 Peg. They are assumed to have formed in situ according to the basic model that a core formed first by accretion of solid particles, then later it captured substantial amounts of gas from the protoplanetary disk. In most of the calculations we prescribe a constant accretion rate for the solid core. The evolution of the gaseous envelope assumes that: (1) it is in quasi-hydrostatic equilibrium, (2) the gas accretion rate is determined by the requirement that the outer radius of the planet is the place at which the thermal velocity of the gas allows it to reach the boundary of the planet's Hill sphere, (3) the gas accretion rate is limited, moreover, by the prescribed maximum rate at which the nebula can supply the gas, and (4) the growth of the planet stops once it obtains approximately the minimum mass determined from radial velocity measurements. Calculations are carried out through an initial phase during which solid accretion dominates, past the point of crossover when the masses of solid and gaseous material are equal, through the phase of rapid gas accretion, and into the final phase of contraction and cooling at constant mass. Alternative calculations are presented for the case of 47 UMa in which the solid accretion rate is calculated, not assumed, and the dissolution of planetesimals within the gaseous envelope is considered. In all cases there is a short phase of high luminosity (1e-3-1e-2 Lsun) associated with rapid gas accretion. The height and duration of this peak depend on uncertain model parameters. The conclusion is reached that in situ formation of all of these companions is possible under some conditions. However, it is more likely that orbital migration was an important component of the evolution, at least for the planets around rho CrB and 51 Peg.


title: Sharp Bounds on Eigenvalues via Spectral Embedding Based on Signless Laplacians

abstract: Using spectral embedding based on the signless Laplacian, we obtain bounds on the spectrum of transition matrices on graphs. As a consequence, we bound return probabilities and the uniform mixing time of simple random walk on graphs. In addition, spectral embedding is used in this article to bound the spectrum of graph adjacency matrices. Our method is adapted from [Lyons and Oveis Gharan, 2017].


title: Proof of some congruence conjectures of Z.-H. Sun involving Ap\'{e}ry-like numbers

abstract: In this paper, we mainly prove the following conjecture of Z.-H. Sun cite{SH20}: Let $p>3$ be a prime. Then $$\sum_{k=0}^{p-1}\binom{2k}k\frac{3k+1}{(-16)^k}f_k\equiv(-1)^{(p-1)/2}p+p^3E_{p-3}\pmod{p^4},$$ where $f_n=\sum_{k=0}^n\binom{n}k^3$ and $E_n$ stand for the $n$th Franel number and $n$th Euler number respectively.


title: Planets or asteroids? A geochemical method to constrain the masses of White Dwarf pollutants

abstract: Polluted white dwarfs that have accreted planetary material provide a unique opportunity to probe the geology of exoplanetary systems. However, the nature of the bodies which pollute white dwarfs is not well understood: are they small asteroids, minor planets, or even terrestrial planets? We present a novel method to infer pollutant masses from detections of Ni, Cr and Si. During core--mantle differentiation, these elements exhibit variable preference for metal and silicate at different pressures (i.e., object masses), affecting their abundances in the core and mantle. We model core--mantle differentiation self-consistently using data from metal--silicate partitioning experiments. We place statistical constraints on the differentiation pressures, and hence masses, of bodies which pollute white dwarfs by incorporating this calculation into a Bayesian framework. We show that Ni observations are best suited to constraining pressure when pollution is mantle-like, while Cr and Si are better for core-like pollution. We find 3 systems (WD0449-259, WD1350-162 and WD2105-820) whose abundances are best explained by the accretion of fragments of small parent bodies ($<0.2M_\oplus$). For 2 systems (GD61 and WD0446-255), the best model suggests the accretion of fragments of Earth-sized bodies, although the observed abundances remain consistent ($<3\sigma$) with the accretion of undifferentiated material. This suggests that polluted white dwarfs potentially accrete planetary bodies of a range of masses. However, our results are subject to inevitable degeneracies and limitations given current data. To constrain pressure more confidently, we require serendipitous observation of (nearly) pure core and/or mantle material.


title: Optimal Oblivious Reconfigurable Networks

abstract: Oblivious routing has a long history in both the theory and practice of networking. In this work we initiate the formal study of oblivious routing in the context of reconfigurable networks, a new architecture that has recently come to the fore in datacenter networking. These networks allow a rapidly changing bounded-degree pattern of interconnections between nodes, but the network topology and the selection of routing paths must both be oblivious to the traffic demand matrix. Our focus is on the trade-off between maximizing throughput and minimizing latency in these networks. For every constant throughput rate, we characterize (up to a constant factor) the minimum latency achievable by an oblivious reconfigurable network design that satisfies the given throughput guarantee. The trade-off between these two objectives turns out to be surprisingly subtle: the curve depicting it has an unexpected scalloped shape reflecting the fact that load-balancing becomes more difficult when the average length of routing paths is not an integer because equalizing all the path lengths is not possible. The proof of our lower bound uses LP duality to verify that Valiant load balancing is the most efficient oblivious routing scheme when used in combination with an optimally-designed reconfigurable network topology. The proof of our upper bound uses an algebraic construction in which the network nodes are identified with vectors over a finite field, the network topology is described by either the elementary basis or a sequence of Vandermonde matrices, and routing paths are constructed by selecting columns of these matrices to yield the appropriate mixture of path lengths within the shortest possible time interval.