title: High efficacy of layered controls for reducing transmission of airborne pathogens

abstract: To optimize strategies for curbing the transmission of airborne pathogens, the efficacy of three key controls -- face masks, ventilation, and physical distancing -- must be well understood. In this study we used the Quadrature-based model of Respiratory Aerosol and Droplets to quantify the reduction in exposure to airborne pathogens from various combinations of controls. For each combination of controls, we simulated thousands of scenarios that represent the tremendous variability in factors governing airborne transmission and the efficacy of mitigation strategies. While the efficacy of any individual control was highly variable among scenarios, combining universal mask-wearing with distancing of 1~m or more reduced the median exposure by more than 99\% relative to a close, unmasked conversation, with further reductions if ventilation is also enhanced. The large reductions in exposure to airborne pathogens translated to large reductions in the risk of initial infection in a new host. These findings suggest that layering controls is highly effective for reducing transmission of airborne pathogens and will be critical for curbing outbreaks of novel viruses in the future.


title: First integral cohomology group of the pure mapping class group of a non-orientable surface of infinite type

abstract: In this work we compute the first integral cohomology of the pure mapping class group of a non-orientable surface of infinite topological type and genus at least 3. To this purpose, we also prove several other results already known for orientable surfaces such as the existence of an Alexander method, the fact that the mapping class group is isomorphic to the automorphism group of the curve graph along with the topological rigidity of the curve graph, and the structure of the pure mapping class group as both a Polish group and a semi-direct product.


title: Find the Gap: Black Hole Population Analysis with an Astrophysically Motivated Mass Function

abstract: We introduce a novel black hole mass function which realistically models the physics of star formation and pair instability supernova with a minimal number of parameters. Applying this to all events in the LIGO-Virgo GWTC-2 catalog, we detect a peak at M_BHMG = 74.8^{+4.3}_{-8.0} MS, followed by a break in the mass function. Repeating the analysis without the black holes from the event GW190521, we find this feature at M_BHMG = 55.4^{+3.0}_{-6.1} MS. The latter result establishes the edge of the anticipated "black hole mass gap" at a value compatible with the expectation from standard stellar structure theory, while the former result is ~ 20MS higher, which would have far-reaching implications if confirmed. Using Bayesian techniques, we establish that our mass function fits a new catalog of black hole masses approximately as well as the pre-existing phenomenological mass functions. We also remark on the implications of these results for constraining or discovering new phenomena in nuclear and particle physics.

categories: astro-ph.CO astro-ph.IM astro-ph.SR gr-qc hep-ph


title: The GALAH Survey and Symbiotic Stars. I. Discovery and follow-up of 33 candidate accreting-only systems

abstract: We have identified a first group of 33 new candidates for symbiotic stars (SySt) of the accreting-only variety among the 600,255 stars so far observed by the GALAH high-resolution spectroscopic survey of the Southern Hemisphere, more than doubling the number of those previously known. GALAH aims to high latitudes and this offers the possibility to sound the Galaxy for new SySt away from the usual Plane and Bulge hunting regions. In this paper we focus on SySt of the M spectral type, showing an Halpha emission with a peak in excess of 0.5 above the adjacent continuum level, and not affected by coherent radial pulsations. These constraints will be relaxed in future studies. The 33 new candidate SySt were subjected to a vast array of follow-up confirmatory observations (X-ray/UV observations with the Swift satellite, search for optical flickering, presence of a near-UV upturn in ground-based photometric and spectroscopic data, radial velocity changes suggestive of orbital motion, variability of the emission line profiles). According to Gaia eDR3 parallaxes, the new SySt are located at the tip of the Giant Branch, sharing the same distribution in M(Ks) of the well established SySt. The accretion luminosities of the new SySt are in the range 1-10 Lsun, corresponding to mass-accretion rates of 0.1-1x10(-9) Msun/yr for WDs of 1 Msun. The M giant of one of the new SySt presents a large Lithium over-abundance.


title: Strumming to the Beat: Audio-Conditioned Contrastive Video Textures

abstract: We introduce a non-parametric approach for infinite video texture synthesis using a representation learned via contrastive learning. We take inspiration from Video Textures, which showed that plausible new videos could be generated from a single one by stitching its frames together in a novel yet consistent order. This classic work, however, was constrained by its use of hand-designed distance metrics, limiting its use to simple, repetitive videos. We draw on recent techniques from self-supervised learning to learn this distance metric, allowing us to compare frames in a manner that scales to more challenging dynamics, and to condition on other data, such as audio. We learn representations for video frames and frame-to-frame transition probabilities by fitting a video-specific model trained using contrastive learning. To synthesize a texture, we randomly sample frames with high transition probabilities to generate diverse temporally smooth videos with novel sequences and transitions. The model naturally extends to an audio-conditioned setting without requiring any finetuning. Our model outperforms baselines on human perceptual scores, can handle a diverse range of input videos, and can combine semantic and audio-visual cues in order to synthesize videos that synchronize well with an audio signal.


title: Pricing without no-arbitrage condition in discrete time

abstract: In a discrete time setting, we study the central problem of giving a fair price to some financial product. For several decades, the no-arbitrage conditions and the martingale measures have played a major role for solving this problem. We propose a new approach for estimating the super-replication cost based on convex duality instead of martingale measures duality: The prices are expressed using Fenchel conjugate and bi-conjugate without using any no-arbitrage condition.The super-hedging problem resolution leads endogenously to a weak no-arbitrage condition called Absence of Instantaneous Profit (AIP) under which prices are finite. We study this condition in details, propose several characterizations and compare it to the no-arbitrage condition.


title: A Student-Teacher Architecture for Dialog Domain Adaptation under the Meta-Learning Setting

abstract: Numerous new dialog domains are being created every day while collecting data for these domains is extremely costly since it involves human interactions. Therefore, it is essential to develop algorithms that can adapt to different domains efficiently when building data-driven dialog models. The most recent researches on domain adaption focus on giving the model a better initialization, rather than optimizing the adaptation process. We propose an efficient domain adaptive task-oriented dialog system model, which incorporates a meta-teacher model to emphasize the different impacts between generated tokens with respect to the context. We first train our base dialog model and meta-teacher model adversarially in a meta-learning setting on rich-resource domains. The meta-teacher learns to quantify the importance of tokens under different contexts across different domains. During adaptation, the meta-teacher guides the dialog model to focus on important tokens in order to achieve better adaptation efficiency. We evaluate our model on two multi-domain datasets, MultiWOZ and Google Schema-Guided Dialogue, and achieve state-of-the-art performance.


title: Constitutive modeling of dilute wormlike micelle solutions: shear-induced structure and transient dynamics

abstract: We present a reformulation of the 'reactive rod model' (RRM) of Dutta and Graham [Dutta, Sarit and Graham, Michael D., JNNFM 251 (2018)], a constitutive model for describing the behavior of dilute wormlike micelle solutions. The RRM treats wormlike micelle solutions as dilute suspensions of rigid Brownian rods undergoing reversible scission and growth in flow. Evolution equations for micelle orientation and stress contribution are coupled to a kinetic reaction equation for a collective micelle length, producing dynamic variations in the length and rotational diffusivity of the rods. This model has previously shown success in capturing many critical steady-state rheological features of dilute wormlike micelle solutions, particularly shear-thickening and -thinning, non-zero normal stress differences, and a reentrant shear stress-shear rate curve, and could fit a variety of steady state experimental data. The present work improves on this framework, which showed difficulty in capturing transient dynamics and high-shear behavior, by reformulating the kinetic equation for micelle growth on a more microstructural (though still highly idealized) basis. In particular, we allow for micelle growth associated with strong alignment of rods and breakage due to tensile stresses along the micelles. This new formulation captures both steady and transient shear rheology in good agreement with experiments. We also find good agreement with available steady state extensional rheology.


title: Localizing Visual Sounds the Hard Way

abstract: The objective of this work is to localize sound sources that are visible in a video without using manual annotations. Our key technical contribution is to show that, by training the network to explicitly discriminate challenging image fragments, even for images that do contain the object emitting the sound, we can significantly boost the localization performance. We do so elegantly by introducing a mechanism to mine hard samples and add them to a contrastive learning formulation automatically. We show that our algorithm achieves state-of-the-art performance on the popular Flickr SoundNet dataset. Furthermore, we introduce the VGG-Sound Source (VGG-SS) benchmark, a new set of annotations for the recently-introduced VGG-Sound dataset, where the sound sources visible in each video clip are explicitly marked with bounding box annotations. This dataset is 20 times larger than analogous existing ones, contains 5K videos spanning over 200 categories, and, differently from Flickr SoundNet, is video-based. On VGG-SS, we also show that our algorithm achieves state-of-the-art performance against several baselines.


title: Exact Limit Theorems for Restricted Integer Partitions

abstract: For a set of positive integers $A$, let $p_A(n)$ denote the number of ways to write $n$ as a sum of integers from $A$, and let $p(n)$ denote the usual partition function. In the early 40s, Erd\H{o}s extended the classical Hardy--Ramanujan formula for $p(n)$ by showing that $A$ has density $\alpha$ if and only if $\log p_A(n) \sim \log p(\alpha n)$. Nathanson asked if Erd\H{o}s's theorem holds also with respect to $A$'s lower density, namely, whether $A$ has lower-density $\alpha$ if and only if $\log p_A(n) / \log p(\alpha n)$ has lower limit $1$. We answer this question negatively by constructing, for every $\alpha > 0$, a set of integers $A$ of lower density $\alpha$, satisfying $$ \liminf_{n \rightarrow \infty} \frac{\log p_A(n)}{\log p(\alpha n)} \geq \left(\frac{\sqrt{6}}{\pi}-o_{\alpha}(1)\right)\log(1/\alpha)\;. $$ We further show that the above bound is best possible (up to the $o_\alpha(1)$ term), thus determining the exact extremal relation between the lower density of a set of integers and the lower limit of its partition function. We also prove an analogous theorem with respect to the upper density of a set of integers, answering another question of Nathanson.


title: Spin/Parity Dependent Level Density

