title: A Deep Learning Study on Osteosarcoma Detection from Histological Images

abstract: In the U.S, 5-10\% of new pediatric cases of cancer are primary bone tumors. The most common type of primary malignant bone tumor is osteosarcoma. The intention of the present work is to improve the detection and diagnosis of osteosarcoma using computer-aided detection (CAD) and diagnosis (CADx). Such tools as convolutional neural networks (CNNs) can significantly decrease the surgeon's workload and make a better prognosis of patient conditions. CNNs need to be trained on a large amount of data in order to achieve a more trustworthy performance. In this study, transfer learning techniques, pre-trained CNNs, are adapted to a public dataset on osteosarcoma histological images to detect necrotic images from non-necrotic and healthy tissues. First, the dataset was preprocessed, and different classifications are applied. Then, Transfer learning models including VGG19 and Inception V3 are used and trained on Whole Slide Images (WSI) with no patches, to improve the accuracy of the outputs. Finally, the models are applied to different classification problems, including binary and multi-class classifiers. Experimental results show that the accuracy of the VGG19 has the highest, 96\%, performance amongst all binary classes and multiclass classification. Our fine-tuned model demonstrates state-of-the-art performance on detecting malignancy of Osteosarcoma based on histologic images.


title: Constrained Serial Rule on the Full Preference Domain

abstract: We study the problem of assigning objects to agents in the presence of arbitrary linear constraints when agents are allowed to be indifferent between objects. Our main contribution is the generalization of the (Extended) Probabilistic Serial mechanism via a new mechanism called the Constrained Serial Rule. This mechanism is computationally efficient and maintains desirable efficiency and fairness properties namely constrained ordinal efficiency and envy-freeness among agents of the same type. Our mechanism is based on a linear programming approach that accounts for all constraints and provides a re-interpretation of the bottleneck set of agents that form a crucial part of the Extended Probabilistic Serial mechanism.


title: Assessing racial inequality in COVID-19 testing with Bayesian threshold tests

abstract: There are racial disparities in the COVID-19 test positivity rate, suggesting that minorities may be under-tested. Here, drawing on the literature on statistically assessing racial disparities in policing, we 1) illuminate a statistical flaw, known as infra-marginality, in using the positivity rate as a metric for assessing racial disparities in under-testing; 2) develop a new type of Bayesian threshold test to measure disparities in COVID-19 testing and 3) apply the test to measure racial disparities in testing thresholds in a real-world COVID-19 dataset.


title: The Harmonic Quantum Szil\'ard Engine

abstract: The Szil\'ard engine is a mechanism (akin to Maxwell's demon) for converting information into energy, which seemingly violates the second law of thermodynamics. Originally a classical thought experiment, it was extended to a quantized treatment by Zurek. Here, we examine a new, elegant model of a quantum Szil\'ard engine by replacing the traditional rigid box with a harmonic potential, extending the scope of the model. Remarkably, almost all calculations are exact. This article is suitable to students, researchers and educators interested in the conceptual links between information, entropy, and quantum measurement.


title: QMUL-SDS @ SardiStance: Leveraging Network Interactions to Boost Performance on Stance Detection using Knowledge Graphs

abstract: This paper presents our submission to the SardiStance 2020 shared task, describing the architecture used for Task A and Task B. While our submission for Task A did not exceed the baseline, retraining our model using all the training tweets, showed promising results leading to (f-avg 0.601) using bidirectional LSTM with BERT multilingual embedding for Task A. For our submission for Task B, we ranked 6th (f-avg 0.709). With further investigation, our best experimented settings increased performance from (f-avg 0.573) to (f-avg 0.733) with same architecture and parameter settings and after only incorporating social interaction features -- highlighting the impact of social interaction on the model's performance.


title: Almost elementariness and fiberwise amenability for \'etale groupoids

abstract: In this paper, we introduce two new approximation properties for \'etale groupoids, almost elementariness and (ubiquitous) fiberwise amenability, inspired by Matui's and Kerr's notions of almost finiteness. In fact, we show that, in their respective scopes of applicability, both notions of almost finiteness are equivalent to the conjunction of our two properties. These new properties stem from viewing \'etale groupoids as coarse geometric objects in the spirit of geometric group theory. Fiberwise amenability is a coarse geometric property of \'etale groupoids that is closely related to the existence of invariant measures on unit spaces and corresponds to the amenability of the acting group in a transformation groupoid. Almost elementariness may be viewed as a better dynamical analogue of the regularity properties of C*-algebras than almost finiteness, since, unlike the latter, the former may also be applied to the purely infinite case. To support this analogy, we show almost elementary minimal groupoids give rise to tracially Z-stable reduced groupoid C*-algebras. In particular, the reduced C*-algebras of second countable amenable almost finite groupoids in Matui's sense are Z-stable.


title: Adversarial Examples in Constrained Domains

abstract: Machine learning algorithms have been shown to be vulnerable to adversarial manipulation through systematic modification of inputs (e.g., adversarial examples) in domains such as image recognition. Under the default threat model, the adversary exploits the unconstrained nature of images; each feature (pixel) is fully under control of the adversary. However, it is not clear how these attacks translate to constrained domains that limit which and how features can be modified by the adversary (e.g., network intrusion detection). In this paper, we explore whether constrained domains are less vulnerable than unconstrained domains to adversarial example generation algorithms. We create an algorithm for generating adversarial sketches: targeted universal perturbation vectors which encode feature saliency within the envelope of domain constraints. To assess how these algorithms perform, we evaluate them in constrained (e.g., network intrusion detection) and unconstrained (e.g., image recognition) domains. The results demonstrate that our approaches generate misclassification rates in constrained domains that were comparable to those of unconstrained domains (greater than 95%). Our investigation shows that the narrow attack surface exposed by constrained domains is still sufficiently large to craft successful adversarial examples; and thus, constraints do not appear to make a domain robust. Indeed, with as little as five randomly selected features, one can still generate adversarial examples.


title: Deep subwavelength thermal switch via resonant coupling in monolayer hexagonal boron nitride

abstract: Unlike the electrical conductance that can be widely modulated within the same material even in deep nanoscale devices, tuning the thermal conductance within a single material system or nanostructure is extremely challenging and requires a large-scale device. This prohibits the realization of robust ON/OFF states in switching the flow of thermal currents. Here, we present the theory of a thermal switch based on resonant coupling of three photonic resonators, in analogy to the field-effect electronic transistor composed of a source, gate, and drain. As a material platform, we capitalize on the extreme tunability and low-loss resonances observed in the dielectric function of monolayer hexagonal boron nitride (hBN) under controlled strain. We derive the dielectric function of hBN from first principles, including the phonon-polariton linewidths computed by considering phonon isotope and anharmonic phonon-phonon scattering. Subsequently, we propose a strain-controlled hBN-based thermal switch that modulates thermal conductance by more than an order of magnitude, corresponding to an ON/OFF contrast ratio of 98%, in a deep subwavelength nanostructure.


title: The very forward CASTOR calorimeter of the CMS experiment

abstract: The physics motivation, detector design, triggers, calibration, alignment, simulation, and overall performance of the very forward CASTOR calorimeter of the CMS experiment are reviewed. The CASTOR Cherenkov sampling calorimeter is located very close to the LHC beam line, at a radial distance of about 1 cm from the beam pipe, and at 14.4 m from the CMS interaction point, covering the pseudorapidity range of $-$6.6 $\lt$ $\eta$ $\lt$ $-$5.2. It was designed to withstand high ambient radiation and strong magnetic fields. The performance of the detector in measurements of forward energy density, jets, and processes characterized by rapidity gaps, is reviewed using data collected in proton and nuclear collisions at the LHC.


title: A positive proportion of cubic fields are not monogenic yet have no local obstruction to being so

abstract: We show that a positive proportion of cubic fields are not monogenic, despite having no local obstruction to being monogenic. Our proof involves the comparison of $2$-descent and $3$-descent in a certain family of Mordell curves $E_k \colon y^2 = x^3 + k$. As a by-product of our methods, we show that, for every $r \geq 0$, a positive proportion of curves $E_k$ have Tate--Shafarevich group with $3$-rank at least $r$.


title: Implementing non-scalar diffraction in Fourier optics via the Braunbek method

abstract: Fourier optics is a powerful and efficient tool for solving many diffraction problems, but relies on the assumption of scalar diffraction theory and ignores the three-dimensional structure and material properties of the diffracting element. Recent experiments of sub-scale starshade external occulters revealed that the inclusion of these physical properties is necessary to explain the observed diffraction at 1e-10 of the incident light intensity. Here, we present a methodology for implementing non-scalar diffraction while maintaining the efficiency and ease of standard Fourier optics techniques. Our methodology is based on that of Braunbek, in which the Kirchhoff boundary values are replaced with the exact field in a narrow seam surrounding the edge of the diffracting element. In this paper, we derive the diffraction equations used to implement non-scalar diffraction and outline the computational implementation used to solve those equations. We also provide experimental results that demonstrate our model can replicate the observational signatures of non-scalar diffraction in sub-scale starshades, in effect validating our model to better than 1e-10 in relative intensity. We believe this method to be an efficient tool for including additional physics to the models of coronagraphs and other optical systems in which a full electromagnetic solution is intractable.


title: RandomForestMLP: An Ensemble-Based Multi-Layer Perceptron Against Curse of Dimensionality

abstract: We present a novel and practical deep learning pipeline termed RandomForestMLP. This core trainable classification engine consists of a convolutional neural network backbone followed by an ensemble-based multi-layer perceptrons core for the classification task. It is designed in the context of self and semi-supervised learning tasks to avoid overfitting while training on very small datasets. The paper details the architecture of the RandomForestMLP and present different strategies for neural network decision aggregation. Then, it assesses its robustness to overfitting when trained on realistic image datasets and compares its classification performance with existing regular classifiers.


title: Exploring the Galaxy's halo and very metal-weak thick disk with SkyMapper and Gaia DR2

