title: High-Speed Imaging Receiver Design for 6G Optical Wireless Communications: A Rate-FOV Trade-Off

abstract: The design of a compact high-speed and wide field of view (FOV) receiver is challenging due to the presence of two well-known trade-offs. The first one is the area-bandwidth trade-off of photodetectors (PDs) and the second one is the gain-FOV trade-off due to the use of optics. The combined effects of these two trade-offs imply that the achievable data rate of an imaging optical receiver is limited by its FOV, i.e., a rate-FOV trade-off. To control the area-bandwidth trade-off, an array of small PDs can be used instead of a single PD. Moreover, in practice, a large-area lens is required to ensure sufficient power collection, which in turn limits the receiver FOV (i.e., gain-FOV trade-off). We propose an imaging receiver design in the form of an array of arrays. To achieve a reasonable receiver FOV, we use individual focusing lens for each PD array rather than a single collection lens for the whole receiver. The proposed array of arrays structure provides an effective method to control both gain-FOV trade-off (via an array of lenses) and area-bandwidth trade-off (via arrays of PDs). We first derive a tractable analytical model for the SNR of an array of PDs where the maximum ratio combining has been employed. Then, we extend the model for the proposed array of arrays structure and the accuracy of the analytical model is verified based on several Optic Studio-based simulations. Next, we formulate an optimization problem to maximize the achievable data rate of the imaging receiver subject to a minimum required FOV. The optimization problem is solved for two commonly used modulation techniques, namely, OOK and direct current biased optical orthogonal frequency division multiplexing with variable rate quadrature amplitude modulation. It is demonstrated that a data rate of ~ 24 Gbps with a FOV of 15 is achievable using OOK with a total receiver size of 2 cm by 2 cm.


title: On Upward-Planar L-Drawings of Graphs

abstract: In an upward-planar L-drawing of a directed acyclic graph (DAG) each edge $e$ is represented as a polyline composed of a vertical segment with its lowest endpoint at the tail of $e$ and of a horizontal segment ending at the head of $e$. Distinct edges may overlap, but not cross. Recently, upward-planar L-drawings have been studied for $st$-graphs, i.e., planar DAGs with a single source $s$ and a single sink $t$ containing an edge directed from $s$ to $t$. It is known that a plane $st$-graph, i.e., an embedded $st$-graph in which the edge $(s,t)$ is incident to the outer face, admits an upward-planar L-drawing if and only if it admits a bitonic $st$-ordering, which can be tested in linear time. We study upward-planar L-drawings of DAGs that are not necessarily $st$-graphs. On the combinatorial side, we show that a plane DAG admits an upward-planar L-drawing if and only if it is a subgraph of a plane $st$-graph admitting a bitonic $st$-ordering. This allows us to show that not every tree with a fixed bimodal embedding admits an upward-planar L-drawing. Moreover, we prove that any acyclic cactus with a single source (or a single sink) admits an upward-planar L-drawing, which respects a given outerplanar embedding if there are no transitive edges. On the algorithmic side, we consider DAGs with a single source (or a single sink). We give linear-time testing algorithms for these DAGs in two cases: (i) when the drawing must respect a prescribed embedding and (ii) when no restriction is given on the embedding, but it is biconnected and series-parallel.


title: Extensible Machine Learning for Encrypted Network Traffic Application Labeling via Uncertainty Quantification

abstract: With the increasing prevalence of encrypted network traffic, cyber security analysts have been turning to machine learning (ML) techniques to elucidate the traffic on their networks. However, ML models can become stale as new traffic emerges that is outside of the distribution of the training set. In order to reliably adapt in this dynamic environment, ML models must additionally provide contextualized uncertainty quantification to their predictions, which has received little attention in the cyber security domain. Uncertainty quantification is necessary both to signal when the model is uncertain about which class to choose in its label assignment and when the traffic is not likely to belong to any pre-trained classes. We present a new, public dataset of network traffic that includes labeled, Virtual Private Network (VPN)-encrypted network traffic generated by 10 applications and corresponding to 5 application categories. We also present an ML framework that is designed to rapidly train with modest data requirements and provide both calibrated, predictive probabilities as well as an interpretable "out-of-distribution" (OOD) score to flag novel traffic samples. We describe calibrating OOD scores using p-values of the relative Mahalanobis distance. We demonstrate that our framework achieves an F1 score of 0.98 on our dataset and that it can extend to an enterprise network by testing the model: (1) on data from similar applications, (2) on dissimilar application traffic from an existing category, and (3) on application traffic from a new category. The model correctly flags uncertain traffic and, upon retraining, accurately incorporates the new data.


title: Tau depolarization at very high energies for neutrino telescopes

abstract: The neutrino interaction length scales with energy, and becomes comparable to Earth's diameter above 10's of TeV energies. Over terrestrial distances, the tau's short lifetime leads to an energetic regenerated tau neutrino flux, tau neutrino to tau to tau neutrino, within the Earth. The next generation of neutrino experiments aim to detect ultra-high energy neutrinos. Many of them rely on detecting either the regenerated tau neutrino, or a tau decay shower. Both of these signatures can be affected by the polarization of the tau through the energy distribution of the secondary particles produced from the tau's decay. While taus produced in weak interactions are nearly 100 percent polarized, it is expected that taus experience some depolarization due to electromagnetic interactions in the Earth. In this paper, for the first time we quantify the depolarization of taus in electromagnetic energy loss. We find that tau depolarization has only small effects on the final energy of tau neutrinos or taus produced by high energy tau neutrinos incident on the Earth. Tau depolarization can be directly implemented in Monte Carlo simulations such as nuPyProp and TauRunner.


title: Benefits of Feedforward for Model Predictive Airpath Control of Diesel Engines

abstract: This paper investigates options to complement a diesel engine airpath feedback controller with a feedforward. The control objective is to track the intake manifold pressure and exhaust gas recirculation (EGR) rate targets by manipulating the EGR valve and variable geometry turbine (VGT) while satisfying state and input constraints. The feedback controller is based on rate-based Model Predictive Control (MPC) that provides integral action for tracking. Two options for the feedforward are considered one based on a look-up table that specifies the feedforward as a function of engine speed and fuel injection rate, and another one based on a (non-rate-based) MPC that generates dynamic feedforward trajectories. The controllers are designed and verified using a high-fidelity engine model in GT-Power and exploit a low-order rate-based linear parameter-varying (LPV) model for prediction which is identified from transient response data generated by the GT-Power model. It is shown that the combination of feedforward and feedback MPC has the potential to improve the performance and robustness of the control design. In particular, the feedback MPC without feedforward can lose stability at low engine speeds, while MPC-based feedforward results in the best transient response. Mechanisms by which feedforward is able to assist in stabilization and improve performance are discussed.


title: Second-Order Asymptotics of Hoeffding-Like Hypothesis Tests

abstract: We consider a binary statistical hypothesis testing problem, where $n$ independent and identically distributed random variables $Z^n$ are either distributed according to the null hypothesis $P$ or the alternate hypothesis $Q$, and only $P$ is known. For this problem, a well-known test is the Hoeffding test, which accepts $P$ if the Kullback-Leibler (KL) divergence between the empirical distribution of $Z^n$ and $P$ is below some threshold. In this paper, we consider Hoeffding-like tests, where the KL divergence is replaced by other divergences, and characterize, for a large class of divergences, the first and second-order terms of the type-II error for a fixed type-I error. Since the considered class includes the KL divergence, we obtain the second-order term of the Hoeffiding test as a special case.


title: On Distributed Adaptive Optimization with Gradient Compression

abstract: We study COMP-AMS, a distributed optimization framework based on gradient averaging and adaptive AMSGrad algorithm. Gradient compression with error feedback is applied to reduce the communication cost in the gradient transmission process. Our convergence analysis of COMP-AMS shows that such compressed gradient averaging strategy yields same convergence rate as standard AMSGrad, and also exhibits the linear speedup effect w.r.t. the number of local workers. Compared with recently proposed protocols on distributed adaptive methods, COMP-AMS is simple and convenient. Numerical experiments are conducted to justify the theoretical findings, and demonstrate that the proposed method can achieve same test accuracy as the full-gradient AMSGrad with substantial communication savings. With its simplicity and efficiency, COMP-AMS can serve as a useful distributed training framework for adaptive gradient methods.


title: Improved decision making with similarity based machine learning: Applications in chemistry

abstract: Despite the fundamental progress in autonomous molecular and materials discovery, data scarcity throughout chemical compound space still severely hampers the use of modern ready-made machine learning models as they rely heavily on the paradigm, 'the bigger the data the better'. Presenting similarity based machine learning (SML), we show an approach to select data and train a model on-the-fly for specific queries, enabling decision making in data scarce scenarios in chemistry. By solely relying on query and training data proximity to choose training points, only a fraction of data is necessary to converge to competitive performance. After introducing SML for the harmonic oscillator and the Rosenbrock function, we describe applications to scarce data scenarios in chemistry which include quantum mechanics based molecular design and organic synthesis planning. Finally, we derive a relationship between the intrinsic dimensionality and volume of feature space, governing the overall model accuracy.


title: Evolution of viral pathogens follows a linear order

abstract: Although lessons have been learned from previous severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) outbreaks, the rapid evolution of the viruses means that future outbreaks of a much larger scale are possible, as shown by the current coronavirus disease 2019 (COVID-19) outbreak. Therefore, it is necessary to better understand the evolution of coronaviruses as well as viruses in general. This study reports a comparative analysis of the amino acid usage within several key viral families and genera that are prone to triggering outbreaks, including coronavirus (SARS-CoV-2, SARS-CoV, MERS-CoV, HCoV-HKU1, HCoV-OC43, HCoV-NL63, HCoV-229E), influenza A (H1N1, H3N2), flavivirus (dengue virus serotypes 1-4, Zika) and ebolavirus (Zaire, Sudan, Bundibugyo ebolavirus). Our analysis reveals that the distribution of amino acid usage in the viral genome is constrained to follow a linear order, and the distribution remains closely related to the viral species within the family or genus. This constraint can be adapted to predict viral mutations and future variants of concern. By studying previous SARS and MERS outbreaks, we have adapted this naturally occurring pattern to determine that although pangolin plays a role in the outbreak of COVID-19, it may not be the sole agent as an intermediate animal. In addition to this study, our findings contribute to the understanding of viral mutations for subsequent development of vaccines and toward developing a model to determine the source of the outbreak.