abstract: It is shown that the Constant Temperature (CT) model of nuclear level density is a direct consequence of a symmetrized Poisson distribution of nuclear level spacings. The standard CT model describing the total level density is shown to be fatally flawed due to discontinuities at the Yrast energies, the onset of new $J^{\pi}$ sequences, that disrupt the exponential formula and cause the back shift parameter to become nonphysically negative. A new CT-JPI level density model is proposed with a constant temperature and separate back shift parameters for each $J^{\pi}$ sequence. The CT-JPI model is also constrained to reproduce the spin distribution predicted by Ericson's spin distribution function at the neutron separation energy. A fitting procedure is described for determining the temperature $T$, back shifts $E_0(J^{\pi})$, and spin cutoff parameters $\sigma_c$ from nuclear structure and resonance data. The CT-JPI model is demonstrated to successfully predict the level densities for a wide range of spins and parities for 46 nuclear with Z=7-92. In variance with earlier predictions the spin cut-off parameters show no mass dependence and instead substantial variation at all mass regions.


title: Merton Investment Problems in Finance and Insurance for the Hawkes-based Models

abstract: We show how to solve Merton optimal investment stochastic control problem for Hawkes-based models in finance and insurance, i.e., for a wealth portfolio X(t) consisting of a bond and a stock price described by general compound Hawkes process (GCHP), and for a capital R(t) of an insurance company with the amount of claims described by the risk model based on GCHP. The novelty of the results consists of the new Hawkes-based models and in the new optimal investment results in finance and insurance for those models.


title: Compactivorous Sets in Banach Spaces

abstract: A set $E$ in a Banach space $X$ is compactivorous if for every compact set $K$ in $X$ there is a nonempty, (relatively) open subset of $K$ which can be translated into $E$. In a separable Banach space, this is a sufficient condition which guarantees the Haar nonnegligibility of Borel subsets. We give some characterisations of this property in both separable and nonseparable Banach spaces and prove an extension of the main theorem to countable products of locally compact Polish groups.


title: Dynamic-Aware Autonomous Exploration in Populated Environments

abstract: Autonomous exploration allows mobile robots to navigate in initially unknown territories in order to build complete representations of the environments. In many real-life applications, environments often contain dynamic obstacles which can compromise the exploration process by temporarily blocking passages, narrow paths, exits or entrances to other areas yet to be explored. In this work, we formulate a novel exploration strategy capable of explicitly handling dynamic obstacles, thus leading to complete and reliable exploration outcomes in populated environments. We introduce the concept of dynamic frontiers to represent unknown regions at the boundaries with dynamic obstacles together with a cost function which allows the robot to make informed decisions about when to revisit such frontiers. We evaluate the proposed strategy in challenging simulated environments and show that it outperforms a state-of-the-art baseline in these populated scenarios.


title: Thermodynamics of vacancies in concentrated solid solutions: From dilute Ni-alloys to the Cantor system

abstract: The vacancy concentration at finite temperatures is studied for a series of (CoCrFeMn)$_{1-x_\mathrm{Ni}}$Ni$_{x_\mathrm{Ni}}$ alloys by grand-canonical Monte-Carlo (MC) simulations. The vacancy formation energies are calculated from a classical interatomic potential and exhibit a distribution due to the different chemical environments of the vacated sites. In dilute alloys, this distribution features multiple discrete peaks, while concentrated alloys exhibit an unimodal distribution as there are many different chemical environments of similar vacancy formation energy. MC simulations using a numerically efficient bond-counting model confirm that the vacancy concentration even in concentrated alloys may be calculated by the established Maxwell-Boltzmann equation weighted by the given distribution of formation energies. We calculate the variation of vacancy concentration as function of Ni content in the (CoCrFeMn)$_{x_\mathrm{Ni}}$Ni$_{1-x_\mathrm{Ni}}$ and prove the excellent agreement of the thermodynamic model and the results from the grand-canonical Monte-Carlo simulations.


title: Constrained Parameterization of Reduced Rank and Co-integrated Vector Autoregression

abstract: The paper provides a parametrization of Vector Autoregression (VAR) that enables one to look at the parameters associated with unit root dynamics and those associated with stable dynamics separately. The task is achieved via a novel factorization of the VAR polynomial that partitions the polynomial spectrum into unit root and stable and zero roots via polynomial factors. The proposed factorization adds to the literature of spectral factorization of matrix polynomials. The main benefit is that using the parameterization, actions could be taken to model the dynamics due to a particular class of roots, e.g. unit roots or zero roots, without changing the properties of the dynamics due to other roots. For example, using the parameterization one is able to estimate cointegrating space with appropriate rank that maintains the root structure of the original VAR processes or one can estimate a reduced rank causal VAR process maintaining the constraints of causality. In essence, this parameterization provides the practitioner an option to perform estimation of VAR processes with constrained root structure (e.g., conintegrated VAR or reduced rank VAR) such that the estimated model maintains the assumed root structure.


title: ReStyle: A Residual-Based StyleGAN Encoder via Iterative Refinement

abstract: Recently, the power of unconditional image synthesis has significantly advanced through the use of Generative Adversarial Networks (GANs). The task of inverting an image into its corresponding latent code of the trained GAN is of utmost importance as it allows for the manipulation of real images, leveraging the rich semantics learned by the network. Recognizing the limitations of current inversion approaches, in this work we present a novel inversion scheme that extends current encoder-based inversion methods by introducing an iterative refinement mechanism. Instead of directly predicting the latent code of a given real image using a single pass, the encoder is tasked with predicting a residual with respect to the current estimate of the inverted latent code in a self-correcting manner. Our residual-based encoder, named ReStyle, attains improved accuracy compared to current state-of-the-art encoder-based methods with a negligible increase in inference time. We analyze the behavior of ReStyle to gain valuable insights into its iterative nature. We then evaluate the performance of our residual encoder and analyze its robustness compared to optimization-based inversion and state-of-the-art encoders.


title: Detecting Bosonic Dark Matter with Neutron Stars

abstract: What if the dark matter-nucleon scattering cross section is too small to be detected by direct detection experiments? It is well known in the literature that some interactions lead to dark matter-nucleon scattering cross sections that can be velocity and momentum suppressed. We show that in the case of bosonic dark matter, neutron star spectroscopy offers a possible detection. Firstly, we discuss the case of scalar dark matter with scalar, pseudoscalar, and vector mediators. Later, we do this exercise for vector dark matter. We show that, depending on the nature of dark matter and the interaction involved, neutron stars can improve the sensitivity on the dark matter-nucleon scattering cross section by orders of magnitude, representing a major step forward in the dark matter siege.


title: Demazure Formulas for Weight Polytopes

abstract: The characters of simple Lie algebras are naturally decomposed into lattice polytope sums. The Brion formula for those polytope sums is remarkably similar to the Weyl character formula. Here we start to investigate if other character formulas have analogs for lattice polytope sums, by focusing on the Demazure character formulas. Using Demazure operators, we write expressions for the lattice sums of the weight polytopes of rank-2 simple Lie algebras, and the rank-3 algebra $A_3$.


title: Earth's carbon deficit caused by early loss through irreversible sublimation

abstract: Carbon is an essential element for life but its behavior during Earth's accretion is not well understood. Carbonaceous grains in meteoritic and cometary materials suggest that irreversible sublimation, and not condensation, governs carbon acquisition by terrestrial worlds. Through astronomical observations and modeling we show that the sublimation front of carbon carriers in the solar nebula, or the soot line, moved inward quickly so that carbon-rich ingredients would be available for accretion at 1 au after the first million years. On the other hand, geological constraints firmly establish a severe carbon deficit in Earth, requiring the destruction of inherited carbonaceous organics in the majority of its building blocks. The carbon-poor nature of the Earth thus implies carbon loss in its precursor material through sublimation within the first million years.


title: Adversarial Robustness under Long-Tailed Distribution

abstract: Adversarial robustness has attracted extensive studies recently by revealing the vulnerability and intrinsic characteristics of deep networks. However, existing works on adversarial robustness mainly focus on balanced datasets, while real-world data usually exhibits a long-tailed distribution. To push adversarial robustness towards more realistic scenarios, in this work we investigate the adversarial vulnerability as well as defense under long-tailed distributions. In particular, we first reveal the negative impacts induced by imbalanced data on both recognition performance and adversarial robustness, uncovering the intrinsic challenges of this problem. We then perform a systematic study on existing long-tailed recognition methods in conjunction with the adversarial training framework. Several valuable observations are obtained: 1) natural accuracy is relatively easy to improve, 2) fake gain of robust accuracy exists under unreliable evaluation, and 3) boundary error limits the promotion of robustness. Inspired by these observations, we propose a clean yet effective framework, RoBal, which consists of two dedicated modules, a scale-invariant classifier and data re-balancing via both margin engineering at training stage and boundary adjustment during inference. Extensive experiments demonstrate the superiority of our approach over other state-of-the-art defense methods. To our best knowledge, we are the first to tackle adversarial robustness under long-tailed distributions, which we believe would be a significant step towards real-world robustness. Our code is available at: https://github.com/wutong16/Adversarial_Long-Tail .


title: Blow the Dog Whistle: A Chinese Dataset for Cant Understanding with Common Sense and World Knowledge

abstract: Cant is important for understanding advertising, comedies and dog-whistle politics. However, computational research on cant is hindered by a lack of available datasets. In this paper, we propose a large and diverse Chinese dataset for creating and understanding cant from a computational linguistics perspective. We formulate a task for cant understanding and provide both quantitative and qualitative analysis for tested word embedding similarity and pretrained language models. Experiments suggest that such a task requires deep language understanding, common sense, and world knowledge and thus can be a good testbed for pretrained language models and help models perform better on other tasks. The code is available at https://github.com/JetRunner/dogwhistle. The data and leaderboard are available at https://competitions.codalab.org/competitions/30451.


title: deepregression: a Flexible Neural Network Framework for Semi-Structured Deep Distributional Regression

