title: The Chemo-Dynamical Groups of Galactic Globular Clusters

abstract: We introduce a multi-component chemo-dynamical method for splitting the Galactic population of Globular Clusters (GCs) into three distinct constituents: bulge, disc, and stellar halo. The latter is further decomposed into the individual large accretion events that built up the Galactic stellar halo: the Gaia-Enceladus-Sausage, Kraken and Sequoia structures, and the Sagittarius and Helmi streams. Our modelling is extensively tested using mock GC samples constructed from the AURIGA suite of hydrodynamical simulations of Milky Way (MW)-like galaxies. We find that, on average, a proportion of the accreted GCs cannot be associated with their true infall group and are left ungrouped, biasing our recovered population numbers to approximately 80 percent of their true value. Furthermore, the identified groups have a completeness and a purity of only 65 percent. This reflects the difficulty of the problem, a result of the large degree of overlap in energy-action space of the debris from past accretion events. We apply the method to the Galactic data to infer, in a statistically robust and easily quantifiable way, the GCs associated with each MW accretion event. The resulting groups' population numbers of GCs, corrected for biases, are then used to infer the halo and stellar masses of the now defunct satellites that built up the halo of the MW.


title: The Uniform Infinite Cubic Planar Graph

abstract: We prove that the random simple cubic planar graph $\mathsf{C}_n$ with an even number $n$ of vertices admits a novel uniform infinite cubic planar graph (UICPG) as quenched local limit. We describe how the limit may be constructed by a series of random blow-up operations applied to the dual map of the type~III Uniform Infinite Planar Triangulation established by Angel and Schramm (Comm. Math. Phys., 2003). Our main technical lemma is a contiguity relation between $\mathsf{C}_n$ and a model where the networks inserted at the links of the largest $3$-connected component of $\mathsf{C}_n$ are replaced by independent copies of a specific Boltzmann network. We prove that the number of vertices of the largest $3$-connected component concentrates at $\kappa n$ for $\kappa \approx 0.85085$, with Airy-type fluctuations of order $n^{2/3}$. The second-largest component is shown to have significantly smaller size $O_p(n^{2/3})$.


title: Majorana bound states in encapsulated bilayer graphene

abstract: The search for robust topological superconductivity and Majorana bound states continues, exploring both one-dimensional (1D) systems such as semiconducting nanowires and two-dimensional (2D) platforms. In this work we study a 2D approach based on graphene bilayers encapsulated in transition metal dichalcogenides that, unlike previous proposals involving the Quantum Hall regime in graphene, requires weaker magnetic fields and does not rely on interactions. The encapsulation induces strong spin-orbit coupling on the graphene bilayer, which opens a sizeable gap and stabilizes fragile pairs of helical edge states. We show that, when subject to an in-plane Zeeman field, armchair edges can be transformed into p-wave one-dimensional topological superconductors by contacting them laterally with conventional superconductors. We demonstrate the emergence of Majorana bound states (MBSs) at the sample corners of crystallographically perfect flakes, belonging either to the D or the BDI symmetry classes depending on parameters. We compute the phase diagram, the resilience of MBSs against imperfections, and their manifestation as a 4$\pi$-periodic effect in Josephson junction geometries, all suggesting the existence of a topological phase within experimental reach.


title: Views, variety and celestial spheres

abstract: This paper describes a continuation of the program of causal views, in which the world consists of nothing but a vast number of partial views of its past. Each view is associated to an event, and is a representation of the immediate causal past of that event. These consists mainly of processes that transfer energy, momentum and other charges to it from its past events. There is fundamentally no space or spacetime, just a large number of events, which are the causes of events to come. This is a development of energetic causal set theories, developed with Marina Cortes. Momentum and energy are fundamental, and are conserved under their transformation from present events to future events. As a result Minkowski spacetime emerges, in a way that preserves causal relations. The locality of events as constructed in the emergent spacetime is a consequence of the conservation of energy-momentum fundamentally. In this paper we propose that the views of events can be represented in terms of degrees of freedom on punctured two surfaces-each puncture corresponding to an immediate past event. This opens up to the use of mathematical tools from twistor theory and the nascent area of celestial spheres. We show that This makes possible versions of the theory that are relativistically invariant.


title: A least squares support vector regression for anisotropic diffusion filtering

abstract: Anisotropic diffusion filtering for signal smoothing as a low-pass filter has the advantage of the edge-preserving, i.e., it does not affect the edges that contain more critical data than the other parts of the signal. In this paper, we present a numerical algorithm based on least squares support vector regression by using Legendre orthogonal kernel with the discretization of the nonlinear diffusion problem in time by the Crank-Nicolson method. This method transforms the signal smoothing process into solving an optimization problem that can be solved by efficient numerical algorithms. In the final analysis, we have reported some numerical experiments to show the effectiveness of the proposed machine learning based approach for signal smoothing.


title: Machine learning based modelling and optimization in hard turning of AISI D6 steel with newly developed AlTiSiN coated carbide tool

abstract: In recent times Mechanical and Production industries are facing increasing challenges related to the shift toward sustainable manufacturing. In this article, machining was performed in dry cutting condition with a newly developed coated insert called AlTiSiN coated carbides coated through scalable pulsed power plasma technique in dry cutting condition and a dataset was generated for different machining parameters and output responses. The machining parameters are speed, feed, depth of cut and the output responses are surface roughness, cutting force, crater wear length, crater wear width, and flank wear. The data collected from the machining operation is used for the development of machine learning (ML) based surrogate models to test, evaluate and optimize various input machining parameters. Different ML approaches such as polynomial regression (PR), random forest (RF) regression, gradient boosted (GB) trees, and adaptive boosting (AB) based regression are used to model different output responses in the hard machining of AISI D6 steel. The surrogate models for different output responses are used to prepare a complex objective function for the germinal center algorithm-based optimization of the machining parameters of the hard turning operation.


title: Positive definite maps on amenable groups

abstract: We describe conditions that characterize amenability for groups in terms of positive definite functions valued in a von Neumann algebra.


title: Combined Pruning for Nested Cross-Validation to Accelerate Automated Hyperparameter Optimization for Embedded Feature Selection in High-Dimensional Data with Very Small Sample Sizes

abstract: Background: Embedded feature selection in high-dimensional data with very small sample sizes requires optimized hyperparameters for the model building process. For this hyperparameter optimization, nested cross-validation must be applied to avoid a biased performance estimation. The resulting repeated training with high-dimensional data leads to very long computation times. Moreover, it is likely to observe a high variance in the individual performance evaluation metrics caused by outliers in tiny validation sets. Therefore, early stopping applying standard pruning algorithms to save time risks discarding promising hyperparameter sets. Result: To speed up feature selection for high-dimensional data with tiny sample size, we adapt the use of a state-of-the-art asynchronous successive halving pruner. In addition, we combine it with two complementary pruning strategies based on domain or prior knowledge. One pruning strategy immediately stops computing trials with semantically meaningless results for the selected hyperparameter combinations. The other is a new extrapolating threshold pruning strategy suitable for nested-cross-validation with a high variance of performance evaluation metrics. In repeated experiments, our combined pruning strategy keeps all promising trials. At the same time, the calculation time is substantially reduced compared to using a state-of-the-art asynchronous successive halving pruner alone. Up to 81.3\% fewer models were trained achieving the same optimization result. Conclusion: The proposed combined pruning strategy accelerates data analysis or enables deeper searches for hyperparameters within the same computation time. This leads to significant savings in time, money and energy consumption, opening the door to advanced, time-consuming analyses.


title: Quantum Machine Learning in Finance: Time Series Forecasting

abstract: We explore the efficacy of the novel use of parametrised quantum circuits (PQCs) as quantum neural networks (QNNs) for forecasting time series signals with simulated quantum forward propagation. The temporal signals consist of several sinusoidal components (deterministic signal), blended together with trends and additive noise. The performance of the PQCs is compared against that of classical bidirectional long short-term memory (BiLSTM) neural networks. Our results show that for time series signals consisting of small amplitude noise variations (up to 40 per cent of the amplitude of the deterministic signal) PQCs, with only a few parameters, perform similar to classical BiLSTM networks, with thousands of parameters, and outperform them for signals with higher amplitude noise variations. Thus, QNNs can be used effectively to model time series having, at the same time, the significant advantage of being trained significantly faster than a classical machine learning model in a quantum computer.


title: Dimension towers of SICs. II. Some constructions

abstract: A SIC is a maximal equiangular tight frame in a finite dimensional Hilbert space. Given a SIC in dimension $d$, there is good evidence that there always exists an aligned SIC in dimension $d(d-2)$, having predictable symmetries and smaller equiangular tight frames embedded in them. We provide a recipe for how to calculate sets of vectors in dimension $d(d-2)$ that share these properties. They consist of maximally entangled vectors in certain subspaces defined by the numbers entering the $d$ dimensional SIC. However, the construction contains free parameters and we have not proven that they can always be chosen so that one of these sets of vectors is a SIC. We give some worked examples that, we hope, may suggest to the reader how our construction can be improved. For simplicity we restrict ourselves to the case of odd dimensions.


title: BEA-Base: A Benchmark for ASR of Spontaneous Hungarian

abstract: Hungarian is spoken by 15 million people, still, easily accessible Automatic Speech Recognition (ASR) benchmark datasets - especially for spontaneous speech - have been practically unavailable. In this paper, we introduce BEA-Base, a subset of the BEA spoken Hungarian database comprising mostly spontaneous speech of 140 speakers. It is built specifically to assess ASR, primarily for conversational AI applications. After defining the speech recognition subsets and task, several baselines - including classic HMM-DNN hybrid and end-to-end approaches augmented by cross-language transfer learning - are developed using open-source toolkits. The best results obtained are based on multilingual self-supervised pretraining, achieving a 45% recognition error rate reduction as compared to the classical approach - without the application of an external language model or additional supervised data. The results show the feasibility of using BEA-Base for training and evaluation of Hungarian speech recognition systems.