title: On Dependent Dirichlet Processes for General Polish Spaces

abstract: We study Dirichlet process-based models for sets of predictor-dependent probability distributions, where the domain and predictor space are general Polish spaces. We generalize the definition of dependent Dirichlet processes, originally constructed on Euclidean spaces, to more general Polish spaces. We provide sufficient conditions under which dependent Dirichlet processes have appealing properties regarding continuity (weak and strong), association structure, and support (under different topologies). We also provide sufficient conditions under which mixture models induced by dependent Dirichlet processes have appealing properties regarding strong continuity, association structure, support, and weak consistency under i.i.d. sampling of both responses and predictors. The results can be easily extended to more general dependent stick-breaking processes.


title: Do you smell something decaying? Updated linear constraints on decaying dark matter scenarios

abstract: The stability of particles in the cosmic soup is an important property that can affect the cosmic evolution. In this work, we update the constraints on the decaying cold dark matter scenario, when the decay products are effectively massless. We assume, as a base case, that all of dark matter is unstable and it can decay on cosmological time scales. We then extend the analysis to include the scenario where only a fraction of dark matter is unstable, while the remaining part is composed of the standard, stable, dark matter. We consider observations of cosmological probes at linear scales, i.e., Planck 2018 cosmic microwave background temperature, polarization, and lensing measurements, along with geometrical information from baryon acoustic oscillation (BAO) measurements from SDSS DR7, BOSS DR12, eBOSS DR16 and 6dFGS, to derive conservative constraints on the dark matter decay rate. We consider these dataset separately, to asses the relative constraining power of each dataset, as well as together to asses the joint constraints. We find the most stringent upper limit on the decay rate of decaying cold dark matter particles to be $\Gammadm < 0.129 \times 10^{-18}\,\mathrm{s}^{-1}$ (or, equivalently, the dark matter lifetime $\taudm > 246$ Gyr) at 95\% C.L. for the combination of Planck primary anisotropies, lensing and BAO. We further explore one-parameter extensions of our baseline DCDM model. Namely, we vary the sum of neutrino masses, the curvature density parameter, and the tensor-to-scalar ratio along with the DCDM parameters. When varying the tensor-to-scalar ratio we also add data from the BICEP/Keck experiment.


title: Variance of the Hellings-Downs Correlation

abstract: Gravitational waves (GWs) create correlations in the arrival times of pulses from different pulsars. The expected correlation $\mu(\gamma)$ as a function of the angle $\gamma$ between the directions to two pulsars was calculated by Hellings and Downs for an isotropic and unpolarized GW background, and several pulsar timing array (PTA) collaborations are working to observe these. We ask: given a set of noise-free observations, are they consistent with that expectation? To answer this, we calculate the expected variance $\sigma^2(\gamma)$ in the correlation for a single GW point source, as pulsar pairs with fixed separation angle $\gamma$ are swept around the sky. We then use this to derive simple analytic expressions for the variance produced by a set of discrete point sources uniformly scattered in space for two cases of interest: (1) point sources radiating GWs at the same frequency, generating confusion noise, and (2) point sources radiating GWs at distinct non-overlapping frequencies. By averaging over all pulsar sky positions at fixed separation angle $\gamma$, we show how this variance may be cleanly split into cosmic variance and pulsar variance, also demonstrating that measurements of the variance can provide information about the nature of GW sources. In a series of technical appendices, we calculate the mean and variance of the Hellings-Downs correlation for an arbitrary (polarized) point source, quantify the impact of neglecting pulsar terms, and calculate the pulsar and cosmic variance for a Gaussian ensemble. The mean and variance of the Gaussian ensemble may also be obtained from the previous discrete-source confusion-noise model in the limit of a high density of weak sources.


title: Few-Shot Parameter-Efficient Fine-Tuning is Better and Cheaper than In-Context Learning

abstract: Few-shot in-context learning (ICL) enables pre-trained language models to perform a previously-unseen task without any gradient-based training by feeding a small number of training examples as part of the input. ICL incurs substantial computational, memory, and storage costs because it involves processing all of the training examples every time a prediction is made. Parameter-efficient fine-tuning (PEFT) (e.g. adapter modules, prompt tuning, sparse update methods, etc.) offers an alternative paradigm where a small set of parameters are trained to enable a model to perform the new task. In this paper, we rigorously compare few-shot ICL and PEFT and demonstrate that the latter offers better accuracy as well as dramatically lower computational costs. Along the way, we introduce a new PEFT method called (IA)$^3$ that scales activations by learned vectors, attaining stronger performance while only introducing a relatively tiny amount of new parameters. We also propose a simple recipe based on the T0 model called T-Few that can be applied to new tasks without task-specific tuning or modifications. We validate the effectiveness of T-Few on completely unseen tasks by applying it to the RAFT benchmark, attaining super-human performance for the first time and outperforming the state-of-the-art by 6% absolute. All of the code used in our experiments is publicly available.


title: Degenerate Distributed Feedback Photonic Structure with Double Grating Exhibiting Degenerate Band Edge

abstract: We propose a degenerate version of the Bragg condition, associated to a degenerate band edge (DBE). A standard Bragg condition can be implemented using a periodic grating that operates at the regular band edge. The structure we propose to realize a DBE is made of two stacked identical gratings that form two coupled periodic waveguides with broken mirror symmetry. The occurrence of the DBE is verified by using both an eigenmode solver and the calculations of scattering parameters, using full-wave simulations. The proposed structure is a good candidate to conceive a degenerate distributed feedback lasers operating at the DBE.


title: Parametrization and Estimation of High-Rank Line-of-Sight MIMO Channels with Reflected Paths

abstract: High-rank line-of-sight (LOS) MIMO systems have attracted considerable attention for millimeter wave and THz communications. The small wavelengths in these frequencies enable spatial multiplexing with massive data rates at long distances. Such systems are also being considered for multi-path non-LOS (NLOS) environments. In these scenarios, standard channel models based on plane waves cannot capture the curvature of each wave front necessary to model spatial multiplexing. This work presents a novel and simple multi-path wireless channel parametrization where each path is replaced by a LOS path with a reflected image source. The model fully is valid for all paths with specular planar reflections and captures the spherical nature of each wave front. Importantly, it is shown that the model uses only two additional parameters relative to the standard plane wave model. Moreover, the parameters can be easily captured in standard ray tracing. The accuracy of the approach is demonstrated on detailed ray tracing simulations at 28 GHz and 140 GHz in a dense urban area.


title: Improved entanglement indicators for optical fields and its application in the event-ready experiment for bright squeezed vacuum with induced non-gaussianity

abstract: Better versions of separability conditions for four mode optical fields, i.e. two beams with two modes of mutually orthogonal polarization are given. Our conditions involve variances. Their meaning is intuitive and their implementation is feasible. Namely, if for a given quantum state the spread of the data around its mean value is smaller than the minimal spread predicted for the set of separable states, then the given state is entangled. Our conditions are formulated for standard quantum Stokes observables and normalized Stokes observables and result to be more efficient that the previous conditions for four mode optical fields involving variances. We test our conditions for bright squeezed vacuum with (and without) induced non-gaussianity obtained by addition or subtraction of photons. Also we propose a practical experimental scheme of how to generate such states for an event-ready experiment.


title: Linking pressure gradients with the stability of molecular clouds in galactic outflows

abstract: The jets launched by actively accreting black holes are capable of launching several of the massive (million or billion solar mass) molecular outflows observed in galaxies. These outflows could suppress or enhance star formation in galaxies. To investigate the stability of clouds capable to form stars in outflows, we modeled CO and HCO+ ALMA data of the galaxy IC5063, in which black-hole jets impact molecular clouds. Using a radiative transfer code that self-consistently performs astrochemical and thermal balance calculations based on the available gas heating sources, we found that mechanical heating and cosmic ray (CR) heating are fully capable of individually reproducing the data. In our best-fit model, CRs provide about 1/3rd of the dense gas heating at the radio lobes, emphasizing the role of this often neglected mechanism in heating the gas and potentially generating outflows. The gas temperature and density indicate that the jet passage leads to an increase of about 1 order of magnitude in the internal pressure Pi of molecular clouds (with Pi/k from 8*10^5 up to 7*10^6 K cm^-3), irrespective of the excitation mechanism. From the fluxes of [S II] and [N II] lines in VLT MUSE data, the external pressure Pe of molecular clouds increases in several regions enough to exceed Pi. This result leads us to conclude that we are observing the expansion of an ionized overpressurized cocoon that compresses molecular clouds and that could lead to their collapse. Some jet-impacted clouds, nonetheless, near pathways that the jet cleared have increased Pi and decreased Pe. They are likely to undergo evaporation of their outer layers. Part of the evaporated layers could mass load the outflow thanks to ram pressure from co-spatial ionized gas flows. The observed pressure changes thus suggest that both star formation enhancement and suppression could simultaneously occur.


title: A New Class of String Transformations for Compressed Text Indexing

abstract: Introduced about thirty years ago in the field of Data Compression, the Burrows-Wheeler Transform (BWT) is a string transformation that, besides being a booster of the performance of memoryless compressors, plays a fundamental role in the design of efficient self-indexing compressed data structures. Finding other string transformations with the same remarkable properties of BWT has been a challenge for many researchers for a long time. Among the known BWT variants, the only one that has been recently shown to be a valid alternative to BWT is the Alternating BWT (ABWT), another invertible string transformation introduced about ten years ago in connection with a generalization of Lyndon words. In this paper, we introduce a whole class of new string transformations, called local orderings-based transformations, which have all the myriad virtues of BWT. We show that this new family is a special case of a much larger class of transformations, based on context adaptive alphabet orderings, that includes BWT and ABWT. Although all transformations support pattern search, we show that, in the general case, the transformations within our larger class may take quadratic time for inversion and pattern search. As a further result, we show that the local orderings-based transformations can be used for the construction of the recently introduced r-index, which makes them suitable also for highly repetitive collections. In this context, we consider the problem of finding, for a given string, the BWT variant that minimizes the number of runs in the transformed string, and we provide an algorithm solving this problem in linear time.