abstract: In this paper we describe the implementation of semi-structured deep distributional regression, a flexible framework to learn conditional distributions based on the combination of additive regression models and deep networks. Our implementation encompasses (1) a modular neural network building system based on the deep learning library \pkg{TensorFlow} for the fusion of various statistical and deep learning approaches, (2) an orthogonalization cell to allow for an interpretable combination of different subnetworks, as well as (3) pre-processing steps necessary to set up such models. The software package allows to define models in a user-friendly manner via a formula interface that is inspired by classical statistical model frameworks such as \pkg{mgcv}. The packages' modular design and functionality provides a unique resource for both scalable estimation of complex statistical models and the combination of approaches from deep learning and statistics. This allows for state-of-the-art predictive performance while simultaneously retaining the indispensable interpretability of classical statistical models.


title: Early volatile depletion on planetesimals inferred from C-S systematics of iron meteorite parent bodies

abstract: During the formation of terrestrial planets, volatile loss may occur through nebular processing, planetesimal differentiation, and planetary accretion. We investigate iron meteorites as an archive of volatile loss during planetesimal processing. The carbon contents of the parent bodies of magmatic iron meteorites are reconstructed by thermodynamic modelling. Calculated solid/molten alloy partitioning of C increases greatly with liquid S concentration and inferred parent body C concentrations range from 0.0004 to 0.11 wt.%. Parent bodies fall into 2 compositional clusters characterized by cores with medium, and low C/S. Both of these require significant planetesimal degassing, as metamorphic devolatilization on chondrite-like precursors is insufficient to account for their C depletions. Planetesimal core formation models, ranging from closed system extraction to degassing of a wholly molten body, show that significant open system silicate melting and volatile loss is required to match medium and low C/S parent body core compositions. Greater depletion in C relative to S is the hallmark of silicate degassing, indicating that parent body core compositions record processes that affect composite silicate/iron planetesimals. Degassing of bare cores stripped of their silicate mantles would deplete S with negligible C loss, and could not account for inferred parent body core compositions. Devolatilization during small-body differentiation is thus a key process in shaping the volatile inventory of terrestrial planets derived from planetesimals and planetary embryos.


title: Macroscale behavior of random lower triangular matrices

abstract: We analyze the macroscale behavior of random lower (and therefore upper) triangular matrices with entries drawn iid from a distribution with nonzero mean and finite variance. We show that such a matrix behaves like a probabilistic version of a Riemann sum and therefore in the limit behaves like the Volterra operator. Specifically, we analyze certain SOT-like and WOT-like modes of convergence for random lower triangular matrices to a scaled Volterra operator. We close with a brief discussion of moments.


title: Realising Einstein's mirror: Optomechanical damping with a thermal photon gas

abstract: In 1909 Einstein described the thermalization of a mirror within a blackbody cavity by collisions with thermal photons. While the time to thermalize the motion of even a microscale or nanoscale object is so long that it is not feasible, we show that it is using the high intensity light from an amplified thermal light source with a well-defined chemical potential. We predict damping of the center-of mass motion due to this effect on times scales of seconds for small optomechanical systems, such as levitated nanoparticles, allowing experimental observation.


title: Adaptive Variants of Optimal Feedback Policies

abstract: The stable combination of optimal feedback policies with online learning is studied in a new control-theoretic framework for uncertain nonlinear systems. The framework can be systematically used in transfer learning and sim-to-real applications, where an optimal policy learned for a nominal system needs to remain effective in the presence of significant variations in parameters. Given unknown parameters within a bounded range, the resulting adaptive control laws guarantee convergence of the closed-loop system to the state of zero cost. Online adjustment of the learning rate is used as a key stability mechanism, and preserves certainty equivalence when designing optimal policies without assuming uncertainty to be within the control range. The approach is illustrated on the familiar mountain car problem, where it yields near-optimal performance despite the presence of parametric model uncertainty.


title: The Multi-Agent Behavior Dataset: Mouse Dyadic Social Interactions

abstract: Multi-agent behavior modeling aims to understand the interactions that occur between agents. We present a multi-agent dataset from behavioral neuroscience, the Caltech Mouse Social Interactions (CalMS21) Dataset. Our dataset consists of trajectory data of social interactions, recorded from videos of freely behaving mice in a standard resident-intruder assay. To help accelerate behavioral studies, the CalMS21 dataset provides benchmarks to evaluate the performance of automated behavior classification methods in three settings: (1) for training on large behavioral datasets all annotated by a single annotator, (2) for style transfer to learn inter-annotator differences in behavior definitions, and (3) for learning of new behaviors of interest given limited training data. The dataset consists of 6 million frames of unlabeled tracked poses of interacting mice, as well as over 1 million frames with tracked poses and corresponding frame-level behavior annotations. The challenge of our dataset is to be able to classify behaviors accurately using both labeled and unlabeled tracking data, as well as being able to generalize to new settings.


title: A Bombieri-Vinogradov theorem for higher rank groups

abstract: We establish a result of Bombieri-Vinogradov type for the Dirichlet coefficients at prime ideals of the standard $L$-function associated to a self-dual cuspidal automorphic representation $\pi$ of $\mathrm{GL}_n$ over a number field $F$ which is not a quadratic twist of itself. Our result does not rely on any unproven progress towards the generalized Ramanujan conjecture or the nonexistence of Landau-Siegel zeros. In particular, when $\pi$ is fixed and not equal to a quadratic twist of itself, we prove the first unconditional Siegel-type lower bound for the twisted $L$-values $|L(1,\pi\otimes\chi)|$ in the $\chi$-aspect, where $\chi$ is a primitive quadratic Hecke character over $F$. Our result improves the levels of distribution in other works that relied on these unproven hypotheses. As applications, when $n=2,3,4$, we prove a $\mathrm{GL}_n$ analogue of the Titchmarsh divisor problem and a nontrivial bound for a certain $\mathrm{GL}_n\times\mathrm{GL}_2$ shifted convolution sum.


title: Hypothesis Formalization: Empirical Findings, Software Limitations, and Design Implications

abstract: Data analysis requires translating higher level questions and hypotheses into computable statistical models. We present a mixed-methods study aimed at identifying the steps, considerations, and challenges involved in operationalizing hypotheses into statistical models, a process we refer to as hypothesis formalization. In a formative content analysis of research papers, we find that researchers highlight decomposing a hypothesis into sub-hypotheses, selecting proxy variables, and formulating statistical models based on data collection design as key steps. In a lab study, we find that analysts fixated on implementation and shaped their analysis to fit familiar approaches, even if sub-optimal. In an analysis of software tools, we find that tools provide inconsistent, low-level abstractions that may limit the statistical models analysts use to formalize hypotheses. Based on these observations, we characterize hypothesis formalization as a dual-search process balancing conceptual and statistical considerations constrained by data and computation, and discuss implications for future tools.


title: The Effects of Process Parameters on Melt-pool Oscillatory Behaviour in Gas Tungsten Arc Welding

abstract: Internal flow behaviour and melt-pool surface oscillations during arc welding are complex and not yet fully understood. In the present work, high-fidelity numerical simulations are employed to describe the effects of welding position, sulphur concentration (60-300 ppm) and travel speed (1.25-5 mm/s) on molten metal flow dynamics in fully-penetrated melt-pools. A wavelet transform is implemented to obtain time-resolved frequency spectra of the oscillation signals, which overcomes the shortcomings of the Fourier transform in rendering time resolution of the frequency spectra. Comparing the results of the present numerical calculations with available analytical and experimental datasets, the robustness of the proposed approach in predicting melt-pool oscillations is demonstrated. The results reveal that changes in the surface morphology of the pool resulting from a change in welding position alter the spatial distribution of arc forces and power-density applied to the molten material, and in turn affect flow patterns in the pool. Under similar welding conditions, changing the sulphur concentration affects the Marangoni flow pattern, and increasing the travel speed decreases the size of the pool and increases the offset between top and bottom melt-pool surfaces, affecting the flow structures (vortex formation) on the surface. Variations in the internal flow pattern affect the evolution of melt-pool shape and its surface oscillations.


title: Lipschitz geometry of operator spaces and Lipschitz-free operator spaces

abstract: We show that there is an operator space notion of Lipschitz embeddability between operator spaces which is strictly weaker than its linear counterpart but which is still strong enough to impose linear restrictions on operator space structures. This shows that there is a nontrivial theory of nonlinear geometry for operator spaces and it answers a question in [BCD21]. For that, we introduce the operator space version of Lipschitz-free Banach spaces and prove several properties of it. In particular, we show that separable operator spaces satisfy a sort of isometric Lipschitz-lifting property in the sense of G. Godefroy and N. Kalton. Gateaux differentiability of Lipschitz maps in the operator space category is also studied.


title: The sum of powers of subtree sizes for conditioned Galton-Watson trees

abstract: We study the additive functional $X_n(\alpha)$ on conditioned Galton-Watson trees given, for arbitrary complex $\alpha$, by summing the $\alpha$th power of all subtree sizes. Allowing complex $\alpha$ is advantageous, even for the study of real $\alpha$, since it allows us to use powerful results from the theory of analytic functions in the proofs. For $\Re\alpha < 0$, we prove that $X_n(\alpha)$, suitably normalized, has a complex normal limiting distribution; moreover, as processes in $\alpha$, the weak convergence holds in the space of analytic functions in the left half-plane. We establish, and prove similar process-convergence extensions of, limiting distribution results for $\alpha$ in various regions of the complex plane. We focus mainly on the case where $\Re\alpha > 0$, for which $X_n(\alpha)$, suitably normalized, has a limiting distribution that is not normal but does not depend on the offspring distribution $\xi$ of the conditioned Galton-Watson tree, assuming only that $E[\xi] = 1$ and $0 < \mathrm{Var} [\xi] < \infty$. Under a weak extra moment assumption on $\xi$, we prove that the convergence extends to moments, ordinary and absolute and mixed, of all orders. At least when $\Re\alpha > \frac12$, the limit random variable $Y(\alpha)$ can be expressed as a function of a normalized Brownian excursion.


title: Second-order PDEs in 3D with Einstein-Weyl conformal structure