title: Meta-Learning Hypothesis Spaces for Sequential Decision-making

abstract: Obtaining reliable, adaptive confidence sets for prediction functions (hypotheses) is a central challenge in sequential decision-making tasks, such as bandits and model-based reinforcement learning. These confidence sets typically rely on prior assumptions on the hypothesis space, e.g., the known kernel of a Reproducing Kernel Hilbert Space (RKHS). Hand-designing such kernels is error prone, and misspecification may lead to poor or unsafe performance. In this work, we propose to meta-learn a kernel from offline data (Meta-KeL). For the case where the unknown kernel is a combination of known base kernels, we develop an estimator based on structured sparsity. Under mild conditions, we guarantee that our estimated RKHS yields valid confidence sets that, with increasing amounts of offline data, become as tight as those given the true unknown kernel. We demonstrate our approach on the kernelized bandit problem (a.k.a.~Bayesian optimization), where we establish regret bounds competitive with those given the true kernel. We also empirically evaluate the effectiveness of our approach on a Bayesian optimization task.


title: Lyapunov Functions and Stability Analysis of Fractional-Order Systems

abstract: This study presents new estimates for fractional derivatives without singular kernels defined by some specific functions. Based on obtained inequalities, we give a useful method to establish the global stability of steady states for fractional-order systems and generalize some works existing in the literature. Finally, we apply our results to prove the global stability of a fractional-order SEIR model with a general incidence rate.


title: Entangling free electrons and optical excitations

abstract: The inelastic interaction between flying particles and optical nanocavities gives rise to entangled states in which some excitations of the latter are paired with changes in the energy or momentum of the former. In particular, entanglement of free electrons and nanocavity modes opens appealing opportunities associated with the strong interaction capabilities of the electrons. However, the degree of entanglement that is currently achievable by electron interaction with optical cavities is severely limited by the lack of external selectivity over the resulting state mixtures. Here, we propose a scheme to generate pure entanglement between designated optical excitations in a cavity and separable free-electron states. Specifically, we shape the electron wave-function profile to dramatically reduce the number of accessible cavity modes and simultaneously associate them with targeted electron scattering directions. We exemplify this concept through a theoretical description of free-electron entanglement with degenerate and nondegenerate plasmon modes in silver nanoparticles as well as atomic vibrations in an inorganic molecule. The generated entanglement can be further propagated through its electron component to extend quantum interactions beyond currently explored protocols.


title: Bayesian Persuasion Meets Mechanism Design: Going Beyond Intractability with Type Reporting

abstract: Bayesian persuasion studies how an informed sender should partially disclose information so as to influence the behavior of self-interested receivers. In the last years, a growing attention has been devoted to relaxing the assumption that the sender perfectly knows receiver's payoffs. The first crucial step towards such an achievement is to study settings where each receiver's payoffs depend on their unknown type, which is randomly determined by a known finite-supported probability distribution. This begets considerable computational challenges, as computing a sender-optimal signaling scheme is inapproximable up to within any constant factor. In this work, we circumvent this issue by leveraging ideas from mechanism design. In particular, we introduce a type reporting step in which the receiver is asked to report their type to the sender, after the latter has committed to a menu defining a signaling scheme for each possible receiver's type. We prove that, with a single receiver, the addition of this type reporting stage makes the sender's computational problem tractable. Then, we extend our framework to settings with multiple receivers, focusing on the case of no inter-agent externalities and binary actions. We show that it is possible to find a sender-optimal solution in polynomial-time by means of the ellipsoid method, given access to a suitable polynomial-time separation oracle. This can be implemented for supermodular and anonymous sender's utility functions. As for the case of submodular sender's utility functions, we first approximately cast the sender's problem into a linearly-constrained mathematical program whose objective function is the multi-linear extension of the sender's utility. Then, we show how to find in polynomial-time an approximate solution to the program by means of a continuous greedy algorithm. This provides a (1 -1/e)-approximation to the problem.


title: Signal Quality Assessment of Photoplethysmogram Signals using Quantum Pattern Recognition and lightweight CNN Architecture

abstract: Photoplethysmography (PPG) signal comprises physiological information related to cardiorespiratory health. However, while recording, these PPG signals are easily corrupted by motion artifacts and body movements, leading to noise enriched, poor quality signals. Therefore ensuring high-quality signals is necessary to extract cardiorespiratory information accurately. Although there exists several rule-based and Machine-Learning (ML) - based approaches for PPG signal quality estimation, those algorithms' efficacy is questionable. Thus, this work proposes a lightweight CNN architecture for signal quality assessment employing a novel Quantum pattern recognition (QPR) technique. The proposed algorithm is validated on manually annotated data obtained from the University of Queensland database. A total of 28366, 5s signal segments are preprocessed and transformed into image files of 20 x 500 pixels. The image files are treated as an input to the 2D CNN architecture. The developed model classifies the PPG signal as `good' or `bad' with an accuracy of 98.3% with 99.3% sensitivity, 94.5% specificity and 98.9% F1-score. Finally, the performance of the proposed framework is validated against the noisy `Welltory app' collected PPG database. Even in a noisy environment, the proposed architecture proved its competence. Experimental analysis concludes that a slim architecture along with a novel Spatio-temporal pattern recognition technique improve the system's performance. Hence, the proposed approach can be useful to classify good and bad PPG signals for a resource-constrained wearable implementation.


title: Associative detachment in anion-atom reactions involving a dipole-bound electron

abstract: Associative electronic detachment (AED) between anions and neutral atoms leads to the detachment of the anion's electron resulting in the formation of a neutral molecule. It plays a key role in chemical reaction networks, like the interstellar medium, the Earth's ionosphere and biochemical processes. Here, a class of AED involving a closed-shell anion (OH^-) and alkali atoms (rubidium) is investigated by precisely controlling the fraction of electronically excited rubidium. Reaction with the ground state atom gives rise to a stable intermediate complex with an electron solely bound via dipolar forces. The stability of the complex is governed by the subtle interplay of diabatic and adiabatic couplings into the autodetachment manifold. The measured rate coefficients are in good agreement with ab initio calculations, revealing pronounced steric effects. For excited state rubidium, however, a lower reaction rate is observed, indicating dynamical stabilization processes suppressing the coupling into the autodetachment region. Our work provides a stringent test of ab initio calculations on anion-neutral collisions and constitutes a generic, conceptual framework for understanding electronic state dependent dynamics in AEDs.


title: Progress on the Kundt conjecture

abstract: The Kundt conjecture states that a Lorentzian manifold of arbitrary dimension which is not characterized by its scalar polynomial curvature invariants (SPIs) allows for a non-twisting, non-shearing and non-expanding (in short, Kundt) null congruence of geodesics. The conjecture has been proven for dimensions 3 and 4. A necessary condition for a spacetime not to be characterized by SPIs is that all covariant derivatives of the Riemann tensor are of aligned type II or more special in the null alignment classification. In arbitrary dimensions, we prove that this property indeed requires the presence of a Kundt null congruence when a certain genericity condition holds, or when the trac-free Ricci orWeyl tensor is of genuine type III or N, thus confirming the validity of the Kundt conjecture in these cases. We also strenghten the results for dimensions 3 and 4 by removing regularity assumptions and showing that only the third covariant derivative is needed to obtain the results. A key tool in our proofs is a new bilinear map acting on tensors of related boost orders relative to a null direction.


title: Semantic of Cloud Computing services for Time Series workflows

abstract: Time series (TS) are present in many fields of knowledge, research, and engineering. The processing and analysis of TS are essential in order to extract knowledge from the data and to tackle forecasting or predictive maintenance tasks among others The modeling of TS is a challenging task, requiring high statistical expertise as well as outstanding knowledge about the application of Data Mining(DM) and Machine Learning (ML) methods. The overall work with TS is not limited to the linear application of several techniques, but is composed of an open workflow of methods and tests. These workflow, developed mainly on programming languages, are complicated to execute and run effectively on different systems, including Cloud Computing (CC) environments. The adoption of CC can facilitate the integration and portability of services allowing to adopt solutions towards services Internet Technologies (IT) industrialization. The definition and description of workflow services for TS open up a new set of possibilities regarding the reduction of complexity in the deployment of this type of issues in CC environments. In this sense, we have designed an effective proposal based on semantic modeling (or vocabulary) that provides the full description of workflow for Time Series modeling as a CC service. Our proposal includes a broad spectrum of the most extended operations, accommodating any workflow applied to classification, regression, or clustering problems for Time Series, as well as including evaluation measures, information, tests, or machine learning algorithms among others.


title: First-order Temporal Logic on Finite Traces: Semantic Properties, Decidable Fragments, and Applications

abstract: Formalisms based on temporal logics interpreted over finite strict linear orders, known in the literature as finite traces, have been used for temporal specification in automated planning, process modelling, (runtime) verification and synthesis of programs, as well as in knowledge representation and reasoning. In this paper, we focus on first-order temporal logic on finite traces. We first investigate preservation of equivalences and satisfiability of formulas between finite and infinite traces, by providing a set of semantic and syntactic conditions to guarantee when the distinction between reasoning in the two cases can be blurred. Moreover, we show that the satisfiability problem on finite traces for several decidable fragments of first-order temporal logic is ExpSpace-complete, as in the infinite trace case, while it decreases to NExpTime when finite traces bounded in the number of instants are considered. This leads also to new complexity results for temporal description logics over finite traces. Finally, we investigate applications to planning and verification, in particular by establishing connections with the notions of insensitivity to infiniteness and safety from the literature.


title: Generalisations of multiple zeta values to rooted forests

abstract: We show that any convergent (shuffle) arborified zeta value admits a series representation. This justifies the introduction of a new generalisation to rooted forests of multiple zeta values, and we study its algebraic properties. As a consequence of the series representation, we derive elementary proofs of some results of Bradley and Zhou for Mordell-Tornheim zeta values and give explicit formulas. The series representation for shuffle arborified zeta values also implies that they are conical zeta values. We characterise which conical zeta values are arborified zeta values and evaluate them as sums of multiple zeta values with rational coefficients.