title: Enhancing gravitational wave anisotropies with peaked scalar sources

abstract: Gravitational wave (GW) backgrounds of cosmological origin are expected to be nearly isotropic, with small anisotropies resembling those of the cosmic microwave background. We analyse the case of a scalar-induced GW background and clarify in the process the relation between two different approaches to calculating GW anisotropies. We focus on GW scenarios sourced by a significantly peaked scalar spectrum, which are frequently considered in the context of primordial black holes production. We show that the resulting GW anisotropies are characterised by a distinct frequency dependence. We explore the observational consequences concentrating on a GW background enhanced in the frequency band of space-based GW detectors. We study the detectability of the signal through both cross-correlations among different space-based GW detectors, and among GW and CMB experiments.


title: CHES: a space-borne astrometric mission for the detection of habitable planets of the nearby solar-type stars

abstract: The Closeby Habitable Exoplanet Survey (CHES) mission is proposed to discover habitable-zone Earth-like planets of the nearby solar-type stars ($\sim 10~\mathrm{pc}$ away from our solar system) via micro-arcsecond relative astrometry. The major scientific objectives of CHES are: to search for Earth Twins or terrestrial planets in habitable zones orbiting 100 FGK nearby stars; further to conduct a comprehensive survey and extensively characterize the nearby planetary systems. The primary payload is a high-quality, low-distortion, high-stability telescope. The optical subsystem is a coaxial three-mirror anastigmat (TMA) with a $1.2 \mathrm{~m}$-aperture, $0.44^{\circ} \times 0.44^{\circ}$ field of view and $500 \mathrm{~nm}-900 \mathrm{~nm}$ working waveband. The camera focal plane is composed of 81 MOSAIC scientific CMOS detectors each with $4 \mathrm{~K} \times 4 \mathrm{~K}$ pixels. The heterodyne laser interferometric calibration technology is employed to ensure micro-arcsecond level (1 $\mu$as) relative astrometry precision to meet the requirements for detection of Earth-like planets. CHES satellite operates at the Sun-Earth L2 point and observes the entire target stars for 5 years. CHES will offer the first direct measurements of true masses and inclinations of Earth Twins and super-Earths orbiting our neighbor stars based on micro-arcsecond astrometry from space. This will definitely enhance our understanding of the formation of diverse nearby planetary systems and the emergence of other worlds for solar-type stars, and finally to reflect the evolution of our own solar system.

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


title: Aggregating Pairwise Semantic Differences for Few-Shot Claim Veracity Classification

abstract: As part of an automated fact-checking pipeline, the claim veracity classification task consists in determining if a claim is supported by an associated piece of evidence. The complexity of gathering labelled claim-evidence pairs leads to a scarcity of datasets, particularly when dealing with new domains. In this paper, we introduce SEED, a novel vector-based method to few-shot claim veracity classification that aggregates pairwise semantic differences for claim-evidence pairs. We build on the hypothesis that we can simulate class representative vectors that capture average semantic differences for claim-evidence pairs in a class, which can then be used for classification of new instances. We compare the performance of our method with competitive baselines including fine-tuned BERT/RoBERTa models, as well as the state-of-the-art few-shot veracity classification method that leverages language model perplexity. Experiments conducted on the FEVER and SCIFACT datasets show consistent improvements over competitive baselines in few-shot settings. Our code is available.


title: Minimal Representations of Tropical Rational Functions

abstract: This paper studies the following question: given a piecewise-linear function, find its minimal algebraic representation as a tropical rational signomial. We put forward two different notions of minimality, one based on monomial length, the other based on factorization length. We show that in dimension one, both notions coincide, but this is not true in dimensions two or more. We prove uniqueness of the minimal representation for dimension one and certain subclasses of piecewise-linear functions in dimension two. As a proof step, we obtain counting formulas and lower bounds for the number of regions in an arrangement of tropical hypersurfaces, giving a small extension for a result by Mont\'{u}far, Ren and Zhang. As an equivalent formulation, it gives a lower bound on the number of vertices in a regular mixed subdivision of a Minkowski sum, giving a small extension for Adiprasito's Lower Bound Theorem for Minkowski sums.


title: A Computationally Governed Log-domain Interior-point Method for Model Predictive Control

abstract: This paper introduces a computationally efficient approach for solving Model Predictive Control (MPC) reference tracking problems with state and control constraints. The approach consists of three key components: First, a log-domain interior-point quadratic programming method that forms the basis of the overall approach; second, a method of warm-starting this optimizer by using the MPC solution from the previous timestep; and third, a computational governor that bounds the suboptimality of the warm-start by altering the reference command provided to the MPC problem. As a result, the closed-loop system is altered in a manner so that MPC solutions can be computed using fewer optimizer iterations per timestep. In a numerical experiment, the computational governor reduces the worst-case computation time of a standard MPC implementation by 90, while maintaining good closed-loop performance.


title: Any-k Algorithms for Enumerating Ranked Answers to Conjunctive Queries

abstract: We study ranked enumeration for Conjunctive Queries (CQs) where the answers are ordered by a given ranking function (e.g., an ORDER BY clause in SQL). We develop "any-k" algorithms, which, without knowing the number k of desired answers, push down the ranking into joins by carefully ordering the computation of intermediate tuples and avoiding materialization of join answers until they are needed. For this to be possible, the ranking function needs to obey a particular type of monotonicity. Supported ranking functions include the common sum-of-weights, where answers are compared by the sum of input-tuple weights, as well as any commutative selective dioid. Our results extend a well-known unranked-enumeration dichotomy, which states that only free-connex CQs are tractable (under certain hardness hypotheses and for CQs without self-joins). For this class of queries and with n denoting the size of the input, the data complexity of our ranked enumeration approach for the time to the kth CQ answer is O(n+klogk), which is only a logarithmic factor slower than the O(n+k) unranked-enumeration guarantee. A core insight of our work is that ranked enumeration for CQs is closely related to Dynamic Programming and the fundamental task of path enumeration in a weighted DAG. We uncover a previously unknown tradeoff: one any-k algorithm has lower complexity when the number of returned answers is small, the other when their number is large. This tradeoff is eliminated under a stricter monotonicity property that we define and exploit for a novel algorithm that asymptotically dominates all previously known alternatives, including Eppstein's algorithm for sum-of-weights path enumeration. We empirically demonstrate the findings of our theoretical analysis in an experimental study that highlights the superiority of our approach over the join-then-rank approach that existing database systems follow.


title: Invariants of knots and links at roots of unity

abstract: We present a comprehensive classification of invariants of knots and links associated with irreducible representations of \uqslN{}, when the parameter of quantization $q$ is a root of unity. We demonstrate that, besides the standard HOMFLY-PT invariants, which are associated with representations with highest and lowest weights, non-trivial invariants can be associated only with nilpotent representations with parameters. We define the corresponding invariants and discuss their relations with standard invariants at particular values of parameters.


title: Joint OAM Radar-Communication Systems: Target Recognition and Beam Optimization

abstract: Orbital angular momentum (OAM) radars are able to estimate the azimuth angle and the rotation velocity of multiple targets without relative motion or beam scanning. Moreover, OAM wireless communications can achieve high spectral efficiency (SE) by utilizing a set of information-bearing modes on the same frequency channel. Benefitting from the above advantages, in this paper, we design a novel radar-centric joint OAM radar-communication (RadCom) scheme based on uniform circular arrays (UCAs), which modulates information signals on the existing OAM radar waveform. In details, we first propose an OAM-based three-dimensional (3-D) super-resolution position estimation and rotation velocity detection method, which can accurately estimate the 3-D position and rotation velocity of multiple targets. Then, we derive the posterior Cramer-Rao bound (PCRB) of the OAM-based estimates and, finally, we analyze the transmission rate of the integrated communication system. To achieve the best trade-off between imaging and communication, the transmitted integrated OAM beams are optimized by means of an exhaustive search method. Both mathematical analysis and simulation results show that the proposed radar-centric joint OAM RadCom scheme can accurately estimate the 3-D position and rotation velocity of multiple targets while ensuring the transmission rate of the communication receiver, which can be regarded as an effective supplement to existing joint RadCom schemes.


title: Linear instability of turbulent channel flow

abstract: Laminar-turbulent pattern formation is a distinctive feature of the intermittency regime in subcritical plane shear flows. By performing extensive numerical simulations of the plane channel flow, we show that the pattern emerges from a spatial modulation of the turbulent flow, due to a linear instability. We sample over many realizations the linear response of the fluctuating turbulent field to a temporal impulse, in the regime where the turbulent flow is stable, just before the onset of the instability. The dispersion relation is constructed from the ensemble-averaged relaxation rates. As the instability threshold is approached, the relaxation rate of the least damped modes eventually reaches zero. The method allows, despite the presence of turbulent fluctuations and without any closure model, for an accurate estimation of the wavevector of the modulation at onset.


title: The First Optimal Algorithm for Smooth and Strongly-Convex-Strongly-Concave Minimax Optimization