abstract: Einstein-Weyl geometry is a triple (D,g,w), where D is a symmetric connection, [g] is a conformal structure and w is a covector such that: (i) connection D preserves the conformal class [g], that is, Dg=wg; (ii) trace-free part of the symmetrised Ricci tensor of D vanishes. Three-dimensional Einstein-Weyl structures arise naturally on solutions of second-order dispersionless integrable PDEs in 3D. In this context, [g] coincides with the characteristic conformal structure and is therefore uniquely determined by the equation. On the contrary, the covector w is a somewhat more mysterious object, recovered from the Einstein-Weyl conditions. We demonstrate that, for generic second-order PDEs (for instance, for all equations not of Monge-Ampere type), the covector w is also expressible in terms of the equation, thus providing an efficient dispersionless integrability test. The knowledge of g and w provides a dispersionless Lax pair by an explicit formula which is apparently new. Some partial classification results of PDEs with Einstein-Weyl characteristic conformal structure are obtained. A rigidity conjecture is proposed according to which for any generic second-order PDE with Einstein-Weyl property, all dependence on the 1-jet variables can be eliminated via a suitable contact transformation.


title: Discussion of "A Gibbs sampler for a class of random convex polytopes"

abstract: An exciting new algorithmic breakthrough has been advanced for how to carry out inferences in a Dempster-Shafer (DS) formulation of a categorical data generating model. The developed sampling mechanism, which draws on theory for directed graphs, is a clever and remarkable achievement, as this has been an open problem for many decades. In this discussion, I comment on important contributions, central questions, and prevailing matters of the article.


title: Attractors with Large Complex Structure for One-Parameter Families of Calabi-Yau Manifolds

abstract: The attractor equations for an arbitrary one-parameter family of Calabi-Yau manifolds are studied in the large complex structure region. These equations are solved iteratively, generating what we term an N-expansion, which is a power series in the Gromov-Witten invariants of the manifold. The coefficients of this series are associated with integer partitions. In important cases we are able to find closed-form expressions for the general term of this expansion. To our knowledge, these are the first generic solutions to attractor equations that incorporate instanton contributions. In particular, we find a simple closed-form formula for the entropy associated to rank two attractor points, including those recently discovered. The applications of our solutions are briefly discussed. Most importantly, we are able to give an expression for the Wald entropy of black holes that includes all genus 0 instanton corrections.


title: NGC 5128 globular cluster candidates out to 150 kpc: a comprehensive catalog from Gaia and ground based data

abstract: We present a new catalog of 40502 globular cluster (GC) candidates in NGC 5128 out to a projected radius of $\sim$150 kpc, based on data from the Panoramic Imaging Survey of Centaurus and Sculptor (PISCeS), Gaia Data Release 2, and the NOAO Source Catalog. Ranking these candidates based on the likelihood that they are true GCs, we find that approximately 1900 belong to our top two ranking categories and should be the highest priority for spectroscopic follow-up for confirmation. Taking into account our new data and a vetting of previous GC catalogs, we estimate a total GC population of $1450 \pm 160$ GCs. We show that a substantial number of sources previously argued to be low-velocity GCs are instead foreground stars, reducing the inferred GC velocity dispersion. This work showcases the power of Gaia to identify slightly extended sources at the $\sim 4$ Mpc distance of NGC 5128, enabling accurate identification of GCs throughout the entire extended halo, not just the inner regions that have been the focus of most previous work.


title: Interactions-disorder duality and critical phenomena in nodal semimetals, dilute gases and other systems

abstract: We investigate classes of interacting systems that allow for a mapping to disordered noninteracting systems. As we show, such a mapping is possible for interacting systems with a suppressed density of states at the chemical potential, leading to suppressed screening, and systems near BCS-type instabilities. The mapping can also be applied qualitatively to other classes of systems that are not exactly dual to each other. The established duality suggests a new approach to analytical and numerical studies of many-body and disorder-driven phenomena in a variety of systems and allows to predict, e.g., new phase transitions dual to the previously known ones. Using the established duality, we predict new disorder-driven transitions in nodal-line semimetals and systems with long-range hopping dual to, respectively, the BCS and BEC-vacuum transitions in interacting systems and new interaction-driven transitions dual to previously known non-Anderson disorder-driven transitions. The established principle can also be used to classify and describe phase transitions in dissipative systems described by non-Hermitian Hamiltonians.

categories: cond-mat.mes-hall cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.quant-gas cond-mat.str-el


title: Hierarchical compressed sensing

abstract: Compressed sensing is a paradigm within signal processing that provides the means for recovering structured signals from linear measurements in a highly efficient manner. Originally devised for the recovery of sparse signals, it has become clear that a similar methodology would also carry over to a wealth of other classes of structured signals. In this work, we provide an overview over the theory of compressed sensing for a particularly rich family of such signals, namely those of hierarchically structured signals. Examples of such signals are constituted by blocked vectors, with only few non-vanishing sparse blocks. We present recovery algorithms based on efficient hierarchical hard-thresholding. The algorithms are guaranteed to converge, in a stable fashion both with respect to measurement noise as well as to model mismatches, to the correct solution provided the measurement map acts isometrically restricted to the signal class. We then provide a series of results establishing the required condition for large classes of measurement ensembles. Building upon this machinery, we sketch practical applications of this framework in machine-type communications and quantum tomography.


title: Multi-spacecraft observations of coronal loops to verify a force-free field reconstruction and infer loop cross sections

abstract: Active region EUV loops are believed to trace a subset of magnetic field lines through the corona. Malanushenko et al. (2009) proposed a method, using loop images and line-of-sight photospheric magnetograms, to infer the three-dimensional shape and field strength along each loop. McCarthy et al. (2019) used this novel method to compute the total magnetic flux interconnecting a pair of active regions observed by SDO/AIA. They adopted the common assumption that each loop had a circular cross section. The accuracy of inferred shape and circularity of cross sections can both be tested using observations of the same loops from additional vantage points as provided by STEREO/EUVI. Here, we use multiple viewing angles to confirm the three-dimensional structure of loops. Of 151 viable cases, 105 (69.5%) matched some form of visible coronal structure when viewed approximately in quadrature. A loop with a circular cross-section should appear of a similar width in different perspectives. In contradiction to this, we find a puzzling lack of correlation between loop diameters seen from different perspectives, even an anti-correlation in some cases. Features identified as monolithic loops in AIA may, in fact, be more complex density enhancements. The 30.5% of reconstructions from AIA which did not match any feature in EUVI might be such enhancements. Others may be genuine loop structures, but with elliptical cross sections. We observe an anti-correlation between diameter and brightness, lending support to the latter hypothesis. Of 13 suitable for width analysis, four loops are consistent with non-circular cross sections, where we find anti-correlation in both comparisons.


title: Parton distributions in the SMEFT from high-energy Drell-Yan tails

abstract: The high-energy tails of charged- and neutral-current Drell-Yan processes provide important constraints on the light quark and anti-quark parton distribution functions (PDFs) in the large-x region. At the same time, short-distance new physics effects such as those encoded by the Standard Model Effective Field Theory (SMEFT) would induce smooth distortions to the same high-energy Drell-Yan tails. In this work, we assess for the first time the interplay between PDFs and EFT effects for high-mass Drell-Yan processes at the LHC and quantify the impact that the consistent joint determination of PDFs and Wilson coefficients has on the bounds derived for the latter. We consider two well-motivated new physics scenarios: $1)$ electroweak oblique corrections $(\hat W, \hat Y)$ and $2)$ four-fermion interactions potentially related to the LHCb anomalies in $R(K^{(*)})$. We account for available Drell-Yan data, both from unfolded cross sections and from searches, and carry out dedicated projections for the High-Luminosity LHC. Our main finding is that, while the interplay between PDFs and EFT effects remains moderate for the current dataset, it will become a significant challenge for EFT analyses at the HL-LHC.


title: Relaxing the Conditional Independence Assumption of CTC-based ASR by Conditioning on Intermediate Predictions

abstract: This paper proposes a method to relax the conditional independence assumption of connectionist temporal classification (CTC)-based automatic speech recognition (ASR) models. We train a CTC-based ASR model with auxiliary CTC losses in intermediate layers in addition to the original CTC loss in the last layer. During both training and inference, each generated prediction in the intermediate layers is summed to the input of the next layer to condition the prediction of the last layer on those intermediate predictions. Our method is easy to implement and retains the merits of CTC-based ASR: a simple model architecture and fast decoding speed. We conduct experiments on three different ASR corpora. Our proposed method improves a standard CTC model significantly (e.g., more than 20 % relative word error rate reduction on the WSJ corpus) with a little computational overhead. Moreover, for the TEDLIUM2 corpus and the AISHELL-1 corpus, it achieves a comparable performance to a strong autoregressive model with beam search, but the decoding speed is at least 30 times faster.


title: Complex charge density waves at Van Hove singularity on hexagonal lattices: Haldane-model phase diagram and potential realization in kagome metals $\text{AV}_3\text{Sb}_5$

abstract: We investigate how the real and imaginary charge density waves interplay at the Van Hove singularity on the hexagonal lattices. A phenomenological analysis indicates the formation of $3Q$ complex orders at all three nesting momenta. Under a total phase condition, unequal phases at the three momenta break the rotation symmetry generally. The $3Q$ complex orders constitute a rich Haldane-model phase diagram. When effective time-reversal symmetries arise under 1-site translations, the Dirac semimetals are protected. The breakdown of these symmetries gaps the Dirac points and leads to the trivial and Chern insulator phases. These phases are deformations of purely real and imaginary orders, which exhibit trivial site and/or bond density and chiral flux orders, respectively. The exotic single-Dirac-point semimetals also appear along the gapless phase boundary. We further show that the theoretical model offers transparent interpretations of experimental observations in the kagome metals $\text{AV}_3\text{Sb}_5$ with $\text{A}=\text{K},\text{Rb},\text{Cs}$. The topological charge density waves may be identified with the complex orders in the Chern insulator phase. Meanwhile, the lower-temperature symmetry-breaking phenomena may be interpreted as the secondary orders from the complex order ground states. Our work sheds light on the nature of the topological charge density waves in the kagome metals $\text{AV}_3\text{Sb}_5$ and may offer useful indications to the experimentally observed charge orders in the future experiments.