abstract: In this work we combine spectroscopic information from the \textit{SkyMapper survey for Extremely Metal-Poor stars} and astrometry from Gaia DR2 to investigate the kinematics of a sample of 475 stars with a metallicity range of $ -6.5 \leq \rm [Fe/H] \leq -2.05$ dex. Exploiting the action map, we identify 16 and 40 stars dynamically consistent with the \textit{Gaia Sausage} and \textit{Gaia Sequoia} accretion events, respectively. The most metal-poor of these candidates have metallicities of $\rm [Fe/H]=-3.31$ and $\rm [Fe/H]=-3.74$, respectively, helping to define the low-metallicity tail of the progenitors involved in the accretion events. We also find, consistent with other studies, that $\sim$21\% of the sample have orbits that remain confined to within 3~kpc of the Galactic plane, i.e., |Z$_{max}$| $\leq$ 3~kpc. Of particular interest is a sub-sample ($\sim$11\% of the total) of low |Z$_{max}$| stars with low eccentricities and prograde motions. The lowest metallicity of these stars has [Fe/H] = --4.30 and the sub-sample is best interpreted as the very low-metallicity tail of the metal-weak thick disk population. The low |Z$_{max}$|, low eccentricity stars with retrograde orbits are likely accreted, while the low |Z$_{max}$|, high eccentricity pro- and retrograde stars are plausibly associated with the \textit{Gaia Sausage} system. We find that a small fraction of our sample ($\sim$4\% of the total) is likely escaping from the Galaxy, and postulate that these stars have gained energy from gravitational interactions that occur when infalling dwarf galaxies are tidally disrupted.


title: Ribbon distance bounds from Bar-Natan Homology and $\alpha$-Homology

abstract: We prove a lower bound on the ribbon distance via Bar-Natan Homology. This lower bound agrees with Alishahi's lower bound on the unknotting number via Bar-Natan Homology, which furthers a pattern first observed by Sarkar: Sarkar's lower bound on the ribbon distance via Lee Homology agreed with Alishahi and Dowlin's lower bound on the unknotting number via Lee Homology. We also prove a lower bound on the ribbon distance via $\alpha$-homology, a variant of Khovanov homology defined recently by Khovanov and Robert.


title: Reducing Neural Network Parameter Initialization Into an SMT Problem

abstract: Training a neural network (NN) depends on multiple factors, including but not limited to the initial weights. In this paper, we focus on initializing deep NN parameters such that it performs better, comparing to random or zero initialization. We do this by reducing the process of initialization into an SMT solver. Previous works consider certain activation functions on small NNs, however the studied NN is a deep network with different activation functions. Our experiments show that the proposed approach for parameter initialization achieves better performance comparing to randomly initialized networks.


title: Budget Sharing for Multi-Analyst Differential Privacy

abstract: Large organizations that collect data about populations (like the US Census Bureau) release summary statistics that are used by multiple stakeholders for resource allocation and policy making problems. These organizations are also legally required to protect the privacy of individuals from whom they collect data. Differential Privacy (DP) provides a solution to release useful summary data while preserving privacy. Most DP mechanisms are designed to answer a single set of queries. In reality, there are often multiple stakeholders that use a given data release and have overlapping but not-identical queries. This introduces a novel joint optimization problem in DP where the privacy budget must be shared among different analysts. We initiate study into the problem of DP query answering across multiple analysts. To capture the competing goals and priorities of multiple analysts, we formulate three desiderata that any mechanism should satisfy in this setting -- The Sharing Incentive, Non-Interference, and Adaptivity -- while still optimizing for overall error. We demonstrate how existing DP query answering mechanisms in the multi-analyst settings fail to satisfy at least one of the desiderata. We present novel DP algorithms that provably satisfy all our desiderata and empirically show that they incur low error on realistic tasks.


title: Pointwise Lineability in Sequence Spaces

abstract: We prove that several results of lineability/spaceability in the framework of sequence spaces are valid in a stricter sense.


title: Effective General Relativistic Description of Jamming in Granular Matter