abstract: In this paper, we revisit the smooth and strongly-convex-strongly-concave minimax optimization problem. Zhang et al. (2021) and Ibrahim et al. (2020) established the lower bound $\Omega\left(\sqrt{\kappa_x\kappa_y} \log \frac{1}{\epsilon}\right)$ on the number of gradient evaluations required to find an $\epsilon$-accurate solution, where $\kappa_x$ and $\kappa_y$ are condition numbers for the strong convexity and strong concavity assumptions. However, the existing state-of-the-art methods do not match this lower bound: algorithms of Lin et al. (2020) and Wang and Li (2020) have gradient evaluation complexity $\mathcal{O}\left( \sqrt{\kappa_x\kappa_y}\log^3\frac{1}{\epsilon}\right)$ and $\mathcal{O}\left( \sqrt{\kappa_x\kappa_y}\log^3 (\kappa_x\kappa_y)\log\frac{1}{\epsilon}\right)$, respectively. We fix this fundamental issue by providing the first algorithm with $\mathcal{O}\left(\sqrt{\kappa_x\kappa_y}\log\frac{1}{\epsilon}\right)$ gradient evaluation complexity. We design our algorithm in three steps: (i) we reformulate the original problem as a minimization problem via the pointwise conjugate function; (ii) we apply a specific variant of the proximal point algorithm to the reformulated problem; (iii) we compute the proximal operator inexactly using the optimal algorithm for operator norm reduction in monotone inclusions.


title: Tuning Parameter Selection for Penalized Estimation via $R^2$

abstract: The tuning parameter selection strategy for penalized estimation is crucial to identify a model that is both interpretable and predictive. However, popular strategies (e.g., minimizing average squared prediction error via cross-validation) tend to select models with more predictors than necessary. This paper proposes a simple, yet powerful cross-validation strategy based on maximizing squared correlations between the observed and predicted values, rather than minimizing squared error loss for the purposes of support recovery. The strategy can be applied to all penalized least-squares estimators and we show that, under certain conditions, the metric implicitly performs a bias adjustment. Specific attention is given to the Lasso estimator, in which our strategy is closely related to the Relaxed Lasso estimator. We demonstrate our approach on a functional variable selection problem to identify optimal placement of surface electromyogram sensors to control a robotic hand prosthesis.


title: Probing Top-quark Operators with Precision Electroweak Measurements

abstract: In the Standard Model Effective Field Theory (SMEFT), operators involving the top quark are generally difficult to probe, and can generate sizable loop contributions to the electroweak precision observables, measured by past and future lepton colliders. Could the high precision of the electroweak measurements compensate the loop suppression and provide competitive reaches on these operators? Would the inclusion of these contributions introduce too many additional parameters for a meaningful global electroweak analysis to be done? In this paper, we perform a detailed phenomenological study to address these two important questions. Focusing on eight dimension-6 operators that generate anomalous couplings between the electroweak gauge bosons and the third-generation quarks, we calculate their one loop contributions to the $e^+e^- \to f\bar{f}$ processes both on and off the Z-pole and the $e^-e^+ \to WW$ process. A global analysis is performed with these eight operators and the ones that contribute to the above processes at tree level, using the measurements at LEP, SLC and several low energy experiments. We find that, while the current electroweak precision measurements are sensitive to the one-loop effects of the top-quark operators, it is difficult to separate them from the operators that contribute at the tree level, making a global analysis rather challenging. Under more assumptions, competitive reaches could be obtained in a global fit. Another important finding of our study is that the two operators that generate dipole interactions of the bottom quark have significant impacts in the Z-pole measurements and should not be omitted. We also discuss the implication of the recently reported W-boson mass measurement at CDF to our results. Finally, we estimate the reaches of future lepton colliders in probing the top-quark operators with precision electroweak measurements.


title: Ontology-Driven and Weakly Supervised Rare Disease Identification from Clinical Notes

abstract: Computational text phenotyping is the practice of identifying patients with certain disorders and traits from clinical notes. Rare diseases are challenging to be identified due to few cases available for machine learning and the need for data annotation from domain experts. We propose a method using ontologies and weak supervision, with recent pre-trained contextual representations from Bi-directional Transformers (e.g. BERT). The ontology-based framework includes two steps: (i) Text-to-UMLS, extracting phenotypes by contextually linking mentions to concepts in Unified Medical Language System (UMLS), with a Named Entity Recognition and Linking (NER+L) tool, SemEHR, and weak supervision with customised rules and contextual mention representation; (ii) UMLS-to-ORDO, matching UMLS concepts to rare diseases in Orphanet Rare Disease Ontology (ORDO). The weakly supervised approach is proposed to learn a phenotype confirmation model to improve Text-to-UMLS linking, without annotated data from domain experts. We evaluated the approach on three clinical datasets, MIMIC-III discharge summaries, MIMIC-III radiology reports, and NHS Tayside brain imaging reports from two institutions in the US and the UK, with annotations. The improvements in the precision were pronounced (by over 30% to 50% absolute score for Text-to-UMLS linking), with almost no loss of recall compared to the existing NER+L tool, SemEHR. Results on radiology reports from MIMIC-III and NHS Tayside were consistent with the discharge summaries. The overall pipeline processing clinical notes can extract rare disease cases, mostly uncaptured in structured data (manually assigned ICD codes). We discuss the usefulness of the weak supervision approach and propose directions for future studies.


title: Fibered Universal Algebra for First-Order Logics

abstract: We extend Lawvere-Pitts prop-categories (aka. hyperdoctrines) to develop a general framework for providing "algebraic" semantics for nonclassical first-order logics. This framework includes a natural notion of substitution, which allows first-order logics to be considered as structural closure operators just as propositional logics are in abstract algebraic logic. We then establish an extension of the homomorphism theorem from universal algebra for generalized prop-categories and characterize two natural closure operators on the prop-categorical semantics. The first closes a class of structures (which are interpreted as morphisms of prop-categories) under the satisfaction of their common first-order theory and the second closes a class of prop-categories under their associated first-order consequence. It turns out these closure operators have characterizations that closely mirror Birkhoff's characterization of the closure of a class of algebras under the satisfaction of their common equational theory and Blok and J\'onsson's characterization of closure under equational consequence, respectively. These "algebraic" characterizations of the first-order closure operators are unique to the prop-categorical semantics and do not have analogs, for example, in the Tarskian semantics for classical first-order logic. The prop-categories we consider are much more general than traditional intuitionistic prop-categories or triposes (i.e., topos representing indexed partially ordered sets). Nonetheless, to our knowledge, our results are still new even when restricted to these special classes of prop-categories.


title: Uhlmann holonomy against Lindblad dynamics of topological systems at finite temperatures

abstract: The Uhlmann phase, which reflects the holonomy as the purified state of a density matrix traverses a loop in the parameter space, has been used to characterize topological properties of several systems at finite temperatures. We test the Uhlmann holonomy against various system-environment couplings in quantum dynamics described by the Lindblad equations of three topological systems, including the Su-Schrieffer-Heeger (SSH) model, Kitaev chain, and Bernevig-Hughes-Zhang (BHZ) model. The Uhlmann phase is shown to remain quantized in all the examples if the initial state is topological and only certain types of the Lindblad jump operators are present. Topological protection at finite temperatures against environmental effects in quantum dynamics is therefore demonstrated albeit only for a restricted class of system-environment couplings.


title: Ranked Prioritization of Groups in Combinatorial Bandit Allocation

abstract: Preventing poaching through ranger patrols protects endangered wildlife, directly contributing to the UN Sustainable Development Goal 15 of life on land. Combinatorial bandits have been used to allocate limited patrol resources, but existing approaches overlook the fact that each location is home to multiple species in varying proportions, so a patrol benefits each species to differing degrees. When some species are more vulnerable, we ought to offer more protection to these animals; unfortunately, existing combinatorial bandit approaches do not offer a way to prioritize important species. To bridge this gap, (1) We propose a novel combinatorial bandit objective that trades off between reward maximization and also accounts for prioritization over species, which we call ranked prioritization. We show this objective can be expressed as a weighted linear sum of Lipschitz-continuous reward functions. (2) We provide RankedCUCB, an algorithm to select combinatorial actions that optimize our prioritization-based objective, and prove that it achieves asymptotic no-regret. (3) We demonstrate empirically that RankedCUCB leads to up to 38% improvement in outcomes for endangered species using real-world wildlife conservation data. Along with adapting to other challenges such as preventing illegal logging and overfishing, our no-regret algorithm addresses the general combinatorial bandit problem with a weighted linear objective.


title: Energy relaxation in a vacuum capacitor-resistor circuit: measurement of multiple decays with divergent time constants

abstract: The decay of the electrical energy in a resistor-vacuum capacitor circuit is shown to involve multiple relaxation processes, with dramatically different time constants. This is measured using a vacuum capacitor to eliminate the effect of a dielectric between the plates (polypropylene capacitors are shown to exhibit similar behavior). A simple phenomenological model accounts for this behavior in spite of the difficulty in applying Maxwell's equations to such a circuit. These results will lead to a revision of our understanding of the physics of circuits, having particular impact on applications that use capacitors as sensors in collecting precision data (such as found in quantum measurements and dielectric spectroscopy).


title: How Platform-User Power Relations Shape Algorithmic Accountability: A Case Study of Instant Loan Platforms and Financially Stressed Users in India

abstract: Accountability, a requisite for responsible AI, can be facilitated through transparency mechanisms such as audits and explainability. However, prior work suggests that the success of these mechanisms may be limited to Global North contexts; understanding the limitations of current interventions in varied socio-political conditions is crucial to help policymakers facilitate wider accountability. To do so, we examined the mediation of accountability in the existing interactions between vulnerable users and a 'high-risk' AI system in a Global South setting. We report on a qualitative study with 29 financially-stressed users of instant loan platforms in India. We found that users experienced intense feelings of indebtedness for the 'boon' of instant loans, and perceived huge obligations towards loan platforms. Users fulfilled obligations by accepting harsh terms and conditions, over-sharing sensitive data, and paying high fees to unknown and unverified lenders. Users demonstrated a dependence on loan platforms by persisting with such behaviors despite risks of harms such as abuse, recurring debts, discrimination, privacy harms, and self-harm to them. Instead of being enraged with loan platforms, users assumed responsibility for their negative experiences, thus releasing the high-powered loan platforms from accountability obligations. We argue that accountability is shaped by platform-user power relations, and urge caution to policymakers in adopting a purely technical approach to fostering algorithmic accountability. Instead, we call for situated interventions that enhance agency of users, enable meaningful transparency, reconfigure designer-user relations, and prompt a critical reflection in practitioners towards wider accountability. We conclude with implications for responsibly deploying AI in FinTech applications in India and beyond.