title: Creativity and Machine Learning: A Survey

abstract: There is a growing interest in the area of machine learning and creativity. This survey presents an overview of the history and the state of the art of computational creativity theories, key machine learning techniques (including generative deep learning), and corresponding automatic evaluation methods. After presenting a critical discussion of the key contributions in this area, we outline the current research challenges and emerging opportunities in this field.


title: The Reservoir Learning Power across Quantum Many-Boby Localization Transition

abstract: Harnessing the quantum computation power of the present noisy-intermediate-size-quantum devices has received tremendous interest in the last few years. Here we study the learning power of a one-dimensional long-range randomly-coupled quantum spin chain, within the framework of reservoir computing. In time sequence learning tasks, we find the system in the quantum many-body localized (MBL) phase holds long-term memory, which can be attributed to the emergent local integrals of motion. On the other hand, MBL phase does not provide sufficient nonlinearity in learning highly-nonlinear time sequences, which we show in a parity check task. This is reversed in the quantum ergodic phase, which provides sufficient nonlinearity but compromises memory capacity. In a complex learning task of Mackey-Glass prediction that requires both sufficient memory capacity and nonlinearity, we find optimal learning performance near the MBL-to-ergodic transition. This leads to a guiding principle of quantum reservoir engineering at the edge of quantum ergodicity reaching optimal learning power for generic complex reservoir learning tasks. Our theoretical finding can be readily tested with present experiments.


title: Unbiased likelihood-free inference of the Hubble constant from light standard sirens

abstract: Multi-messenger observations of binary neutron star mergers offer a promising path towards resolution of the Hubble constant ($H_0$) tension, provided their constraints are shown to be free from systematics such as the Malmquist bias. In the traditional Bayesian framework, accounting for selection effects in the likelihood requires calculation of the expected number (or fraction) of detections as a function of the parameters describing the population and cosmology; a potentially costly and/or inaccurate process. This calculation can, however, be bypassed completely by performing the inference in a framework in which the likelihood is never explicitly calculated, but instead fit using forward simulations of the data, which naturally include the selection. This is Likelihood-Free Inference (LFI). Here, we use density-estimation LFI, coupled to neural-network-based data compression, to infer $H_0$ from mock catalogues of binary neutron star mergers, given noisy redshift, distance and peculiar velocity estimates for each object. We demonstrate that LFI yields statistically unbiased estimates of $H_0$ in the presence of selection effects, with precision matching that of sampling the full Bayesian hierarchical model. Marginalizing over the bias increases the $H_0$ uncertainty by only $6\%$ for training sets consisting of $O(10^4)$ populations. The resulting LFI framework is applicable to population-level inference problems with selection effects across astrophysics.


title: Configuration spaces of clusters as $E_d$-algebras

abstract: It is a classical result that configuration spaces of labelled particles in $\mathbb{R}^d$ are free $E_d$-algebras and that their $d$-fold bar construction is equivalent to the $d$-fold suspension of the labelling space. In this paper, we study a variation of these spaces, namely configuration spaces of labelled clusters of particles. These configuration spaces are again $E_d$-algebras, and we give geometric models for their iterated bar construction in two different ways: one establishes a description of these configuration spaces of clusters as cellular $E_1$-algebras, and the other one uses an additional verticality constraint. In the last section, we apply these results in order to calculate the stable homology of certain vertical configuration spaces.


title: Fermion masses and mixings, dark matter, leptogenesis and $g-2$ muon anomaly in an extended 2HDM with inverse seesaw

abstract: We propose a predictive $Q_4$ flavored 2HDM model, where the scalar sector is enlarged by the inclusion of several gauge singlet scalars and the fermion sector by the inclusion of right handed Majorana neutrinos. In our model, the $Q_4$ family symmetry is supplemented by several auxiliary cyclic symmetries, whose spontaneous breaking produces the observed pattern of SM charged fermion masses and quark mixing angles. The light active neutrino masses are generated from an inverse seesaw mechanism at one loop level thanks to a remnant preserved $Z_2$ symmetry. Our model succesfully reproduces the measured dark matter relic abundance and is consistent with direct detection constraints for masses of the DM candidate around $\sim$ 6.3 TeV. Furthermore, our model is also consistent with the lepton and baryon asymmetries of the Universe as well as with the muon anomalous magnetic moment


title: Structural Organization of Space Polymers

abstract: Extra-terrestrial polymers of glycine with iron have been characterized by mass spectrometry to have a core mass of 1494Da with dominant rod-like variants at m/z 1567 and m/z 1639 [1]. Several principal macro-structural morphologies are observed in solvent extracts from CV3 class meteoritic material. The first is an extended sheet of linked triskelia containing the 1494Da core entity that encloses gas bubbles in the solvent [1]. A second is of fiber-like crystals found here, via X-ray diffraction, to be multiple-walled nanotubes made from a square lattice of the 1494Da polymer. A third is a dispersion of floating phantom-like short tubes of up to 100micron length [1] with characteristic angled bends that suggest the influence of a specific underlying protein structure. Here it is proposed that the angled tubes are the observable result of a space-filling super-polymerization of 1638Da polymer subunits guided by the tetragonal symmetry of linking silicon bonds. Distorted hexagonal sheets are linked by perpendicular subunits in a three-dimensional hexagonal diamond structure to fill the largest possible volume. This extended very low-density structure is conjectured to have dominated in a process of chemical selection because it captured a maximum amount of molecular raw material in the ultra-low density of molecular clouds or of the proto-solar nebula. This could have led ultimately to the accretion of the earliest planetary bodies.


title: Dirac-like Hamiltonians associated to Schr\"odinger factorizations

abstract: In this work, we have extended the factorization method of scalar shape-invariant Schr\"o\-din\-ger Hamiltonians to a class of Dirac-like matrix Hamiltonians. The intertwining operators of the Schr\"odinger equations have been implemented in the Dirac-like shape invariant equations. We have considered also another kind of anti-intertwining operators changing the sign of energy. The Dirac-like Hamiltonians can be obtained from reduction of higher dimensional spin systems. Two examples have been worked out, one obtained from the sphere ${\cal S}^2$ and a second one, having a non-Hermitian character, from the hyperbolic space ${\cal H}^2$.


title: On the Quenched Free Energy of JT Gravity and Supergravity

abstract: The quenched free energy, $F_Q(T){=}{-}T\langle \ln Z(T)\rangle$, of various JT gravity and supergravity theories is explored, taking into account the key non-perturbative physics that is accessible using their matrix model formulations. The leading low energy physics of these systems can be modelled by the Airy and (a family of) Bessel models, which arise from scaling limits of matrix ensembles. The $F_Q(T)$s of these models are directly computed by explicit sampling of the matrix ensembles, and how their properties are connected to the statistical mechanics of the underlying discrete spectrum of the ensembles is elucidated. Some of the low temperature ($T$) features of the results confirm recent observations by Jassen and Mirbabayi. The results are then used as benchmarks for exploring an intriguing formula proposed by Okuyama for computing $F_Q(T)$ in terms of the connected correlators of its partition function, the wormholes of the gravity theory. A low $T$ truncation of the correlators helps render the formula practical, but it is shown that this is at the expense of much of its accuracy. The significance of the statistical interpretation of $F_Q(T)$ for black hole microphysics is discussed.


title: Online change-point detection for a transient change

abstract: We consider a popular online change-point problem of detecting a transient change in distributions of i.i.d. random variables. For this change-point problem, several change-point procedures are formulated and some advanced results for a particular procedure are surveyed. Some new approximations for the average run length to false alarm are offered and the power of these procedures for detecting a transient change in mean of a sequence of normal random variables is compared.


title: Visual Vibration Tomography: Estimating Interior Material Properties from Monocular Video

abstract: An object's interior material properties, while invisible to the human eye, determine motion observed on its surface. We propose an approach that estimates heterogeneous material properties of an object from a monocular video of its surface vibrations. Specifically, we show how to estimate Young's modulus and density throughout a 3D object with known geometry. Knowledge of how these values change across the object is useful for simulating its motion and characterizing any defects. Traditional non-destructive testing approaches, which often require expensive instruments, generally estimate only homogenized material properties or simply identify the presence of defects. In contrast, our approach leverages monocular video to (1) identify image-space modes from an object's sub-pixel motion, and (2) directly infer spatially-varying Young's modulus and density values from the observed modes. We demonstrate our approach on both simulated and real videos.


title: Diagnosing collisions in the interior of a wormhole

abstract: Two distant black holes can be connected in the interior through a wormhole. Such a wormhole has been interpreted as an entangled state shared between two exterior regions. If Alice and Bob send signals into each of the black holes, they can meet in the interior. In this letter, we interpret this meeting in terms of the quantum circuit that prepares the entangled state: Alice and Bob sending signals creates growing perturbations in the circuit, whose overlap represents their meeting inside the wormhole. We argue that such overlap in the circuit is quantified by a particular six-point correlation function. Therefore, exterior observers in possession of the entangled qubits can use this correlation function to diagnose the collision in the interior without having to jump in themselves.


title: Neural Network-based Control for Multi-Agent Systems from Spatio-Temporal Specifications

abstract: We propose a framework for solving control synthesis problems for multi-agent networked systems required to satisfy spatio-temporal specifications. We use Spatio-Temporal Reach and Escape Logic (STREL) as a specification language. For this logic, we define smooth quantitative semantics, which captures the degree of satisfaction of a formula by a multi-agent team. We use the novel quantitative semantics to map control synthesis problems with STREL specifications to optimization problems and propose a combination of heuristic and gradient-based methods to solve such problems. As this method might not meet the requirements of a real-time implementation, we develop a machine learning technique that uses the results of the off-line optimizations to train a neural network that gives the control inputs at current states. We illustrate the effectiveness of the proposed framework by applying it to a model of a robotic team required to satisfy a spatial-temporal specification under communication constraints.


title: Strain control of the metal-insulator transition in epitaxial SrCrO$_{3}$ thin films