title: Similarity Learning based Few Shot Learning for ECG Time Series Classification

abstract: Using deep learning models to classify time series data generated from the Internet of Things (IoT) devices requires a large amount of labeled data. However, due to constrained resources available in IoT devices, it is often difficult to accommodate training using large data sets. This paper proposes and demonstrates a Similarity Learning-based Few Shot Learning for ECG arrhythmia classification using Siamese Convolutional Neural Networks. Few shot learning resolves the data scarcity issue by identifying novel classes from very few labeled examples. Few Shot Learning relies first on pretraining the model on a related relatively large database, and then the learning is used for further adaptation towards few examples available per class. Our experiments evaluate the performance accuracy with respect to K (number of instances per class) for ECG time series data classification. The accuracy with 5- shot learning is 92.25% which marginally improves with further increase in K. We also compare the performance of our method against other well-established similarity learning techniques such as Dynamic Time Warping (DTW), Euclidean Distance (ED), and a deep learning model - Long Short Term Memory Fully Convolutional Network (LSTM-FCN) with the same amount of data and conclude that our method outperforms them for a limited dataset size. For K=5, the accuracies obtained are 57%, 54%, 33%, and 92% approximately for ED, DTW, LSTM-FCN, and SCNN, respectively.


title: CREAM: a Package to Compute [Auto, Endo, Iso, Mono, Epi]-morphisms, Congruences, Divisors and More for Algebras of Type $(2^n,1^n)$

abstract: The CREAM GAP package computes automorphisms, congruences, endomorphisms and subalgebras of algebras with an arbitrary number of binary and unary operations; it also decides if between two such algebras there exists a monomorphism, an epimorphism, an isomorphism or if one is a divisor of the other. Thus it finds those objects for almost all algebras used in practice (groups, quasigroups in their various signatures, semigroups possibly with many unary operations, fields, semi-rings, quandles, logic algebras, etc). As a one-size-fits-all package, it only relies on universal algebra theorems, without taking advantage of specific theorems about, eg, groups or semigroups to reduce the search space. Canon and Holt produced very fast code to compute automorphisms of groups that outperform CREAM on orders larger than 128. Similarly, Mitchell et al. take advantage of deep theorems to compute automorphisms and congruences of completely 0-simple semigroups in a very efficient manner. However these domains (groups of order above 128 and completely 0-simple semigroups) are among the very few examples of GAP code faster than our general purpose package CREAM. For the overwhelming majority of other classes of algebras, either ours is the first code computing the above mentioned objects, or the existing algorithms are outperformed by CREAM, in some cases by several orders of magnitude. To get this performance, CREAM uses a mixture of universal algebra algorithms together with GAP coupled with artificial intelligence theorem proving tools (AITP) and very delicate C implementations. As an example of the latter, we re-implement Freese's very clever algorithm for computing congruences in universal algebras, in a way that outperforms all other known implementations.


title: $C$-normality of rank-one perturbations of normal operators

abstract: For a separable complex Hilbert space $H$, we say that a bounded linear operator $T$ acting on $H$ is $C$-normal, where $C$ is a conjugation on $H$, if it satisfies $CT^*TC=TT^*$. For a normal operator, we give geometric conditions which guarantee that its rank-one perturbation is a $C$-normal for some conjugation $C$.


title: Electron Scattering Emission in the Light Curves of Stars with Centrifugal Magnetospheres

abstract: Strongly magnetic, rapidly rotating B-type stars with relatively weak winds form centrifugal magnetospheres (CMs), as the stellar wind becomes magnetically confined above the Kepler co-rotation radius. Approximating the magnetic field as a dipole tilted by an angle $\beta$ with respect to the rotation axis, the CM plasma is concentrated in clouds at and above the Kepler radius along the intersection of the rotational and magnetic equatorial planes. Stellar rotation can bring such clouds in front of the stellar disk, leading to absorption of order 0.1 magnitude ($\sim 10 \%$ of continuum flux). However some stars with prominent CMs, such as $\sigma$ Ori E, show an emission bump in addition to absorption dips, which has been so far unexplained. We show that emission can occur from electron scattering toward the observer when CM clouds are projected off the stellar limb. Using the Rigidly Rotating Magnetosphere model, modified with a centrifugal breakout density scaling, we present a model grid of photometric light curves spanning parameter space in observer inclination angle $i$, magnetic obliquity angle $\beta$, critical rotation fraction $W$, and optical depth at the Kepler radius $\tau_{\text{K}}$. We show that $\tau_{\text{K}}$ of order unity can produce emission bumps of the magnitude $\sim 0.05$ seen in $\sigma$ Ori E. We discuss the implications for modeling the light curves of CM stars, as well as future work for applying the radiative transfer model developed here to 3D MHD simulations of CMs.


title: Non-existence of cusps for degenerate Alt-Caffarelli functionals

abstract: We eliminate the existence of cusps in a class of \textit{degenerate} free-boundary problems for the Alt-Caffarelli functional $J_{Q}(v, \Omega):= \int_{\Omega}|\nabla v|^2 + Q^2(x)\chi_{\{v>0\}}dx,$ so-called because $Q(x) = \text{dist}(x, \Gamma)^{\gamma}$ for $\Gamma$ an affine $k$-plane and $0< \gamma$. This problem is inspired by a generalization of the variational formulation of the Stokes Wave by Arama and Leoni. The elimination of cusps implies that the results of [Mccurdy20] in fact describe the entire free-boundary as it intersects $\Gamma$.


title: A General, Evolution-Inspired Reward Function for Social Robotics

abstract: The field of social robotics will likely need to depart from a paradigm of designed behaviours and imitation learning and adopt modern reinforcement learning (RL) methods to enable robots to interact fluidly and efficaciously with humans. In this paper, we present the Social Reward Function as a mechanism to provide (1) a real-time, dense reward function necessary for the deployment of RL agents in social robotics, and (2) a standardised objective metric for comparing the efficacy of different social robots. The Social Reward Function is designed to closely mimic those genetically endowed social perception capabilities of humans in an effort to provide a simple, stable and culture-agnostic reward function. Presently, datasets used in social robotics are either small or significantly out-of-domain with respect to social robotics. The use of the Social Reward Function will allow larger in-domain datasets to be collected close to the behaviour policy of social robots, which will allow both further improvements to reward functions and to the behaviour policies of social robots. We believe this will be the key enabler to developing efficacious social robots in the future.


title: Penalized Estimation of Frailty-Based Illness-Death Models for Semi-Competing Risks