title: Deep Architecture Connectivity Matters for Its Convergence: A Fine-Grained Analysis

abstract: Advanced deep neural networks (DNNs), designed by either human or AutoML algorithms, are growing increasingly complex. Diverse operations are connected by complicated connectivity patterns, e.g., various types of skip connections. Those topological compositions are empirically effective and observed to smooth the loss landscape and facilitate the gradient flow in general. However, it remains elusive to derive any principled understanding of their effects on the DNN capacity or trainability, and to understand why or in which aspect one specific connectivity pattern is better than another. In this work, we theoretically characterize the impact of connectivity patterns on the convergence of DNNs under gradient descent training in fine granularity. By analyzing a wide network's Neural Network Gaussian Process (NNGP), we are able to depict how the spectrum of an NNGP kernel propagates through a particular connectivity pattern, and how that affects the bound of convergence rates. As one practical implication of our results, we show that by a simple filtration on "unpromising" connectivity patterns, we can trim down the number of models to evaluate, and significantly accelerate the large-scale neural architecture search without any overhead. Code is available at: https://github.com/VITA-Group/architecture_convergence.


title: Classical and quantum butterfly effect in nonlinear vector mechanics

abstract: We establish the correspondence between the classical and quantum butterfly effects in nonlinear vector mechanics with the broken $O(N)$ symmetry. On one hand, we analytically calculate the out-of-time ordered correlation functions and the quantum Lyapunov exponent using the augmented Schwinger-Keldysh technique in the large-$N$ limit. On the other hand, we numerically estimate the classical Lyapunov exponent in the high-temperature limit, where the classical chaotic behavior emerges. In both cases, Lyapunov exponents approximately coincide and scale as $\kappa \approx 1.3 \sqrt[4]{\lambda T}/N$ with temperature $T$, number of degrees of freedom $N$, and coupling constant $\lambda$.

categories: hep-th cond-mat.stat-mech nlin.CD quant-ph


title: Process, Bias and Temperature Scalable CMOS Analog Computing Circuits for Machine Learning

abstract: Analog computing is attractive compared to digital computing due to its potential for achieving higher computational density and higher energy efficiency. However, unlike digital circuits, conventional analog computing circuits cannot be easily mapped across different process nodes due to differences in transistor biasing regimes, temperature variations and limited dynamic range. In this work, we generalize the previously reported margin-propagation-based analog computing framework for designing novel \textit{shape-based analog computing} (S-AC) circuits that can be easily cross-mapped across different process nodes. Similar to digital designs S-AC designs can also be scaled for precision, speed, and power. As a proof-of-concept, we show several examples of S-AC circuits implementing mathematical functions that are commonly used in machine learning (ML) architectures. Using circuit simulations we demonstrate that the circuit input/output characteristics remain robust when mapped from a planar CMOS 180nm process to a FinFET 7nm process. Also, using benchmark datasets we demonstrate that the classification accuracy of a S-AC based neural network remains robust when mapped across the two processes and to changes in temperature.


title: AdS$_3$/AdS$_2$ degression of Fronsdal fields

abstract: We analyze the Kaluza-Klein type procedure in AdS$_3$ space called the dimensional degression. The topological theory of the Fronsdal field in AdS$_3$ is reformulated in terms of the fields propagating in AdS$_2$. We find that the Fronsdal field in AdS$_3$ leads to finitely many Kaluza-Klein modes. Namely, the obtained spectrum is the massive Klein-Gordon and Proca fields in AdS$_2$. The result is derived by using the specific mode expansion, the gauge fixing, and 2-dimensional Schouten identities.


title: Identifying concept libraries from language about object structure

abstract: Our understanding of the visual world goes beyond naming objects, encompassing our ability to parse objects into meaningful parts, attributes, and relations. In this work, we leverage natural language descriptions for a diverse set of 2K procedurally generated objects to identify the parts people use and the principles leading these parts to be favored over others. We formalize our problem as search over a space of program libraries that contain different part concepts, using tools from machine translation to evaluate how well programs expressed in each library align to human language. By combining naturalistic language at scale with structured program representations, we discover a fundamental information-theoretic tradeoff governing the part concepts people name: people favor a lexicon that allows concise descriptions of each object, while also minimizing the size of the lexicon itself.


title: Bubble dynamics in fluids with N-body simulations

abstract: We present a new approach to studies of bubble dynamics in fluids. Relying on particle-based simulations, this method is general and suitable for cases where the commonly used perfect fluid description fails. We study expanding true vacuum bubbles surrounded by free or self-interacting particles and quantify how self-interactions affect the terminal bubble wall velocity. We find that, for sufficiently strongly self-interacting fluids, local thermal equilibrium is maintained around the bubble wall and the fluid profile is similar to that obtained with the perfect fluid description.


title: Short Proofs of Linear Growth of Quantum Circuit Complexity

abstract: The complexity of a quantum gate, defined as the minimal number of elementary gates to build it, is an important concept in quantum information and computation. It is shown recently that the complexity of quantum gates built from random quantum circuits almost surely grows linearly with the number of building blocks. In this article, we provide two short proofs of this fact. We also discuss a discrete version of quantum circuit complexity growth.


title: Ocean signatures in the total flux and polarization spectra of Earth-like exoplanets

abstract: Numerical simulations of starlight that is reflected by Earth-like exoplanets predict habitability signatures that can be searched for with future telescopes. We explore signatures of water oceans in the flux and polarization spectra of this reflected light. With an adding-doubling algorithm, we compute the total flux F, polarized flux Q and degree of polarization P of starlight reflected by dry and ocean model planets with Earth-like atmospheres and patchy clouds. The oceans consist of Fresnel reflecting surfaces with wind-ruffled waves, foam and wave shadows, above natural blue seawater. Our results are presented as functions of wavelength (from 300 to 2500 nm with 1 nm resolution) and as functions of the planetary phase angle from 90 to 170 degrees. The ocean glint increases F, |Q| and P with increasing phase angle at non-absorbing wavelengths, and causes the spectra of F and |Q| for the various phase angles to intersect. In the near-infrared, Q is negative, i.e. the direction of polarization is perpendicular to the plane through the star, planet, and observer. In the P-spectra, the glint leaves dips (instead of peaks) in gaseous absorption bands. All those signatures are missing in the spectra of dry planets. The dips in P, and the negative Q in the near-infrared, can be searched for at a phase angle of 90 degrees, where the planet-star separation is largest. Those ocean signatures in polarized light do not suffer from false positive glint signals that could be due to clouds or reflecting dry surfaces. For heavily cloudy planets, ocean detection is possible when the glint is (partially) cloud-free. When modelling signals of planets with oceans, using horizontally inhomogeneous cloud covers is thus crucial. Observations spread over time would increase the probability of catching a cloud-free glint and detecting an ocean.


title: Short and long-time path tightness of the continuum directed random polymer

abstract: We consider the point-to-point continuum directed random polymer ($\mathsf{CDRP}$) model that arises as a scaling limit from $1+1$ dimensional directed polymers in the intermediate disorder regime. We show that the annealed law of a point-to-point $\mathsf{CDRP}$ of length $t$ converges to the Brownian bridge under diffusive scaling when $t \downarrow 0$. In case that $t$ is large, we show that the transversal fluctuations of point-to-point $\mathsf{CDRP}$ are governed by the $2/3$ exponent. More precisely, as $t$ tends to infinity, we prove tightness of the annealed path measures of point-to-point $\mathsf{CDRP}$ of length $t$ upon scaling the length by $t$ and fluctuations of paths by $t^{2/3}$. The $2/3$ exponent is tight such that the one-point distribution of the rescaled paths converges to the geodesics of the directed landscape. This point-wise convergence can be enhanced to process-level modulo a conjecture. Our short and long-time tightness results also extend to point-to-line $\mathsf{CDRP}$. In the course of proving our main results, we establish quantitative versions of quenched modulus of continuity estimates for long-time $\mathsf{CDRP}$ which are of independent interest.


title: RepSR: Training Efficient VGG-style Super-Resolution Networks with Structural Re-Parameterization and Batch Normalization

abstract: This paper explores training efficient VGG-style super-resolution (SR) networks with the structural re-parameterization technique. The general pipeline of re-parameterization is to train networks with multi-branch topology first, and then merge them into standard 3x3 convolutions for efficient inference. In this work, we revisit those primary designs and investigate essential components for re-parameterizing SR networks. First of all, we find that batch normalization (BN) is important to bring training non-linearity and improve the final performance. However, BN is typically ignored in SR, as it usually degrades the performance and introduces unpleasant artifacts. We carefully analyze the cause of BN issue and then propose a straightforward yet effective solution. In particular, we first train SR networks with mini-batch statistics as usual, and then switch to using population statistics at the later training period. While we have successfully re-introduced BN into SR, we further design a new re-parameterizable block tailored for SR, namely RepSR. It consists of a clean residual path and two expand-and-squeeze convolution paths with the modified BN. Extensive experiments demonstrate that our simple RepSR is capable of achieving superior performance to previous SR re-parameterization methods among different model sizes. In addition, our RepSR can achieve a better trade-off between performance and actual running time (throughput) than previous SR methods. Codes will be available at https://github.com/TencentARC/RepSR.


title: A Discrete Droplet Method for Modelling Thin Film Flows

abstract: In this paper, we present a new model to simulate the formation, evolution, and break up of a thin film of fluid flowing over a curved surface. Referred to as the discrete droplet method (DDM), the model captures the evolution of thin fluid films by tracking individual moving fluid droplets. In contrast to existing thin film models that solve a PDE to determine the film height, here, we compute the film height by numerical integration based on the aggregation of droplets. The novelty of this approach in using droplets makes it suitable for simulating the formation of fluid films, and modelling thin film flows on partially wetted surfaces. The DDM is a Lagrangian approach, with a force balance on each droplet governing the motion, and derivatives approximated using a smoothed particle hydrodynamics (SPH) like approach. The proposed model is thoroughly validated by comparing results against analytical solutions, against the results of the shallow-water equations for thin film flow, and also against results from a full 3D resolved Navier Stokes model. We also present the use of the DDM on an industrial test case. The results highlight the effectiveness of the model for simulations of flows with thin films.