abstract: In the perovskite oxide SrCrO$_{3}$ the interplay between crystal structure, strain and orbital ordering enables a transition from a metallic to an insulating electronic structure under certain conditions. We identified a narrow window of oxygen partial pressure in which highly strained SrCrO$_{3}$ thin films can be grown using radio-frequency (RF) off-axis magnetron sputtering on three different substrates, (LaAlO$_{3}$)$_{0.3}$-(Sr$_{2}$TaAlO$_{6}$)$_{0.7}$ (LSAT), SrTiO$_{3}$ (STO) and DyScO$_{3}$ (DSO). X-ray diffraction and atomic force microscopy confirmed the quality of the films and a metal-insulator transition driven by the substrate induced strain was demonstrated.


title: Differentially Private Histograms in the Shuffle Model from Fake Users

abstract: There has been much recent work in the shuffle model of differential privacy, particularly for approximate $d$-bin histograms. While these protocols achieve low error, the number of messages sent by each user -- the message complexity -- has so far scaled with $d$ or the privacy parameters. The message complexity is an informative predictor of a shuffle protocol's resource consumption. We present a protocol whose message complexity is two when there are sufficiently many users. The protocol essentially pairs each row in the dataset with a fake row and performs a simple randomization on all rows. We show that the error introduced by the protocol is small, using rigorous analysis as well as experiments on real-world data. We also prove that corrupt users have a relatively low impact on our protocol's estimates.


title: The Varchenko-Gel'fand Ring of a Cone

abstract: For a hyperplane arrangement in a real vector space, the coefficients of its Poincar\'{e} polynomial have many interpretations. An interesting one is provided by the Varchenko-Gel'fand ring, which is the ring of functions from the chambers of the arrangement to the integers with pointwise addition and multiplication. Varchenko and Gel'fand gave a simple presentation for this ring, along with a filtration and associated graded ring whose Hilbert series is the Poincar\'{e} polynomial. We generalize these results to cones defined by intersections of halfspaces of some of the hyperplanes and prove a novel result for the Varchenko-Gel'fand ring of an arrangement: when the arrangement is supersolvable the associated graded ring of the arrangement is Koszul.


title: Conditional physics informed neural networks

abstract: We introduce conditional PINNs (physics informed neural networks) for estimating the solution of classes of eigenvalue problems. The concept of PINNs is expanded to learn not only the solution of one particular differential equation but the solutions to a class of problems. We demonstrate this idea by estimating the coercive field of permanent magnets which depends on the width and strength of local defects. When the neural network incorporates the physics of magnetization reversal, training can be achieved in an unsupervised way. There is no need to generate labeled training data. The presented test cases have been rigorously studied in the past. Thus, a detailed and easy comparison with analytical solutions is made. We show that a single deep neural network can learn the solution of partial differential equations for an entire class of problems.


title: $\texttt{nrCascadeSim}$ -- A simulation tool for nuclear recoil cascades resulting from neutron capture

abstract: Neutron capture-induced nuclear recoils have emerged as an important tool for detector calibrations in direct dark matter detection and coherent elastic neutrino-nucleus scattering (CE${\nu}$NS). $\texttt{nrCascadeSim}$ is a command-line tool for generating simulation data for energy deposits resulting from neutron capture on pure materials. Presently, silicon, germanium, neon, and argon are supported. While the software was developed for solid state detector calibration, it can be used for any application which requires simulated neutron capture-induced nuclear recoil data.


title: Approximate Robust NMPC using Reinforcement Learning

abstract: We present a Reinforcement Learning-based Robust Nonlinear Model Predictive Control (RL-RNMPC) framework for controlling nonlinear systems in the presence of disturbances and uncertainties. An approximate Robust Nonlinear Model Predictive Control (RNMPC) of low computational complexity is used in which the state trajectory uncertainty is modelled via ellipsoids. Reinforcement Learning is then used in order to handle the ellipsoidal approximation and improve the closed-loop performance of the scheme by adjusting the MPC parameters generating the ellipsoids. The approach is tested on a simulated Wheeled Mobile Robot (WMR) tracking a desired trajectory while avoiding static obstacles.


title: Incompleteness of the large-$N$ analysis of the $O(N)$ models: Nonperturbative cuspy fixed points and their nontrivial homotopy at finite $N$

abstract: We summarize the usual implementations of the large $N$ limit of $O(N)$ models and show in detail why and how they can miss some physically important fixed points when they become singular in the limit $N\to\infty$. Using Wilson's renormalization group in its functional nonperturbative versions, we show how the singularities build up as $N$ increases. In the Wilson-Polchinski version of the nonperturbative renormalization group, we show that the singularities are cusps, which become boundary layers for finite but large values of $N$. The corresponding fixed points being never close to the Gaussian, are out of reach of the usual perturbative approaches. We find four new fixed points and study them in all dimensions and for all $N>0$ and show that they play an important role for the tricritical physics of $O(N)$ models. Finally, we show that some of these fixed points are bi-valued when they are considered as functions of $d$ and $N$ thus revealing important and nontrivial homotopy structures. The Bardeen-Moshe-Bander phenomenon that occurs at $N=\infty$ and $d=3$ is shown to play a crucial role for the internal consistency of all our results.


title: InverseForm: A Loss Function for Structured Boundary-Aware Segmentation

abstract: We present a novel boundary-aware loss term for semantic segmentation using an inverse-transformation network, which efficiently learns the degree of parametric transformations between estimated and target boundaries. This plug-in loss term complements the cross-entropy loss in capturing boundary transformations and allows consistent and significant performance improvement on segmentation backbone models without increasing their size and computational complexity. We analyze the quantitative and qualitative effects of our loss function on three indoor and outdoor segmentation benchmarks, including Cityscapes, NYU-Depth-v2, and PASCAL, integrating it into the training phase of several backbone networks in both single-task and multi-task settings. Our extensive experiments show that the proposed method consistently outperforms baselines, and even sets the new state-of-the-art on two datasets.


title: Proof of the Theory-to-Practice Gap in Deep Learning via Sampling Complexity bounds for Neural Network Approximation Spaces

abstract: We study the computational complexity of (deterministic or randomized) algorithms based on point samples for approximating or integrating functions that can be well approximated by neural networks. Such algorithms (most prominently stochastic gradient descent and its variants) are used extensively in the field of deep learning. One of the most important problems in this field concerns the question of whether it is possible to realize theoretically provable neural network approximation rates by such algorithms. We answer this question in the negative by proving hardness results for the problems of approximation and integration on a novel class of neural network approximation spaces. In particular, our results confirm a conjectured and empirically observed theory-to-practice gap in deep learning. We complement our hardness results by showing that approximation rates of a comparable order of convergence are (at least theoretically) achievable.


title: Dissecting the Quadruple Binary Hyad vA 351 -- Masses for three M Dwarfs and a White Dwarf

abstract: We extend results first announced by Franz et al. (1998), that identified vA 351 = H346 in the Hyades as a multiple star system containing a white dwarf. With Hubble Space Telescope Fine Guidance Sensor fringe tracking and scanning, and more recent speckle observations, all spanning 20.7 years, we establish a parallax, relative orbit, and mass fraction for two components, with a period, $P=2.70$y and total mass 2.1 Msun. With ground-based radial velocities from the McDonald Observatory Otto Struve 2.1m telescope Sandiford Spectrograph, and Center for Astrophysics Digital Speedometers, spanning 37 years, we find that component B consists of BC, two M dwarf stars orbiting with a very short period (P_ BC=0.749 days), having a mass ratio M_C/M_B=0.95. We confirm that the total mass of the system can only be reconciled with the distance and component photometry by including a fainter, higher mass component. The quadruple system consists of three M dwarfs (A,B,C) and one white dwarf (D). We determine individual M dwarf masses M_A=0.53+/-0.10 Msun, M_B=0.43+/-0.04Msun, and M_C=0.41+/-0.04Msun. The WD mass, 0.54+/-0.04Msun, comes from cooling models, an assumed Hyades age of 670My, and consistency with all previous and derived astrometric, photometric, and RV results. Velocities from H-alpha and He I emission lines confirm the BC period derived from absorption lines, with similar (He I) and higher (H-alpha) velocity amplitudes. We ascribe the larger H-alpha amplitude to emission from a region each component shadows from the other, depending on the line of sight.


title: Communication-Efficient Agnostic Federated Averaging

abstract: In distributed learning settings such as federated learning, the training algorithm can be potentially biased towards different clients. Mohri et al. (2019) proposed a domain-agnostic learning algorithm, where the model is optimized for any target distribution formed by a mixture of the client distributions in order to overcome this bias. They further proposed an algorithm for the cross-silo federated learning setting, where the number of clients is small. We consider this problem in the cross-device setting, where the number of clients is much larger. We propose a communication-efficient distributed algorithm called Agnostic Federated Averaging (or AgnosticFedAvg) to minimize the domain-agnostic objective proposed in Mohri et al. (2019), which is amenable to other private mechanisms such as secure aggregation. We highlight two types of naturally occurring domains in federated learning and argue that AgnosticFedAvg performs well on both. To demonstrate the practical effectiveness of AgnosticFedAvg, we report positive results for large-scale language modeling tasks in both simulation and live experiments, where the latter involves training language models for Spanish virtual keyboard for millions of user devices.


title: Heuristics2Annotate: Efficient Annotation of Large-Scale Marathon Dataset For Bounding Box Regression