abstract: We propose here that certain observational features of granular matter in the infrared limit, exhibiting the phenomenon of {\it jamming}, arise from an underlying effective general relativistic description. The proposal stems from the assumption (which we justify on physical grounds) that grains in granular matter move freely in an {\it effective} curved Riemannian space. The termination of their trajectories at the onset of jamming is obtained from the focussing of a converging congruence of geodesics in such a space, as a solution of the Raychaudhuri equation for such congruences. This may happen irrespective of whether or not the curvature is sourced by external stresses (via an effective Einstein equation), although the properties of the resultant jammed state solution do differ in the two cases. A definite prediction of this geometrical approach is the negative role played by those trajectories which twist about each other, in reaching the jammed state. The local symmetries of granular interaction, translational and rotational invariance (corresponding to `force balance' and `torque balance' in standard force-based approaches to jamming) are inherent in the effective general relativity framework. A recently-proposed effective elasticity model of the jammed state, based on a tensorial variant of standard electrostatics (Vector Charge Theory), is seen to be entirely subsumed within the linearized version of the effective general relativistic description.


title: The solutions of Dirac equation on the hyperboloid under perpendicular magnetic fields

abstract: In this study, firstly it is reviewed how the solutions of the Dirac-Weyl equation for a massless charge on the hyperboloid under perpendicular magnetic fields are obtained by using supersymmetric (SUSY) quantum mechanics methods. Then, the solutions of the Dirac equation for a massive charge under magnetic fields have been computed in terms of the solutions which were found before for the Dirac-Weyl equation. As an example, the case of a constant magnetic field on the hyperbolic surface for massless and massive charges has been worked out.


title: The Devil is in the Details: Evaluating Limitations of Transformer-based Methods for Granular Tasks

abstract: Contextual embeddings derived from transformer-based neural language models have shown state-of-the-art performance for various tasks such as question answering, sentiment analysis, and textual similarity in recent years. Extensive work shows how accurately such models can represent abstract, semantic information present in text. In this expository work, we explore a tangent direction and analyze such models' performance on tasks that require a more granular level of representation. We focus on the problem of textual similarity from two perspectives: matching documents on a granular level (requiring embeddings to capture fine-grained attributes in the text), and an abstract level (requiring embeddings to capture overall textual semantics). We empirically demonstrate, across two datasets from different domains, that despite high performance in abstract document matching as expected, contextual embeddings are consistently (and at times, vastly) outperformed by simple baselines like TF-IDF for more granular tasks. We then propose a simple but effective method to incorporate TF-IDF into models that use contextual embeddings, achieving relative improvements of up to 36% on granular tasks.


title: Cross-entropy method in application to SIRC model

abstract: The study considers the usage of a probabilistic optimization method called Cross-Entropy (CE). This is the version of the Monte Carlo method created by Reuven Rubinstein (1997). It was developed in the context of determining rare events. Here we will present the way in which the CE method can be used for problems of optimization of epidemiological models, and more specifically the optimization of the SIRC (Susceptible - Infectious - Recovered - Cross-immune) model based on the functions supervising the care of specific groups in the model. With the help of weighted sampling, an attempt was made to find the fastest and most accurate version of the algorithm.


title: Tuning valleys and wave functions of van der Waals heterostructures by varying the number of layers: A first-principles study

abstract: In van der Waals heterostructures of two-dimensional transition-metal dichalcogenides (2D TMDCs) electron and hole states are spatially localized in different layers forming long-lived interlayer excitons. Here, we have investigated, from first principles, the influence of additional electron or hole layers on the electronic properties of a MoS2/WSe2 heterobilayer (HBL), which is a direct band gap material. Additional layers modify the interlayer hybridization, mostly affecting the quasiparticle energy and real-space extend of hole states at the G and electron states at the Q valleys. For a sufficient number of additional layers, the band edges move from K to Q or G, respectively. Adding electron layers to the HBL leads to more delocalized Q states, while G states do not extend much beyond the HBL, even when more hole layers are added. These results suggest a simple and yet powerful way to tune band edges and the real-space extend of the electron and hole wave function in TMDC heterostructures, strongly affecting the lifetime and dynamics of interlayer excitons.


title: High dimensional regimes of non-stationary Gaussian correlated Wishart matrices

abstract: We study the high-dimensional asymptotic regimes of correlated Wishart matrices $d^{-1}\mathcal{Y}\mathcal{Y}^T$, where $\mathcal{Y}$ is a $n\times d$ Gaussian random matrix with correlated and non-stationary entries. We prove that under different normalizations, two distinct regimes emerge as both $n$ and $d$ grow to infinity. The first regime is the one of central convergence, where the law of the properly renormalized Wishart matrices becomes close in Wasserstein distance to that of a Gaussian orthogonal ensemble matrix. In the second regime, a non-central convergence happens, and the law of the normalized Wishart matrices becomes close in Wasserstein distance to that of the so-called Rosenblatt-Wishart matrix recently introduced by Nourdin and Zheng. We then proceed to show that the convergences stated above also hold in a functional setting, namely as weak convergence in $C([a,b];M_n(\mathbb{R}))$. As an application of our main result (in the central convergence regime), we show that it can be used to prove convergence in expectation of the empirical spectral distributions of the Wishart matrices to the semicircular law. Our findings complement and extend a rich collection of results on the study of the fluctuations of Gaussian Wishart matrices, and we provide explicit examples based on Gaussian entries given by normalized increments of a bi-fractional or a sub-fractional Brownian motion.


title: Mathematical simulation of package delivery optimization using a combination of carriers

abstract: A variety of goods and services in the contemporary world requires permanent improvement of services e-commerce platform performance. Modern society is so deeply integrated with mail deliveries, purchasing of goods and services online, that makes competition between service and good providers a key selection factor. As long as logistic, timely, and cost-effective delivery plays important part authors decided to analyze possible ways of improvements in the current field, especially for regions distantly located from popular distribution centers. Considering both: fast and lazy delivery the factor of costs is playing an important role for each end-user. Given work proposes a simulation that analyses the current cost of delivery for e-commerce orders in the context of delivery by the Supplier Fleet, World-Wide delivery service fleet, and possible vendor drop-ship and checks of the alternative ways can be used to minimize the costs. The main object of investigation is focused around mid and small businesses living far from big distribution centers (except edge cases like lighthouses, edge rocks with very limited accessibility) but actively using e-commerce solutions for daily activities fulfillment. Authors analyzed and proposed a solution for the problem of cost optimization for packages delivery for long-distance deliveries using a combination of paths delivered by supplier fleets, worldwide and local carriers. Data models and Add-ons of contemporary Enterprise Resource Planning systems were used, and additional development is proposed in the perspective of the flow selection change. The experiment is based on data sources of the United States companies using a wide range of carriers for delivery services and uses the data sources of the real companies; however, it applies repetitive simulations to analyze variances in obtained solutions.


title: A Newton-Type Active Set Method for Nonlinear Optimization with Polyhedral Constraints

abstract: A Newton-type active set algorithm for large-scale minimization subject to polyhedral constraints is proposed. The algorithm consists of a gradient projection step, a second-order Newton-type step in the null space of the constraint matrix, and a set of rules for branching between the two steps. We show that the proposed method asymptotically takes the Newton step when the active constraints are linearly independent and a strong second-order sufficient optimality condition holds. We also show that the method has a quadratic rate of convergence under standard conditions. Numerical experiments are presented illustrating the performance of the algorithm on the CUTEst and on a specific class of problems for which finding second-order stationary points is critical.


title: Some finiteness results on triangular automorphisms

abstract: In this note we prove that every finite collection of connected algebraic subgroups of the group of triangular automorphisms of the affine space generates a connected solvable algebraic subgroup.


title: K-theoretic Hall algebras, quantum groups and super quantum groups

abstract: We first prove that the K-theoretic Hall algebra of a preprojective algebra of affine type is isomorphic to the positive half of a quantum toroidal quantum group. An essential step consists to deform the K-theoretic Hall algebra so that the deformation is torsion free over some polynomial subalgebra. Next, we compare super toroidal quantum groups of type A with K-theoretic Hall algebras of quivers with potential, which are defined using the Grothendieck groups of categories of singularities of some Landau-Ginzburg models.


title: Recent Advances in Two-Dimensional Metal Monochalcogenides

abstract: The family of emerging low-symmetry and structural in-plane anisotropic 2D materials have been expanding rapidly in recent years. As an important emerging anisotropic 2D material, the black phosphorene (BP) analog group IVA-VI metal monochalcogenides (MMCs) have been surged recently due to their distinctive crystalline symmetries, exotic in-plane anisotropic electronic and optical response, earth abundant and environmentally friendly characteristics. In this article, we review the recent research advancements in the field of anisotropic 2D MMCs. At first, the unique wavy crystal structures together with the optical and electronic properties of such materials are discussed. The review continues with the various methods adopted for the synthesis of layered MMCs including micromechanical and liquid phase exfoliation as well as physical vapor deposition. The last part of the article focuses on the application of the structural anisotropic response of 2D MMCs in field effect transistors, photovoltaic cells nonlinear optics and valleytronic devices. Besides presenting the significant research in the field of this emerging class of 2D materials, this review also delineates the existing limitations and discusses emerging possibilities and future prospects.


title: A Learning Theoretic Perspective on Local Explainability

abstract: In this paper, we explore connections between interpretable machine learning and learning theory through the lens of local approximation explanations. First, we tackle the traditional problem of performance generalization and bound the test-time accuracy of a model using a notion of how locally explainable it is. Second, we explore the novel problem of explanation generalization which is an important concern for a growing class of finite sample-based local approximation explanations. Finally, we validate our theoretical results empirically and show that they reflect what can be seen in practice.


title: Development of the Poincare cross-section method: Visualization the three-dimensional sections of four-dimensional flows

abstract: The theme of the article is the application of the Poincare section method for visual classification of attractors in the four-dimensional phase space; the purpose of the study is to introduce consideration of three-dimensional Poincare sections, and develop an algorithm for their use for classification of four-dimensional flows by the type of attractors. The article contains the obtained in three-dimensional Poincare sections traces of 4D attractors, allowing you to get an idea of their topology, as well as to conduct their primary classification. It is important that the model discussed in the article is formulated as universal, which allows one to study the evolution of systems (if they can be described iterative equations) consisting of subsystems interacting along arbitrary laws.


title: Tinker-HP : Accelerating Molecular Dynamics Simulations of Large Complex Systems with Advanced Point Dipole Polarizable Force Fields using GPUs and Multi-GPUs systems

abstract: We present the extension of the Tinker-HP package (Lagard\`ere et al., Chem. Sci., 2018,9, 956-972) to the use of Graphics Processing Unit (GPU) cards to accelerate molecular dynamics simulations using polarizable many-body force fields. The new high-performance module allows for an efficient use of single- and multi-GPU architectures ranging from research laboratories to modern supercomputer centers. After detailing an analysis of our general scalable strategy that relies on OpenACC and CUDA, we discuss the various capabilities of the package. Among them, the multi-precision possibilities of the code are discussed. If an efficient double precision implementation is provided to preserve the possibility of fast reference computations, we show that a lower precision arithmetic is preferred providing a similar accuracy for molecular dynamics while exhibiting superior performances. As Tinker-HP is mainly dedicated to accelerate simulations using new generation point dipole polarizable force field, we focus our study on the implementation of the AMOEBA model. Testing various NVIDIA platforms including 2080Ti, 3090, V100 and A100 cards, we provide illustrative benchmarks of the code for single- and multi-cards simulations on large biosystems encompassing up to millions of atoms. The new code strongly reduces time to solution and offers the best performances to date obtained using the AMOEBA polarizable force field. Perspectives toward the strong-scaling performance of our multi-node massive parallelization strategy, unsupervised adaptive sampling and large scale applicability of the Tinker-HP code in biophysics are discussed. The present software has been released in phase advance on GitHub in link with the High Performance Computing community COVID-19 research efforts and is free for Academics (see https://github.com/TinkerTools/tinker-hp).

categories: physics.comp-ph cs.DC cs.MS physics.chem-ph


title: Extended field theories are local and have classifying spaces

abstract: We show that all extended functorial field theories, both topological and nontopological, are local. We define the smooth (infinity,d)-category of bordisms with geometric data, such as Riemannian metrics or geometric string structures, and prove that it satisfies codescent with respect to the target S, which implies the locality property. We apply this result to construct a classifying space for concordance classes of functorial field theories with geometric data, solving a conjecture of Stolz and Teichner about the existence of such a space. We use our classifying space construction to develop a geometric theory of power operations, following the recent work of Barthel, Berwick-Evans, and Stapleton.


title: Increasing the achievable contrast of infrared interferometry with an error correlation model

abstract: Interferometric observables are strongly correlated, yet it is common practice to ignore these correlations in the data analysis process. We develop an empirical model for the correlations present in Very Large Telescope Interferometer GRAVITY data and show that properly accounting for them yields fainter detection limits and increases the reliability of potential detections. We extracted the correlations of the (squared) visibility amplitudes and the closure phases directly from intermediate products of the GRAVITY data reduction pipeline and fitted our empirical models to them. Then, we performed model fitting and companion injection and recovery tests with both simulated and real GRAVITY data, which are affected by correlated noise, and compared the results when ignoring the correlations and when properly accounting for them with our empirical models. When accounting for the correlations, the faint source detection limits improve by a factor of up to $\sim 2$ at angular separations $> 20~\rm{mas}$. For commonly used detection criteria based on $\chi^2$ statistics, this mostly results in claimed detections being more reliable. Ignoring the correlations present in interferometric data is a dangerous assumption which might lead to a large number of false detections. The commonly used detection criteria (e.g. in the model fitting pipeline CANDID) are only reliable when properly accounting for the correlations; furthermore, instrument teams should work on providing full covariance matrices instead of statistically independent error bars as part of the official data reduction pipelines.


title: Focus on the present: a regularization method for the ASR source-target attention layer

abstract: This paper introduces a novel method to diagnose the source-target attention in state-of-the-art end-to-end speech recognition models with joint connectionist temporal classification (CTC) and attention training. Our method is based on the fact that both, CTC and source-target attention, are acting on the same encoder representations. To understand the functionality of the attention, CTC is applied to compute the token posteriors given the attention outputs. We found that the source-target attention heads are able to predict several tokens ahead of the current one. Inspired by the observation, a new regularization method is proposed which leverages CTC to make source-target attention more focused on the frames corresponding to the output token being predicted by the decoder. Experiments reveal stable improvements up to 7\% and 13\% relatively with the proposed regularization on TED-LIUM 2 and LibriSpeech.


title: A Program for Multi-Messenger Standard Siren Cosmology in the Era of LIGO A+, Rubin Observatory, and Beyond

abstract: The most promising variation of the standard siren technique combines gravitational-wave (GW) data for binary neutron star (BNS) mergers with redshift measurements enabled by their electromagnetic (EM) counterparts, to constrain cosmological parameters such as $H_0$, $\Omega_m$, and $w_0$. Here we evaluate the near- and long-term prospects of multi-messenger cosmology in the era of future GW observatories: Advanced LIGO Plus (A+, 2025), Voyager-like detectors (2030s), and Cosmic Explorer-like detectors (CE, 2035 and beyond). We show that the BNS horizon distance of $\approx 700$ Mpc for A+ is well-matched to the sensitivity of the Vera C. Rubin Observatory (VRO) for kilonova detections. We find that one year of joint A+ and VRO observations will constrain the value of $H_0$ to percent-level precision, given a small investment of VRO time dedicated to target-of-opportunity GW follow-up. In the Voyager era, the BNS-kilonova observations begin to constrain $\Omega_m$ with an investment of a few percent of VRO time. With the larger BNS horizon distance in the Cosmic Explorer era, on-axis short gamma-ray bursts (SGRBs) and their afterglows (though accompanying only some of the GW-detected mergers) supplant kilonovae as the most promising counterparts for redshift identification. We show that five years of joint observations with Cosmic Explorer-like facilities and a next-generation gamma-ray satellite with localization capabilities similar to that presently possible with Swift could constrain both $\Omega_m$ and $w_0$ to $15-20\%$. We therefore advocate for a robust target-of-opportunity (ToO) program with VRO, and a wide-field gamma-ray satellite with improved sensitivity in the 2030s, to enable standard siren cosmology with next-generation gravitational wave facilities.


title: Lessons from the $B^{0,+}\to K^{*0,+}\mu^+\mu^-$ angular analyses

abstract: We perform an analysis within the Standard Model of $B^{0,+} \to K^{*0,+} \mu^+ \mu^-$ decays in light of the recent measurements from the LHCb experiment, showing that new data strengthen the need for sizable hadronic contributions and correlations among them. We then extend our analysis to New Physics via the Standard Model Effective Theory, and carry out a state-of-the-art fit of available $b \to s \ell^+ \ell^-$ data, including possible hadronic contributions. We find the case of a fully left-handed operator standing out as the simplest scenario with a significance of almost $6\sigma$.


title: Spin precession as a new window into disformal scalar fields

abstract: We launch a first investigation into how a light scalar field coupled both conformally and disformally to matter influences the evolution of spinning point-like bodies. Working directly at the level of the equations of motion, we derive novel spin-orbit and spin-spin effects accurate to leading order in a nonrelativistic and weak-field expansion. Crucially, unlike the spin-independent effects induced by the disformal coupling, which have been shown to vanish in circular binaries due to rotational symmetry, the spin-dependent effects we study here persist even in the limit of zero eccentricity, and so provide a new and qualitatively distinct way of probing these kinds of interactions. To illustrate their potential, we confront our predictions with spin-precession measurements from the Gravity Probe B experiment and find that the resulting constraint improves upon existing bounds from perihelion precession by over 5 orders of magnitude. Our results therefore establish spin effects as a promising window into the disformally coupled dark sector.


title: Employing Partial Least Squares Regression with Discriminant Analysis for Bug Prediction

abstract: Forecasting defect proneness of source code has long been a major research concern. Having an estimation of those parts of a software system that most likely contain bugs may help focus testing efforts, reduce costs, and improve product quality. Many prediction models and approaches have been introduced during the past decades that try to forecast bugged code elements based on static source code metrics, change and history metrics, or both. However, there is still no universal best solution to this problem, as most suitable features and models vary from dataset to dataset and depend on the context in which we use them. Therefore, novel approaches and further studies on this topic are highly necessary. In this paper, we employ a chemometric approach - Partial Least Squares with Discriminant Analysis (PLS-DA) - for predicting bug prone Classes in Java programs using static source code metrics. To our best knowledge, PLS-DA has never been used before as a statistical approach in the software maintenance domain for predicting software errors. In addition, we have used rigorous statistical treatments including bootstrap resampling and randomization (permutation) test, and evaluation for representing the software engineering results. We show that our PLS-DA based prediction model achieves superior performances compared to the state-of-the-art approaches (i.e. F-measure of 0.44-0.47 at 90% confidence level) when no data re-sampling applied and comparable to others when applying up-sampling on the largest open bug dataset, while training the model is significantly faster, thus finding optimal parameters is much easier. In terms of completeness, which measures the amount of bugs contained in the Java Classes predicted to be defective, PLS-DA outperforms every other algorithm: it found 69.3% and 79.4% of the total bugs with no re-sampling and up-sampling, respectively.


title: Continuous and Diverse Image-to-Image Translation via Signed Attribute Vectors

abstract: Recent image-to-image (I2I) translation algorithms focus on learning the mapping from a source to a target domain. However, the continuous translation problem that synthesizes intermediate results between two domains has not been well-studied in the literature. Generating a smooth sequence of intermediate results bridges the gap of two different domains, facilitating the morphing effect across domains. Existing I2I approaches are limited to either intra-domain or deterministic inter-domain continuous translation. In this work, we present an effectively signed attribute vector, which enables continuous translation on diverse mapping paths across various domains. In particular, we introduce a unified attribute space shared by all domains that utilize the sign operation to encode the domain information, thereby allowing the interpolation on attribute vectors of different domains. To enhance the visual quality of continuous translation results, we generate a trajectory between two sign-symmetrical attribute vectors and leverage the domain information of the interpolated results along the trajectory for adversarial training. We evaluate the proposed method on a wide range of I2I translation tasks. Both qualitative and quantitative results demonstrate that the proposed framework generates more high-quality continuous translation results against the state-of-the-art methods.


title: On cosmological expansion and local physics

abstract: We find an exact convergence in the local dynamics described by two supposedly antagonistic approaches in modern cosmology: one starting from an expanding universe perspective such as FLRW, the other based on a local model ignoring any notion of expansion, such as Schwarzschild dS. Both models are in complete agreement when the local effects of the expansion are circumscribed to the presence of the cosmological constant. We elaborate in the relevant role of static backgrounds like the Schwarzschild-dS metric in standard form as the most proper coordinatizations to describe physics at the local scale. Finally, making use of an old and too often forgotten relativistic kinematical invariant, we clarify some widespread misunderstandings on space expansion, cosmological and gravitational redshifts. As a byproduct we propose a {\sl unique and unambiguous} prescription to match the local and cosmological expression of a specific observable.


title: Prediction against a limited adversary

abstract: We study the problem of prediction with expert advice with adversarial corruption where the adversary can at most corrupt one expert. Using tools from viscosity theory, we characterize the long-time behavior of the value function of the game between the forecaster and the adversary. We provide lower and upper bounds for the growth rate of regret without relying on a comparison result. We show that depending on the description of regret, the limiting behavior of the game can significantly differ.


title: Asymptotic Theory of Expectile Neural Networks

abstract: Neural networks are becoming an increasingly important tool in applications. However, neural networks are not widely used in statistical genetics. In this paper, we propose a new neural networks method called expectile neural networks. When the size of parameter is too large, the standard maximum likelihood procedures may not work. We use sieve method to constrain parameter space. And we prove its consistency and normality under nonparametric regression framework.


title: Estimating County-Level COVID-19 Exponential Growth Rates Using Generalized Random Forests

abstract: Rapid and accurate detection of community outbreaks is critical to address the threat of resurgent waves of COVID-19. A practical challenge in outbreak detection is balancing accuracy vs. speed. In particular, while estimation accuracy improves with longer fitting windows, speed degrades. This paper presents a machine learning framework to balance this tradeoff using generalized random forests (GRF), and applies it to detect county level COVID-19 outbreaks. This algorithm chooses an adaptive fitting window size for each county based on relevant features affecting the disease spread, such as changes in social distancing policies. Experiment results show that our method outperforms any non-adaptive window size choices in 7-day ahead COVID-19 outbreak case number predictions.


title: Data taking network for COMET Phase-I

abstract: An experiment to search for mu-e conversion named COMET is being constructed at J-PARC. The experiment will be carried out using a two-stage approach of Phase-I and Phase-II. The data taking system of Phase-I is developed based on common network technology. The data taking system consists of two kinds of networks. One is a front-end network. Its network bundles around twenty front-end devices that have a 1-Gb optical network port. And a front-end computer accepts data from the devices via its network. The other is a back-end network that collects all event fragments from the front-end computers using a 10-Gb network. We used a low price 1Gb/10Gb optical network switch for the front-end network. And direct connection between the front-end PC and an event building PC using 10-Gb optical network devices was used for the back-end network. The event building PC has ten 10-Gb network ports. And each network port of the event building PC is connected to the front-end PC's port without using a network switch. We evaluated data taking performance with an event building on these two kinds of networks. The event building throughput of the front-end network achieved 337 MiB/s. And the event building throughput of the back-end networks achieved 1.2 GiB/s. It means that we could reduce the construction cost of the data taking network using this structure without deteriorating performance. Moreover, we evaluated the writing speed of the local storage RAID disk system connected to a back-end PC by a SAS interface, and a long-distance network copy from the experiment location to the lasting storage.


title: Interference traps waves in open system: Bound states in the continuum

abstract: I review the four mechanisms of bound states in the continuum (BICs) in application to microwave and acoustic cavities open to directional waveguides. The most simple are the symmetry protected BICs which are localized inside the cavity because of the orthogonality of the eigenmodes to the propagating modes of waveguides. However, the most general and interesting is the Friedrich-Wintgen mechanism when the BICs are result of full destructive interference of outgoing resonant modes. The third type of the BICs, the Fabry-Perot BICs, occur in a double resonator system when each resonator can serve as an ideal mirror. At last, the accidental BICs can be realized in the open cavities with no symmetry like the open Sinai billiard in which the eigenmode of the resonator can become orthogonal to the continuum of the waveguide accidentally by a smooth deformation of the eigenmode. We also review the one-dimensional systems in which the BICs occur owing to full destructive interference of two waves separated by spin or polarization or by paths in the Aharonov-Bohm rings. We widely use the method of effective non-Hermitian Hamiltonian equivalent to the coupled mode theory which detects bound states in the continuum (BICs) by finding zero widths resonances.


title: A Low-Power Time-to-Digital Converter for the CMS Endcap Timing Layer (ETL) Upgrade

abstract: We present the design and test results of a Time-to-Digital-Converter (TDC). The TDC will be a part of the readout ASIC, called ETROC, to read out Low-Gain Avalanche Detectors (LGADs) for the CMS Endcap Timing Layer (ETL) of High-Luminosity LHC upgrade. One of the challenges of the ETROC design is that the TDC is required to consume less than 200 W for each pixel at the nominal hit occupancy of 1%. To meet the low-power requirement, we use a single delay line for both the Time of Arrival (TOA) and the Time over Threshold (TOT) measurements without delay control. A double-strobe self-calibration scheme is used to compensate for process variation, temperature, and power supply voltage. The TDC is fabricated in a 65 nm CMOS technology. The overall performances of the TDC have been evaluated. The TOA has a bin size of 17.8 ps within its effective dynamic range of 11.6 ns. The effective measurement precision of the TOA is 5.6 ps and 9.9 ps with and without the nonlinearity correction, respectively. The TDC block consumes 97 W at the hit occupancy of 1%. Over a temperature range from 23 C to 78 C and a power supply voltage range from 1.05 V to 1.35 V (the nominal value of 1.20 V), the self-calibrated bin size of the TOA varies within 0.4%. The measured TDC performances meet the requirements except that more tests will be performed in the future to verify that the TDC complies with the radiation-tolerance specifications.


title: Comprehensible Counterfactual Explanation on Kolmogorov-Smirnov Test

abstract: The Kolmogorov-Smirnov (KS) test is popularly used in many applications, such as anomaly detection, astronomy, database security and AI systems. One challenge remained untouched is how we can obtain an explanation on why a test set fails the KS test. In this paper, we tackle the problem of producing counterfactual explanations for test data failing the KS test. Concept-wise, we propose the notion of most comprehensible counterfactual explanations, which accommodates both the KS test data and the user domain knowledge in producing explanations. Computation-wise, we develop an efficient algorithm MOCHE (for MOst CompreHensible Explanation) that avoids enumerating and checking an exponential number of subsets of the test set failing the KS test. MOCHE not only guarantees to produce the most comprehensible counterfactual explanations, but also is orders of magnitudes faster than the baselines. Experiment-wise, we present a systematic empirical study on a series of benchmark real datasets to verify the effectiveness, efficiency and scalability of most comprehensible counterfactual explanations and MOCHE.


title: A Lane-Changing Prediction Method Based on Temporal Convolution Network

abstract: Lane-changing is an important driving behavior and unreasonable lane changes can result in potentially dangerous traffic collisions. Advanced Driver Assistance System (ADAS) can assist drivers to change lanes safely and efficiently. To capture the stochastic time series of lane-changing behavior, this study proposes a temporal convolutional network (TCN) to predict the long-term lane-changing trajectory and behavior. In addition, the convolutional neural network (CNN) and recurrent neural network (RNN) methods are considered as the benchmark models to demonstrate the learning ability of the TCN. The lane-changing dataset was collected by the driving simulator. The prediction performance of TCN is demonstrated from three aspects: different input variables, different input dimensions and different driving scenarios. Prediction results show that the TCN can accurately predict the long-term lane-changing trajectory and driving behavior with shorter computational time compared with two benchmark models. The TCN can provide accurate lane-changing prediction, which is one key information for the development of accurate ADAS.


title: Beyond the acceptance limit of DRAGON: the case of the $\mathrm{^6Li(\alpha,\gamma)^{10}B}$ reaction

abstract: Radiative capture reactions play a pivotal role for our understanding of the origin of the elements in the cosmos. Recoil separators provide an effective way to study these reactions, in inverse kinematics, and take advantage of the use of radioactive ion beams. However, a limiting factor in the study of radiative capture reactions in inverse kinematics is the momentum spread of the product nuclei, which can result in an angular spread larger than the geometric acceptance of the separator. The DRAGON facility at TRIUMF is a versatile recoil separator, designed to study radiative capture reactions relevant to astrophysics in the $\mathrm{A \sim 10-30}$ region. In this work we present the first attempt to study with DRAGON a reaction, $\mathrm{^6Li(\alpha,\gamma)^{10}B}$ , for which the recoil angular spread exceeds DRAGON's acceptance. Our result is in good agreement with the literature value, showing that DRAGON can measure resonance strengths of astrophysically important reactions even when not all the recoils enter the separator.


title: Sample-efficient reinforcement learning using deep Gaussian processes

abstract: Reinforcement learning provides a framework for learning to control which actions to take towards completing a task through trial-and-error. In many applications observing interactions is costly, necessitating sample-efficient learning. In model-based reinforcement learning efficiency is improved by learning to simulate the world dynamics. The challenge is that model inaccuracies rapidly accumulate over planned trajectories. We introduce deep Gaussian processes where the depth of the compositions introduces model complexity while incorporating prior knowledge on the dynamics brings smoothness and structure. Our approach is able to sample a Bayesian posterior over trajectories. We demonstrate highly improved early sample-efficiency over competing methods. This is shown across a number of continuous control tasks, including the half-cheetah whose contact dynamics have previously posed an insurmountable problem for earlier sample-efficient Gaussian process based models.


title: Reducing the number of redundant pair-wise interactions in hydrodynamic meshless methods

abstract: Widely used Lagrangian numerical codes that compute the physical interaction with neighbouring resolution elements (particles), duplicate the calculation of the interaction between pairs of particles. We developed an algorithm that reduces the number of redundant calculations. The algorithm makes use of a hash function to flag already computed interactions and eventual collisions. The result of the hashing is stored in two caches. Without limiting the cache memory usage, all duplicated calculations can be avoided, achieving the speed-up of a factor on two. We show that, limiting the cache size (in bits) to double the typical number of neighbouring particles, 70 per cent of the redundant calculations can be avoided, yielding a speed-up of almost 35 per cent.


title: Dynamics of Droplet Generation from Corneal Tear Film during Non-contact Eye Procedure in the Context of COVID-19

abstract: Non-invasive medical diagnostics demonstrate a propensity for droplet generation and should be studied to devise risk mitigation strategies against the spread of the SARS-CoV-2 virus. We investigate the air-puff tonometry, which uses a short-timed air-puff to applanate the human eye in a bid to detect the early onset of glaucoma by measuring the intraocular pressure. The air-puff consists of a vortex trailed by a high-speed jet. High-speed imaging of the eye during a typical tonometry measurement reveals a sequence of events starting with the interaction between the tear layer and the air puff leading to an initial sheet ejection. It is immediately followed by the trailing jet applanating the central corneal section, causing capillary waves to form and interact with the highly 3D transient expanding sheet. Such interaction with the capillary waves and the surrounding airfield due to the trailing jet causes the expanding sheet to undergo bag breakup, finger formation by Rayleigh Taylor instability and further break up into subsequent droplets by Rayleigh Plateau instability, which eventually splashes onto nearby objects, potentially forming fomites or aerosols which can lead to infections. The complex spatiotemporal phenomenon is carefully documented by rigorous experiments and corroborated using comprehensive theoretical analyses.


title: Evolution of the chemical enrichment and the Mass-Metallicity relation in CALIFA galaxies

abstract: We use fossil record techniques on the CALIFA sample to study how galaxies in the local universe have evolved in terms of their chemical content. We show how the metallicity and the mass-metallicity relation (MZR) evolve through time for the galaxies in our sample and how this evolution varies when we divide them based on their mass, morphology and star-forming status. We also check the impact of measuring the metallicity at the centre or the outskirts. We find the expected results that the most massive galaxies got enriched faster, with the MZR getting steeper at higher redshifts. However, once we separate the galaxies into morphology bins this behaviour is not as clear, which suggests that morphology is a primary factor to determine how fast a galaxy gets enriched, with mass determining the amount of enrichment. We also find that star-forming galaxies appear to be converging in their chemical evolution, that is, the metallicity of star-forming galaxies of different mass is very similar at recent times compared to several Gyr ago.


title: Messung der Masse des Schwarzen Loches in LMC X-3 (Measuring the mass of the black hole in LMC X-3)

abstract: Because of its high mass, the X-ray binary system LMC X-3 is a candidate for a binary system hosting a black hole. The goal of this work is to determine the mass of this black hole and the orbital inclination of the system. Therefore, data from the GROND instrument at the MPI/ESO 2.2m telescope at La Silla, Chile that was recorded between Oktober 20, 2007 and December 13, 2013 was analyzed. The program 'XRbinary' from Edward L. Robinson was used to describe the data and invert for the black hole mass and the orbital inclination of the binary system. The mass of the black hole was determined to ($7.15\pm0.65$) solar masses and the inclination to ($68^{+2}_{-3}$) degrees. From these two quantities, estimations about the rotation of the black hole in LMC X-3 can be derived and the mass distribution of black holes in X-ray binaries can be expanded.


title: Parameter Estimation for RANS Models Using Approximate Bayesian Computation

abstract: We use approximate Bayesian computation (ABC) to estimate unknown parameter values, as well as their uncertainties, in Reynolds-averaged Navier-Stokes (RANS) simulations of turbulent flows. The ABC method approximates posterior distributions of model parameters, but does not require the direct computation, or estimation, of a likelihood function. Compared to full Bayesian analyses, ABC thus provides a faster and more flexible parameter estimation for complex models and a wide range of reference data. In this paper, we describe the ABC approach, including the use of a calibration step, adaptive proposal, and Markov chain Monte Carlo (MCMC) technique to accelerate the parameter estimation, resulting in an improved ABC approach, denoted ABC-IMCMC. As a test of the classic ABC rejection algorithm, we estimate parameters in a nonequilibrium RANS model using reference data from direct numerical simulations of periodically sheared homogeneous turbulence. We then demonstrate the use of ABC-IMCMC to estimate parameters in the Menter shear-stress-transport (SST) model using experimental reference data for an axisymmetric transonic bump. We show that the accuracy of the SST model for this test case can be improved using ABC-IMCMC, indicating that ABC-IMCMC is a promising method for the calibration of RANS models using a wide range of reference data.


title: Investigating Clumpy Galaxies in the Sloan Digital Sky Survey Stripe 82 using the Galaxy Zoo

abstract: Giant, star-forming clumps are a common feature prevalent amongst high-redshift star-forming galaxies and play a critical role in shaping their chaotic morphologies and yet, their nature and role in galaxy evolution remains to be fully understood. A majority of the effort to study clumps has been focused at high redshifts, and local clump studies have often suffered from small sample sizes. In this work, we present an analysis of clump properties in the local universe, and for the first time, performed with a statistically significant sample. With the help of the citizen science-powered Galaxy Zoo: Hubble project, we select a sample of 92 $z<0.06$ clumpy galaxies in Sloan Digital Sky Survey Stripe 82 galaxies. Within this sample, we identify 543 clumps using a contrast-based image analysis algorithm and perform photometry as well as estimate their stellar population properties. The overall properties of our $z<0.06$ clump sample are comparable to the high-redshift clumps. However, contrary to the high-redshift studies, we find no evidence of a gradient in clump ages or masses as a function of their galactocentric distances. Our results challenge the inward migration scenario for clump evolution for the local universe, potentially suggesting a larger contribution of ex-situ clumps and/or longer clump migration timescales.


title: Identifying resonances of the Galactic bar in Gaia DR2: II. Clues from angle space

abstract: The Milky Way disk exhibits intricate orbit substructure of still-debated dynamical origin. The angle variables $(\theta_\phi,\theta_R)$ -- which are conjugates to the actions $(L_z,J_R)$, and describe a star's location along its orbit -- are a powerful diagnostic to identify $l$:$m$ resonances via the orbit shape relation $\Delta \theta_R / \Delta \theta_\phi = -m/l$. In the past, angle signatures have been hidden by survey selection effects (SEs). Using test particle simulations of a barred galaxy, we demonstrate that \emph{Gaia} should allow us to identify the Galactic bar's Outer Lindblad Resonance ($l=+1, m=2$, OLR) in angle space. We investigate strategies to overcome SEs. In the angle data of the \emph{Gaia} DR2 RVS sample, we independently identify four candidates for the OLR and therefore for the pattern speed $\Omega_\text{bar}$. The strongest candidate, $\Omega_\text{bar}\sim1.4\Omega_0$, positions the OLR above the `Sirius' moving group, agrees with measurements from the Galactic center, and might be supported by higher-order resonances around the `Hercules/Horn'. But it misses the classic orbit orientation flip, as discussed in the companion study on actions. The candidate $\Omega_\text{bar}\sim1.2\Omega_0$ was also suggested by the action-based study, has the OLR at the `Hat', is consistent with \emph{slow bar} models, but still affected by SEs. Weaker candidates are $\Omega_\text{bar}=1.6\Omega_0$ and $1.74\Omega_0$. In addition, we show that the stellar angles do not support the `Hercules/Horn' being created by the OLR of a \emph{fast bar}. We conclude that -- to resolve if `Sirius' or `Hat' is related to the bar's OLR -- more complex dynamical explanations and more extended data with well-behaved SEs are required.


title: Constraining the Halo Mass of Damped Ly$\alpha$ Absorption Systems (DLAs) at $z=2-3.5$ using the Quasar-CMB Lensing Cross-correlation

abstract: We study the cross correlation of damped Ly$\alpha$ systems (DLAs) and their background quasars, using the most updated DLA catalog and the Planck 2018 CMB lensing convergence field. Our measurement suggests that the DLA bias $b_{\rm DLA}$ is smaller than $3.1$, corresponding to $\log(M/M_\odot h^{-1})\leq 12.3$ at a confidence of $90\%$. These constraints are broadly consistent with Alonso et al. (2018) and previous measurements by cross-correlation between DLAs and the Ly$\alpha$ forest (e.g. Font-Ribera et al. 2012; Perez-Rafols et al. 2018). Further, our results demonstrate the potential of obtaining a more precise measurement of the halo mass of high-redshift sources using next generation CMB experiments with a higher angular resolution. The python-based codes and data products of our analysis are available at https://github.com/LittleLin1999/CMB-lensingxDLA.


title: The Elephant in the Bathtub: when the physics of star formation regulate the baryon cycle of galaxies

abstract: In simple models of galaxy formation and evolution, star formation is solely regulated by the amount of gas present in the galaxy. However, it has recently been shown that star formation can be suppressed by galactic dynamics in galaxies that contain a dominant spheroidal component and a low gas fraction. This 'dynamical suppression' is hypothesised to also contribute to quenching gas-rich galaxies at high redshift, but its impact on the galaxy population at large remains unclear. In this paper, we assess the importance of dynamical suppression in the context of gas regulator models of galaxy evolution through hydrodynamic simulations of isolated galaxies, with gas-to-stellar mass ratios of 0.01-0.20 and a range of galactic gravitational potentials from disc-dominated to spheroidal. Star formation is modelled using a dynamics-dependent efficiency per free-fall time, which depends on the virial parameter of the gas. We find that dynamical suppression becomes more effective at lower gas fractions and quantify its impact on the star formation rate as a function of gas fraction and stellar spheroid mass surface density. We combine the results of our simulations with observed scaling relations that describe the change of galaxy properties across cosmic time, and determine the galaxy mass and redshift range where dynamical suppression may affect the baryon cycle. We predict that the physics of star formation can limit and regulate the baryon cycle at low redshifts ($z \lesssim 1.4$) and high galaxy masses ($M_{\ast} \gtrsim 3 \times 10^{10}~M_{\odot}$), where dynamical suppression can drive galaxies off the star formation main sequence.


title: Hot Jupiter and ultra-cold Saturn formation in dense star clusters

abstract: The discovery of high incidence of hot Jupiters in dense clusters challenges the field-based hot Jupiter formation theory. In dense clusters, interactions between planetary systems and flyby stars are relatively common. This has a significant impact on planetary systems, dominating hot Jupiter formation. In this paper, we perform high precision, few-body simulations of stellar flybys and subsequent planet migration in clusters. A large parameter space exploration demonstrates that close flybys that change the architecture of the planetary system can activate high eccentricity migration mechanisms: Lidov-Kozai and planet-planet scattering, leading to high hot Jupiter formation rate in dense clusters. Our simulations predict that many of the hot Jupiters are accompanied by "ultra-cold Saturns", expelled to apastra of thousands of AU. This increase is particularly remarkable for planetary systems originally hosting two giant planets with semi-major axis ratios $\sim$ 4 and the flyby star approaching nearly perpendicular to the planetary orbital plane. The estimated lower limit to the hot Jupiter formation rate of a virialized cluster is $\sim 1.6\times10^{-4}({\sigma}/{\rm 1kms^{-1}})^5({a_{\rm p}}/{\rm 20 AU})({M_{\rm c}}/{\rm 1000M_\odot})^{-2}$Gyr$^{-1}$ per star, where $\sigma$ is the cluster velocity dispersion, $a_{\rm p}$ is the size of the planetary system and $M_{\rm c}$ is the mass of the cluster. Our simulations yield a hot Jupiter abundance which is $\sim$ 50 times smaller than that observed in the old open cluster M67. We expect that interactions involving binary stars, as well as a third or more giant planets, will close the discrepancy.


title: Neutrino Masses and Mixing in Models with Large Extra Dimensions and Localized Fermions

abstract: Using a low-energy effective field theory approach, we study some properties of models with large extra dimensions, in which quarks and leptons have localized wave functions in the extra dimensions. We consider models with two types of gauge groups: (i) the Standard-Model gauge group, and (ii) the left-right symmetric (LRS) gauge group. Our main focus is on the lepton sector of models with $n=2$ extra dimensions, in particular, neutrino masses and mixing. We analyze the requisite conditions that the models must satisfy to be in accord with data and present a solution for lepton wave functions in the extra dimensions that fulfills these conditions. As part of our work, we also present a new solution for quark wave function centers. Issues with flavor-changing neutral current effects are assessed. Finally, we remark on baryogenesis and dark matter in these models.


title: Observational insights on the origin of giant low surface brightness galaxies

abstract: Giant low surface brightness galaxies (gLSBGs) with dynamically cold stellar discs reaching the radius of 130 kpc challenge currently considered galaxy formation mechanisms. We analyse new deep long-slit optical spectroscopic observations, archival optical images and published HI and optical spectroscopic data for a sample of seven gLSBGs, for which we performed mass modelling and estimated the parameters of dark matter haloes assuming the Burkert dark matter density profile. Our sample is not homogeneous by morphology, parameters of stellar populations and total mass, however, six of seven galaxies sit on the high-mass extension of the baryonic Tully-Fisher relation. In UGC 1382 we detected a global counterrotation of the stellar high surface brightness (HSB) disc with respect to the extended LSB disc. In UGC 1922 with signatures of a possible merger, the gas counterrotation is seen in the inner disc. Six galaxies host active galactic nuclei, three of which have the estimated black hole masses substantially below those expected for their (pseudo-)bulge properties suggesting poor merger histories. Overall, the morphology, internal dynamics, and low star formation efficiency in the outer discs indicate that the three formation scenarios shape gLSBGs: (i) a two-stage formation when an HSB galaxy is formed first and then grows an LSB disc by accreting gas from an external supply; (ii) an unusual shallow and extended dark matter halo; (iii) a major merger with fine-tuned orbital parameters and morphologies of the merging galaxies.


title: Supersymmetry and Quantum Computation

abstract: The interplay between supersymmetry and classical and quantum computation is discussed. First, it is shown that the problem of computing the Witten index of $\mathcal N \leq 2$ quantum mechanical systems is $\#P$-complete and therefore intractable. Then, the notions of supersymmetry in the space of qubits and supersymmetric quantum circuits are introduced and some of their properties discussed. In particular, it is shown that these define a nontrivial subclass of quantum algorithms with robustness properties typical of supersymmetric systems. Concrete examples, including the supersymmetric SYK model and fermion hard-core models are discussed. Some applications and open questions are suggested.


title: New Pathways to the Relic Abundance of Vector-Portal Dark Matter

abstract: We fully explore the thermal freezeout histories of a vector-portal dark matter model, in the region of parameter space in which the ratio of masses of the dark photon $A^{\prime}$ and dark matter $\chi$ is in the range $1 \lesssim m_{A^{\prime}}/m_{\chi} \lesssim 2$. In this region $2 \rightarrow 2$ and $3 \rightarrow 2$ annihilation processes within the dark sector, as well as processes that transfer energy between the dark sector and the Standard Model, play important roles in controlling the thermal freezeout of the dark matter. We carefully track the temperatures of all species, relaxing the assumption of previous studies that the dark and Standard Model sectors remain in thermal equilibrium throughout dark matter freezeout. Our calculations reveal a rich set of novel pathways which lead to the observed relic density of dark matter, and we develop a simple analytic understanding of these different regimes. The viable parameter space in our model provides a target for future experiments searching for light (MeV-GeV) dark matter, and includes regions where the dark matter self-interaction cross section is large enough to affect the small-scale structure of galaxies.


title: The outer disc in shambles: blind detection of Monoceros and ACS with Gaia's astrometric sample

abstract: The astrometric sample of Gaia allows us to study the outermost Galactic disc, the halo and their interface. It is precisely at the very edge of the disc where the effects of external perturbations are expected to be the most noticeable. Our goal is to detect the kinematic substructure present in the halo and at the edge of the Milky Way (MW) disc, and provide observational constraints on their phase-space distribution. We download, one HEALpix at a time, the proper motion histogram of distant stars, to which we apply a Wavelet Transformation to reveal the significant overdensities. We then analyse the large coherent structures that appear in the sky. We reveal a sharp yet complex anticentre dominated by Monoceros (MNC) and the Anticentre Stream (ACS) in the north, which we find with an intensity comparable to the Magellanic clouds and the Sagittarius stream, and by MNC south and TriAnd at negative latitudes. Our method allows us to perform a morphological analysis of MNC and ACS, both spanning more than 100$^\circ$ in longitude, and to provide a high purity sample of giants with which we track MNC down to latitudes as low as $\sim$5$^\circ$. Their colour-magnitude diagram is consistent with extended structures at a distance of $\sim$10-11 kpc originated in the disc, with a very low ratio of RR Lyrae over M giants, and kinematics compatible with the rotation curve at those distances or only slightly slower. We present a precise characterisation of MNC and ACS, two previously known structures that our method reveals naturally, allowing us to detect them without limiting ourselves to a particular stellar type and, for the first time, using only kinematics. Our results allow future studies to model their chemo-dynamics and evolution, thus constraining some of the most influential processes that shaped the MW.


title: Ensuring Robustness in Training Set Based Global 21-cm Cosmology Analysis

abstract: We present a methodology for ensuring the robustness of our analysis pipeline in separating the global 21-cm hydrogen cosmology signal from large systematics based on singular value decomposition (SVD) of training sets. We show how traditional goodness-of-fit metrics such as the $\chi^2$ statistic that assess the fit to the full data may not be able to detect a suboptimal extraction of the 21-cm signal when it is fit alongside one or more additional components due to significant covariance between them. However, we find that comparing the number of SVD eigenmodes for each component chosen by the pipeline for a given fit to the distribution of eigenmodes chosen for synthetic data realizations created from training set curves can detect when one or more of the training sets is insufficient to optimally extract the signal. Furthermore, this test can distinguish which training set (e.g. foreground, 21-cm signal) needs to be modified in order to better describe the data and improve the quality of the 21-cm signal extraction. We also extend this goodness-of-fit testing to cases where a prior distribution derived from the training sets is applied and find that, in this case, the $\chi^2$ statistic as well as the recently introduced $\psi^2$ statistic are able to detect inadequacies in the training sets due to the increased restrictions imposed by the prior. Crucially, the tests described in this paper can be performed when analyzing any type of observations with our pipeline.


title: Classification of Periodic Variable Stars with Novel Cyclic-Permutation Invariant Neural Networks

abstract: Neural networks (NNs) have been shown to be competitive against state-of-the-art feature engineering and random forest (RF) classification of periodic variable stars. Although previous work utilising NNs has made use of periodicity by period folding multiple-cycle time-series into a single cycle -- from time-space to phase-space -- no approach to date has taken advantage of the fact that network predictions should be invariant to the initial phase of the period-folded sequence. Initial phase is exogenous to the physical origin of the variability and should thus be factored out. Here, we present cyclic-permutation invariant networks, a novel class of NNs for which invariance to phase shifts is guaranteed through polar coordinate convolutions, which we implement by means of "Symmetry Padding." Across three different datasets of variable star light curves, we show that two implementations of the cyclic-permutation invariant network: the iTCN and the iResNet, consistently outperform non-invariant baselines and reduce overall error rates by between 4% to 22%. Over a 10-class OGLE-III sample, the iTCN/iResNet achieves an average per-class accuracy of 93.4%/93.3%, compared to RNN/RF accuracies of 70.5%/89.5% in a recent study using the same data. Finding improvement on a non-astronomy benchmark, we suggest that the methodology introduced here should also be applicable to a wide range of science domains where periodic data abounds due to physical symmetries.

categories: astro-ph.IM astro-ph.SR cs.LG physics.data-an


title: Nonuniversal interstellar density spectra probed by pulsars

abstract: The Galactic interstellar turbulence affects the density distribution and star formation. We introduce a new method of measuring interstellar turbulent density spectra by using the dispersion measures (DMs) of a large sample of pulsars. Without the need of invoking multiple tracers, we obtain nonuniversal density spectra in the multi-phase interstellar medium over different ranges of length scales. By comparing the analytical structure function of DMs with the observationally measured one in different areas of sky, we find a shallow density spectrum arising from the supersonic turbulence in cold interstellar phases, and a Kolmogorov-like density spectrum in the diffuse warm ionized medium (WIM). Both spectra extend up to hundreds of pc. On larger scales, we for the first time identify a steep density spectrum in the diffuse WIM extending up to several kpc. Our results show that the DMs of pulsars can provide unique new information on the interstellar turbulence.


title: Atmospheric Rossiter-McLaughlin effect and transmission spectroscopy of WASP-121b with ESPRESSO

abstract: WASP-121b is one of the most studied Ultra-hot Jupiters: many recent analyses of its atmosphere report interesting features at different wavelength ranges. In this paper we analyze one transit of WASP-121b acquired with the high-resolution spectrograph ESPRESSO at VLT in 1-telescope mode, and one partial transit taken during the commissioning of the instrument in 4-telescope mode. We investigate the anomalous in-transit radial velocity curve and study the transmission spectrum of the planet. By analysing the in-transit radial velocities we were able to infer the presence of the atmospheric Rossiter-McLaughlin effect. We measured the height of the planetary atmospheric layer that correlates with the stellar mask (mainly Fe) to be 1.052$\pm$0.015 Rp and we also confirmed the blueshift of the planetary atmosphere. By examining the planetary absorption signal on the stellar cross-correlation functions we confirmed the presence of a temporal variation of its blueshift during transit, which could be investigated spectrum-by-spectrum. We detected significant absorption in the transmission spectrum for Na, H, K, Li, CaII, Mg, and we certified their planetary nature by using the 2D tomographic technique. Particularly remarkable is the detection of Li, with a line contrast of $\sim$0.2% detected at the 6$\sigma$ level. With the cross-correlation technique we confirmed the presence of FeI, FeII, CrI and VI. H$\alpha$ and CaII are present up to very high altitudes in the atmosphere ($\sim$1.44 Rp and $\sim$2 Rp, respectively), and also extend beyond the transit-equivalent Roche lobe radius of the planet. These layers of the atmosphere have a large line broadening that is not compatible with being caused by the tidally-locked rotation of the planet alone, and could arise from vertical winds or high-altitude jets in the evaporating atmosphere.


title: Massive black hole binary systems and the NANOGrav 12.5 year results

abstract: The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has recently reported evidence for the presence of a common stochastic signal across their array of pulsars. The origin of this signal is still unclear. One of the possibilities is that it is due to a stochastic gravitational wave background (SGWB) in the $\sim 1-10\,{\rm nHz}$ frequency region. Taking the NANOGrav observational result at face value, we show that this signal would be fully consistent with a SGWB produced by an unresolved population of in-spiralling massive black hole binaries (MBHBs) predicted by current theoretical models. Considering an astrophysically agnostic model we find that the MBHB merger rate is loosely constrained to the range $10^{-11} - 2$ $\mathrm{Mpc}^{-3}\,\mathrm{Gyr}^{-1}$. Including additional constraints from galaxy pairing fractions and MBH-bulge scaling relations, we find that the MBHB merger rate is $10^{-5} - 5\times10^{-4}$ $\mathrm{Mpc}^{-3}\,\mathrm{Gyr}^{-1}$, the MBHB merger time-scale is $\le 3\,\mathrm{Gyr}$ and the norm of the $M_\mathrm{BH}-M_\mathrm{bulge}$ relation $\ge 1.2\times 10^{8}\,M_\odot$ (all intervals quoted at 90\% confidence). Regardless of the astrophysical details of MBHB assembly, this result would imply that a sufficiently large population of massive black holes pair up, form binaries and merge within a Hubble time.

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


title: Entanglement of formation of mixed many-body quantum states via Tree Tensor Operators

abstract: We present a numerical strategy to efficiently estimate bipartite entanglement measures, and in particular the Entanglement of Formation, for many-body quantum systems on a lattice. Our approach exploits the Tree Tensor Operator tensor network ansatz, a positive loopless representation for density matrices which, as we demonstrate, efficiently encodes information on bipartite entanglement, enabling the up-scaling of entanglement estimation. Employing this technique, we observe a finite-size scaling law for the entanglement of formation in 1D critical lattice models at finite temperature for up to 128 spins, extending to mixed states the scaling law for the entanglement entropy.


title: Bouncing Cosmology in f(Q) Symmetric Teleparallel Gravity

abstract: We consider f(Q) extended symmetric teleparallel cosmologies, where Q is the non-metricity scalar, and constrain its functional form through the order reduction method. By using this technique, we are able to reduce and integrate the field equations and thus to select the corresponding models giving rise to bouncing cosmology. The selected Lagrangian is then used to develop the Hamiltonian formalism and to obtain the Wave Function of the Universe which suggests that classical observable universes can be recovered according to the Hartle Criterion.


title: Universal AdS black holes in theories with sixteen supercharges and their microstates

abstract: We provide a universal microscopic counting for the microstates of the asymptotically AdS black holes and black strings that arise as solutions of the half-maximal gauged supergravity in four and five dimensions. These solutions can be embedded in all M-theory and type II string backgrounds with an AdS vacuum and sixteen supercharges and provide an infinite set of examples dual to $\mathcal{N}=2$ and $\mathcal{N}=4$ conformal field theories in four and three dimensions, respectively. The counting is universal and it is performed by either studying the large $N$ limit of the relevant supersymmetric index of the dual field theory or by using the charged Cardy formula.


title: Highly turbulent gas on GMC-scales in NGC 3256, the nearest luminous infrared galaxy

abstract: We present the highest resolution CO (2-1) observations obtained to date (0.25") of NGC 3256 and use them to determine the detailed properties of the molecular interstellar medium in the central 6 kpc of this merger. Distributions of physical quantities are reported from pixel-by-pixel measurements at 55 and 120 pc scales and compared to disc galaxies observed by PHANGS-ALMA. Mass surface densities range from 8 to 5500 M$_{\odot}$ pc$^{-2}$ and velocity dispersions from 10 to 200 km s$^{-1}$. Peak brightness temperatures as large as 37 K are measured, indicating the gas in NGC 3256 may be hotter than all regions in nearby disc galaxies measured by PHANGS-ALMA. Brightness temperatures even surpass those in the overlap region of NGC 4038/9 at the same scales. The majority of the gas appears unbound with median virial parameters of 7 to 19, although external pressure may bind some of the gas. High internal turbulent pressures of 10$^{5}$ to 10$^{10}$ K cm$^{-3}$ are found. Given the lack of significant trends in surface density, brightness temperature, and velocity dispersion with physical scale we argue the molecular gas is made up of a smooth medium down to 55 pc scales, unlike the more structured medium found in the PHANGS-ALMA disc galaxies.


title: NSF Convergence Approach to Transition Basic Research into Practice

abstract: The National Science Foundation Convergence Accelerator addresses national-scale societal challenges through use-inspired convergence research. Leveraging a convergence approach the Convergence Accelerator builds upon basic research and discovery to make timely investments to strengthen the Nations innovation ecosystem associated with several key R&D priority areas and practices to include the coronavirus disease 2019, harnessing the data revolution, the future of work, and quantum technology. Artificial Intelligence is a key underlying theme across all of these areas.


title: Aspects of quantum information in finite density field theory

abstract: We study different aspects of quantum field theory at finite density using methods from quantum information theory. For simplicity we focus on massive Dirac fermions with nonzero chemical potential, and work in $1+1$ space-time dimensions. Using the entanglement entropy on an interval, we construct an entropic $c$-function that is finite. Unlike what happens in Lorentz-invariant theories, this $c$-function exhibits a strong violation of monotonicity; it also encodes the creation of long-range entanglement from the Fermi surface. Motivated by previous works on lattice models, we next calculate numerically the Renyi entropies and find Friedel-type oscillations; these are understood in terms of a defect operator product expansion. Furthermore, we consider the mutual information as a measure of correlation functions between different regions. Using a long-distance expansion previously developed by Cardy, we argue that the mutual information detects Fermi surface correlations already at leading order in the expansion. We also analyze the relative entropy and its Renyi generalizations in order to distinguish states with different charge and/or mass. In particular, we show that states in different superselection sectors give rise to a super-extensive behavior in the relative entropy. Finally, we discuss possible extensions to interacting theories, and argue for the relevance of some of these measures for probing non-Fermi liquids.

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


title: Transmission in a Fano-Anderson chain with a topological defect

abstract: The Fano-Anderson chain consists of a linear lattice with a discrete side-unit, and exhibits Fano-resonant scattering due to coupling between the discrete states of the side-unit with the tight-binding continuum. We study Fano-resonance-assisted transport for the case of a topologically non-trivial side unit. We find that the topology of the side unit influences the transmission characteristics which thus can be an effective detection tool of the topological phases of the side unit. Furthermore, we explore the role of dual links between the linear tight-binding chain and the side unit. The secondary connection between the main chain and the side unit can modify the position or width of the Fano resonance dip in the transmission probability, and thus yield additional control.


title: Uncovering Majorana nature through a precision measurement of $CP$ phase

abstract: We show the possibility to discover the neutrino nature by measuring the Majorana CP phase at the DUNE experiment. This phase is turned on by a decoherence environment, possibly originated by physics at the Planck scale. A sizable distortion in the measurement of the Dirac CP violation phase $\delta_{\mathrm{CP}}$ is observed at DUNE when compared with T2HK measurement due to decoherence and non-null Majorana phase. Being that, when the measurement of the Majorana phase is performed at DUNE, it reaches a precision of 23 (21) $\%$ for a decoherence parameter $\Gamma=4.5(5.5)\times 10^{-24} \mathrm{GeV}$ and a Majorana phase equal to $1.5 \pi$. The latter precision is similar to the one obtained at the T2K experiment at its current Dirac CP violation phase measurement.


title: Quantum gravitational decoherence from fluctuating minimal length and deformation parameter at the Planck scale

abstract: Schemes of gravitationally induced decoherence are being actively investigated as possible mechanisms for the quantum-to-classical transition. Here, we introduce a decoherence process due to quantum gravity effects. We assume a foamy quantum spacetime with a fluctuating minimal length coinciding on average with the Planck scale. Considering deformed canonical commutation relations with a fluctuating deformation parameter, we derive a Lindblad master equation that yields localization in energy space and decoherence times consistent with the currently available observational evidence. Compared to other schemes of gravitational decoherence, we find that the decoherence rate predicted by our model is extremal, being minimal in the deep quantum regime below the Planck scale and maximal in the mesoscopic regime beyond it. We discuss possible experimental tests of our model based on cavity optomechanics setups with ultracold massive molecular oscillators and we provide preliminary estimates on the values of the physical parameters needed for actual laboratory implementations.

categories: gr-qc cond-mat.stat-mech hep-th quant-ph


title: Regluing graphs of Free Groups

abstract: Answering a question due to Min, we prove that a finite graph of roses admits a regluing such that the resulting graph of roses has hyperbolic fundamental group.


title: Approximating the long time average of the density operator: Diagonal ensemble

abstract: For an isolated generic quantum system out of equilibrium, the long time average of observables is given by the diagonal ensemble, i.e. the mixed state with the same probability for energy eigenstates as the initial state but without coherences between different energies. In this work we present a method to approximate the diagonal ensemble using tensor networks. Instead of simulating the real time evolution, we adapt a filtering scheme introduced earlier in [Phys. Rev. B 101, 144305 (2020)] to this problem. We analyze the performance of the method on a non-integrable spin chain, for which we observe that local observables converge towards thermal values polynomially with the inverse width of the filter.


title: Site-to-Site Internet Traffic Control

abstract: Queues allow network operators to control traffic: where queues build, they can enforce scheduling and shaping policies. In the Internet today, however, there is a mismatch between where queues build and where control is most effectively enforced; queues build at bottleneck links that are often not under the control of the data sender. To resolve this mismatch, we propose a new kind of middlebox, called Bundler. Bundler uses a novel inner control loop between a sendbox (in the sender's site) and a receivebox (in the receiver's site) to determine the aggregate rate for the bundle, leaving the end-to-end connections and their control loops intact. Enforcing this sending rate ensures that bottleneck queues that would have built up from the bundle's packets now shift from the bottleneck to the sendbox. The sendbox then exercises control over its traffic by scheduling packets to achieve higher-level objectives. We have implemented Bundler in Linux and evaluated it with real-world and emulation experiments. We find that Bundler allows the sender-chosen policy to be effective: when configured to implement Stochastic Fairness Queueing (SFQ), it improves median flow completion time (FCT) by between 28% and 97% across various scenarios.


title: Optimal Measurement of Field Properties with Quantum Sensor Networks

abstract: We consider a quantum sensor network of qubit sensors coupled to a field $f(\vec{x};\vec{\theta})$ analytically parameterized by the vector of parameters $\vec\theta$. The qubit sensors are fixed at positions $\vec{x}_1,\dots,\vec{x}_d$. While the functional form of $f(\vec{x};\vec{\theta})$ is known, the parameters $\vec{\theta}$ are not. We derive saturable bounds on the precision of measuring an arbitrary analytic function $q(\vec{\theta})$ of these parameters and construct the optimal protocols that achieve these bounds. Our results are obtained from a combination of techniques from quantum information theory and duality theorems for linear programming. They can be applied to many problems, including optimal placement of quantum sensors, field interpolation, and the measurement of functionals of parametrized fields.


title: The Choquard logarithmic equation involving a nonlinearity with exponential growth

abstract: In the present work we are concerned with the Choquard Logarithmic equation $-\Delta u + au + \lambda (\ln|\cdot|\ast |u|^{2})u = f(u)$ in $\mathbb{R}^2$, for $ a>0 $, $ \lambda >0 $ and a nonlinearity $f$ with exponential critical growth. We prove the existence of a nontrivial solution at the mountain pass level and a nontrivial ground state solution. Also, we provide these results under a symmetric setting, taking into account subgroups of $ O(2) $.


title: Realization of high-fidelity CZ and ZZ-free iSWAP gates with a tunable coupler

abstract: High-fidelity two-qubit gates at scale are a key requirement to realize the full promise of quantum computation and simulation. The advent and use of coupler elements to tunably control two-qubit interactions has improved operational fidelity in many-qubit systems by reducing parasitic coupling and frequency crowding issues. Nonetheless, two-qubit gate errors still limit the capability of near-term quantum applications. The reason, in part, is the existing framework for tunable couplers based on the dispersive approximation does not fully incorporate three-body multi-level dynamics, which is essential for addressing coherent leakage to the coupler and parasitic longitudinal ($ZZ$) interactions during two-qubit gates. Here, we present a systematic approach that goes beyond the dispersive approximation to exploit the engineered level structure of the coupler and optimize its control. Using this approach, we experimentally demonstrate CZ and $ZZ$-free iSWAP gates with two-qubit interaction fidelities of $99.76 \pm 0.07$% and $99.87 \pm 0.23$%, respectively, which are close to their $T_1$ limits.


title: Structural Signatures for Thermodynamic Stability in Vitreous Silica: Insight from Machine Learning and Molecular Dynamics Simulations

abstract: The structure-thermodynamic stability relationship in vitreous silica is investigated using machine learning and a library of 24,157 inherent structures generated from melt-quenching and replica exchange molecular dynamics simulations. We find the thermodynamic stability, i.e., enthalpy of the inherent structure ($e_{\mathrm{IS}}$), can be accurately predicted by both linear and nonlinear machine learning models from numeric structural descriptors commonly used to characterize disordered structures. We find short-range features become less indicative of thermodynamic stability below the fragile-to-strong transition. On the other hand, medium-range features, especially those between 2.8-~6 $\unicode{x212B}$;, show consistent correlations with $e_{\mathrm{IS}}$ across the liquid and glass regions, and are found to be the most critical to stability prediction among features from different length scales. Based on the machine learning models, a set of five structural features that are the most predictive of the silica glass stability is identified.