title: Can we actually constrain $f_{\rm NL}$ using the scale-dependent bias effect? An illustration of the impact of galaxy bias uncertainties using the BOSS DR12 galaxy power spectrum

abstract: The scale-dependent bias effect on the galaxy power spectrum is a very promising probe of the local primordial non-Gaussianity (PNG) parameter $f_{\rm NL}$, but the amplitude of the effect is proportional to $f_{\rm NL}b_{\phi}$, where $b_{\phi}$ is the linear PNG galaxy bias parameter. Our knowledge of $b_{\phi}$ is currently very limited, yet nearly all existing $f_{\rm NL}$ constraints and forecasts assume precise knowledge for it. Here, we use the BOSS DR12 galaxy power spectrum to illustrate how our uncertain knowledge of $b_{\phi}$ currently prevents us from constraining $f_{\rm NL}$ with a given statistical precision $\sigma_{f_{\rm NL}}$. Assuming different fixed choices for the relation between $b_{\phi}$ and the linear density bias $b_1$, we find that $\sigma_{f_{\rm NL}}$ can vary by as much as an order of magnitude. Our strongest bound is $f_{\rm NL} = 16 \pm 16\ (1\sigma)$, while the loosest is $f_{\rm NL} = 230 \pm 226\ (1\sigma)$ for the same BOSS data. The impact of $b_{\phi}$ can be especially pronounced because it can be close to zero. We also show how marginalizing over $b_{\phi}$ with wide priors is not conservative, and leads in fact to biased constraints through parameter space projection effects. Independently of galaxy bias assumptions, the scale-dependent bias effect can only be used to detect $f_{\rm NL} \neq 0$ by constraining the product $f_{\rm NL}b_{\phi}$, but the error bar $\sigma_{f_{\rm NL}}$ remains undetermined and the results cannot be compared with the CMB; we find $f_{\rm NL}b_{\phi} \neq 0$ with $1.6\sigma$ significance. We also comment on why these issues are important for analyses with the galaxy bispectrum. Our results strongly motivate simulation-based research programs aimed at robust theoretical priors for the $b_{\phi}$ parameter, without which we may never be able to competitively constrain $f_{\rm NL}$ using galaxy data.


title: Vector stability in quadratic metric-affine theories

abstract: In this work we study the stability of the four vector irreducible pieces of the torsion and the nonmetricity tensors in the general quadratic metric-affine Lagrangian in 4 dimensions. The goal will be to elucidate under which conditions the spin-1 modes associated to such vectors can propagate in a safe way, together with the graviton. This highly constrains the theory reducing the parameter space of the quadratic curvature part from 16 to 5 parameters. We also study the sub-case of Weyl-Cartan gravity, proving that the stability of the vector sector is only compatible with an Einstein-Proca theory for the Weyl vector.


title: NTIRE 2022 Challenge on Efficient Super-Resolution: Methods and Results

abstract: This paper reviews the NTIRE 2022 challenge on efficient single image super-resolution with focus on the proposed solutions and results. The task of the challenge was to super-resolve an input image with a magnification factor of $\times$4 based on pairs of low and corresponding high resolution images. The aim was to design a network for single image super-resolution that achieved improvement of efficiency measured according to several metrics including runtime, parameters, FLOPs, activations, and memory consumption while at least maintaining the PSNR of 29.00dB on DIV2K validation set. IMDN is set as the baseline for efficiency measurement. The challenge had 3 tracks including the main track (runtime), sub-track one (model complexity), and sub-track two (overall performance). In the main track, the practical runtime performance of the submissions was evaluated. The rank of the teams were determined directly by the absolute value of the average runtime on the validation set and test set. In sub-track one, the number of parameters and FLOPs were considered. And the individual rankings of the two metrics were summed up to determine a final ranking in this track. In sub-track two, all of the five metrics mentioned in the description of the challenge including runtime, parameter count, FLOPs, activations, and memory consumption were considered. Similar to sub-track one, the rankings of five metrics were summed up to determine a final ranking. The challenge had 303 registered participants, and 43 teams made valid submissions. They gauge the state-of-the-art in efficient single image super-resolution.


title: Revisiting Random Channel Pruning for Neural Network Compression

abstract: Channel (or 3D filter) pruning serves as an effective way to accelerate the inference of neural networks. There has been a flurry of algorithms that try to solve this practical problem, each being claimed effective in some ways. Yet, a benchmark to compare those algorithms directly is lacking, mainly due to the complexity of the algorithms and some custom settings such as the particular network configuration or training procedure. A fair benchmark is important for the further development of channel pruning. Meanwhile, recent investigations reveal that the channel configurations discovered by pruning algorithms are at least as important as the pre-trained weights. This gives channel pruning a new role, namely searching the optimal channel configuration. In this paper, we try to determine the channel configuration of the pruned models by random search. The proposed approach provides a new way to compare different methods, namely how well they behave compared with random pruning. We show that this simple strategy works quite well compared with other channel pruning methods. We also show that under this setting, there are surprisingly no clear winners among different channel importance evaluation methods, which then may tilt the research efforts into advanced channel configuration searching methods.


title: HULC: 3D Human Motion Capture with Pose Manifold Sampling and Dense Contact Guidance

abstract: Marker-less monocular 3D human motion capture (MoCap) with scene interactions is a challenging research topic relevant for extended reality, robotics and virtual avatar generation. Due to the inherent depth ambiguity of monocular settings, 3D motions captured with existing methods often contain severe artefacts such as incorrect body-scene inter-penetrations, jitter and body floating. To tackle these issues, we propose HULC, a new approach for 3D human MoCap which is aware of the scene geometry. HULC estimates 3D poses and dense body-environment surface contacts for improved 3D localisations, as well as the absolute scale of the subject. Furthermore, we introduce a 3D pose trajectory optimisation based on a novel pose manifold sampling that resolves erroneous body-environment inter-penetrations. Although the proposed method requires less structured inputs compared to existing scene-aware monocular MoCap algorithms, it produces more physically-plausible poses: HULC significantly and consistently outperforms the existing approaches in various experiments and on different metrics. Project page: https://vcai.mpi-inf.mpg.de/projects/HULC/.


title: RISP: Rendering-Invariant State Predictor with Differentiable Simulation and Rendering for Cross-Domain Parameter Estimation

abstract: This work considers identifying parameters characterizing a physical system's dynamic motion directly from a video whose rendering configurations are inaccessible. Existing solutions require massive training data or lack generalizability to unknown rendering configurations. We propose a novel approach that marries domain randomization and differentiable rendering gradients to address this problem. Our core idea is to train a rendering-invariant state-prediction (RISP) network that transforms image differences into state differences independent of rendering configurations, e.g., lighting, shadows, or material reflectance. To train this predictor, we formulate a new loss on rendering variances using gradients from differentiable rendering. Moreover, we present an efficient, second-order method to compute the gradients of this loss, allowing it to be integrated seamlessly into modern deep learning frameworks. We evaluate our method in rigid-body and deformable-body simulation environments using four tasks: state estimation, system identification, imitation learning, and visuomotor control. We further demonstrate the efficacy of our approach on a real-world example: inferring the state and action sequences of a quadrotor from a video of its motion sequences. Compared with existing methods, our approach achieves significantly lower reconstruction errors and has better generalizability among unknown rendering configurations.


title: Linking stellar populations to HII regions across nearby galaxies: I. Constraining pre-supernova feedback from young clusters in NGC1672

abstract: One of the fundamental factors regulating the evolution of galaxies is stellar feedback. However, we still do not have strong observational constraints on the relative importance of the different feedback mechanisms (e.g. radiation, ionised gas pressure, stellar winds) in driving HII region evolution and molecular cloud disruption. In this letter, we constrain the relative importance of the various feedback mechanisms from young massive star populations by resolving HII regions across the disk of the nearby star-forming galaxy NGC 1672. We combine measurements of ionised gas nebular lines obtained by PHANGS-MUSE, with high-resolution imaging from the HST in both the narrow-band H{\alpha} and broad-band filters. We identify a sample of 40 isolated, compact HII regions in the HST H{\alpha} image, for which we measure the sizes that were previously unresolved in seeing-limited ground-based observations. Additionally, we identify the ionisation source(s) for each HII region from catalogues produced as part of the PHANGS-HST survey. We find that the HII regions investigated are mildly dominated by thermal or wind pressure, yet their elevation above the radiation pressure is within the expected uncertainty range. We see that radiation pressure provides a substantially higher contribution to the total pressure than previously found in the literature over similar size scales. In general, we find higher pressures within more compact HII regions, which is driven by the inherent size scaling relations of each pressure term, albeit with significant scatter introduced by the variation in the stellar population properties (e.g. luminosity, mass, age, metallicity). For nearby galaxies, here we provide a promising approach that could yield the statistics required to map out how the importance of different stellar feedback mechanisms evolve over the lifetime of an HII region.


title: Ion beam analysis and big data: How data science can support next-generation instrumentation

abstract: With a growing demand for accurate ion beam analysis on a large number of samples, it becomes an issue of how to ensure the quality standards and consistency over hundreds or thousands of samples. In this sense, a virtual assistant that checks the data quality, emitting certificates of quality, is highly desired. Even the processing of a massive number of spectra is a problem regarding the consistency of the analysis. In this work, we report the design and first results of a virtual layer under implementation in our laboratory. It consists of a series of systems running in the cloud that perform the mentioned tasks and serves as a virtual assistant for member staff and users. We aim to bring the concept of the Internet of Things and artificial intelligence closer to the laboratory to support a new generation of instrumentation.