abstract: Annotating a large-scale in-the-wild person re-identification dataset especially of marathon runners is a challenging task. The variations in the scenarios such as camera viewpoints, resolution, occlusion, and illumination make the problem non-trivial. Manually annotating bounding boxes in such large-scale datasets is cost-inefficient. Additionally, due to crowdedness and occlusion in the videos, aligning the identity of runners across multiple disjoint cameras is a challenge. We collected a novel large-scale in-the-wild video dataset of marathon runners. The dataset consists of hours of recording of thousands of runners captured using 42 hand-held smartphone cameras and covering real-world scenarios. Due to the presence of crowdedness and occlusion in the videos, the annotation of runners becomes a challenging task. We propose a new scheme for tackling the challenges in the annotation of such large dataset. Our technique reduces the overall cost of annotation in terms of time as well as budget. We demonstrate performing fps analysis to reduce the effort and time of annotation. We investigate several annotation methods for efficiently generating tight bounding boxes. Our results prove that interpolating bounding boxes between keyframes is the most efficient method of bounding box generation amongst several other methods and is 3x times faster than the naive baseline method. We introduce a novel way of aligning the identity of runners in disjoint cameras. Our inter-camera alignment tool integrated with the state-of-the-art person re-id system proves to be sufficient and effective in the alignment of the runners across multiple cameras with non-overlapping views. Our proposed framework of annotation reduces the annotation cost of the dataset by a factor of 16x, also effectively aligning 93.64% of the runners in the cross-camera setting.


title: Joint Learning for Spatial Context-based Seismic Inversion of Multiple Datasets for Improved Generalizability and Robustness

abstract: Seismic inversion plays a very useful role in detailed stratigraphic interpretation of seismic data. Seismic inversion enables estimation of rock properties over the complete seismic section. Traditional and machine learning-based seismic inversion workflows are limited to inverting each seismic trace independently of other traces to estimate impedance profiles, leading to lateral discontinuities in the presence of noise and large geological variations in the seismic data. In addition, machine learning-based approaches suffer the problem of overfitting if there is a small number of wells on which the model is trained. We propose a two-pronged strategy to overcome these problems. We present a Temporal Convolutional Network that models seismic traces temporally. We further inject spatial context for each trace into its estimations of the impedance profile. To counter the problem of limited labeled data, we also present a joint learning scheme whereby multiple datasets are simultaneously used for training, sharing beneficial information among each other. This results in the improvement in generalization performance on all datasets. We present a case study of acoustic impedance inversion using the open source SEAM and Marmousi 2 datasets. Our evaluations show that our proposed approach is able to perform robust and laterally consistent estimations of impedance in the presence of noisy seismic data and limited labels. We compare and contrast our approach to other learning-based seismic inversion methodologies in the literature. On SEAM, we are able to obtain an average MSE of 0.0966, the lowest among all other methodologies.


title: An ideal rapid-cycle Thouless pump

abstract: Thouless pumping is a fundamental instance of quantized transport, which is topologically protected. Although its theoretical importance, the adiabaticity condition is an obstacle for further practical applications. Here, focusing on the Rice-Mele model, we provide a family of finite-frequency examples that ensure both the absence of excitations and the perfect quantization of the pumped charge at the end of each cycle. This family, which contains an adiabatic protocol as a limiting case, is obtained through a mapping onto the zero curvature representation of the Euclidean sinh-Gordon equation.


title: Multiscale Governance

abstract: Future societal systems will be characterized by heterogeneous human behaviors and also collective action. The interaction between local systems and global systems will be complex. Humandemics will propagate because of the pathways that connect the different systems and several invariant behaviors and patterns that have emerged globally. On the contrary, infodemics of misinformation can be a risk as it has occurred in the COVID-19 pandemic. The emerging fragility or robustness of the system will depend on how this complex network of systems is governed. Future societal systems will not be only multiscale in terms of the social dimension, but also in the temporality. Necessary and proper prevention and response systems based on complexity, ethic and multi-scale governance will be required. Real-time response systems are the basis for resilience to be the foundation of robust societies. A top-down approach led by Governmental organs for managing humandemics is not sufficient and may be only effective if policies are very restrictive and their efficacy depends not only in the measures implemented but also on the dynamics of the policies and the population perception and compliance. This top-down approach is even weaker if there is not national and international coordination. Coordinating top-down agencies with bottom-up constructs will be the design principle. Multi-scale governance integrates decision-making processes with signaling, sensing and leadership mechanisms to drive thriving societal systems with real-time sensitivity.


title: Fiscal Stimulus of Last Resort

abstract: I examine global dynamics in a monetary model with overlapping generations of finite-horizon agents and a binding lower bound on nominal interest rates. Debt targeting rules exacerbate the possibility of self-fulfilling liquidity traps, for agents expect austerity following deflationary slumps. Conversely, activist but sustainable fiscal policy regimes - implementing intertemporally balanced tax cuts and/or transfer increases in response to disinflationary trajectories - are capable of escaping liquidity traps and embarking inflation into a globally stable path that converges to the target. Should fiscal stimulus of last resort be overly aggressive, however, spiral dynamics around the liquidity-trap steady state exist, causing global indeterminacy.


title: Machine Learning-Driven Virtual Bidding with Electricity Market Efficiency Analysis

abstract: This paper develops a machine learning-driven portfolio optimization framework for virtual bidding in electricity markets considering both risk constraint and price sensitivity. The algorithmic trading strategy is developed from the perspective of a proprietary trading firm to maximize profit. A recurrent neural network-based Locational Marginal Price (LMP) spread forecast model is developed by leveraging the inter-hour dependencies of the market clearing algorithm. The LMP spread sensitivity with respect to net virtual bids is modeled as a monotonic function with the proposed constrained gradient boosting tree. We leverage the proposed algorithmic virtual bid trading strategy to evaluate both the profitability of the virtual bid portfolio and the efficiency of U.S. wholesale electricity markets. The comprehensive empirical analysis on PJM, ISO-NE, and CAISO indicates that the proposed virtual bid portfolio optimization strategy considering the price sensitivity explicitly outperforms the one that neglects the price sensitivity. The Sharpe ratio of virtual bid portfolios for all three electricity markets are much higher than that of the S&P 500 index. It was also shown that the efficiency of CAISO's two-settlement system is lower than that of PJM and ISO-NE.


title: A Quadratic Programming Approach to Manipulation in Real-Time Using Modular Robots

abstract: Motion planning in high-dimensional space is a challenging task. In order to perform dexterous manipulation in an unstructured environment, a robot with many degrees of freedom is usually necessary, which also complicates its motion planning problem. Real-time control brings about more difficulties in which robots have to maintain the stability while moving towards the target. Redundant systems are common in modular robots that consist of multiple modules and are able to transformed into different configurations with respect to different needs. Different from robots with fixed geometry or configurations, the kinematics model of a modular robotic system can alter as the robot reconfigures itself, and developing a generic control and motion planning approach for such systems is difficult, especially when multiple motion goals are coupled. A new manipulation planning framework is developed in this paper. The problem is formulated as a sequential linearly constrained quadratic program (QP) that can be solved efficiently. Some constraints can be incorporated into this QP, including a novel way to approximate environment obstacles. This solution can be used directly for real-time applications or as an off-line planning tool, and it is validated and demonstrated on the CKBot and SMORES-EP modular robot platforms.


title: Efficient transfer learning for NLP with ELECTRA

abstract: Clark et al. [2020] claims that the ELECTRA approach is highly efficient in NLP performances relative to computation budget. As such, this reproducibility study focus on this claim, summarized by the following question: Can we use ELECTRA to achieve close to SOTA performances for NLP in low-resource settings, in term of compute cost?


title: Exploring Targeted Universal Adversarial Perturbations to End-to-end ASR Models

abstract: Although end-to-end automatic speech recognition (e2e ASR) models are widely deployed in many applications, there have been very few studies to understand models' robustness against adversarial perturbations. In this paper, we explore whether a targeted universal perturbation vector exists for e2e ASR models. Our goal is to find perturbations that can mislead the models to predict the given targeted transcript such as "thank you" or empty string on any input utterance. We study two different attacks, namely additive and prepending perturbations, and their performances on the state-of-the-art LAS, CTC and RNN-T models. We find that LAS is the most vulnerable to perturbations among the three models. RNN-T is more robust against additive perturbations, especially on long utterances. And CTC is robust against both additive and prepending perturbations. To attack RNN-T, we find prepending perturbation is more effective than the additive perturbation, and can mislead the models to predict the same short target on utterances of arbitrary length.


title: Hydrodynamic interactions and extreme particle clustering in turbulence

abstract: From new detailed experimental data, we found that the Radial Distribution Function (RDF) of inertial particles in turbulence grows explosively with $r^{-6}$ scaling as the collision radius is approached. We corrected a theory by Yavuz et al. (Phys. Rev. Lett. 120, 244504 (2018)) based on hydrodynamic interactions between pairs of weakly inertial particles, and demonstrate that even this corrected theory cannot explain the observed RDF behavior. We explore several alternative mechanisms for the discrepancy that were not included in the theory and show that none of them are likely the explanation, suggesting new, yet to be identified physical mechanisms are at play.


title: Lower Bounds Implementing Mediators in Asynchronous Systems

abstract: Abraham, Dolev, Geffner, and Halpern proved that, in asynchronous systems, a $(k,t)$-robust equilibrium for $n$ players and a trusted mediator can be implemented without the mediator as long as $n > 4(k+t)$, where an equilibrium is $(k,t)$-robust if, roughly speaking, no coalition of $t$ players can decrease the payoff of any of the other players, and no coalition of $k$ players can increase their payoff by deviating. We prove that this bound is tight, in the sense that if $n \le 4(k+t)$ there exist $(k,t)$-robust equilibria with a mediator that cannot be implemented by the players alone. Even though implementing $(k,t)$-robust mediators seems closely related to implementing asynchronous multiparty $(k+t)$-secure computation \cite{BCG93}, to the best of our knowledge there is no known straightforward reduction from one problem to another. Nevertheless, we show that there is a non-trivial reduction from a slightly weaker notion of $(k+t)$-secure computation, which we call $(k+t)$-strict secure computation, to implementing $(k,t)$-robust mediators. We prove the desired lower bound by showing that there are functions on $n$ variables that cannot be $(k+t)$-strictly securely computed if $n \le 4(k+t)$. This also provides a simple alternative proof for the well-known lower bound of $4t+1$ on asynchronous secure computation in the presence of up to $t$ malicious agents.


title: Generalized pentagonal geometries

abstract: A pentagonal geometry PENT($k$, $r$) is a partial linear space, where every line is incident with $k$ points, every point is incident with $r$ lines, and for each point $x$, there is a line incident with precisely those points that are not collinear with $x$. Here we generalize the concept by allowing the points not collinear with $x$ to form the point set of a Steiner system $S(2,k,w)$ whose blocks are lines of the geometry.