abstract: Semi-competing risks refers to the survival analysis setting where the occurrence of a non-terminal event is subject to whether a terminal event has occurred, but not vice versa. Semi-competing risks arise in a broad range of clinical contexts, with a novel example being the pregnancy condition preeclampsia, which can only occur before the `terminal' event of giving birth. Models that acknowledge semi-competing risks enable investigation of relationships between covariates and the joint timing of the outcomes, but methods for model selection and prediction of semi-competing risks in high dimensions are lacking. Instead, researchers commonly analyze only a single or composite outcome, losing valuable information and limiting clinical utility -- in the obstetric setting, this means ignoring valuable insight into timing of delivery after preeclampsia has onset. To address this gap we propose a novel penalized estimation framework for frailty-based illness-death multi-state modeling of semi-competing risks. Our approach combines non-convex and structured fusion penalization, inducing global sparsity as well as parsimony across submodels. We perform estimation and model selection via a pathwise routine for non-convex optimization, and prove the first statistical error bound results in this setting. We present a simulation study investigating estimation error and model selection performance, and a comprehensive application of the method to joint risk modeling of preeclampsia and timing of delivery using pregnancy data from an electronic health record.


title: New Characterizations of Core Imputations of Matching and $b$-Matching Games

abstract: We give new characterizations of core imputations for the following games: * The assignment game. * Concurrent games, i.e., general graph matching games having non-empty core. * The unconstrained bipartite $b$-matching game (edges can be matched multiple times). * The constrained bipartite $b$-matching game (edges can be matched at most once). The classic paper of Shapley and Shubik \cite{Shapley1971assignment} showed that core imputations of the assignment game are precisely optimal solutions to the dual of the LP-relaxation of the game. Building on this, Deng et al. \cite{Deng1999algorithms} gave a general framework which yields analogous characterizations for several fundamental combinatorial games. Interestingly enough, their framework does not apply to the last two games stated above. In turn, we show that some of the core imputations of these games correspond to optimal dual solutions and others do not. This leads to the tantalizing question of understanding the origins of the latter. We also present new characterizations of the profits accrued by agents and teams in core imputations of the first two games. Our characterization for the first game is stronger than that for the second; the underlying reason is that the characterization of vertices of the Birkhoff polytope is stronger than that of the Balinski polytope.


title: Double Copy for Massive Scalar Field Theories

abstract: We explore extensions of the double copy to massive theories and find a new cubic theory with a local double copy. We consider the nonlinear sigma model and the special galileon theory, massless versions of which are known to be related through the double copy. We show that by performing a Kaluza-Klein reduction of these theories from five dimensions down to four, a double copy relation exists between the resulting massive four-dimensional scalar field theories. This requires the vanishing contribution of new galileon terms arising in high dimensions. We further explore if other interactions that do not arise from a dimensional reduction of the nonlinear sigma model could be double copied and find a new cubic interaction which satisfies the BCJ relations up to 5-point amplitudes.


title: Observation of an isoscalar resonance with exotic $J^{PC}=1^{-+}$ quantum numbers in $J/\psi\rightarrow\gamma\eta\eta'$

abstract: Using a sample of (10.09$\pm$0.04)$\times$10$^{9}$ $J/\psi$ events collected with the BESIII detector operating at the BEPCII storage ring, a partial wave analysis of the decay $J/\psi \rightarrow \gamma\eta\eta'$ is performed. The first observation of an isoscalar state with exotic quantum numbers $J^{PC}=1^{-+}$, denoted as $\eta_1(1855)$, is reported in the process $J/\psi \rightarrow \gamma\eta_1(1855)$ with $\eta_1(1855)\rightarrow\eta\eta'$. Its mass and width are measured to be (1855$\pm$9$_{-1}^{+6}$)MeV/$c^{2}$ and (188$\pm$18$_{-8}^{+3}$)MeV, respectively, where the first uncertainties are statistical and the second are systematic, and its statistical significance is estimated to be larger than 19$\sigma$.


title: Datamodels: Predicting Predictions from Training Data

abstract: We present a conceptual framework, datamodeling, for analyzing the behavior of a model class in terms of the training data. For any fixed "target" example $x$, training set $S$, and learning algorithm, a datamodel is a parameterized function $2^S \to \mathbb{R}$ that for any subset of $S' \subset S$ -- using only information about which examples of $S$ are contained in $S'$ -- predicts the outcome of training a model on $S'$ and evaluating on $x$. Despite the potential complexity of the underlying process being approximated (e.g., end-to-end training and evaluation of deep neural networks), we show that even simple linear datamodels can successfully predict model outputs. We then demonstrate that datamodels give rise to a variety of applications, such as: accurately predicting the effect of dataset counterfactuals; identifying brittle predictions; finding semantically similar examples; quantifying train-test leakage; and embedding data into a well-behaved and feature-rich representation space. Data for this paper (including pre-computed datamodels as well as raw predictions from four million trained deep neural networks) is available at https://github.com/MadryLab/datamodels-data .


title: Partial wave analysis of $J/\psi\rightarrow\gamma\eta\eta'$

abstract: Based on a sample of (10.09$\pm$0.04)$\times$10$^{9}$ $J/\psi$ events collected with the BESIII detector operating at the BEPCII storage ring, a partial wave analysis of the decay $J/\psi \rightarrow \gamma\eta\eta'$ is performed. An isoscalar state with exotic quantum numbers $J^{PC}=1^{-+}$, denoted as $\eta_1(1855)$, has been observed for the first time with statistical significance larger than 19$\sigma$. Its mass and width are measured to be (1855$\pm$9$_{-1}^{+6}$)~MeV/$c^{2}$ and (188$\pm$18$_{-8}^{+3}$)~MeV, respectively. The product branching fraction ${\cal B}(J/\psi$$\rightarrow$$ \gamma\eta_1(1855)$$\rightarrow$$\gamma\eta\eta')$ is measured to be (2.70$\pm 0.41 _{-0.35}^{+0.16}) \times$10$^{-6}$. The first uncertainties are statistical and the second are systematic. In addition, an upper limit on the branching ratio ${\cal B}(f_0(1710)$$\rightarrow$$\eta\eta')$/${\cal B}(f_0(1710)$$\rightarrow$$\pi\pi)$ is determined to be $1.61 \times 10^{-3}$ at 90\% confidence level, which lends support to the hypothesis that the $f_0(1710)$ has a large glueball component.


title: Magnetohydrodynamic wave mode identification in circular and elliptical sunspot umbrae: evidence for high order modes

abstract: In this paper we provide clear direct evidence of multiple concurrent higher order magnetohydrodynamic (MHD) modes in circular and elliptical sunspots by applying both Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD) techniques on solar observational data. These techniques are well documented and validated in the areas of fluid mechanics, hydraulics, and granular flows, yet are relatively new to the field of solar physics. While POD identifies modes based on orthogonality in space and it provides a clear ranking of modes in terms of their contribution to the variance of the signal, DMD resolves modes that are orthogonal in time. The clear presence of the fundamental slow sausage and kink body modes, as well as higher order slow sausage and kink body modes have been identified using POD and DMD analysis of the chromospheric H$\alpha$ line at 6562.808~{\AA} for both the circular and elliptical sunspots. Additionally, to the various slow body modes, evidence for the presence of the fast surface kink mode was found in the circular sunspot. All the MHD modes patterns were cross-correlated with their theoretically predicted counterparts and we demonstrated that ellipticity cannot be neglected when interpreting MHD wave modes. The higher-order MHD wave modes are even more sensitive to irregularities in umbral cross-sectional shapes, hence this must be taken into account for more accurate modelling of the modes in sunspots and pores.


title: Black-box Bayesian inference for economic agent-based models

abstract: Simulation models, in particular agent-based models, are gaining popularity in economics. The considerable flexibility they offer, as well as their capacity to reproduce a variety of empirically observed behaviours of complex systems, give them broad appeal, and the increasing availability of cheap computing power has made their use feasible. Yet a widespread adoption in real-world modelling and decision-making scenarios has been hindered by the difficulty of performing parameter estimation for such models. In general, simulation models lack a tractable likelihood function, which precludes a straightforward application of standard statistical inference techniques. Several recent works have sought to address this problem through the application of likelihood-free inference techniques, in which parameter estimates are determined by performing some form of comparison between the observed data and simulation output. However, these approaches are (a) founded on restrictive assumptions, and/or (b) typically require many hundreds of thousands of simulations. These qualities make them unsuitable for large-scale simulations in economics and can cast doubt on the validity of these inference methods in such scenarios. In this paper, we investigate the efficacy of two classes of black-box approximate Bayesian inference methods that have recently drawn significant attention within the probabilistic machine learning community: neural posterior estimation and neural density ratio estimation. We present benchmarking experiments in which we demonstrate that neural network based black-box methods provide state of the art parameter inference for economic simulation models, and crucially are compatible with generic multivariate time-series data. In addition, we suggest appropriate assessment criteria for future benchmarking of approximate Bayesian inference procedures for economic simulation models.


title: Production of Fock Mixtures in Trapped Ions for Motional Metrology

abstract: We present a protocol to produce a class of non-thermal Fock state mixtures in trapped ions. This class of states features a clear metrological advantage with respect to the ground state, thus overcoming the standard quantum limit without the need for full sideband cooling and Fock-state preparation on a narrow electronic transition. The protocol consists in the cyclic repetition of red-sideband, measurement and preparation laser pulses. By means of the Kraus map representation of the protocol, it is possible to relate the length of the red sideband pulses to the specific class of states that can be generated. With the help of numerical simulations, we analyze the parametric regime where these states can be reliably reproduced.


title: Log-Concavity of Infinite Product Generating Functions

abstract: In the $1970$s Nicolas proved that the coefficients $p_d(n)$ defined by the generating function \begin{equation*} \sum_{n=0}^{\infty} p_d(n) \, q^n = \prod_{n=1}^{\infty} \left( 1- q^n\right)^{-n^{d-1}} \end{equation*} are log-concave for $d=1$. Recently, Ono, Pujahari, and Rolen have extended the result to $d=2$. Note that $p_1(n)=p(n)$ is the partition function and $p_2(n)=\func{pp}\left( n\right) $ is the number of plane partitions. In this paper, we invest in properties for $p_d(n)$ for general $d$. Let $n \geq 6$. Then $p_d(n)$ is almost log-concave for $n$ divisible by $3$ and almost strictly log-convex otherwise.


title: Regret Minimization with Performative Feedback

abstract: In performative prediction, the deployment of a predictive model triggers a shift in the data distribution. As these shifts are typically unknown ahead of time, the learner needs to deploy a model to get feedback about the distribution it induces. We study the problem of finding near-optimal models under performativity while maintaining low regret. On the surface, this problem might seem equivalent to a bandit problem. However, it exhibits a fundamentally richer feedback structure that we refer to as performative feedback: after every deployment, the learner receives samples from the shifted distribution rather than only bandit feedback about the reward. Our main contribution is an algorithm that achieves regret bounds scaling only with the complexity of the distribution shifts and not that of the reward function. The algorithm only relies on smoothness of the shifts and does not assume convexity. Moreover, its final iterate is guaranteed to be near-optimal. The key algorithmic idea is careful exploration of the distribution shifts that informs a novel construction of confidence bounds on the risk of unexplored models. More broadly, our work establishes a conceptual approach for leveraging tools from the bandits literature for the purpose of regret minimization with performative feedback.


title: Bayesian inference and prediction for mean-mixtures of normal distributions

abstract: We study frequentist risk properties of predictive density estimators for mean mixtures of multivariate normal distributions, involving an unknown location parameter $\theta \in \mathbb{R}^d$, and which include multivariate skew normal distributions. We provide explicit representations for Bayesian posterior and predictive densities, including the benchmark minimum risk equivariant (MRE) density, which is minimax and generalized Bayes with respect to an improper uniform density for $\theta$. For four dimensions or more, we obtain Bayesian densities that improve uniformly on the MRE density under Kullback-Leibler loss. We also provide plug-in type improvements, investigate implications for certain type of parametric restrictions on $\theta$, and illustrate and comment the findings based on numerical evaluations.


title: Covid-19 vaccine hesitancy and mega-influencers

abstract: Covid-19 vaccines are widely available in the United States, yet our Covid-19 vaccination rates have remained far below 100%. Not only that, but CDC data shows that even in places where vaccine acceptance was proportionally high at the outset of the Covid-19 vaccination effort, that willingness has not necessarily translated into high rates of vaccination over the subsequent months. We model how such a shift could have arisen, using parameters in agreement with data from the state of Alabama. The simulations suggest that in Alabama, local interactions would have favored the emergence of tight consensus around the initial majority view, which was to accept the Covid-19 vaccine. Yet this is not what happened. We therefore add to our model the impact of mega-influencers such as mass media, the governor of the state, etc. Our simulations show that a single vaccine-hesitant mega-influencer, reaching a large fraction of the population, can indeed cause the consensus to shift radically, from acceptance to hesitancy. Surprisingly this is true even when the mega-influencer only reaches individuals who are already somewhat inclined to agree with them, and under the conservative assumption that individuals give no more weight to the mega-influencer than they would give to a single one of their friends or neighbors. Our simulations also suggest that a competing mega-influencer with the opposite view can shift the mean population opinion back, but cannot restore the tightness of consensus around that view. Our code and data are distributed in the ODyN (Opinion Dynamic Networks) library available at https://github.com/annahaensch/ODyN.


title: FiNCAT: Financial Numeral Claim Analysis Tool

abstract: While making investment decisions by reading financial documents, investors need to differentiate between in-claim and outof-claim numerals. In this paper, we present a tool which does it automatically. It extracts context embeddings of the numerals using one of the transformer based pre-trained language model called BERT. After this, it uses a Logistic Regression based model to detect whether the numerals is in-claim or out-of-claim. We use FinNum-3 (English) dataset to train our model. After conducting rigorous experiments we achieve a Macro F1 score of 0.8223 on the validation set. We have open-sourced this tool and it can be accessed from https://github.com/sohomghosh/FiNCAT_Financial_Numeral_Claim_Analysis_Tool


title: Reducing the Amount of Real World Data for Object Detector Training with Synthetic Data

abstract: A number of studies have investigated the training of neural networks with synthetic data for applications in the real world. The aim of this study is to quantify how much real world data can be saved when using a mixed dataset of synthetic and real world data. By modeling the relationship between the number of training examples and detection performance by a simple power law, we find that the need for real world data can be reduced by up to 70% without sacrificing detection performance. The training of object detection networks is especially enhanced by enriching the mixed dataset with classes underrepresented in the real world dataset. The results indicate that mixed datasets with real world data ratios between 5% and 20% reduce the need for real world data the most without reducing the detection performance.


title: Efficient Policy Space Response Oracles

abstract: Policy Space Response Oracle methods (PSRO) provide a general solution to learn Nash equilibrium in two-player zero-sum games but suffer from two drawbacks: (1) the computation inefficiency due to the need for consistent meta-game evaluation via simulations, and (2) the exploration inefficiency due to finding the best response against a fixed meta-strategy at every epoch. In this work, we propose Efficient PSRO (EPSRO) that largely improves the efficiency of the above two steps. Central to our development is the newly-introduced subroutine of no-regret optimization on the unrestricted-restricted (URR) game. By solving URR at each epoch, one can evaluate the current game and compute the best response in one forward pass without the need for meta-game simulations. Theoretically, we prove that the solution procedures of EPSRO offer a monotonic improvement on the exploitability, which none of existing PSRO methods possess. Furthermore, we prove that the no-regret optimization has a regret bound of $\mathcal{O}(\sqrt{T\log{[(k^2+k)/2]}})$, where $k$ is the size of restricted policy set. Most importantly, a desirable property of EPSRO is that it is parallelizable, this allows for highly efficient exploration in the policy space that induces behavioral diversity. We test EPSRO on three classes of games, and report a 50x speedup in wall-time and 10x data efficiency while maintaining similar exploitability as existing PSRO methods on Kuhn and Leduc Poker games.


title: Experimental demonstration of Gaussian boson sampling with displacement

abstract: Gaussian boson sampling (GBS) is quantum sampling task in which one has to draw samples from the photon-number distribution of a large-dimensional nonclassical squeezed state of light. In an effort to make this task intractable for a classical computer, experiments building GBS machines have mainly focused on increasing the dimensionality and squeezing strength of the nonclassical light. However, no experiment has yet demonstrated the ability to displace the squeezed state in phase-space, which is generally required for practical applications of GBS. In this work, we build a GBS machine which achieves the displacement by injecting a laser beam alongside a two-mode squeezed vacuum state into a 15-mode interferometer. We focus on two new capabilities. Firstly, we use the displacement to reconstruct the multimode Gaussian state at the output of the interferometer. Our reconstruction technique is in situ and requires only three measurements settings regardless of the state dimension. Secondly, we study how the addition of classical laser light in our GBS machine affects the complexity of sampling its output photon statistics. We introduce and validate approximate semi-classical models which reduce the computational cost when a significant fraction of the detected light is classical.


title: A Simple Proof of Siegel's Theorem Using Mellin Transform

abstract: In this paper, we present a simple analytic proof of Siegel's theorem that concerns the lower bound of $L(1,\chi)$ for primitive quadratic $\chi$. Our new method compares an elementary lower bound with an analytic upper bound obtained by the inverse Mellin transform of $\Gamma(s)$.


title: Differentially Private Community Detection for Stochastic Block Models

abstract: The goal of community detection over graphs is to recover underlying labels/attributes of users (e.g., political affiliation) given the connectivity between users (represented by adjacency matrix of a graph). There has been significant recent progress on understanding the fundamental limits of community detection when the graph is generated from a stochastic block model (SBM). Specifically, sharp information theoretic limits and efficient algorithms have been obtained for SBMs as a function of $p$ and $q$, which represent the intra-community and inter-community connection probabilities. In this paper, we study the community detection problem while preserving the privacy of the individual connections (edges) between the vertices. Focusing on the notion of $(\epsilon, \delta)$-edge differential privacy (DP), we seek to understand the fundamental tradeoffs between $(p, q)$, DP budget $(\epsilon, \delta)$, and computational efficiency for exact recovery of the community labels. To this end, we present and analyze the associated information-theoretic tradeoffs for three broad classes of differentially private community recovery mechanisms: a) stability based mechanism; b) sampling based mechanisms; and c) graph perturbation mechanisms. Our main findings are that stability and sampling based mechanisms lead to a superior tradeoff between $(p,q)$ and the privacy budget $(\epsilon, \delta)$; however this comes at the expense of higher computational complexity. On the other hand, albeit low complexity, graph perturbation mechanisms require the privacy budget $\epsilon$ to scale as $\Omega(\log(n))$ for exact recovery. To the best of our knowledge, this is the first work to study the impact of privacy constraints on the fundamental limits for community detection.


abstract: In 1916, Riesz proved that the Riemann hypothesis is equivalent to the bound $\sum_{n=1}^\infty \frac{\mu(n)}{n^2} \exp\left( - \frac{x}{n^2} \right) = O_{\epsilon} \left( x^{-\frac{3}{4} + \epsilon} \right)$, as $x \rightarrow\infty$, for any $\epsilon >0$. Around the same time, Hardy and Littlewood gave another equivalent criteria for the Riemann hypothesis while correcting an identity of Ramanujan. In the present paper, we establish a one-variable generalization of the identity of Hardy and Littlewood and as an application, we provide Riesz-type criteria for the Riemann hypothesis. In particular, we obtain the bound given by Riesz as well as the bound of Hardy and Littlewood.


title: Effect of dispersants on bacterial colonization of oil droplets: a microfluidic approach

abstract: Bacteria biodegradation of immiscible oil requires cell-droplet encounters, surface attachment, and hydrocarbon metabolism. Chemical dispersants are applied to oil spills to reduce the mean dispersed droplet size, thereby increasing the available surface area for attachment, in attempts to facilitate bacterial biodegradation. However, their effectiveness remains contentious as studies have shown that dispersants can inhibit, enhance, or have no effect on biodegradation. Therefore, questions remain on whether dispersants affect surface attachment or cell viability. Here, using microfluidics and time-lapse microscopy, we directly observe the attachment and growth of the marine bacterium, \emph{Alcanivorax borkumensis}, on stationary crude oil droplets ($5$ \textmu m $< R < 150$ \textmu m) in the presence of Corexit 9500. We show that the average colonization time, or the time comprised of encounters, attachment, and growth, is dependent on droplet size and primarily driven by diffusive encounters. Our results suggest that dispersants do not inhibit or enhance these biophysical processes.


title: Doubly stochastic arrays with small support

abstract: An $n \times m$ non-negative matrix with row sum $m$ and column sum $n$ is called doubly stochastic. We answer the problem of finding doubly stochastic matrices of smallest posible support for every $1 <n \leq m$. Any matrix of minimum support is extremal in the sence of convexity, while examples of extremal matrices that are not of minimum support are given. But when $n,m$ are coprime integers extremal matrices are precisely those of minimum support.


title: Rewiring What-to-Watch-Next Recommendations to Reduce Radicalization Pathways

abstract: Recommender systems typically suggest to users content similar to what they consumed in the past. If a user happens to be exposed to strongly polarized content, she might subsequently receive recommendations which may steer her towards more and more radicalized content, eventually being trapped in what we call a "radicalization pathway". In this paper, we study the problem of mitigating radicalization pathways using a graph-based approach. Specifically, we model the set of recommendations of a "what-to-watch-next" recommender as a d-regular directed graph where nodes correspond to content items, links to recommendations, and paths to possible user sessions. We measure the "segregation" score of a node representing radicalized content as the expected length of a random walk from that node to any node representing non-radicalized content. High segregation scores are associated to larger chances to get users trapped in radicalization pathways. Hence, we define the problem of reducing the prevalence of radicalization pathways by selecting a small number of edges to "rewire", so to minimize the maximum of segregation scores among all radicalized nodes, while maintaining the relevance of the recommendations. We prove that the problem of finding the optimal set of recommendations to rewire is NP-hard and NP-hard to approximate within any factor. Therefore, we turn our attention to heuristics, and propose an efficient yet effective greedy algorithm based on the absorbing random walk theory. Our experiments on real-world datasets in the context of video and news recommendations confirm the effectiveness of our proposal.


title: The Cosmological Constant as Event Horizon

abstract: General Relativity allows for a cosmological constant ($\Lambda$) which has inspired models of cosmic Inflation and Dark Energy. We show instead that $r_\Lambda = \sqrt{3/\Lambda}$ corresponds to an event horizon: a causal boundary term in the action. Our Universe is expanding inside its Schwarzschild radius $r_S=\rL=2GM$, which could have originated from a uniform free falling cloud of mass $M$ that collapsed as a Black Hole (BH) 25 Gyrs ago. Such a BH Universe allows for large-scale structure formation without the need of Inflation or Dark Energy.


title: Parameter estimation for the FOU(p) process with the same lambda

abstract: The FOU(p) processes can be considered as an alternative to ARMA (or ARFIMA) processes to model time series. Also, there is no substantial loss when we model a time series using FOU(p) processes with the same lambda, than using differents values of lambda. In this work we propose a new method to estimate the unique value of lambda in a FOU(p) process. Under certain conditions, we will prove consistency and asymptotic normality. We will show that this new method is more easy and fast to compute. By simulations, we show that the new procedure work well and is more efficient than the general method. Also, we include an application to real data, and we show that the new method work well too and outperforms the family of ARMA(p, q).


title: Topological invariants for words of linear factor complexity

abstract: Given a finite alphabet $\Sigma$ and a right-infinite word $w$ over the alphabet $\Sigma$, we construct a topological space ${\rm Rec}(w)$ consisting of all right-infinite recurrent words whose factors are all factors of $w$, where we work up to an equivalence in which two words are equivalent if they have the exact same set of factors (finite contiguous subwords). We show that ${\rm Rec}(w)$ can be endowed with a natural topology and we show that if $w$ is word of linear factor complexity then ${\rm Rec}(w)$ is a finite topological space. In addition, we note that there are examples which show that if $f:\mathbb{N}\to \mathbb{N}$ is a function that tends to infinity as $n\to \infty$ then there is a word whose factor complexity function is ${\rm O}(nf(n))$ such that ${\rm Rec}(w)$ is an infinite set. Finally, we pose a realization problem: which finite topological spaces can arise as ${\rm Rec}(w)$ for a word of linear factor complexity?


title: Effective medium theory for second-gradient elasticity with chirality

abstract: We derive effective, parsimonious models from a heterogeneous second-gradient nonlinear elastic material taking into account chiral scale-size effects. Our classification of the effective equations depends on the hierarchy of four characteristic lengths: The size of the heterogeneities $\ell$, the intrinsic lengths of the constituents $\ell_{\rm SG}$ and $\ell_{\rm chiral}$, and the overall characteristic length of the domain ${\rm L}$. Depending on the different scale interactions between $\ell_{\rm SG}$, $\ell_{\rm chiral}$, $\ell$, and ${\rm L}$ we obtain either an effective Cauchy continuum or an effective second-gradient continuum. The working technique combines scaling arguments with the periodic homogenization asymptotic procedure. Both the passage to the homogenization limit and the unveiling of the correctors' structure rely on a suitable use of the periodic unfolding and related operators.


title: Generalization Analysis of Message Passing Neural Networks on Large Random Graphs

abstract: Message passing neural networks (MPNN) have seen a steep rise in popularity since their introduction as generalizations of convolutional neural networks to graph-structured data, and are now considered state-of-the-art tools for solving a large variety of graph-focused problems. We study the generalization error of MPNNs in graph classification and regression. We assume that graphs of different classes are sampled from different random graph models. We show that, when training a MPNN on a dataset sampled from such a distribution, the generalization gap increases in the complexity of the MPNN, and decreases, not only with respect to the number of training samples, but also with the average number of nodes in the graphs. This shows how a MPNN with high complexity can generalize from a small dataset of graphs, as long as the graphs are large. The generalization bound is derived from a uniform convergence result, that shows that any MPNN, applied on a graph, approximates the MPNN applied on the geometric model that the graph discretizes.


title: Spontaneous flow created by active topological defects

abstract: Topological defects are at the root of the large-scale organization of liquid crystals. In two-dimensional active nematics, two classes of topological defects of charges $\pm 1/2$ are known to play a major role due to active stresses. Despite this importance, few analytical results have been obtained on the flow-field and active-stress patterns around active topological defects. Using the generic hydrodynamic theory of active systems, we investigate the flow and stress patterns around these topological defects in unbounded, two-dimensional active nematics. Under generic assumptions, we derive analytically the spontaneous velocity and stall force of self-advected defects in the presence of both shear and rotational viscosities. Applying our formalism to the dynamics of monolayers of elongated cells at confluence, we show that the non-conservation of cell number generically increases the self-advection velocity and could provide an explanation for their observed role in cellular extrusion and multilayering. We finally investigate numerically the influence of the Ericksen stress. Our work paves the way to a generic study of the role of topological defects in active nematics, and in particular in monolayers of elongated cells.


title: Radial Kohn-Sham problem via integral-equation approach

abstract: We present a numerical tool for solving the non-relativistic Kohn-Sham problem for spherically-symmetric atoms. It treats the Schr\"{o}dinger equation as an integral equation relying heavily on convolutions. The solver supports different types of exchange-correlation functionals including screened and long-range corrected hybrids. We implement a new method for treating range separation based on the complementary error function kernel. The present tool is applied in non-relativistic total energy calculations of atoms. A comparison with ultra-precise reference data[Cinal, JOMC 58, 1571 (2020)] shows a 14-digit agreement for Hartree-Fock results. We provide further benchmark data obtained with 5 different exchange-correlation functionals.


title: Numerical Evidence for Exponential Speed-up of QAOA over Unstructured Search for Approximate Constrained Optimization

abstract: Despite much recent work, the true promise and limitations of the Quantum Alternating Operator Ansatz (QAOA) are unclear. A critical question regarding QAOA is to what extent its performance scales with the input size of the problem instance, in particular the necessary growth in the number of QAOA rounds to reach a high approximation ratio. We present numerical evidence for an exponential speed-up of QAOA over Grover-style unstructured search in finding approximate solutions to constrained optimization problems. Our result provides a strong hint that QAOA is able to exploit the structure of an optimization problem and thus overcome the lower bound for unstructured search. To this end, we conduct a comprehensive numerical study on several Hamming-weight constrained optimization problems for which we include combinations of all standardly studied mixer and phase separator Hamiltonians (Ring mixer, Clique mixer, Objective Value phase separator) as well as quantum minimum-finding inspired Hamiltonians (Grover mixer, Threshold-based phase separator). We identify Clique-Objective-QAOA with an exponential speed-up over Grover-Threshold-QAOA and tie the latter's scaling to that of unstructured search, with all other QAOA combinations coming in at a distant third. Our result suggests that maximizing QAOA performance requires a judicious choice of mixer and phase separator, and should trigger further research into other QAOA variations.


title: Timely Gossiping with File Slicing and Network Coding

abstract: We consider a system consisting of a large network of $n$ users and a library of files, wherein inter-user communication is established based upon gossip mechanisms. Each file is initially present at exactly one node, which is designated as the file \emph{source}. The source gets updated with newer versions of the file according to an arbitrary distribution in real time, and the other users in the network wish to acquire the latest possible version of the file. We present a class of gossip protocols that achieve $O(1)$ age at a typical node in a single-file system and $O(n)$ age at a typical node for a given file in an $n$-file system. We show that file slicing and network coding based protocols fall under the presented class of protocols. Numerical evaluation results are presented to confirm the aforementioned bounds.


title: AI-based Medical e-Diagnosis for Fast and Automatic Ventricular Volume Measurement in the Patients with Normal Pressure Hydrocephalus

abstract: Based on CT and MRI images acquired from normal pressure hydrocephalus (NPH) patients, using machine learning methods, we aim to establish a multi-modal and high-performance automatic ventricle segmentation method to achieve efficient and accurate automatic measurement of the ventricular volume. First, we extract the brain CT and MRI images of 143 definite NPH patients. Second, we manually label the ventricular volume (VV) and intracranial volume (ICV). Then, we use machine learning method to extract features and establish automatic ventricle segmentation model. Finally, we verify the reliability of the model and achieved automatic measurement of VV and ICV. In CT images, the Dice similarity coefficient (DSC), Intraclass Correlation Coefficient (ICC), Pearson correlation, and Bland-Altman analysis of the automatic and manual segmentation result of the VV were 0.95, 0.99, 0.99, and 4.2$\pm$2.6 respectively. The results of ICV were 0.96, 0.99, 0.99, and 6.0$\pm$3.8 respectively. The whole process takes 3.4$\pm$0.3 seconds. In MRI images, the DSC, ICC, Pearson correlation, and Bland-Altman analysis of the automatic and manual segmentation result of the VV were 0.94, 0.99, 0.99, and 2.0$\pm$0.6 respectively. The results of ICV were 0.93, 0.99, 0.99, and 7.9$\pm$3.8 respectively. The whole process took 1.9$\pm$0.1 seconds. We have established a multi-modal and high-performance automatic ventricle segmentation method to achieve efficient and accurate automatic measurement of the ventricular volume of NPH patients. This can help clinicians quickly and accurately understand the situation of NPH patient's ventricles.


title: Hydrodynamics and multiscale order in confluent epithelia

abstract: We formulate a hydrodynamic theory of confluent epithelia: i.e. monolayers of epithelial cells adhering to each other without gaps. Taking advantage of recent progresses toward establishing a general hydrodynamic theory of p-atic liquid crystals, we demonstrate that collectively migrating epithelia feature both nematic (i.e. p=2) and hexatic (i.e. p=6) order, with the former being dominant at large and the latter at small length scales. Such a remarkable multiscale liquid crystal order leaves a distinct signature in the system's structure factor, which exhibits two different power law scaling regimes, reflecting both the hexagonal geometry of small cells clusters, as well as the uniaxial structure of the global cellular flow. We support these analytical predictions with two different cell-resolved models of epithelia -- i.e. the self-propelled Voronoi model and the multiphase field model -- and highlight how momentum dissipation and noise influence the range of fluctuations at small length scales, thereby affecting the degree of cooperativity between cells. Our construction provides a theoretical framework to conceptualize the recent observation of multiscale order in layers of Madin-Darby canine kidney cells and pave the way for further theoretical developments.


title: Sexaquark dilemma in neutron stars and its solution by quark deconfinement

abstract: Following the idea that a stable sexaquark state with quark content (uuddss) would have gone unnoticed by experiment so far and that such a particle would be a good dark matter candidate, we investigate the possible role of a stable sexaquark in the physics of compact stars given the stringent constraints on the equation of state that stem from observations of high mass pulsars and GW170817 bounds on the compactness of intermediate mass stars. We find that there is a "sexaquark dilemma" (analogous to the hyperon dilemma) for which the dissociation of the sexaquark in quark matter is a viable solution fulfilling all present constraints from multi-messenger astronomy. The parameters needed to model the hybrid star including sexaquarks are in line with parameters of pre-existing quark- and hadronic-matter models. We find that current constraints -- tidal deformability in accordance with GW170817 and maximum mass above the lower limit from PSR J0740+6620 -- can be satisfied two ways: with early quark deconfinement such that neither sexaquarks nor hyperons are present in any NS interiors, or with later deconfinement such that a neutron-sexaquark shell surrounds the inner quark matter core.


title: The threshold model with anticonformity under random sequential updating

abstract: We study an asymmetric version of the threshold model with anticonformity under asynchronous update mode that mimics continuous time. We study this model on a complete graph using three different approaches: mean-field approximation, Monte Carlo simulation, and the Markov chain approach. The latter approach yields analytical results for arbitrarily small systems, in contrast to the mean-field approach, which is strictly correct only for an infinite system. We show that for sufficiently large systems, all three approaches produce the same results, as expected. We consider two cases: (1) homogeneous, in which all agents have the same tolerance threshold, and (2) heterogeneous, in which the thresholds are given by a beta distribution parametrized by two positive shape parameters $\alpha$ and $\beta$. The heterogeneous case can be treated as a generalized model that reduces to a homogeneous model in special cases. We show that particularly interesting behaviors, including social hysteresis and critical mass, arise only for values of $\alpha$ and $\beta$ that yield the shape of the distribution observed in real social systems.


title: Using the non-hydrodynamic mode to study the onset of hydrodynamic behavior in ultraperipheral symmetric nuclear collisions

abstract: With the attempts of extending the hydrodynamic framework of heavy-ion collision to proton-proton and other small and low energy systems, we are confronted with the question of how small the system can get and still be safely modelled as a fluid. One of the transport coefficients required in the $2^{nd}$ order relativistic viscous hydrodynamics is the shear relaxation time, inclusion of which solves the causality violation problem in the Navier-Stokes equation. In phenomenological studies this coefficient has been taken as a constant and much attention has gone into finding and fixing the shear viscosity to entropy density ratio, $\eta/s$. This transport coefficient also happens to control the non-hydrodynamic mode of the out-of-equilibrium hydrodynamics theory. It has been predicted that for decreasing system size, observables become sensitive to variation in shear relaxation time as a result of increasing dominance of non-hydrodynamic mode, which could potentially indicate breakdown of hydrodynamics. In this study, we try to test this prediction in the peripheral Pb-Pb collisions at $2.76$ TeV and Au-Au collisions at $200$ GeV, with IPGlasma initial condition and $(2+1)-$Dimensional viscous hydrodynamics. We find that elliptic flow does show adequate sensitivity to variation in relaxation time for decreasing system size. The multiplicity rapidity density limit for applicability of hydrodynamics is found to be around $dN/dy\approx10$, with the possibility of refinement in this value given a way to improve the centrality resolution in experimental data for referencing in peripheral collisions.


title: Topological Black Holes with curvature induced scalarization in the extended scalar-tensor theories

abstract: We study the perturbative behaviour of topological black holes in the presence of a cosmological constant and a scalar field coupled to the Gauss-Bonnet term. We calculate both analytically and numerically the quasi-normal modes of scalar perturbations in the extended scalar-tensor-Gauss-Bonnet gravity. In the case of small black holes we find a phase transition of the topological black hole to a hairy configuration.


title: Quasi-integrable modules, a class of non-highest weight modules over twisted affine Lie superalgebras

abstract: In this paper, we characterize quasi-integrable modules, of nonzero level, over twisted affine Lie superalgebras. We show that quasi-integrable modules are not necessarily highest weight modules. We prove that each quasi-integrable module is parabolically induced from a cuspidal module, over a finite dimensional Lie superalgebra having a Cartan subalgebra whose corresponding root system just contain real roots; in particular, the classification of quasi-integrable modules is reduced to the known classification of cuspidal modules over such Lie superalgebras.


title: Discovering Invariant Spatial Features in Electron Energy Loss Spectroscopy Images on the Mesoscopic and Atomic Levels

abstract: Over the last two decades, Electron Energy Loss Spectroscopy (EELS) imaging with a scanning transmission electron microscope (STEM) has emerged as a technique of choice for visualizing complex chemical, electronic, plasmonic, and phononic phenomena in complex materials and structures. The availability of the EELS data necessitates the development of methods to analyze multidimensional datasets with complex spatial and energy structures. Traditionally, the analysis of these data sets has been based on analysis of individual spectra, one at a time, whereas the spatial structure and correlations between individual spatial pixels containing the relevant information of the physics of underpinning processes have generally been ignored and analyzed only via the visualization as 2D maps. Here we develop a machine learning-based approach and workflows for the analysis of spatial structures in 3D EELS data sets using a combination of dimensionality reduction and multichannel rotationally-invariant variational autoencoders. This approach is illustrated for the analysis of both the plasmonic phenomena in a system of nanowires and in the core excitations in functional oxides using low loss and core loss EELS, respectively. The code developed in this manuscript is open sourced and freely available and provided as a Jupyter notebook for the interested reader here.


title: Scalable Fragment-Based 3D Molecular Design with Reinforcement Learning

abstract: Machine learning has the potential to automate molecular design and drastically accelerate the discovery of new functional compounds. Towards this goal, generative models and reinforcement learning (RL) using string and graph representations have been successfully used to search for novel molecules. However, these approaches are limited since their representations ignore the three-dimensional (3D) structure of molecules. In fact, geometry plays an important role in many applications in inverse molecular design, especially in drug discovery. Thus, it is important to build models that can generate molecular structures in 3D space based on property-oriented geometric constraints. To address this, one approach is to generate molecules as 3D point clouds by sequentially placing atoms at locations in space -- this allows the process to be guided by physical quantities such as energy or other properties. However, this approach is inefficient as placing individual atoms makes the exploration unnecessarily deep, limiting the complexity of molecules that can be generated. Moreover, when optimizing a molecule, organic and medicinal chemists use known fragments and functional groups, not single atoms. We introduce a novel RL framework for scalable 3D design that uses a hierarchical agent to build molecules by placing molecular substructures sequentially in 3D space, thus attempting to build on the existing human knowledge in the field of molecular design. In a variety of experiments with different substructures, we show that our agent, guided only by energy considerations, can efficiently learn to produce molecules with over 100 atoms from many distributions including drug-like molecules, organic LED molecules, and biomolecules.


title: Stay Positive: Non-Negative Image Synthesis for Augmented Reality

abstract: In applications such as optical see-through and projector augmented reality, producing images amounts to solving non-negative image generation, where one can only add light to an existing image. Most image generation methods, however, are ill-suited to this problem setting, as they make the assumption that one can assign arbitrary color to each pixel. In fact, naive application of existing methods fails even in simple domains such as MNIST digits, since one cannot create darker pixels by adding light. We know, however, that the human visual system can be fooled by optical illusions involving certain spatial configurations of brightness and contrast. Our key insight is that one can leverage this behavior to produce high quality images with negligible artifacts. For example, we can create the illusion of darker patches by brightening surrounding pixels. We propose a novel optimization procedure to produce images that satisfy both semantic and non-negativity constraints. Our approach can incorporate existing state-of-the-art methods, and exhibits strong performance in a variety of tasks including image-to-image translation and style transfer.


title: Interactron: Embodied Adaptive Object Detection

abstract: Over the years various methods have been proposed for the problem of object detection. Recently, we have witnessed great strides in this domain owing to the emergence of powerful deep neural networks. However, there are typically two main assumptions common among these approaches. First, the model is trained on a fixed training set and is evaluated on a pre-recorded test set. Second, the model is kept frozen after the training phase, so no further updates are performed after the training is finished. These two assumptions limit the applicability of these methods to real-world settings. In this paper, we propose Interactron, a method for adaptive object detection in an interactive setting, where the goal is to perform object detection in images observed by an embodied agent navigating in different environments. Our idea is to continue training during inference and adapt the model at test time without any explicit supervision via interacting with the environment. Our adaptive object detection model provides a 7.2 point improvement in AP (and 12.7 points in AP50) over DETR, a recent, high-performance object detector. Moreover, we show that our object detection model adapts to environments with completely different appearance characteristics, and performs well in them. The code is available at: https://github.com/allenai/interactron .


title: When Do Flat Minima Optimizers Work?

abstract: Recently, flat-minima optimizers, which seek to find parameters in low-loss neighborhoods, have been shown to improve a neural network's generalization performance over stochastic and adaptive gradient-based optimizers. Two methods have received significant attention due to their scalability: 1. Stochastic Weight Averaging (SWA), and 2. Sharpness-Aware Minimization (SAM). However, there has been limited investigation into their properties and no systematic benchmarking of them across different domains. We fill this gap here by comparing the loss surfaces of the models trained with each method and through broad benchmarking across computer vision, natural language processing, and graph representation learning tasks. We discover several surprising findings from these results, which we hope will help researchers further improve deep learning optimizers, and practitioners identify the right optimizer for their problem.


title: Systemic Risk Models for Disjoint and Overlapping Groups with Equilibrium Strategies

abstract: We analyze the systemic risk for disjoint and overlapping groups (e.g., central clearing counterparties (CCP)) by proposing new models with realistic game features. Specifically, we generalize the systemic risk measure proposed in [F. Biagini, J.-P. Fouque, M. Frittelli, and T. Meyer-Brandis, Finance and Stochastics, 24(2020), 513--564] by allowing individual banks to choose their preferred groups instead of being assigned to certain groups. We introduce the concept of Nash equilibrium for these new models, and analyze the optimal solution under Gaussian distribution of the risk factor. We also provide an explicit solution for the risk allocation of the individual banks, and study the existence and uniqueness of Nash equilibrium both theoretically and numerically. The developed numerical algorithm can simulate scenarios of equilibrium, and we apply it to study the bank-CCP structure with real data and show the validity of the proposed model.


title: Measuring the Accessibility of Domain Name Encryption and Its Impact on Internet Filtering

abstract: Most online communications rely on DNS to map domain names to their hosting IP address(es). Previous work has shown that DNS-based network interference is widespread due to the unencrypted and unauthenticated nature of the original DNS protocol. In addition to DNS, accessed domain names can also be monitored by on-path observers during the TLS handshake when the SNI extension is used. These lingering issues with exposed plaintext domain names have led to the development of a new generation of protocols that keep accessed domain names hidden. DNS-over-TLS (DoT) and DNS-over-HTTPS (DoH) hide the domain names of DNS queries, while Encrypted Server Name Indication (ESNI) encrypts the domain name in the SNI extension. We present DNEye, a measurement system built on top of a network of distributed vantage points, which we used to study the accessibility of DoT/DoH and ESNI, and to investigate whether these protocols are tampered with by network providers (e.g., for censorship). Moreover, we evaluate the efficacy of these protocols in circumventing network interference when accessing content blocked by traditional DNS manipulation. We find evidence of blocking efforts against domain name encryption technologies in several countries, including China, Russia, and Saudi Arabia. At the same time, we discover that domain name encryption can help with unblocking more than 55% and 95% of censored domains in China and other countries where DNS-based filtering is heavily employed.


title: Controller confidentiality for nonlinear systems under sensor attacks

abstract: Controller confidentiality under sensor attacks refers to whether the internal states of the controller can be estimated when the adversary knows the model of the plant and controller, while only having access to sensors, but not the actuators. We show that the controller's state can be estimated accurately when the nonlinear closed-loop system is detectable. In the absence of detectability, controller confidentiality can still be breached with a periodic probing scheme via the sensors under a robust observability assumption, which allows for the controller's state to be estimated with arbitrary accuracy during the probing period, and with bounded error during the non-probing period. Further, stealth can be maintained by choosing an appropriate probing duration. This study shows that the controller confidentiality for nonlinear systems can be breached by balancing the estimation precision and the stealthiness of the adversary.


title: Tutorial on amortized optimization

abstract: Optimization is a ubiquitous modeling tool and is often deployed in settings which repeatedly solve similar instances of the same problem. Amortized optimization methods use learning to predict the solutions to problems in these settings, exploiting the shared structure between similar problem instances. These methods have been crucial in variational inference and reinforcement learning and are capable of solving optimization problems many orders of magnitudes times faster than traditional optimization methods that do not use amortization. This tutorial presents an introduction to the amortized optimization foundations behind these advancements and overviews their applications in variational inference, sparse coding, gradient-based meta-learning, control, reinforcement learning, convex optimization, optimal transport, and deep equilibrium networks. The source code for this tutorial is available at https://github.com/facebookresearch/amortized-optimization-tutorial.


title: Locally Typical Sampling

abstract: Today's probabilistic language generators fall short when it comes to producing coherent and fluent text despite the fact that the underlying models perform well under standard metrics, e.g., perplexity. This discrepancy has puzzled the language generation community for the last few years. In this work, we posit that the abstraction of natural language generation as a discrete stochastic process--which allows for an information-theoretic analysis--can provide new insights into the behavior of probabilistic language generators, e.g., why high-probability texts can be dull or repetitive. Humans use language as a means of communicating information, aiming to do so in a simultaneously efficient and error-minimizing manner; in fact, psycholinguistics research suggests humans choose each word in a string with this subconscious goal in mind. We formally define the set of strings that meet this criterion: those for which each word has an information content close to the expected information content, i.e., the conditional entropy of our model. We then propose a simple and efficient procedure for enforcing this criterion when generating from probabilistic models, which we call locally typical sampling. Automatic and human evaluations show that, in comparison to nucleus and top-k sampling, locally typical sampling offers competitive performance (in both abstractive summarization and story generation) in terms of quality while consistently reducing degenerate repetitions.


title: DKM: Dense Kernelized Feature Matching for Geometry Estimation

abstract: Feature matching is a challenging computer vision task that involves finding correspondences between two images of a 3D scene. In this paper we consider the dense approach instead of the more common sparse paradigm, thus striving to find all correspondences. Perhaps counter-intuitively, dense methods have previously shown inferior performance to their sparse and semi-sparse counterparts for estimation of two-view geometry. This changes with our novel dense method, which outperforms both dense and sparse methods on geometry estimation. The novelty is threefold: First, we propose a kernel regression global matcher. Secondly, we propose warp refinement through stacked feature maps and depthwise convolution kernels. Thirdly, we propose learning dense confidence through consistent depth and a balanced sampling approach for dense confidence maps. Through extensive experiments we confirm that our proposed dense method, \textbf{D}ense \textbf{K}ernelized Feature \textbf{M}atching, sets a new state-of-the-art on multiple geometry estimation benchmarks. In particular, we achieve an improvement on MegaDepth-1500 of +4.9 and +8.9 AUC$@5^{\circ}$ compared to the best previous sparse method and dense method respectively. Our code is provided at https://github.com/Parskatt/dkm


title: Epithelia are multiscale active liquid crystals

abstract: Biological processes such as embryogenesis, wound healing and cancer progression, crucially rely on the ability of epithelial cells to coordinate their mechanical activity over length scales order of magnitudes larger than the typical cellular size. While regulated by various signalling pathways, it has recently become evident that this behavior can additionally hinge on a minimal toolkit of physical mechanisms, of which liquid crystal order is the most prominent example. Yet, experimental and theoretical studies have given so far inconsistent results in this respect: whereas nematic order is often invoked in the interpretation of experimental data, computational models have instead suggested that hexatic order could serve as a linchpin for collective migration in confluent cell layers. In this article we resolve this dilemma. Using a combination of in vitro experiments on Madin-Darby canine kidney cells (MDCK), numerical simulations and analytical work, we demonstrate that both nematic and hexatic order are, in fact, present in epithelial layers, with the former being dominant at large length scales and the latter at small length scales. In MDCK GII cells on uncoated glass, these different types of liquid crystal order crossover at a length scale of the order of ten cell sizes. Our work sheds light on the emergent organization of living matter, provides a new framework for deciphering the structure of epithelia, and paves the way toward a comprehensive and predictive mesoscopic theory of tissues.


title: Superfluidity in the 1D Bose-Hubbard Model

abstract: We study superfluidity in the 1D Bose-Hubbard model using a variational matrix product state technique. We determine the superfluid density as a function of the Hubbard parameters by calculating the energy cost of phase twists in the thermodynamic limit. As the system is critical, correlation functions decay as power laws and the entanglement entropy grows with the bond dimension of our variational state. We relate the resulting scaling laws to the superfluid density. We compare two different algorithms for optimizing the infinite matrix product state and develop a physical explanation why one of them (VUMPS) is more efficient than the other (iDMRG). Finally, we comment on finite-temperature superfluidity in one dimension and how our results can be realized in cold atom experiments.


title: Higher-order initial conditions with massive neutrinos

abstract: The discovery that neutrinos have mass has important consequences for cosmology. The main effect of massive neutrinos is to suppress the growth of cosmic structure on small scales. Such growth can be accurately modelled using cosmological $N$-body simulations, but doing so requires accurate initial conditions (ICs). There is a trade-off, especially with first-order ICs, between truncation errors for late starts and discreteness and relativistic errors for early starts. Errors can be minimized by starting simulations at late times using higher-order ICs. In this paper, we show that neutrino effects can be absorbed into scale-independent coefficients in higher-order Lagrangian perturbation theory (LPT). This clears the way for the use of higher-order ICs for massive neutrino simulations. We demonstrate that going to higher order substantially improves the accuracy of simulations. To match the sensitivity of surveys like DESI and Euclid, errors in the matter power spectrum should be well below 1%. However, we find that first-order Zel'dovich ICs lead to much larger errors, even when starting as early as $z=127$, exceeding 1% at $z=0$ for $k>0.5\text{ Mpc}^{-1}$ for the power spectrum and $k>0.1\text{ Mpc}^{-1}$ for the equilateral bispectrum in our simulations. Ratios of power spectra with different neutrino masses are more robust than absolute statistics, but still depend on the choice of ICs. For all statistics considered, we obtain 1% agreement between 2LPT and 3LPT at $z=0$.


title: Equivalence principle and HBAR entropy of an atom falling into a quantum corrected black hole

abstract: In this work, we have investigated the phenomenon of acceleration radiation exhibited by an atom falling into a quantum corrected Schwarzschild black hole. We observe that the excitation-probability of the atom with simultaneous emission of a photon satisfies the equivalence principle when we compare it to the excitation probability of a mirror accelerating with respect to an atom. We also demonstrate the validity of the equivalence principle for a generic black hole geometry. Then we calculate the horizon brightened acceleration radiation (HBAR) entropy for this quantum corrected black hole geometry. We observed that the HBAR entropy has the form identical to that of Bekenstein-Hawking black hole entropy along with universal quantum gravity corrections.


title: Relation Between the Partial Derivatives of the Kinetic Energy in the Lagrangian and Hamiltonian Formalisms of Dynamics

abstract: The partial derivative of the kinetic energy of a dynamical system with respect to a generalized coordinate as it appears in the Lagrangian formalism is not equal to the derivative of the kinetic energy with respect to the same coordinate in the Hamiltonian formalism but differs by a sign. We find another exact relation between the two partial derivatives in the case of a conservative system. We also identify another form of kinetic energy whose partial derivative with respect to a generalized coordinate vanishes identically.


title: Visualizing Automatic Speech Recognition -- Means for a Better Understanding?

abstract: Automatic speech recognition (ASR) is improving ever more at mimicking human speech processing. The functioning of ASR, however, remains to a large extent obfuscated by the complex structure of the deep neural networks (DNNs) they are based on. In this paper, we show how so-called attribution methods, that we import from image recognition and suitably adapt to handle audio data, can help to clarify the working of ASR. Taking DeepSpeech, an end-to-end model for ASR, as a case study, we show how these techniques help to visualize which features of the input are the most influential in determining the output. We focus on three visualization techniques: Layer-wise Relevance Propagation (LRP), Saliency Maps, and Shapley Additive Explanations (SHAP). We compare these methods and discuss potential further applications, such as in the detection of adversarial examples.