title: Accelerating reactive-flow simulations using vectorized chemistry integration

abstract: The high cost of chemistry integration is a significant computational bottleneck for realistic reactive-flow simulations using operator splitting. Here we present a methodology to accelerate the solution of the chemical kinetic ordinary differential equations using single-instruction, multiple-data vector processing on CPUs using the OpenCL framework. First, we compared several vectorized integration algorithms using chemical kinetic source terms and analytical Jacobians from the pyJac software against a widely used integration code, CVODEs. Next, we extended the OpenFOAM computational fluid dynamics library to incorporate the vectorized solvers, and we compared the accuracy of a fourth-order linearly implicit integrator -- both in vectorized form and a corresponding method native to OpenFOAM -- with the community standard chemical kinetics library Cantera. We then applied our methodology to a variety of chemical kinetic models, turbulent intensities, and simulation scales to examine a range of engineering and scientific scale problems, including (pseudo) steady-state as well as time-dependent Reynolds-averaged Navier--Stokes simulations of the Sandia flame D and the Volvo Flygmotor bluff-body stabilized, premixed flame. Subsequently, we compared the performance of the vectorized and native OpenFOAM integrators over the studied models and simulations and found that our vectorized approach performs up to 33--35x faster than the native OpenFOAM solver with high accuracy.


title: Explainable Computational Creativity

abstract: Human collaboration with systems within the Computational Creativity (CC) field is often restricted to shallow interactions, where the creative processes, of systems and humans alike, are carried out in isolation, without any (or little) intervention from the user, and without any discussion about how the unfolding decisions are taking place. Fruitful co-creation requires a sustained ongoing interaction that can include discussions of ideas, comparisons to previous/other works, incremental improvements and revisions, etc. For these interactions, communication is an intrinsic factor. This means giving a voice to CC systems and enabling two-way communication channels between them and their users so that they can: explain their processes and decisions, support their ideas so that these are given serious consideration by their creative collaborators, and learn from these discussions to further improve their creative processes. For this, we propose a set of design principles for CC systems that aim at supporting greater co-creation and collaboration with their human collaborators.


title: Impact of Charge on the Complexity of Static Sphere in $f(R,\textbf{T}^{2})$ Gravity

abstract: This paper investigates the complexity of a charged static sphere filled with anisotropic matter in the background of energy-momentum squared gravity. For this purpose, we evaluate the modified field and conservation equations to determine the structure of celestial system. The mass function is calculated through Misner-Sharp as well as Tolman mass definitions. The complexity of a self-gravitating system depends on different factors such as anisotropic pressure, electromagnetic field, energy density inhomogeneity, etc. We formulate the structure scalars by the orthogonal decomposition of the Riemann tensor to develop a complexity factor containing all vital features of the stellar structure. The vanishing complexity condition is achieved by setting the complexity factor equal to zero. Finally, we construct two static solutions by utilizing the energy density of Gokhroo-Mehra solution as well as the polytropic equation of state along with the zero complexity condition. It is found that electromagnetic field decreases the complexity of stellar structure.


title: A Closer Look at Audio-Visual Multi-Person Speech Recognition and Active Speaker Selection

abstract: Audio-visual automatic speech recognition is a promising approach to robust ASR under noisy conditions. However, up until recently it had been traditionally studied in isolation assuming the video of a single speaking face matches the audio, and selecting the active speaker at inference time when multiple people are on screen was put aside as a separate problem. As an alternative, recent work has proposed to address the two problems simultaneously with an attention mechanism, baking the speaker selection problem directly into a fully differentiable model. One interesting finding was that the attention indirectly learns the association between the audio and the speaking face even though this correspondence is never explicitly provided at training time. In the present work we further investigate this connection and examine the interplay between the two problems. With experiments involving over 50 thousand hours of public YouTube videos as training data, we first evaluate the accuracy of the attention layer on an active speaker selection task. Secondly, we show under closer scrutiny that an end-to-end model performs at least as well as a considerably larger two-step system that utilizes a hard decision boundary under various noise conditions and number of parallel face tracks.


title: Search for baryon junctions in photonuclear processes and isobar collisions at RHIC

abstract: During the early development of Quantum Chromodynamics, it was proposed that baryon number could be carried by a non-perturbative Y-shaped topology of gluon fields, called the gluon junction, rather than by the valence quarks as in the QCD standard model. A puzzling feature of ultra-relativistic nucleus-nucleus collisions is the apparent substantial baryon excess in the midrapidity region that could not be adequately accounted for in most conventional models of quark and diquark transport. The transport of baryonic gluon junctions is predicted to lead to a characteristic exponential distribution of net-baryon density with rapidity and could resolve the puzzle. In this context we point out that the rapidity density of net-baryons near midrapidity indeed follows an exponential distribution with a slope of $-0.61\pm0.03$ as a function of beam rapidity in the existing global data from A+A collisions at AGS, SPS and RHIC energies. To further test if quarks or gluon junctions carry the baryon quantum number, we propose to study the absolute magnitude of the baryon vs. charge stopping in isobar collisions at RHIC. We also argue that semi-inclusive photon-induced processes ($\gamma+p$/A) at RHIC kinematics provide an opportunity to search for the signatures of the baryon junction and to shed light onto the mechanisms of observed baryon excess in the mid-rapidity region in ultra-relativistic nucleus-nucleus collisions. Such measurements can be further validated in A+A collisions at the LHC and $e+p$/A collisions at the EIC.


title: Diatomic-py: A python module for calculating the rotational and hyperfine structure of $^1\Sigma$ molecules

abstract: We present a computer program to calculate the quantised rotational and hyperfine energy levels of $^{1}\Sigma $ diatomic molecules in the presence of dc electric, dc magnetic, and off-resonant optical fields. Our program is applicable to the bialkali molecules used in ongoing state-of-the-art experiments with ultracold molecular gases. We include functions for the calculation of space-fixed electric dipole moments, magnetic moments and transition dipole moments.


title: The Value-Added Catalog of ASAS-SN Eclipsing Binaries: Parameters of Thirty Thousand Detached Systems

abstract: Detached eclipsing binaries are a fundamental tool for measuring the physical parameters of stars that are effectively evolving in isolation. Starting from more than 40,000 eclipsing binary candidates identified by the All-Sky Automated Survey for Supernovae (ASAS-SN), we use PHOEBE to determine the sum of the fractional radii, the ratio of effective temperatures, the inclinations, and the eccentricities for 35,464 systems. We visually inspect all the light curve models to verify the model fits and examine the TESS light curves, when available, to select systems with evidence for additional physics, such as spots, mass transfer, and hierarchical triples. We examine the distributions of the eclipsing binary model parameters and the orbital parameters. We identify two groups in the sum of the fractional radii and effective temperature ratio parameter space that may distinguish systems approaching the semidetached limit. Combining Gaia EDR3 with extinction estimates from 3-dimensional dust maps, we examine the properties of the systems as a function of their absolute magnitude and evolutionary state. Finally, we present light curves of selected eclipsing binaries that may be of interest for follow-up studies.


title: Quantum complexity and topological phases of matter

abstract: In this work, we find that the complexity of quantum many-body states, defined as a spread in the Krylov basis, may serve as a new probe that distinguishes topological phases of matter. We illustrate this analytically in one of the representative examples, the Su-Schrieffer-Heeger model, finding that spread complexity becomes constant in the topological phase. Moreover, in the same setup, we analyze exactly solvable quench protocols where the evolution of the spread complexity shows distinct dynamical features depending on the topological vs non-topological phase of the initial state as well as the quench Hamiltonian.


title: Greedy optimization of the geometry of Majorana Josephson junctions

abstract: Josephson junctions in a two-dimensional electron gas with spin-orbit coupling are a promising candidate to realize topological superconductivity. While it is known that the geometry of the junction strongly influences the size of the topological gap, the question of how to construct optimal geometries remains unexplored. We introduce a greedy numerical algorithm to optimize the shape of Majorana junctions. The core of the algorithm relies on perturbation theory and is embarrassingly parallel, which allows it to explore the design space efficiently. By introducing stochastic variations in the junction Hamiltonian, we avoid overfitting geometries to specific system parameters. Furthermore, we constrain the optimizer to produce smooth geometries by applying image filtering and fabrication resolution constraints. We run the algorithm in various setups and find that it reliably produces geometries with increased topological gaps over large parameter ranges. The results are robust to variations in the optimization starting point and the presence of disorder, which suggests the optimizer is capable of finding global maxima.


title: Wide twin binaries are extremely eccentric: evidence of twin binary formation in circumbinary disks

abstract: The Gaia mission recently revealed an excess population of equal-mass "twin" wide binaries, with mass ratio $q\gtrsim 0.95$, extending to separations of at least 1000 AU. The origin of this population is an enigma: twin binaries are thought to form via correlated accretion in circumbinary disks, but the typical observed protostellar disks have radii of $\sim100$ AU, far smaller than the separations of the widest twins. Here, we infer the eccentricity distribution of wide twins from the distribution of their $v$-$r$ angles, i.e., the angle between the components' separation and relative velocity vectors. We find that wide twins must be on extremely eccentric orbits. For the excess-twin population at 400-1000 AU, we infer a near-delta function excess of high-eccentricity system, with eccentricity $0.95 \lesssim e \leq 1$. These high eccentricities for wide twins imply pericenter distances of order $10$ AU and suggest that their orbits were scattered via dynamical interactions in their birth environments, consistent with a scenario in which twins are born in circumbinary disks and subsequently widened. These results further establish twin wide binaries as a distinct population and imply that wide twins can be used as a probe of the dynamical history of stellar populations.


title: Quark formation and phenomenology in binary neutron-star mergers using V-QCD