title: Lidar-Monocular Surface Reconstruction Using Line Segments

abstract: Structure from Motion (SfM) often fails to estimate accurate poses in environments that lack suitable visual features. In such cases, the quality of the final 3D mesh, which is contingent on the accuracy of those estimates, is reduced. One way to overcome this problem is to combine data from a monocular camera with that of a LIDAR. This allows fine details and texture to be captured while still accurately representing featureless subjects. However, fusing these two sensor modalities is challenging due to their fundamentally different characteristics. Rather than directly fusing image features and LIDAR points, we propose to leverage common geometric features that are detected in both the LIDAR scans and image data, allowing data from the two sensors to be processed in a higher-level space. In particular, we propose to find correspondences between 3D lines extracted from LIDAR scans and 2D lines detected in images before performing a bundle adjustment to refine poses. We also exploit the detected and optimized line segments to improve the quality of the final mesh. We test our approach on the recently published dataset, Newer College Dataset. We compare the accuracy and the completeness of the 3D mesh to a ground truth obtained with a survey-grade 3D scanner. We show that our method delivers results that are comparable to a state-of-the-art LIDAR survey while not requiring highly accurate ground truth pose estimates.


title: Efficient slave-boson approach for multiorbital two-particle response functions and superconductivity

abstract: We develop an efficient approach for computing two-particle response functions and interaction vertices for multiorbital strongly correlated systems based on fluctuation around rotationally-invariant slave-boson saddle-point. The method is applied to the degenerate three-orbital Hubbard-Kanamori model for investigating the origin of the s-wave orbital antisymmetric spin-triplet superconductivity in the Hund's metal regime, previously found in the dynamical mean-field theory studies. By computing the pairing interaction considering the particle-particle and the particle-hole scattering channels, we identify the mechanism leading to the pairing instability around Hund's metal crossover arises from the particle-particle channel, containing the local electron pair fluctuation between different particle-number sectors of the atomic Hilbert space. On the other hand, the particle-hole spin fluctuations induce the s-wave pairing instability before entering the Hund's regime. Our approach paves the way for investigating the pairing mechanism in realistic correlated materials.


title: Reduced etch lag and high aspect ratios by deep reactive ion etching (DRIE)

abstract: Deep reactive ion etching (DRIE) with the Bosch process is one of the key procedures used to manufacture micron-sized structures for MEMS and microfluidic applications in silicon and, hence, of increasing importance for miniaturization in biomedical research. While guaranteeing high aspect ratio structures and providing high design flexibility, the etching procedure suffers from reactive ion etching lag and often relies on complex oxide masks to enable deep etching. In this work, we introduce an optimized Bosch process that reduces the etch lag to below 1.5 %. Furthermore, we improved a three-step Bosch process, allowing the fabrication of structures with 6 micrometre thickness at depths up to 180 micrometre while maintaining their stability.


title: PageRank centrality and algorithms for weighted, directed networks with applications to World Input-Output Tables

abstract: PageRank (PR) is a fundamental tool for assessing the relative importance of the nodes in a network. In this paper, we propose a measure, weighted PageRank (WPR), extended from the classical PR for weighted, directed networks with possible non-uniform node-specific information that is dependent or independent of network structure. A tuning parameter leveraging node degree and strength is introduced. An efficient algorithm based on R program has been developed for computing WPR in large-scale networks. We have tested the proposed WPR on widely used simulated network models, and found it outperformed other competing measures in the literature. By applying the proposed WPR to the real network data generated from World Input-Output Tables, we have seen the results that are consistent with the global economic trends, which renders it a preferred measure in the analysis.


title: Selective game version of q-points

abstract: This work presents the selection principle $S_1^*(\tau_x,CD)$ that characterizes $q$-points. We also discuss the induced topological game $G_1^*(\tau_x,CD)$ and its relations with $W$-points and $\widetilde{W}$-points, as well as with the game $G_1(\Omega_x,\Omega_x)$.


title: Exactly solvable models for 2+1D topological phases derived from crossed modules of semisimple Hopf algebras

abstract: We define an exactly solvable model for 2+1D topological phases of matter on a triangulated surface derived from a crossed module of semisimple finite-dimensional Hopf algebras, the `Hopf-algebraic higher Kitaev model'. This model generalizes both the Kitaev quantum double model for a semisimple Hopf algebra and the full higher Kitaev model derived from a 2-group, and can hence be interpreted as a Hopf-algebraic discrete higher gauge theory. We construct a family of crossed modules of semisimple Hopf algebras, $(\mathscr{F}_{\mathbb{C}}(X) \otimes \mathbb{C}E \xrightarrow{\partial} \mathscr{F}_{\mathbb{C}}(Y) \rtimes \mathbb{C} G, \triangleright)$, that depends on four finite groups, $E,G,X$ and $Y$. We calculate the ground-state spaces of the resulting model on a triangulated surface when $G=E=\{1\}$ and when $Y=\{1\}$, prove that those ground-state spaces are canonically independent of the triangulations, and so depend only on the underlying surface; and moreover we find a 2+1D TQFT whose state spaces on surfaces give the ground-state spaces. These TQFTs are particular cases of Quinn's finite total homotopy TQFT and hence the state spaces assigned to surfaces are free vector spaces on sets of homotopy classes of maps from a surface to homotopy finite spaces, in this case obtained as classifying spaces of finite groupoids and finite crossed modules of groupoids. We leave it as an open problem whether the ground-state space of the Hopf-algebraic higher Kitaev model on a triangulated surface is independent of the triangulation for general crossed modules of semisimple Hopf algebras, whether a TQFT always exists whose state space on a surface gives the ground-state space of the model, and whether the ground-state space of the model obtained from $E,G,X,Y$ can always be given a homotopical explanation.


title: Profunctors between posets and Alexander duality

abstract: We consider profunctors $f : P \promap Q$ between posets and introduce their {\em graph} and {\em ascent}. The profunctors $\Pro(P,Q)$ form themselves a poset, and we consider a partition $\cI \sqcup \cF$ of this into a down-set $\cI$ and up-set $\cF$, called a {\it cut}. To elements of $\cF$ we associate their graphs, and to elements of $\cI$ we associate their ascents. Our basic result is that this, suitable refined, preserves being a cut: We get a cut in the Boolean lattice of subsets of the underlying set of $Q \times P$. Cuts in finite Booleans lattices correspond precisely to finite simplicial complexes. We apply this in commutative algebra where these give classes of Alexander dual square-free monomial ideals giving the full and natural generalized setting of isotonian ideals and letterplace ideals for posets. We study $\Pro(\NN, \NN)$. Such profunctors identify as order preserving maps $f : \NN \pil \NN \cup \{\infty \}$. For our applications when $P$ and $Q$ are infinite, we also introduce a topology on $\Pro(P,Q)$, in particular on profunctors $\Pro(\NN,\NN)$.


title: Robust Semantic Interpretability: Revisiting Concept Activation Vectors

abstract: Interpretability methods for image classification assess model trustworthiness by attempting to expose whether the model is systematically biased or attending to the same cues as a human would. Saliency methods for feature attribution dominate the interpretability literature, but these methods do not address semantic concepts such as the textures, colors, or genders of objects within an image. Our proposed Robust Concept Activation Vectors (RCAV) quantifies the effects of semantic concepts on individual model predictions and on model behavior as a whole. RCAV calculates a concept gradient and takes a gradient ascent step to assess model sensitivity to the given concept. By generalizing previous work on concept activation vectors to account for model non-linearity, and by introducing stricter hypothesis testing, we show that RCAV yields interpretations which are both more accurate at the image level and robust at the dataset level. RCAV, like saliency methods, supports the interpretation of individual predictions. To evaluate the practical use of interpretability methods as debugging tools, and the scientific use of interpretability methods for identifying inductive biases (e.g. texture over shape), we construct two datasets and accompanying metrics for realistic benchmarking of semantic interpretability methods. Our benchmarks expose the importance of counterfactual augmentation and negative controls for quantifying the practical usability of interpretability methods.


title: Variable selection with missing data in both covariates and outcomes: Imputation and machine learning

abstract: The missing data issue is ubiquitous in health studies. Variable selection in the presence of both missing covariates and outcomes is an important statistical research topic but has been less studied. Existing literature focuses on parametric regression techniques that provide direct parameter estimates of the regression model. Flexible nonparametric machine learning methods considerably mitigate the reliance on the parametric assumptions, but do not provide as naturally defined variable importance measure as the covariate effect native to parametric models. We investigate a general variable selection approach when both the covariates and outcomes can be missing at random and have general missing data patterns. This approach exploits the flexibility of machine learning modeling techniques and bootstrap imputation, which is amenable to nonparametric methods in which the covariate effects are not directly available. We conduct expansive simulations investigating the practical operating characteristics of the proposed variable selection approach, when combined with four tree-based machine learning methods, XGBoost, Random Forests, Bayesian Additive Regression Trees (BART) and Conditional Random Forests, and two commonly used parametric methods, lasso and backward stepwise selection. Numeric results suggest that when combined with bootstrap imputation, XGBoost and BART have the overall best variable selection performance with respect to the $F_1$ score and Type I error across various settings. In general, there is no significant difference in the variable selection performance due to imputation methods. We further demonstrate the methods via a case study of risk factors for 3-year incidence of metabolic syndrome with data from the Study of Women's Health Across the Nation.


title: 3d conformal fields with manifest $sl(2,\mathbb{C})$

abstract: In the present paper we construct all short representation of $so(3,2)$ with the $sl(2,\mathbb{C})$ symmetry made manifest due to the use of $sl(2,\mathbb{C})$ spinors. This construction has a natural connection to the spinor-helicity formalism for massless fields in AdS${}_4$ suggested earlier. We then study unitarity of the resulting representations, identify them as the lowest-weight modules and as conformal fields in the three-dimensional Minkowski space. Finally, we compare these results with the existing literature and discuss the properties of these representations under contraction of $so(3,2)$ to the Poincare algebra.