abstract: Using full 3+1 dimensional general-relativistic hydrodynamic simulations of equal- and unequal-mass neutron-star binaries with properties that are consistent with those inferred from the inspiral of GW170817, we perform a detailed study of the quark-formation processes that could take place after merger. We use three equations of state consistent with current pulsar observations derived from a novel finite-temperature framework based on V-QCD, a non-perturbative gauge/gravity model for Quantum Chromodynamics. In this way, we identify three different post-merger stages at which mixed baryonic and quark matter, as well as pure quark matter, are generated. A phase transition triggered collapse already $\lesssim 10\,\rm{ms}$ after the merger reveals that the softest version of our equations of state is actually inconsistent with the expected second-long post-merger lifetime of GW170817. Our results underline the impact that multi-messenger observations of binary neutron-star mergers can have in constraining the equation of state of nuclear matter, especially in its most extreme regimes.


title: Probing sign structure using measurement-induced entanglement

abstract: The sign structure of quantum states is closely connected to quantum phases of matter, yet detecting such fine-grained properties of amplitudes is subtle. Here we employ as a diagnostic measurement-induced entanglement (MIE): the average entanglement generated between two parties after measuring the rest of the system. We propose that for a sign-free state, the MIE upon measuring in the sign-free basis decays no slower than correlations in the state before measurement. Concretely, we prove that MIE is upper bounded by mutual information for sign-free stabilizer states (essentially CSS codes), which establishes a bound between scaling dimensions of conformal field theories describing measurement-induced critical points in stabilizer systems. We also show that for sign-free qubit wavefunctions, MIE between two qubits is upper bounded by a simple two-point correlation function, and we verify our proposal in several critical ground states of one-dimensional systems, including the transverse field and tri-critical Ising models. In contrast, for states with sign structure, such bounds can be violated, as we illustrate in critical hybrid circuits involving both Haar or Clifford random unitaries and measurements, and gapless symmetry-protected topological states.


title: Towards Explicit Discrete Holography: Aperiodic Spin Chains from Hyperbolic Tilings

abstract: We propose a new example of discrete holography that provides a new step towards establishing the AdS/CFT duality for discrete spaces. A class of boundary Hamiltonians is obtained in a natural way from regular tilings of the hyperbolic Poincar\'e disk, via an inflation rule that allows to construct the tiling using concentric layers of tiles. The models in this class are aperiodic spin chains, whose sequences of couplings are obtained from the bulk inflation rule. We explicitly choose the aperiodic XXZ spin chain with spin 1/2 degrees of freedom as an example. The properties of this model are studied by using strong disorder renormalization group techniques, which provide a tensor network construction for the ground state of this spin chain. This can be regarded as discrete bulk reconstruction. Moreover we compute the entanglement entropy in this setup in two different ways: a discretization of the Ryu-Takayanagi formula and a generalization of the standard computation for the boundary aperiodic Hamiltonian. For both approaches, a logarithmic growth of the entanglement entropy in the subsystem size is identified. The coefficients, i.e. the effective central charges, depend on the bulk discretization parameters in both cases, albeit in a different way.


title: SkyHopper mission science case I: Identification of high redshift Gamma-Ray Bursts through space-based near-infrared afterglow observations

abstract: Long-duration gamma-ray burst (GRB) afterglow observations offer cutting-edge opportunities to characterise the star formation history of the Universe back to the epoch of reionisation, and to measure the chemical composition of interstellar and intergalactic gas through absorption spectroscopy. The main barrier to progress is the low efficiency in rapidly and confidently identifying which bursts are high redshift ($z > 5$) candidates before they fade, as this requires low-latency follow-up observations at near-infrared wavelengths (or longer) to determine a reliable photometric redshift estimate. So far this task has been performed by instruments on the ground, but sky visibility and weather constraints limit the number of GRB targets that can be observed and the speed at which follow-up is possible. In this work we develop a Monte Carlo simulation framework to investigate an alternative approach based on the use of a rapid-response near-infrared nano-satellite, capable of simultaneous imaging in four bands from $0.8$ to $1.7\mu$m (a mission concept called SkyHopper). We find that such a nano-satellite is capable of detecting in the H band (1.6 $\mu$m) $72.5\% \pm 3.1\%$ of GRBs concurrently observable with the Swift satellite via its UVOT instrument (and $44.1\% \pm 12.3\%$ of high redshift ($z>5$) GRBs) within 60 minutes of the GRB prompt emission. This corresponds to detecting $\sim 55$ GRB afterglows per year, of which 1-3 have $z > 5$. These rates represent a substantial contribution to the field of high-$z$ GRB science, as only 23 $z > 5$ GRBs have been collectively discovered by the entire astronomical community over the last $\sim 24$ years. Additionally, we find that launching a mini-constellation of 3 near-infrared nano-satellites would increase the detection fraction of afterglows to $\sim 83\%$ and substantially reduce the latency in the photometric redshift determination.


title: Fracton magnetohydrodynamics

abstract: We extend recent work on hydrodynamics with global multipolar symmetries -- known as "fracton hydrodynamics" -- to systems in which the multipolar symmetries are gauged. We refer to the latter as "fracton magnetohydrodynamics", in analogy to conventional magnetohydrodynamics (MHD), which governs systems with gauged charge conservation. We show that fracton MHD arises naturally from higher-rank Maxwell's equations and in systems with one-form symmetries obeying certain constraints; while we focus on "minimal" higher-rank generalizations of MHD that realize diffusion, our methods may also be used to identify other, more exotic hydrodynamic theories (e.g., with magnetic subdiffusion). In contrast to semi-microscopic derivations of MHD, our approach elucidates the origin of the hydrodynamic modes by identifying the corresponding higher-form symmetries. Being rooted in symmetries, the hydrodynamic modes may persist even when the semi-microscopic equations no longer provide an accurate description of the system.


title: Direct Imaging and Spectroscopy of Extrasolar Planets

abstract: Direct imaging and spectroscopy is the likely means by which we will someday identify, confirm, and characterize an Earth-like planet around a nearby Sun-like star. This Chapter summarizes the current state of knowledge regarding discovering and characterizing exoplanets by direct imaging and spectroscopy. We detail instruments and software needed for direct imaging detections and summarize the current inventory of confirmed and candidate directly-imaged exoplanets. Direct imaging and spectroscopy in the past decade has provided key insights into jovian planet atmospheres, probed the demographics of the outskirts of planetary systems, and shed light on gas giant planet formation. We forecast the new tools and future facilities on the ground and in space that will enhance our capabilities for exoplanet imaging and will likely image habitable zone rocky planets around the nearest stars.


title: A Well-Posed UV Completion for Simulating Scalar Galileons

abstract: The Galileon scalar field theory is a prototypical example of an effective field theory that exhibits the Vainshtein screening mechanism, which is incorporated into many extensions to Einstein gravity. The Galileon describes the helicity zero mode of gravitational radiation, the presence of which has significant implications for predictions of gravitational waves from orbiting objects, and for tests of gravity sensitive to additional polarizations. Because of the derivative nature of their interactions, Galileons are superficially not well-posed as effective field theories. Although this property is properly understood merely as an artifact of the effective field theory truncation, and is not theoretically worrisome, at the practical level it nevertheless renders numerical simulation highly problematic. Notwithstanding, previous numerical approaches have successfully evolved the system for reasonable initial data by slowly turning on the interactions. We present here two alternative approaches to improving numerical stability in Galileon numerical simulations. One of these is a minor modification of previous approaches, which introduces a low pass filter that amounts to imposing a UV cutoff together with a relaxation method of turning on interactions. The second approach amounts to constructing a (numerical) UV completion for which the dynamics of the high momentum modes is under control, and for which it is unnecessary to slowly turn on nonlinear interactions. We show that numerical simulations of the UV theory successfully reproduce the correct Galileon dynamics at low energies, consistent with the low-pass filter method and with previous numerical simulations.


title: VERTICO II: effects of HI-identified environmental mechanisms on molecular gas

abstract: In this VERTICO early science paper we explore in detail how environmental mechanisms, identified in HI, affect the resolved properties of molecular gas reservoirs in cluster galaxies. The molecular gas is probed using ALMA ACA (+TP) observations of 12CO(2-1) in 51 spiral galaxies in the Virgo cluster (of which 49 are detected), all of which are included in the VIVA HI survey. The sample spans a stellar mass range of 9 < log M*/Msol < 11. We study molecular gas radial profiles, isodensity radii, and surface densities as a function of galaxy HI deficiency and morphology. There is a weak correlation between global HI and H2 deficiencies, and resolved properties of molecular gas correlate with HI deficiency: galaxies that have large HI deficiencies have relatively steep and truncated molecular gas radial profiles, which is due to the removal of low-surface density molecular gas on the outskirts. Therefore, while the environmental mechanisms observed in HI also affect molecular gas reservoirs, there is only a moderate reduction of the total amount of molecular gas.


title: Proper Motions, Orbits, and Tidal Influences of Milky Way Dwarf Spheroidal Galaxies

abstract: We combine Gaia EDR3 astrometry with accurate photometry and utilize a probabilistic mixture model to measure the systemic proper motion of 52 dwarf spheroidal (dSph) satellite galaxies of the Milky Way (MW). For the 46 dSphs with literature line-of-sight velocities we compute orbits in both a MW and a combined MW + Large Magellanic Cloud (LMC) potential and identify Car II, Car III, Hor I, Hyi I, Phx II, and Ret II as likely LMC satellites. 40% of our dSph sample has a >25% change in pericenter and/or apocenter with the MW + LMC potential. For these orbits, we Monte Carlo sample over the observational uncertainties for each dSph and the uncertainties in the MW and LMC potentials. We predict that Ant II, Boo III, Cra II, Gru II, and Tuc III should be be tidally disrupting by comparing each dSph's average density relative to the MW density at its pericenter. dSphs with large ellipticity (CVn I, Her, Tuc V, UMa I, UMa II, UMi, Wil 1) show a preference for their orbital direction to align with their major axis even for dSphs with large pericenters. We compare the dSph radial orbital phase to subhalos in MW-like N-body simulations and infer that there is not an excess of satellites near their pericenter. With projections of future Gaia data releases, we find dSph orbital precision will be limited by uncertainties in the distance and/or MW potential rather than proper motion precision. Finally, we provide our membership catalogs to enable community follow-up.