title: The Music Streaming Sessions Dataset

abstract: At the core of many important machine learning problems faced by online streaming services is a need to model how users interact with the content they are served. Unfortunately, there are no public datasets currently available that enable researchers to explore this topic. In order to spur that research, we release the Music Streaming Sessions Dataset (MSSD), which consists of 160 million listening sessions and associated user actions. Furthermore, we provide audio features and metadata for the approximately 3.7 million unique tracks referred to in the logs. This is the largest collection of such track metadata currently available to the public. This dataset enables research on important problems including how to model user listening and interaction behaviour in streaming, as well as Music Information Retrieval (MIR), and session-based sequential recommendations. Additionally, a subset of sessions were collected using a uniformly random recommendation setting, enabling their use for counterfactual evaluation of such sequential recommendations. Finally, we provide an analysis of user behavior and suggest further research problems which can be addressed using the dataset.


title: The Mott-semiconducting state in the magic angle bilayer graphene

abstract: Using non-perturbative theoretical method, we address the problem of strong correlations in twisted bilayer-layer graphene at the magic angle. We concentrate on the solution without symmetry breaking, where conventional Mott insulating state is expected for all integer fillings. At Coulomb repulsion corresponding to dielectric constant $\varepsilon\approx 5$ and several integer fillings we find a Mott-semiconducting state, which simultaneously hosts the Mott state, and inside the Mott gap, a second much smaller semiconduting gap. The presence of these Mott-ingap states, which are located at the $\Gamma$ point, makes the Mott state strongly temperature dependent and leads to a bad-metal phase at elevated temperatures. The system is insulating at the charge neutrality point and at even fillings away from it.


title: Universal four-dimensional representation of $H \to \gamma \gamma$ at two loops through the Loop-Tree Duality

abstract: We extend useful properties of the $H\to\gamma\gamma$ unintegrated dual amplitudes from one- to two-loop level, using the Loop-Tree Duality formalism. In particular, we show that the universality of the functional form -- regardless of the nature of the internal particle -- still holds at this order. We also present an algorithmic way to renormalise two-loop amplitudes, by locally cancelling the ultraviolet singularities at integrand level, thus allowing a full four-dimensional numerical implementation of the method. Our results are compared with analytic expressions already available in the literature, finding a perfect numerical agreement. The success of this computation plays a crucial role for the development of a fully local four-dimensional framework to compute physical observables at Next-to-Next-to Leading order and beyond.


title: Multi-modal dialog for browsing large visual catalogs using exploration-exploitation paradigm in a joint embedding space

abstract: We present a multi-modal dialog system to assist online shoppers in visually browsing through large catalogs. Visual browsing is different from visual search in that it allows the user to explore the wide range of products in a catalog, beyond the exact search matches. We focus on a slightly asymmetric version of the complete multi-modal dialog where the system can understand both text and image queries but responds only in images. We formulate our problem of "showing $k$ best images to a user" based on the dialog context so far, as sampling from a Gaussian Mixture Model in a high dimensional joint multi-modal embedding space, that embed both the text and the image queries. Our system remembers the context of the dialog and uses an exploration-exploitation paradigm to assist in visual browsing. We train and evaluate the system on a multi-modal dialog dataset that we generate from large catalog data. Our experiments are promising and show that the agent is capable of learning and can display relevant results with an average cosine similarity of 0.85 to the ground truth. Our preliminary human evaluation also corroborates the fact that such a multi-modal dialog system for visual browsing is well-received and is capable of engaging human users.


title: Generalized Minors and Tensor Invariants

abstract: Berenstein, Fomin and Zelevinsky defined functions on double Bruhat cells which they called generalized minors. By relating certain double Bruhat cells to configuration spaces of flags, we give formulas for these generalized minors as tensor invariants. This allows us to verify certain weight identities. The weights of the tensor invariants can then be used to construct the quiver for the cluster structure on the configuration space of three flags. We also show a converse statement--that the weights of tensor invariants can by computed from the structure of the quiver. The weight identities are important because they are necessary for the existence of cluster structures on the moduli space of framed local systems.


title: Optimal verification of general bipartite pure states

abstract: The efficient and reliable verification of quantum states plays a crucial role in various quantum information processing tasks. We consider the task of verifying entangled states using one-way and two-way classical communication and completely characterize the optimal strategies via convex optimization. We solve these optimization problems using both analytical and numerical methods, and the optimal strategies can be constructed for any bipartite pure state. Compared with the nonadaptive approach, our adaptive strategies significantly improve the efficiency of quantum state verification. Moreover, these strategies are experimentally feasible, as only few local projective measurements are required.


title: Unambiguous Phase Spaces for Subregions

abstract: The covariant phase space technique is a powerful formalism for understanding the Hamiltonian description of covariant field theories. However, applications of this technique to problems involving subregions, such as the exterior of a black hole, have heretofore been plagued by boundary ambiguities. We provide a resolution of these ambiguities by directly computing the symplectic structure from the path integral, showing that it may be written as a contour integral around a partial Cauchy surface. This result has implications for gauge symmetry and entanglement.


title: Utility Preserving Secure Private Data Release

abstract: Differential privacy mechanisms that also make reconstruction of the data impossible come at a cost - a decrease in utility. In this paper, we tackle this problem by designing a private data release mechanism that makes reconstruction of the original data impossible and also preserves utility for a wide range of machine learning algorithms. We do so by combining the Johnson-Lindenstrauss (JL) transform with noise generated from a Laplace distribution. While the JL transform can itself provide privacy guarantees \cite{blocki2012johnson} and make reconstruction impossible, we do not rely on its differential privacy properties and only utilize its ability to make reconstruction impossible. We present novel proofs to show that our mechanism is differentially private under single element changes as well as single row changes to any database. In order to show utility, we prove that our mechanism maintains pairwise distances between points in expectation and also show that its variance is proportional to the dimensionality of the subspace we project the data into. Finally, we experimentally show the utility of our mechanism by deploying it on the task of clustering.


title: Graphs with a unique maximum open packing

abstract: A set $S$ of vertices in a graph is an open packing if (open) neighborhoods of any two distinct vertices in $S$ are disjoint. In this paper, we consider the graphs that have a unique maximum open packing. We characterize the trees with this property by using four local operations such that any nontrivial tree with a unique maximum open packing can be obtained by a sequence of these operations starting from $P_2$. We also prove that the decision version of the open packing number is NP-complete even when restricted to graphs of girth at least $6$. Finally, we show that the recognition of the graphs with a unique maximum open packing is polynomially equivalent to the recognition of the graphs with a unique maximum independent set, and we prove that the complexity of both problems is not polynomial, unless P=NP.


title: Measurement of ultra-low heating rates of a single antiproton in a cryogenic Penning trap

abstract: We report on the first detailed study of motional heating in a cryogenic Penning trap using a single antiproton. Employing the continuous Stern-Gerlach effect we observe cyclotron quantum transition rates of 6(1) quanta/h and an electric field noise spectral density below $7.5(3.4)\times 10^{-20}\,\text{V}^{2}\text{m}^{-2} \text{Hz}^{-1}$, which corresponds to a scaled noise spectral density below $8.8(4.0)\times 10^{-12}\,\text{V}^{2}\text{m}^{-2}$, results which are more than two orders of magnitude smaller than those reported by other ion trap experiments.


title: Tumor Connectomics: Mapping the intra-tumoral complex interaction network

abstract: Tumors are extremely heterogeneous and comprise of a number of intratumor microenvironments or sub-regions. These tumor microenvironments may interact with eac based on complex high-level relationships, which could provide important insight into the organizational structure of the tumor network. To that end, we developed a tumor connectomics framework (TCF) to understand and model the complex functional and morphological interactions within the tumor. Then, we demonstrate the TCF's potential in predicting treatment response in breast cancer patients being treated with neoadjuvant chemotherapy. The TCF was implemented on a breast cancer patient cohort of thirty-four patients with dynamic contrast enhanced (DCE) magnetic resonance imaging (MRI) undergoing neodjuvant chemotherapy treatment. The intra-tumor network connections (tumor connectome) before and after treatment were modeled using advanced graph theoretic centrality, path length and clustering metrics from the DCE-MRI. The percentage change of the graph metrics between two time-points (Baseline and 1st cycle) was computed to predict the patient's final response to treatment. The TCF visualized the inter-voxel network connections across multiple time-points and was able to evaluate specific changes in the tumor connectome with treatment. Degree centrality was identified as the most significant predictor of treatment response with an AUC of 0.83 for classifying responders from non-responders. In conclusion, the TCF graph metrics produced excellent biomarkers for prediction of breast cancer treatment response with improved visualization and interpretability of changes both locally and globally in the tumor.


title: Line shape and $D^{(\ast)}\bar D^{(\ast)}$ probabilities of $\psi(3770)$ from the $e^+e^-\to D\bar D$ reaction

abstract: We have performed a calculation of the $D\bar D$, $D\bar D^\ast$, $D^\ast\bar D$, $D^\ast\bar D^\ast$ components in the wave function of the $\psi(3770)$. For this we make use of the $^3P_0$ model to find the coupling of $\psi(3770)$ to these components, that with an elaborate angular momentum algebra can be obtained with only one parameter. Then we use data for the $e^+e^-\to D\bar D$ reaction, from where we determine a form factor needed in the theoretical frame work, as well as other parameters needed to evaluate the meson-meson selfenergy of the $\psi(3770)$. Once this is done we determine the $Z$ probability to still have a vector core and the probability to have the different meson components. We find $Z$ about $80\sim85\%$, and the individual meson-meson components are rather small, providing new empirical information to support the largely $q\bar q$ component of vector mesons, and the $\psi(3770)$ in particular.


title: Efficient Multiparty Interactive Coding for Insertions, Deletions and Substitutions

abstract: In the field of interactive coding, two or more parties wish to carry out a distributed computation over a communication network that may be noisy. The ultimate goal is to develop efficient coding schemes that can tolerate a high level of noise while increasing the communication by only a constant factor (i.e., constant rate). In this work we consider synchronous communication networks over an arbitrary topology, in the powerful adversarial insertion-deletion noise model. Namely, the noisy channel may adversarially alter the content of any transmitted symbol, as well as completely remove a transmitted symbol or inject a new symbol into the channel. We provide efficient, constant rate schemes that successfully conduct any computation with high probability as long as the adversary corrupts at most $\varepsilon /m$ fraction of the total communication, where $m$ is the number of links in the network and $\varepsilon$ is a small constant. This scheme assumes the parties share a random string to which the adversarial noise is oblivious. We can remove this assumption at the price of being resilient to $\varepsilon / (m\log m)$ adversarial error. While previous work considered the insertion-deletion noise model in the two-party setting, to the best of our knowledge, our scheme is the first multiparty scheme that is resilient to insertions and deletions. Furthermore, our scheme is the first computationally efficient scheme in the multiparty setting that is resilient to adversarial noise.


title: Computing electrostatic potentials using regularization based on the range-separated tensor format

abstract: In this paper, we apply the range-separated (RS) tensor format [6] for the construction of new regularization scheme for the Poisson-Boltzmann equation (PBE) describing the electrostatic potential in biomolecules. In our approach, we use the RS tensor representation to the discretized Dirac delta [21] to construct an efficient RS splitting of the PBE solution in the solute (molecular) region. The PBE then needs to be solved with a regularized source term, and thus black-box solvers can be applied. The main computational benefits are due to the localization of the modified right-hand side within the molecular region and automatic maintaining of the continuity in the Cauchy data on the interface. Moreover, this computational scheme only includes solving a single system of FDM/FEM equations for the smooth long-range (i.e., regularized) part of the collective potential represented by a low-rank RS-tensor with a controllable precision. The total potential is obtained by adding this solution to the directly precomputed rank-structured tensor representation for the short-range contribution. Enabling finer grids in PBE computations is another advantage of the proposed techniques. In the numerical experiments, we consider only the free space electrostatic potential for proof of concept. We illustrate that the classical Poisson equation (PE) model does not accurately capture the solution singularities in the numerical approximation as compared to the new approach by the RS tensor format.


title: Asynchronous Accelerated Proximal Stochastic Gradient for Strongly Convex Distributed Finite Sums

abstract: In this work, we study the problem of minimizing the sum of strongly convex functions split over a network of $n$ nodes. We propose the decentralized and asynchronous algorithm ADFS to tackle the case when local functions are themselves finite sums with $m$ components. ADFS converges linearly when local functions are smooth, and matches the rates of the best known finite sum algorithms when executed on a single machine. On several machines, ADFS enjoys a $O (\sqrt{n})$ or $O(n)$ speed-up depending on the leading complexity term as long as the diameter of the network is not too big with respect to $m$. This also leads to a $\sqrt{m}$ speed-up over state-of-the-art distributed batch methods, which is the expected speed-up for finite sum algorithms. In terms of communication times and network parameters, ADFS scales as well as optimal distributed batch algorithms. As a side contribution, we give a generalized version of the accelerated proximal coordinate gradient algorithm using arbitrary sampling that we apply to a well-chosen dual problem to derive ADFS. Yet, ADFS uses primal proximal updates that only require solving one-dimensional problems for many standard machine learning applications. Finally, ADFS can be formulated for non-smooth objectives with equally good scaling properties. We illustrate the improvement of ADFS over state-of-the-art approaches with simulations.


title: Connecting Cycles for Concentric Circles

abstract: We study perimeters of connecting cycles for concentric circles. More precisely, we are interested in characterization of those connecting cycles which are critical points of perimeter considered as a function on the product of given circles. Specifically, we aim at showing that, generically, perimeter is a Morse function on the configuration space, and computing Morse indices of critical configurations. In particular, we prove that the diametrically aligned configurations are critical and their indices can be calculated from an explicitly given tridiagonal matrix. For four concentric circles, we give examples of non-generic collections of radii and describe a pitchfork type bifurcation of stationary connecting cycles.


title: It could be worse, it could be raining: reliable automatic meteorological forecasting

abstract: Meteorological forecasting provides reliable prediction about the future weather within a given interval of time. Meteorological forecasting can be viewed as a form of hybrid diagnostic reasoning and can be mapped onto an integrated conceptual framework. The automation of the forecasting process would be helpful in a number of contexts, in particular: when the amount of data is too wide to be dealt with manually; to support forecasters education; when forecasting about underpopulated geographic areas is not interesting for everyday life (and then is out from human forecasters' tasks) but is central for tourism sponsorship. We present logic MeteoLOG, a framework that models the main steps of the reasoner the forecaster adopts to provide a bulletin. MeteoLOG rests on several traditions, mainly on fuzzy, temporal and probabilistic logics. On this basis, we also introduce the algorithm Tournament, that transforms a set of MeteoLOG rules into a defeasible theory, that can be implemented into an automatic reasoner. We finally propose an example that models a real world forecasting scenario.


title: Boundary integral formula for harmonic functions on Riemann surfaces

abstract: We construct a boundary integral formula for harmonic functions on open, smoothly-bordered subdomains of Riemann surfaces embeddable into $\C\P^2$. The formula may be considered as an analogue of the Green's formula for domains in $\C$.


title: Space and complexities of territorial systems

abstract: The spatial character of territorial systems plays a crucial role in the emergence of their complexities. This contribution aims at illustrating to what extent different types of complexities can be exhibited in models of such systems. We develop from a theoretical viewpoint some arguments illustrating ontological complexity, in the sense of the diversity and multidimensionality of possible representations, and then complexity in the sense of emergence, i.e. the necessity of the existence of several autonomous levels. We then propose numerical experiments to explore properties of complexity (dynamical complexity and co-evolution) within two simple models of urban morphogenesis. We finally suggest other dimensions of complexity which could be typical of territorial systems.


title: Self-Assembly of Geometric Space from Random Graphs

abstract: We present a Euclidean quantum gravity model in which random graphs dynamically self-assemble into discrete manifold structures. Concretely, we consider a statistical model driven by a discretisation of the Euclidean Einstein-Hilbert action; contrary to previous approaches based on simplicial complexes and Regge calculus our discretisation is based on the Ollivier curvature, a coarse analogue of the manifold Ricci curvature defined for generic graphs. The Ollivier curvature is generally difficult to evaluate due to its definition in terms of optimal transport theory, but we present a new exact expression for the Ollivier curvature in a wide class of relevant graphs purely in terms of the numbers of short cycles at an edge. This result should be of independent intrinsic interest to network theorists. Action minimising configurations prove to be cubic complexes up to defects; there are indications that such defects are dynamically suppressed in the macroscopic limit. Closer examination of a defect free model shows that certain classical configurations have a geometric interpretation and discretely approximate vacuum solutions to the Euclidean Einstein-Hilbert action. Working in a configuration space where the geometric configurations are stable vacua of the theory, we obtain direct numerical evidence for the existence of a continuous phase transition; this makes the model a UV completion of Euclidean Einstein gravity. Notably, this phase transition implies an area-law for the entropy of emerging geometric space. Certain vacua of the theory can be interpreted as baby universes; we find that these configurations appear as stable vacua in a mean field approximation of our model, but are excluded dynamically whenever the action is exact indicating the dynamical stability of geometric space. The model is intended as a setting for subsequent studies of emergent time mechanisms.


title: The Brown-Erd\H{o}s-S\'os Conjecture in finite abelian groups

abstract: The Brown-Erd\H{o}s-S\'{o}s conjecture, one of the central conjectures in extremal combinatorics, states that for any integer $m\geq 6,$ if a 3-uniform hypergraph on $n$ vertices contains no $m$ vertices spanning at least $m-3$ edges, then the number of edges is $o(n^2).$ We prove the conjecture for triple systems coming from finite abelian groups.


title: Dynamical Generation of Elementary Fermion Mass: First Lattice Evidence

abstract: Using lattice simulations we demonstrate from first principles the existence of a non-perturbative mechanism for elementary particle mass generation in models with gauge fields, fermions and scalars, if an exact invariance forbids power divergent fermion masses and fermionic chiral symmetries broken at UV scale are maximally restored. We show that in the Nambu-Goldstone phase a fermion mass term, unrelated to the Yukawa operator, is dynamically generated. In models with electro-weak interactions weak boson masses are also generated opening new scenarios for beyond the Standard Model physics.


title: Physical Access Control Management System Based on Permissioned Blockchain

abstract: Using blockchain as a decentralized backend infrastructure has grabbed the attention of many startups entrepreneurs and developers. Blockchain records transactions permanently and protects them from undesirable tampering. It provides a reliable tamper-proof database which can be considered as a trustable source for tracking the previous system state. In this paper, we present our access control application based on Hyperledger Fabric Blockchain and Hyperledger Composer to control access to physical places. The system components and modular architecture are illustrated, and we have extracted metadata include historian transactions details arising from our demo test. Finally, the performance metrics and resources consumption are provided using Hyperledger Caliper, a benchmark framework for measuring Hyperledger blockchains performance.


title: Confronting gravitational-wave observations with modern nuclear physics constraints

abstract: Multi-messenger observations of neutron star (NS) mergers have the potential to revolutionize nuclear astrophysics. They will improve our understanding of nucleosynthesis, provide insights about the equation of state (EOS) of strongly-interacting matter at high densities, and enable tests of the theory of gravity and of dark matter. Here, we focus on the EOS, where both gravitational waves (GWs) from neutron-star mergers and X-ray observations from space-based detectors such as NICER will provide more stringent constraints on the structure of neutron stars. Furthermore, recent advances in nuclear theory have enabled reliable calculations of the EOS at low densities using effective field theory based Hamiltonians and advanced techniques to solve the quantum many-body problem. In this paper, we address how the first observation of GWs from GW170817 can be combined with modern calculations of the EOS to extract useful insights about the EOS of matter encountered inside neutron stars. We analyze the impact of various uncertainties, the role of phase transitions in the NS core, and discuss how future observations will improve our understanding of dense matter.


title: Eruptions from Quiet Sun Coronal Bright Points. II. Non-Potential Modeling

abstract: Our recent observational study shows that the majority of coronal bright points (CBPs) in the quiet Sun are sources of one or more eruptions during their lifetime. Here, we investigate the non-potential time dependent structure of the magnetic field of the CBP regions with special emphasis on the time evolving magnetic structure at the spatial locations where the eruptions are initiated. The magnetic structure is evolved in time using a Non-Linear Force Free Field (NLFFF) relaxation approach, based on a time series of Helioseismic and Magnetic Imager (HMI) longitudinal magnetograms. This results in a continuous time series of NLFFFs. The time series is initiated with a potential field extrapolation based on a magnetogram taken well before the time of the eruptions. This initial field is then evolved in time in response to the observed changes in the magnetic field distribution at the photosphere. The local and global magnetic field structures from the time series of NLFFF field solutions are analysed in the vicinity of the eruption sites at the approximate times of the eruptions. The analysis shows that many of the CBP eruptions reported in Mou etl. (2018) contain twisted flux tube located at the sites of eruptions. The presence of flux ropes at these locations provides in many cases a direct link between the magnetic field structure, their eruption and the observation of mini coronal mass ejections (mini-CMEs). It is found that all repetitive eruptions are homologous. The NLFFF simulations show that twisted magnetic field structures are created at the locations hosting eruptions in CBPs. These twisted structures are produced by footpoint motions imposed by changes in the photospheric magnetic field observations. The true nature of the micro-flares remains unknown. Further 3D data-driven MHD modelling is required to show how these twisted regions become unstable and erupt.


title: Nonnegative $C^2(\mathbb{R}^2)$ interpolation

abstract: In this paper, we prove two improved versions of the Finiteness Principle for nonnegative $ C^2(\mathbb{R}^2) $ interpolation, previously proven by Fefferman, Israel, and Luli. The first version sharpens the finiteness constant to $ 64 $, and the second version carries better computational practicality. Along the way, we also provide detailed construction of nonnegative $ C^2 $ interpolants in one-dimension, and prove the nonexistence of a bounded linear $ C^2 $-extension operator that preserves nonnegativity.


title: Dominating Sets and Connected Dominating Sets in Dynamic Graphs

abstract: In this paper we study the dynamic versions of two basic graph problems: Minimum Dominating Set and its variant Minimum Connected Dominating Set. For those two problems, we present algorithms that maintain a solution under edge insertions and edge deletions in time $O(\Delta\cdot \text{polylog}~n)$ per update, where $\Delta$ is the maximum vertex degree in the graph. In both cases, we achieve an approximation ratio of $O(\log n)$, which is optimal up to a constant factor (under the assumption that $P \ne NP$). Although those two problems have been widely studied in the static and in the distributed settings, to the best of our knowledge we are the first to present efficient algorithms in the dynamic setting. As a further application of our approach, we also present an algorithm that maintains a Minimal Dominating Set in $O(min(\Delta, \sqrt{m}))$ per update.


title: CapsAttacks: Robust and Imperceptible Adversarial Attacks on Capsule Networks

abstract: Capsule Networks preserve the hierarchical spatial relationships between objects, and thereby bears a potential to surpass the performance of traditional Convolutional Neural Networks (CNNs) in performing tasks like image classification. A large body of work has explored adversarial examples for CNNs, but their effectiveness on Capsule Networks has not yet been well studied. In our work, we perform an analysis to study the vulnerabilities in Capsule Networks to adversarial attacks. These perturbations, added to the test inputs, are small and imperceptible to humans, but can fool the network to mispredict. We propose a greedy algorithm to automatically generate targeted imperceptible adversarial examples in a black-box attack scenario. We show that this kind of attacks, when applied to the German Traffic Sign Recognition Benchmark (GTSRB), mislead Capsule Networks. Moreover, we apply the same kind of adversarial attacks to a 5-layer CNN and a 9-layer CNN, and analyze the outcome, compared to the Capsule Networks to study differences in their behavior.


title: Chinese syzygies by insertions

abstract: We construct a finite convergent semi-quadratic presentation for the Chinese monoid by adding column generators and using combinatorial properties of insertion algorithms on Chinese staircases. We extend this presentation into a coherent one whose generators are columns, rewriting rules are defined by insertion algorithms, and whose syzygies are defined as relations among insertion algorithms. Such a coherent presentation is used for representations of Chinese monoids, in particular, it is a way to describe actions of Chinese monoids on categories.


title: Quantum simulation of spontaneous pair creation in 2D optical lattices

abstract: One of the fundamental predictions of Quantum Electrodynamics (QED) is the spontaneous creation of particle--antiparticle pairs from vacuum in presence of a very strong electric field. Under these extreme conditions a strongly bound state can fetch an otherwise unobservable electron from the Dirac sea, leaving behind a hole representing a positron. Although generally known for many decades, the effect has not yet been demonstrated experimentally. We propose an analogue model of the quantum Dirac field, realized by ultra--cold fermionic atoms in an optical lattice, aiming at an experimental simulation of this intriguing non--perturbative phenomenon. Numerical simulations demonstrate the effect of spontaneous pair creation in the optical analogue system, in qualitative agreement with QED: in the adiabatic regime the vacuum can be destabilized only by supercritical fields exceeding a critical threshold.


title: Scalable Metropolis-Hastings for Exact Bayesian Inference with Large Datasets

abstract: Bayesian inference via standard Markov Chain Monte Carlo (MCMC) methods is too computationally intensive to handle large datasets, since the cost per step usually scales like $\Theta(n)$ in the number of data points $n$. We propose the Scalable Metropolis-Hastings (SMH) kernel that exploits Gaussian concentration of the posterior to require processing on average only $O(1)$ or even $O(1/\sqrt{n})$ data points per step. This scheme is based on a combination of factorized acceptance probabilities, procedures for fast simulation of Bernoulli processes, and control variate ideas. Contrary to many MCMC subsampling schemes such as fixed step-size Stochastic Gradient Langevin Dynamics, our approach is exact insofar as the invariant distribution is the true posterior and not an approximation to it. We characterise the performance of our algorithm theoretically, and give realistic and verifiable conditions under which it is geometrically ergodic. This theory is borne out by empirical results that demonstrate overall performance benefits over standard Metropolis-Hastings and various subsampling algorithms.


title: Laser Cooling of Radium Ions

abstract: The unstable radium nucleus is appealing for probing new physics due to its high mass, octupole deformation and energy level structure. Ion traps, with long hold times and low particle numbers, are excellent for work with radioactive species, such as radium and radium-based molecular ions, where low activity, and hence low total numbers, is desirable. We address the challenges associated with the lack of stable isotopes in a tabletop experiment with a low-activity ($\sim 10 \ \mu\mathrm{Ci}$) source where we laser-cool trapped radium ions. With a laser-cooled radium ion we measured the $7p\ ^2P_{1/2}^o$ state's branching fractions to the ground state, $7s\ ^2S_{1/2}$, and a metastable excited state, $6d\ ^2D_{3/2}$, to be $p=0.9104(7)$ and $0.0896(7)$, respectively. With a nearby tellurium reference line we measured the $7s\ ^2S_{1/2} \rightarrow 7p\ ^2P_{1/2}^o$ transition frequency, 640.09663(6) THz.


title: Propensity Score Matching underestimates Real Treatment Effect, in a simulated theoretical multivariate model

abstract: Propensity Score Matching (PSM) is an useful method to reduce the impact ofTreatment - Selection Bias in the estimation of causal effects in observational studies. After matching, the PSM significantly reduces the sample under investigation, which may lead to other possible biases. In this sense, we want to analyse the behaviour of this PSM compared with other widely used method to deal with non-comparable groups, as is the Multivariate Regression Model (MRM). Monte Carlo Simulations are made to construct groups with different effects in order to compare the behaviour of PSM and MRM estimating this effects. Also the Treatment - Selection Bias reduction for the PSM is calculated. With the PSM a reduction in the Treatment - Selection Bias is achieved, with a reduction in the Relative Real Treatment Effect Estimation Error, but despite of this bias reduction and estimation error reduction, the MRM significantly reduces more this estimation error compared with the PSM. Also the PSM leads to a not insignificant reduction of the sample. This loss of information derived from the matching process may lead to another not known bias, and thus, to the inaccurate of the effect estimation compared with the MRM.


title: Sub-threshold binary neutron star search in Advanced LIGO's first observing run

abstract: We present a search for gravitational waves from double neutron star binaries inspirals in Advanced LIGO's first observing run. The search considers a narrow range of binary chirp masses motivated by the population of known double neutron star binaries in the nearby universe. This search differs from previously published results by providing the most sensitive published survey of neutron stars in Advanced LIGO's first observing run within this narrow mass range and including times when only one of the two LIGO detectors was in operation in the analysis. The search was sensitive to binary neutron star inspirals to an average distance of ~85 Mpc over 93.2 days. We do not identify any unambiguous gravitational wave signals in our sample of 103 sub-threshold candidates with false-alarm-rates of less than one per day. However, given the expected binary neutron star merger rate of R = 100 - 4000 Gpc^(-3) yr^(-1), we expect O(1) gravitational wave events within our candidate list. This suggests the possibility that one or more of these candidates is in fact a binary neutron star merger. Although the contamination fraction in our candidate list is ~99%, it might be possible to correlate these events with other messengers to identify a potential multi-messenger signal. We provide an online candidate list with the times and sky locations for all events in order to enable multi-messenger searches.


title: Towards an Extremal Network Theory -- Robust GDoF Gain of Transmitter Cooperation over TIN

abstract: Significant progress has been made recently in Generalized Degrees of Freedom (GDoF) characterizations of wireless interference channels (IC) and broadcast channels (BC) under the assumption of finite precision channel state information at the transmitters (CSIT), especially for smaller or highly symmetric network settings. A critical barrier in extending these results to larger and asymmetric networks is the inherent combinatorial complexity of such networks. Motivated by other fields such as extremal combinatorics and extremal graph theory, we explore the possibility of an extremal network theory, i.e., a study of extremal networks within particular regimes of interest. As our test application, we study the GDoF benefits of transmitter cooperation in a $K$ user IC over the simple scheme of power control and treating interference as Gaussian noise (TIN) for three regimes of interest -- a TIN regime where TIN was shown to be GDoF optimal for the $K$ user interference channel, a CTIN regime where the GDoF region achievable by TIN is convex without time-sharing, and an SLS regime where a simple layered superposition (SLS) scheme is shown to be optimal in the $K$ user MISO BC, albeit only for $K\leq 3$. As our first result, we show that under finite precision CSIT, TIN is GDoF optimal for the $K$ user IC throughout the CTIN regime. Furthermore, under finite precision CSIT, appealing to extremal network theory we obtain the following results. In both TIN and CTIN regimes, we show that the extremal GDoF gain from transmitter cooperation over TIN is bounded regardless of the number of users: the gain is exactly a factor of $3/2$ in the TIN regime, and $2-1/K$ in the CTIN regime, for arbitrary number of users $K>1$. However, in the SLS regime, the gain is $\Theta(\log_2(K))$, i.e., it scales logarithmically with the number of users.


title: CoCoNet: A Collaborative Convolutional Network

abstract: We present an end-to-end deep network for fine-grained visual categorization called Collaborative Convolutional Network (CoCoNet). The network uses a collaborative layer after the convolutional layers to represent an image as an optimal weighted collaboration of features learned from training samples as a whole rather than one at a time. This gives CoCoNet more power to encode the fine-grained nature of the data with limited samples. We perform a detailed study of the performance with 1-stage and 2-stage transfer learning. The ablation study shows that the proposed method outperforms its constituent parts consistently. CoCoNet also outperforms few state-of-the-art competing methods. Experiments have been performed on the fine-grained bird species classification problem as a representative example, but the method may be applied to other similar tasks. We also introduce a new public dataset for fine-grained species recognition, that of Indian endemic birds and have reported initial results on it.


title: On the Units of GANs (Extended Abstract)

abstract: Generative Adversarial Networks (GANs) have achieved impressive results for many real-world applications. As an active research topic, many GAN variants have emerged with improvements in sample quality and training stability. However, visualization and understanding of GANs is largely missing. How does a GAN represent our visual world internally? What causes the artifacts in GAN results? How do architectural choices affect GAN learning? Answering such questions could enable us to develop new insights and better models. In this work, we present an analytic framework to visualize and understand GANs at the unit-, object-, and scene-level. We first identify a group of interpretable units that are closely related to concepts with a segmentation-based network dissection method. We quantify the causal effect of interpretable units by measuring the ability of interventions to control objects in the output. Finally, we examine the contextual relationship between these units and their surrounding by inserting the discovered concepts into new images. We show several practical applications enabled by our framework, from comparing internal representations across different layers, models, and datasets, to improving GANs by locating and removing artifact-causing units, to interactively manipulating objects in the scene. We will open source our interactive tools to help researchers and practitioners better understand their models.


title: Federated Collaborative Filtering for Privacy-Preserving Personalized Recommendation System

abstract: The increasing interest in user privacy is leading to new privacy preserving machine learning paradigms. In the Federated Learning paradigm, a master machine learning model is distributed to user clients, the clients use their locally stored data and model for both inference and calculating model updates. The model updates are sent back and aggregated on the server to update the master model then redistributed to the clients. In this paradigm, the user data never leaves the client, greatly enhancing the user' privacy, in contrast to the traditional paradigm of collecting, storing and processing user data on a backend server beyond the user's control. In this paper we introduce, as far as we are aware, the first federated implementation of a Collaborative Filter. The federated updates to the model are based on a stochastic gradient approach. As a classical case study in machine learning, we explore a personalized recommendation system based on users' implicit feedback and demonstrate the method's applicability to both the MovieLens and an in-house dataset. Empirical validation confirms a collaborative filter can be federated without a loss of accuracy compared to a standard implementation, hence enhancing the user's privacy in a widely used recommender application while maintaining recommender performance.


title: Numerical and Exact Analyses of Bures and Hilbert-Schmidt Separability and PPT-Probabilities

abstract: We employ a quasirandom methodology, recently developed by Martin Roberts, to estimate the separability probabilities, with respect to the Bures (minimal monotone/statistical distinguishability) measure, of generic two-qubit and two-rebit states. This procedure, based on generalized properties of the golden ratio, yielded, in the course of almost seventeen billion iterations (recorded at intervals of five million), two-qubit estimates repeatedly close to nine decimal places to $\frac{25}{341} =\frac{5^2}{11 \cdot 31} \approx 0.073313783$. Howeer, despite the use of over twenty-three billion iterations, we do not presently perceive an exact value (rational or otherwise) for an estimate of 0.15709623 for the Bures two-rebit separability probability. The Bures qubit-qutrit case--for which Khvedelidze and Rogojin gave an estimate of 0.0014--is analyzed too. The value of $\frac{1}{715}=\frac{1}{5 \cdot 11 \cdot 13} \approx 0.00139860$ is a well-fitting value to an estimate of 0.00139884. Interesting values ($\frac{16}{12375} =\frac{4^2}{3^2 \cdot 5^3 \cdot 11}$ and $\frac{625}{109531136}=\frac{5^4}{2^{12} \cdot 11^2 \cdot 13 \cdot 17}$) are conjectured for the Hilbert-Schmidt (HS) and Bures qubit-qudit ($2 \times 4$) positive-partial-transpose (PPT)-probabilities. We re-examine, strongly supporting, conjectures that the HS qubit-{\it qutrit} and rebit-{\it retrit} separability probabilities are $\frac{27}{1000}=\frac{3^3}{2^3 \cdot 5^3}$ and $\frac{860}{6561}= \frac{2^2 \cdot 5 \cdot 43}{3^8}$, respectively. Prior studies have demonstrated that the HS two-rebit separability probability is $\frac{29}{64}$ and strongly pointed to the HS two-qubit counterpart being $\frac{8}{33}$, and a certain operator monotone one (other than the Bures) being $1 -\frac{256}{27 \pi^2}$.

categories: quant-ph math-ph math.MP math.PR physics.comp-ph


title: Few-shot Learning with Meta Metric Learners

abstract: Few-shot Learning aims to learn classifiers for new classes with only a few training examples per class. Existing meta-learning or metric-learning based few-shot learning approaches are limited in handling diverse domains with various number of labels. The meta-learning approaches train a meta learner to predict weights of homogeneous-structured task-specific networks, requiring a uniform number of classes across tasks. The metric-learning approaches learn one task-invariant metric for all the tasks, and they fail if the tasks diverge. We propose to deal with these limitations with meta metric learning. Our meta metric learning approach consists of task-specific learners, that exploit metric learning to handle flexible labels, and a meta learner, that discovers good parameters and gradient decent to specify the metrics in task-specific learners. Thus the proposed model is able to handle unbalanced classes as well as to generate task-specific metrics. We test our approach in the `$k$-shot $N$-way' few-shot learning setting used in previous work and new realistic few-shot setting with diverse multi-domain tasks and flexible label numbers. Experiments show that our approach attains superior performances in both settings.


title: See Better Before Looking Closer: Weakly Supervised Data Augmentation Network for Fine-Grained Visual Classification

abstract: Data augmentation is usually adopted to increase the amount of training data, prevent overfitting and improve the performance of deep models. However, in practice, random data augmentation, such as random image cropping, is low-efficiency and might introduce many uncontrolled background noises. In this paper, we propose Weakly Supervised Data Augmentation Network (WS-DAN) to explore the potential of data augmentation. Specifically, for each training image, we first generate attention maps to represent the object's discriminative parts by weakly supervised learning. Next, we augment the image guided by these attention maps, including attention cropping and attention dropping. The proposed WS-DAN improves the classification accuracy in two folds. In the first stage, images can be seen better since more discriminative parts' features will be extracted. In the second stage, attention regions provide accurate location of object, which ensures our model to look at the object closer and further improve the performance. Comprehensive experiments in common fine-grained visual classification datasets show that our WS-DAN surpasses the state-of-the-art methods, which demonstrates its effectiveness.


title: A Black-box Attack on Neural Networks Based on Swarm Evolutionary Algorithm

abstract: Neural networks play an increasingly important role in the field of machine learning and are included in many applications in society. Unfortunately, neural networks suffer from adversarial samples generated to attack them. However, most of the generation approaches either assume that the attacker has full knowledge of the neural network model or are limited by the type of attacked model. In this paper, we propose a new approach that generates a black-box attack to neural networks based on the swarm evolutionary algorithm. Benefiting from the improvements in the technology and theoretical characteristics of evolutionary algorithms, our approach has the advantages of effectiveness, black-box attack, generality, and randomness. Our experimental results show that both the MNIST images and the CIFAR-10 images can be perturbed to successful generate a black-box attack with 100\% probability on average. In addition, the proposed attack, which is successful on distilled neural networks with almost 100\% probability, is resistant to defensive distillation. The experimental results also indicate that the robustness of the artificial intelligence algorithm is related to the complexity of the model and the data set. In addition, we find that the adversarial samples to some extent reproduce the characteristics of the sample data learned by the neural network model.


title: A simple electronic device to experiment with the Hopf bifurcation

abstract: We present a simple low-cost electronic circuit that is able to show two different dynamical regimens with oscillations of voltages and with constant values of them. This device is designed as a negative feedback three-node network inspired in the genetic repressilator. The circuit's behavior is modeled by a system of differential equations which is studied in several different ways by applying the dynamical system formalism, making numerical simulations and constructing and measuring it experimentally. We find that the most important characteristics of the Hopf bifurcation can be found and controlled. Particularly, a resistor value plays the role of the bifurcation parameter, which can be easily varied experimentally. As a result, this system can be employed to introduce many aspects of a research in a real physical system and it enables us to study one of the most important kinds of bifurcation.


title: Composing Distributed Data-intensive Web Services Using a Flexible Memetic Algorithm

abstract: Web Service Composition (WSC) is a particularly promising application of Web services, where multiple individual services with specific functionalities are composed to accomplish a more complex task, which must fulfil functional requirements and optimise Quality of Service (QoS) attributes, simultaneously. Additionally, large quantities of data, produced by technological advances, need to be exchanged between services. Data-intensive Web services, which manipulate and deal with those data, are of great interest to implement data-intensive processes, such as distributed Data-intensive Web Service Composition (DWSC). Researchers have proposed Evolutionary Computing (EC) fully-automated WSC techniques that meet all the above factors. Some of these works employed Memetic Algorithms (MAs) to enhance the performance of EC through increasing its exploitation ability of in searching neighbourhood area of a solution. However, those works are not efficient or effective. This paper proposes an MA-based approach to solving the problem of distributed DWSC in an effective and efficient manner. In particular, we develop an MA that hybridises EC with a flexible local search technique incorporating distance of services. An evaluation using benchmark datasets is carried out, comparing existing state-of-the-art methods. Results show that our proposed method has the highest quality and an acceptable execution time overall.


title: Modularization of End-to-End Learning: Case Study in Arcade Games

abstract: Complex environments and tasks pose a difficult problem for holistic end-to-end learning approaches. Decomposition of an environment into interacting controllable and non-controllable objects allows supervised learning for non-controllable objects and universal value function approximator learning for controllable objects. Such decomposition should lead to a shorter learning time and better generalisation capability. Here, we consider arcade-game environments as sets of interacting objects (controllable, non-controllable) and propose a set of functional modules that are specialized on mastering different types of interactions in a broad range of environments. The modules utilize regression, supervised learning, and reinforcement learning algorithms. Results of this case study in different Atari games suggest that human-level performance can be achieved by a learning agent within a human amount of game experience (10-15 minutes game time) when a proper decomposition of an environment or a task is provided. However, automatization of such decomposition remains a challenging problem. This case study shows how a model of a causal structure underlying an environment or a task can benefit learning time and generalization capability of the agent, and argues in favor of exploiting modular structure in contrast to using pure end-to-end learning approaches.


title: Quadratic torsion subgroups of modular Jacobian varieties

abstract: Let $D$ be an odd square-free positive integer and $C$ a divisor of $D$. For any quadratic character $\chi$ modulo $C$, we prove that the $\chi$-part of the group $J_0(DC)_\text{tor}$ of torsion points of $J_0(DC)$ coincides with the $\chi$-part of its cuspidal subgroup, away from those primes of bad reduction or where possible congruences between oldforms and newforms occur.


title: Sobolev inequalities for the symmetric gradient in arbitrary domains

abstract: A form of Sobolev inequalities for the symmetric gradient of vector-valued functions is proposed, which allows for arbitrary ground domains in $\mathbb R ^n$. In the relevant inequalities, boundary regularity of domains is replaced with information on boundary traces of trial functions. The inequalities so obtained exhibit the same exponents as in classical inequalities for the full gradient of Sobolev functions, in regular domains. Furthermore, they involve constants independent of the geometry of the domain, and hence yield novel results yet for smooth domains. Our approach relies upon a pointwise estimate for the functions in question via a Riesz potential of their symmetric gradient and an unconventional potential depending on their boundary trace.


title: About an extension of the Davenport-Rado result to the Herzog-Schonheim conjecture for free groups

abstract: Let $G$ be a group and $H_1$,...,$H_s$ be subgroups of $G$ of indices $d_1$,...,$d_s$ respectively. In 1974, M. Herzog and J. Sch\"onheim conjectured that if $\{H_i\alpha_i\}_{i=1}^{i=s}$, $\alpha_i\in G$, is a coset partition of $G$, then $d_1$,..,$d_s$ cannot be distinct. We consider the Herzog-Sch\"onheim conjecture for free groups of finite rank and propose a new approach, based on an extension of the Davenport-Rado result for $G=\mathbb{Z}$.


title: Non-geodesic incompleteness in Poincar\'e gauge gravity

abstract: In this work we review the study of singularities in Poincar\'e gauge theories of gravity. Since one of the most recent studies uses the appearance of black hole regions of arbitrary dimension as an indicator of singular behaviour, we also give some explicit examples of these structures and study how particles behave around them.


title: Pinning of Diffusional Patterns by Non-Uniform Curvature

abstract: Diffusion-driven patterns appear on curved surfaces in many settings, initiated by unstable modes of an underlying Laplacian operator. On a flat surface or perfect sphere, the patterns are degenerate, reflecting translational/rotational symmetry. Deformations, e.g. by a bulge or indentation, break symmetry and can pin a pattern. We adapt methods of conformal mapping and perturbation theory to examine how curvature inhomogeneities select and pin patterns, and confirm the results numerically. The theory provides an analogy to quantum mechanics in a geometry-dependent potential and yields intuitive implications for cell membranes, tissues, thin films, and noise-induced quasipatterns.


title: Asymptotic Performance Analysis of Generalized User Selection for Interference-Limited Multiuser Secondary Networks

abstract: We analyze the asymptotic performance of a generalized multiuser diversity scheme for an interference-limited secondary multiuser network of underlay cognitive radio systems. Assuming a large number of secondary users and that the noise at each secondary user's receiver is negligible compared to the interference from the primary transmitter, the secondary transmitter transmits information to the $k$-th best secondary user, namely, the one with the $k$-th highest signal-to-interference ratio (SIR). We use extreme value theory to show that the $k$-th highest SIR converges uniformly in distribution to an inverse gamma random variable for a fixed $k$ and large number of secondary users. We use this result to derive asymptotic expressions for the average throughput, effective throughput, average bit error rate and outage probability of the $k$-th best secondary user under continuous power adaptation at the secondary transmitter, which ensures satisfaction of the instantaneous interference constraint at the primary receiver caused by the secondary transmitter. Numerical simulations show that our derived asymptotic expressions are accurate for different values of system parameters.


title: The CM Algorithm for the Maximum Mutual Information Classifications of Unseen Instances

abstract: The Maximum Mutual Information (MMI) criterion is different from the Least Error Rate (LER) criterion. It can reduce failing to report small probability events. This paper introduces the Channels Matching (CM) algorithm for the MMI classifications of unseen instances. It also introduces some semantic information methods, which base the CM algorithm. In the CM algorithm, label learning is to let the semantic channel match the Shannon channel (Matching I) whereas classifying is to let the Shannon channel match the semantic channel (Matching II). We can achieve the MMI classifications by repeating Matching I and II. For low-dimensional feature spaces, we only use parameters to construct n likelihood functions for n different classes (rather than to construct partitioning boundaries as gradient descent) and expresses the boundaries by numerical values. Without searching in parameter spaces, the computation of the CM algorithm for low-dimensional feature spaces is very simple and fast. Using a two-dimensional example, we test the speed and reliability of the CM algorithm by different initial partitions. For most initial partitions, two iterations can make the mutual information surpass 99% of the convergent MMI. The analysis indicates that for high-dimensional feature spaces, we may combine the CM algorithm with neural networks to improve the MMI classifications for faster and more reliable convergence.


title: Ranks for families of all theories of given languages

abstract: For families of all theories of arbitrary given languages we describe ranks and degrees. In particular, we characterize (non-)totally transcendental families. We apply these characterizations for the families of all theories of given languages, with models of given finite or infinite cardinality.


title: FPSA: A Full System Stack Solution for Reconfigurable ReRAM-based NN Accelerator Architecture

abstract: Neural Network (NN) accelerators with emerging ReRAM (resistive random access memory) technologies have been investigated as one of the promising solutions to address the \textit{memory wall} challenge, due to the unique capability of \textit{processing-in-memory} within ReRAM-crossbar-based processing elements (PEs). However, the high efficiency and high density advantages of ReRAM have not been fully utilized due to the huge communication demands among PEs and the overhead of peripheral circuits. In this paper, we propose a full system stack solution, composed of a reconfigurable architecture design, Field Programmable Synapse Array (FPSA) and its software system including neural synthesizer, temporal-to-spatial mapper, and placement & routing. We highly leverage the software system to make the hardware design compact and efficient. To satisfy the high-performance communication demand, we optimize it with a reconfigurable routing architecture and the placement & routing tool. To improve the computational density, we greatly simplify the PE circuit with the spiking schema and then adopt neural synthesizer to enable the high density computation-resources to support different kinds of NN operations. In addition, we provide spiking memory blocks (SMBs) and configurable logic blocks (CLBs) in hardware and leverage the temporal-to-spatial mapper to utilize them to balance the storage and computation requirements of NN. Owing to the end-to-end software system, we can efficiently deploy existing deep neural networks to FPSA. Evaluations show that, compared to one of state-of-the-art ReRAM-based NN accelerators, PRIME, the computational density of FPSA improves by 31x; for representative NNs, its inference performance can achieve up to 1000x speedup.


title: Dirac star in the presence of Maxwell and Proca fields

abstract: We consider configurations consisting of a gravitating nonlinear spinor field $\psi$, with a nonlinearity of the type $\lambda\left(\bar\psi\psi\right)^2$, minimally coupled to Maxwell and Proca fields through the coupling constants $Q_M$ [U(1) electric charge] and $Q_P$, respectively. In order to ensure spherical symmetry of the configurations, we use two spin-$1/2$ fields having opposite spins. By means of numerical computations, we find families of equilibrium configurations with a positive Arnowitt-Deser-Misner (ADM) mass described by regular zero-node asymptotically flat solutions for static Maxwell and Proca fields and for stationary spinor fields. For the case of the Maxwell field, it is shown that, with increasing charge $Q_M$, the masses of the objects increase and diverge as the charge tends to a critical value. For negative values of the coupling constant $\lambda$, we demonstrate that, by choosing physically reasonable values of this constant, it is possible to obtain configurations with masses comparable to the Chandrasekhar mass and with effective radii of the order of kilometers. It enables us to speak of an astrophysical interpretation of such systems, regarding them as charged Dirac stars. In turn, for the system with the Proca field, it is shown that the mass of the configurations also grows with increasing both $|\lambda|$ and the coupling constant $Q_P$. Although in this case the numerical calculations do not allow us to make a definite conclusion about the possibility of obtaining masses comparable to the Chandrasekhar mass for physically reasonable values of $\lambda$, one may expect that such masses can be obtained for certain values of free parameters of the system under consideration.


title: Hierarchically Clustered Representation Learning

abstract: The joint optimization of representation learning and clustering in the embedding space has experienced a breakthrough in recent years. In spite of the advance, clustering with representation learning has been limited to flat-level categories, which often involves cohesive clustering with a focus on instance relations. To overcome the limitations of flat clustering, we introduce hierarchically-clustered representation learning (HCRL), which simultaneously optimizes representation learning and hierarchical clustering in the embedding space. Compared with a few prior works, HCRL firstly attempts to consider a generation of deep embeddings from every component of the hierarchy, not just leaf components. In addition to obtaining hierarchically clustered embeddings, we can reconstruct data by the various abstraction levels, infer the intrinsic hierarchical structure, and learn the level-proportion features. We conducted evaluations with image and text domains, and our quantitative analyses showed competent likelihoods and the best accuracies compared with the baselines.


title: Symmetrized p-convexity and Related Some Integral Inequalities

abstract: In this paper, the author introduces the concept of the symmetrized p-convex function, gives Hermite-Hadamard type inequalities for symmetrized p-convex functions.


title: Heat engine model exhibit super-universal feature and capture the efficiencies of different power plants

abstract: We propose a generalized model of a heat engine and calculate the minimum and maximum bounds on the efficiency at maximum power. We obtain a universal form of generalized extreme bounds on the efficiency at maximum power. Our model unifies the bounds on the efficiency and the universality features are observed for various heat engine models. Even though our model is a direct generalization of low-dissipation heat engines, the bounds on the efficiency obtained at a single target function capture those observed in the actual power plants working at different dissipation levels.


title: Effect of Bi Substitution on Thermoelectric Properties of SbSe2-based Layered Compounds NdO$_{0.8}$F$_{0.2}$Sb$_{1-x}$Bi$_x$Se$_2$

abstract: Although SbSe2-based layered compounds have been predicted to be high-performance thermoelectric materials and topological materials, most of these compounds obtained experimentally have been insulators so far. Here, we present the effect of Bi substitution on the thermoelectric properties of SbSe2-based layered compounds NdO0.8F0.2Sb1-xBixSe2 (x = 0-0.4). The room temperature electrical resistivity is decreased to 8.0 * 10^-5 ohmm for x = 0.4. The electrical power factor is calculated to be 1.4 * 10^-4 W/mK^2 at 660 K, which is in reasonable agreement with combined Jonker and Ioffe analysis. The room-temperature lattice thermal conductivity of less than 1 W/mK is almost independent of x, in contrast to the point-defect scattering model for conventional alloys. The present work provides an avenue for exploring SbSe2-based insulating and BiSe2-based conducting systems.


title: Charmonium Excitation functions in $\bar p$A Collisions

abstract: We study the excitation function of the low-lying charmonium state: $\Psi$(3686) in $\bar p$ Au collisions taking into account their in-medium propagation. The time evolution of the spectral functions of the charmonium state is studied with a BUU type transport model. We calculated the excitation function of $\Psi$(3686) production and show that it is strongly effected by the medium. The energy regime will be available for the PANDA experiment.


title: A conditional Berry-Esseen inequality

abstract: As an extension of a central limit theorem established by Svante Janson, we prove a Berry-Esseen inequality for a sum of independent and identically distributed random variables conditioned by a sum of independent and identically distributed integer-valued random variables.


title: Further Spectral Properties of the Weighted Finite Transform Operator and Approximation in Weighted Sobolev Spaces

abstract: In this work, we first give some mathematical preliminairies concerning the generelized prolate spheroidal wave functions(GPSWFs). This set of special functions have been introduced in [21]and [13] and they are defined as the infinite and countable set of the eigenfunctions of a weighted finite Fourier transform operator. Then, we show that the set of the singular values of this operator has a super-exponential decay rate. We also give some local estimates and bounds of these GPSWFs. As an application of the spectral properties of the GPSWFs and their associated eigenvalues, we give their quality of approximation in a weighted Sobolev space.Finally, we provide the reader with some numerical examples that illustrate the different results of this work.


title: Topological structures of velocity and electric field in the vicinity of a cusp-type magnetic null point

abstract: Topological characteristics reveal important physical properties of plasma structures and astrophysical processes. Physical parameters and constraints are linked with topological invariants, which are important for describing magnetic reconnection scenarios. We analyze stationary non-ideal Ohm's law concerning the Poincare classes of all involved physical fields in 2D by calculating the corresponding topological invariants of their Jacobian (here: particularly the eigenvalues) or Hessian matrices. The magnetic field is assumed to have a cusp structure, and the stagnation point of the plasma flow coincides with the cusp. We find that the stagnation point must be hyperbolic. Furthermore, the functions describing both the resistivity and the Ohmic heating have a saddle point structure, being displaced with respect to the cusp point. These results imply that there is no monotonous relation between current density and anomalous resistivity in the case of a 2D standard magnetic cusp.

categories: physics.plasm-ph astro-ph.SR physics.flu-dyn physics.space-ph


title: Analytical study of the parametric instability of an oscillating scalar field in an expanding universe

abstract: We investigate the dynamics of the perturbations of the inflaton scalar field oscillating around a minimum of its effective potential in an expanding universe. With the assumption of smallness of the ratio of the Hubble parameter to the oscillation frequency we apply the technique of separation of fast and slow motions. Considering the oscillation phase and the energy density as fast and slow variables we derive the Hill equation for the fluctuation modes in which the energy density is treated as a slowly varying parameter. We develop a general perturbative approach to solving the equations of this type, which is based on the Floquet theory and asymptotic expansions in the vicinity of the solutions with the "frozen" parameters. As an example, we consider the $\phi ^{2}-\phi ^{4}$ potential and construct the approximate solutions of the corresponding Lam\'{e} equation. The obtained solutions are found to be in a good agreement with the results of the direct numerical integration.


title: Vortex cusps

abstract: We consider pairs of self-similar 2d vortex sheets forming cusps, equivalently single sheets merging into slip condition walls, as in classical Mach reflection at wedges. We derive from the Birkhoff-Rott equation a reduced model yielding formulas for cusp exponents and other quantities as functions of similarity exponent and strain coefficient. Comparison to numerics shows that piecewise quadratic and higher approximation of vortex sheets agree with each other and with the model. In contrast piecewise linear schemes produce spurious results and violate conservation of mass, a problem that may have been undetected in prior work for other vortical flows where even point vortices were sufficient. We find that vortex cusps only exist if the similarity exponent is sufficiently large and if the circulation on the sheet is counterclockwise (for a sheet above the wall with cusp opening to the right), unless a sufficiently positive strain coefficient compensates. Whenever a cusp cannot exist a spiral-ends jet forms instead; we find many jets are so narrow that they appear as false cusps.


title: Clustered Millimeter Wave Networks with Non-Orthogonal Multiple Access

abstract: We introduce clustered millimeter wave networks with invoking non-orthogonal multiple access~(NOMA) techniques, where the NOMA users are modeled as Poisson cluster processes and each cluster contains a base station (BS) located at the center. To provide realistic directional beamforming, an actual antenna array pattern is deployed at all BSs. We propose three distance-dependent user selection strategies to appraise the path loss impact on the performance of our considered networks. With the aid of such strategies, we derive tractable analytical expressions for the coverage probability and system throughput. Specifically, closed-form expressions are deduced under a sparse network assumption to improve the calculation efficiency. It theoretically demonstrates that the large antenna scale benefits the near user, while such influence for the far user is fluctuant due to the randomness of the beamforming. Moreover, the numerical results illustrate that: 1) the proposed system outperforms traditional orthogonal multiple access techniques and the commonly considered NOMA-mmWave scenarios with the random beamforming; 2) the coverage probability has a negative correlation with the variance of intra-cluster receivers; 3) 73 GHz is the best carrier frequency for near user and 28 GHz is the best choice for far user; 4) an optimal number of the antenna elements exists for maximizing the system throughput.


title: Testing Conditional Independence in Supervised Learning Algorithms

abstract: We propose the conditional predictive impact (CPI), a consistent and unbiased estimator of the association between one or several features and a given outcome, conditional on a reduced feature set. Building on the knockoff framework of Cand\`es et al. (2018), we develop a novel testing procedure that works in conjunction with any valid knockoff sampler, supervised learning algorithm, and loss function. The CPI can be efficiently computed for high-dimensional data without any sparsity constraints. We demonstrate convergence criteria for the CPI and develop statistical inference procedures for evaluating its magnitude, significance, and precision. These tests aid in feature and model selection, extending traditional frequentist and Bayesian techniques to general supervised learning tasks. The CPI may also be applied in causal discovery to identify underlying multivariate graph structures. We test our method using various algorithms, including linear regression, neural networks, random forests, and support vector machines. Empirical results show that the CPI compares favorably to alternative variable importance measures and other nonparametric tests of conditional independence on a diverse array of real and simulated datasets. Simulations confirm that our inference procedures successfully control Type I error and achieve nominal coverage probability. Our method has been implemented in an R package, cpi, which can be downloaded from https://github.com/dswatson/cpi.


title: A distant OB association around RAFGL 5475

abstract: Observations of the galactic disk at mid-infrared and longer wavelengths reveal a wealth of structures indicating the existence of complexes of recent massive star formation. However, little or nothing is known about the stellar component of those complexes. We have carried out observations aiming at the identification of early-type stars in the direction of the bright infrared source RAFGL~5475, around which several interstellar medium structures usually associated with the presence of massive stars have been identified. Our observations have the potential of revealing the suspected but thus far unknown stellar component of the region around RAFGL~5475. We have carried out near-infrared imaging observations ($JHK_S$ bands) designed to reveal the presence of early-type stars based on their positions in color-color and color-magnitude diagrams centered on the location of RAFGL~5475. We took into account the possibility that candidates found might belong to a foreground population physically related either to M16 or M17, two giant HII regions lying midway between the Sun and RAFGL~5475. The near-infrared color-color diagram shows clear evidence for the presence of a moderately obscured population of early-type stars in the region imaged. By studying the distribution of extinction in their direction and basic characteristics of the interstellar medium we show that these new early-type stars are most likely associated with RAFGL~5475. By investigating the possible existence of massive early-type stars in the direction of RAFGL~5475 we have discovered the existence of a new OB association. A very preliminary assessment of its contents suggests the presence of several O-type stars, some of them likely to be associated with structures in the interstellar medium. The new association is located at 4 kpc from the Sun in the Scutum-Centaurus arm.


title: Inferring Heterogeneous Causal Effects in Presence of Spatial Confounding

abstract: We address the problem of inferring the causal effect of an exposure on an outcome across space, using observational data. The data is possibly subject to unmeasured confounding variables which, in a standard approach, must be adjusted for by estimating a nuisance function. Here we develop a method that eliminates the nuisance function, while mitigating the resulting errors-in-variables. The result is a robust and accurate inference method for spatially varying heterogeneous causal effects. The properties of the method are demonstrated on synthetic as well as real data from Germany and the US.


title: Boundary Continuity of Nonparametric Prescribed Mean Curvature Surfaces

abstract: We investigate the boundary behavior of variational solutions of Dirichlet problems for prescribed mean curvature equations at smooth boundary points where certain boundary curvature conditions are satisfied (which preclude the existence of local barrier functions). We prove that if the Dirichlet boundary data $\phi$ is continuous at such a point (and possibly nowhere else), then the solution of the variational problem is continuous at this point.


title: Generic magnetic field intensity profiles of interplanetary coronal mass ejections at Mercury, Venus and Earth from superposed epoch analyses

abstract: We study interplanetary coronal mass ejections (ICMEs) measured by probes at different heliocentric distances (0.3-1AU) to investigate the propagation of ICMEs in the inner heliosphere and determine how the generic features of ICMEs change with heliospheric distance. Using data from the MESSENGER, Venus Express and ACE spacecraft, we analyze with the superposed epoch technique the profiles of ICME substructures, namely the sheath and the magnetic ejecta. We determine that the median magnetic field magnitude in the sheath correlates well with ICME speeds at 1 AU and we use this proxy to order the ICMEs at all spacecraft. We then investigate the typical ICME profiles for three categories equivalent to slow, intermediate and fast ICMEs. Contrary to fast ICMEs, slow ICMEs have a weaker solar wind field at the front and a more symmetric magnetic field profile. We find the asymmetry to be less pronounced at Earth than at Mercury, indicating a relaxation taking place as ICMEs propagate. We also find that the magnetic field intensities in the wake region of the ICMEs do not go back to the pre-ICME solar wind intensities, suggesting that the effects of ICMEs on the ambient solar wind last longer than the duration of the transient event. Such results provide an indication of physical processes that need to be reproduced by numerical simulations of ICME propagation. The samples studied here will be greatly improved by future missions dedicated to the exploration of the inner heliosphere, such as Parker Solar Probe and Solar Orbiter.


title: Definable retractions over Henselian valued fields with analytic structure

abstract: Let $K$ be a Henselian, non-trivially valued field with separated analytic structure. We prove the existence of definable retractions onto an arbitrary closed definable subset of $K^{n}$. Hence directly follow definable non-Archimedean versions of the extension theorems by Tietze--Urysohn and Dugundji. This generalizes our previous paper dealing with complete non-Archimedean fields with separated power series and remains true for Henselian valued fields with strictly convergent analytic structure, because every such a structure can be extended in a definitional way to a separated analytic structure. Our proof uses a variant of the one from that paper, based on canonical resolution of singularities, and a model-theoretic compactness argument.


title: Orthogonal Chirp Division Multiplexing for Power Line Sensing via Time-Domain Reflectometry

abstract: In this study, a time-domain reflectometry (TDR) system based on a baseband version of the orthogonal chirp-division multiplexing (OCDM) scheme, which is relies on a modified discrete Fresnel transform (DFnT), is proposed for power line sensing. After a detailed description of the system model, a multiple access scheme that exploits the contolution theorem of the modified DFnT for enabling distributed reflectometric and transferometric sensing of the monitored power line. Considering typical European underground low-voltage and US overhead medium-voltage (MV) power distribution networks, range resolution and maximum unambiguous range for the measurements are assessed for the proposed scheme. Next, a comparison with multiple access schemes based on the Hermitian symmetric orthogonal frequency-division multiplexing (HS-OFDM) discussed in a previous work is performed considering a Brazilian MV overhead scenario, being the number of measurements obtained over time, as well as signal-to-interference-plus-noise ratio (SINR) used as performance metrics. The attained results show that the proposed multiple access scheme results in the same range resolution as the others. Also, the highest number of measurements over time is obtained by the proposed scheme, which produces orthogonality among signals transmitted by different power line modems (PLMs) neither and time nor frequency domains, but rather in the Fresnel domain. Meanwhile, although its yielded SINR values are fair among the PLMs consisting the distributed sensing system, the proposed scheme is slightly outperformed by the HS-OFDM based on code-division multiple access and by the HS-OFDM based on frequency-division multiple access at some PLMs.


title: Guaranteed lower bounds for cost functionals of time-periodic parabolic optimization problems

abstract: In this paper, a new technique is shown for deriving computable, guaranteed lower bounds of functional type (minorants) for two different cost functionals subject to a parabolic time-periodic boundary value problem. Together with previous results on upper bounds (majorants) for one of the cost functionals, both minorants and majorants lead to two-sided estimates of functional type for the optimal control problem. Both upper and lower bounds are derived for the second new cost functional subject to the same parabolic PDE-constraints, but where the target is a desired gradient. The time-periodic optimal control problems are discretized by the multiharmonic finite element method leading to large systems of linear equations having a saddle point structure. The derivation of preconditioners for the minimal residual method for the new optimization problem is discussed in more detail. Finally, several numerical experiments for both optimal control problems are presented confirming the theoretical results obtained. This work provides the basis for an adaptive scheme for time-periodic optimization problems.


title: Testing the Impact of Satellite Anisotropy on Large and Small Scale Intrinsic Alignments using Hydrodynamical Simulations

abstract: Galaxy intrinsic alignments (IAs) have long been recognised as a significant contaminant to weak lensing-based cosmological inference. In this paper we seek to quantify the impact of a common modelling assumption in analytic descriptions of intrinsic alignments: that of spherically symmetric dark matter halos. Understanding such effects is important as the current generation of intrinsic alignment models are known to be limited, particularly on small scales, and building an accurate theoretical description will be essential for fully exploiting the information in future lensing data. Our analysis is based on a catalogue of 113,560 galaxies between $z=0.06-1.00$ from MassiveBlack-II, a hydrodynamical simulation of box length $100 h^{-1}$ Mpc. We find satellite anisotropy contributes at the level of $\geq 30-40\%$ to the small scale alignment correlation functions. At separations larger than $1 h^{-1}$ Mpc the impact is roughly scale-independent, inducing a shift in the amplitude of the IA power spectra of $\sim20\%$. These conclusions are consistent across the redshift range and between the MassiveBlack-II and Illustris simulations. The cosmological implications of these results are tested using a simulated likelihood analysis. Synthetic cosmic shear data is constructed with the expected characteristics (depth, area and number density) of a future LSST-like survey. Our results suggest that modelling alignments using a halo model based upon spherical symmetry could potentially induce cosmological parameter biases at the $\sim 1.5\sigma$ level for $S_8$ and $w$.


title: Hear the Sound of Weyl Fermions

abstract: Quasiparticles and collective modes are two fundamental aspects that characterize a quantum matter in addition to its ground state features. For example, the low energy physics for Fermi liquid phase in He-III was featured not only by Fermionic quasiparticles near the chemical potential but also by fruitful collective modes in the long-wave limit, including several different sound waves that can propagate through it under different circumstances. On the other hand, it is very difficult for sound waves to be carried by the electron liquid in the ordinary metals, due to the fact that long-range Coulomb interaction among electrons will generate plasmon gap for ordinary electron density fluctuation and thus prohibits the propagation of sound waves through it. In the present paper, we propose a unique type of acoustic collective modes in Weyl semimetals under the magnetic field called chiral zero sound. The chiral zero sound can be stabilized under so-called "chiral limit", where the intra-valley scattering time is much shorter than the inter-valley one, and only propagates along an external magnetic field for Weyl semimetals with multiple-pairs of Weyl points. The sound velocity of the chiral zero sound is proportional to the field strength in the weak field limit, whereas it oscillates dramatically in the strong field limit, generating an entirely new mechanism for quantum oscillations through the dynamics of neutral bosonic excitation, which may manifest itself in the thermal conductivity measurements under magnetic field.


title: First study of reionization in tilted flat and untilted non-flat dynamical dark energy inflation models

abstract: We examine the effects of dark energy dynamics and spatial curvature on cosmic reionization by studying reionization in tilted spatially-flat and untilted non-flat XCDM and $\phi$CDM dynamical dark energy inflation models that best fit the Planck 2015 cosmic microwave background (CMB) anisotropy and a large compilation of non-CMB data. We carry out a detailed statistical study, based on a principal component analysis and a Markov chain Monte Carlo analysis of a compilation of lower-redshift reionization data, to estimate the uncertainties in the cosmological model reionization histories. We find that, irrespective of the nature of dark energy, there are significant differences between the reionization histories of the spatially-flat and non-flat models. Although both the flat and non-flat models can accurately match the low-redshift ($z\lesssim6$) reionization observations, there is a clear discrepancy between high-redshift ($z>7$) Lyman-$\alpha$ emitter data and the predictions from non-flat models. This is solely due to the fact that the non-flat models have a significantly larger electron scattering optical depth, $\tau_{\rm el}$, compared to the flat models, which requires an extended and much earlier reionization scenario supported by more high-redshift ionizing sources in the non-flat models. Non-flat models also require strong redshift evolution in the photon escape fraction, that can become unrealistically high ($\gtrsim1$) at some redshifts. However, $\tau_{\rm el}$ is about 0.9-$\sigma$ lower in the tilted flat $\Lambda$CDM model when the new Planck 2018 data are used and this reduction will partially alleviate the tension between the non-flat model predictions and the data.


title: Soft anomalous dimensions for single-top production at three loops

abstract: I present results for soft anomalous dimensions in single top-quark production processes through three loops. I first discuss the cusp anomalous dimension with one massive and one massless line at three loops. I then calculate the three-loop soft anomalous dimension for $tW$ production and for related processes with new physics such as $tZ$, $t\gamma$, and $tH^-$ production. This is followed by results at two and three loops for the soft anomalous dimension matrix in $s$-channel single-top production and, finally, in $t$-channel single-top production. These new results are needed for soft-gluon resummation at N$^3$LL accuracy and for calculations of soft-gluon corrections at N$^3$LO.


title: On the BPS Spectrum of the Rank-1 Minahan-Nemeschansky Theories

abstract: The F-theory realization of the rank-1 Minahan Nemeschansky (MN) $E_6$, $E_7$ and $E_8$ theories leads to a description of the BPS states on the Coulomb branch in terms of Type IIB $(p,q)$-string networks. Subject to a simple ansatz for the types of networks which can occur, we study the representations of the flavor symmetry group which occur in the BPS spectrum. The results we find for the $E_6$ and $E_7$ theories are in perfect agreement with previous calculations by other methods (in particular, we find that arbitrarily large representations occur), but our scheme is easier to implement and more computationally efficient. The string network picture also gives a possible explanation of the experimental observation that in rank-1 MN theories, BPS states whose charge is $n$ times a primitive charge occur with BPS index divisible by $(-1)^{n+1} n$.


title: Short Term Variability of Evolved Massive Stars with TESS

abstract: We present the first results from a study of TESS Sector 1 and 2 light curves for eight evolved massive stars in the LMC: six yellow supergiants (YSGs) and two luminous blue variables (LBVs), including S Doradus. We use an iterative prewhitening procedure to characterize the short-timescale variability in all eight stars. The periodogram of one of the YSGs, HD 269953, displays multiple strong peaks at higher frequencies than its fellows. While the field surrounding HD 269953 is quite crowded, it is the brightest star in the region, and has infrared colors indicating it is dusty. We suggest HD 269953 may be in a post-red supergiant evolutionary phase. We find a signal with a period of $\sim5$ days for the LBV HD 269582. The periodogram of S Doradus shows a complicated structure, with peaks below frequencies of 1.5 cycles per day. We fit the shape of the background noise of all eight light curves, and find a red noise component in all of them. However, the power law slope of the red noise and the timescale over which coherent structures arise changes from star to star. Our results highlight the potential for studying evolved massive stars with TESS.


abstract: We prove an optimal Alexandrov-Bakelman-Pucci type estimate for plurisubharmonic functions without assuming their continuity. This generalizes a result of Y. Wang. As a corollary we generalize an estimate from \cite{DD18}. We also address a problem posed in \cite{Wan12}.


title: Effects of a fully 3D atmospheric structure on exoplanet transmission spectra: retrieval biases due to day-night temperature gradients

abstract: Transmission spectroscopy provides us with information on the atmospheric properties at the limb, which is often intuitively assumed to be a narrow annulus around the planet. Consequently, studies have focused on the effect of atmospheric horizontal heterogeneities along the limb. Here we demonstrate that the region probed in transmission -- the limb -- extends significantly toward the day and night sides of the planet. We show that the strong day-night thermal and compositional gradients expected on synchronous exoplanets create sufficient heterogeneities across the limb to result in important systematic effects on the spectrum and bias its interpretation. To quantify these effects, we developed a 3D radiative transfer model able to generate transmission spectra of atmospheres based on 3D atmospheric structures. We first produce synthetic JWST observations from a simulation of GJ 1214b and demonstrate the necessity of a real 3D approach to model data for such precise observatories. Second, we investigate how day-night temperature gradients cause a systematic bias in retrieval analysis performed with 1D forward models. For that purpose we synthesize a large set of forward spectra for prototypical HD 209458 b and GJ 1214 b type planets varying the temperatures of the day and night sides as well as the width of the transition region. We then perform typical retrieval analyses and compare the retrieved parameters to the ground truth of the input model. This study reveals systematic biases on the retrieved temperature (found to be higher than the terminator temperature) and abundances. This is due to the fact that the hotter dayside is more extended vertically and screens the nightside. These biases are difficult to detect as the 1D profiles used in the retrieval are found to provide an excellent match to the observed spectra. This needs to be kept in mind when interpreting real data.


title: Finite temperature effects on Majorana bound states in chiral $p$-wave superconductors

abstract: We study Majorana fermions bound to vortex cores in a chiral $p$-wave superconductor at temperatures non-negligible compared to the superconducting gap. Thermal occupation of Caroli de Gennes-Matricon states, below the full gap, causes the free energy difference between the two fermionic parity sectors to decay algebraically with increasing temperature. The power law acquires an additional factor of $T^{-1}$ for each bound state thermally excited. The zero-temperature result is exponentially recovered well below the minigap (lowest-lying CdGM level). Our results suggest that temperatures larger than the minigap may not be disastrous for topological quantum computation. We discuss the prospect of precision measurements of pinning forces on vortices as a readout scheme for Majorana qubits.


title: Planetary Giant Impacts: Convergence of High-Resolution Simulations using Efficient Spherical Initial Conditions and SWIFT

abstract: We perform simulations of giant impacts onto the young Uranus using smoothed particle hydrodynamics (SPH) with over 100 million particles. This 100--1000$\times$ improvement in particle number reveals that simulations with below 10^7 particles fail to converge on even bulk properties like the post-impact rotation period, or on the detailed erosion of the atmosphere. Higher resolutions appear to determine these large-scale results reliably, but even 10^8 particles may not be sufficient to study the detailed composition of the debris -- finding that almost an order of magnitude more rock is ejected beyond the Roche radius than with 10^5 particles. We present two software developments that enable this increase in the feasible number of particles. First, we present an algorithm to place any number of particles in a spherical shell such that they all have an SPH density within 1% of the desired value. Particles in model planets built from these nested shells have a root-mean-squared velocity below 1% of the escape speed, which avoids the need for long precursor simulations to produce relaxed initial conditions. Second, we develop the hydrodynamics code SWIFT for planetary simulations. SWIFT uses task-based parallelism and other modern algorithmic approaches to take full advantage of contemporary supercomputer architectures. Both the particle placement code and SWIFT are publicly released.


title: Dressed states from gauge invariance

abstract: The dressed state formalism enables us to define the infrared finite $S$-matrix for QED. In the formalism, asymptotic charged states are dressed by clouds of photons. The dressed asymptotic states are originally obtained by solving the dynamics of the asymptotic Hamiltonian in the far past or future region. However, there was an argument that the obtained dressed states are not gauge invariant. We resolve the problem by imposing a correct gauge invariant condition. We show that the dressed states can be obtained just by requiring the gauge invariance of asymptotic states. In other words, Gauss's law naturally leads to proper asymptotic states for the infrared finite $S$-matrix. We also discuss the relation between the dressed state formalism and the asymptotic symmetry for QED.


title: Detector-size Upper Bounds on Dark Hadron Lifetime from Cosmology

abstract: We show that in a confining hidden valley model where the lightest hidden particles are dark hadrons that have mass splittings larger than $\mathcal{O}(0.1)$ GeV, if the lightest dark hadron is either stable or decays into Standard Model (SM) hadrons/charged leptons during the big-bang nucleosynthesis (BBN), at least one of the heavier dark hadrons needs to decay into SM particles within $\mathcal{O}(10)$ nanosec. Once being produced at collider experiments, this heavier dark hadron is likely to decay within $\mathcal{O}(1)$ meter distance, which strengthens the motivation of searching for long-lived particles with sub-meter scale decay lengths at colliders. To illustrate the idea, we study the lifetime constraint in scenarios where the lightest dark particle is a pseudo-scalar meson, and dark hadrons couple to SM particles either through kinetic mixing between the SM and dark photons or by mixing between the SM and dark Higgs. We study the annihilation and decay of dark hadrons in a thermal bath and calculate upper bounds on the lightest vector meson (scalar hadron) lifetime in the kinetic mixing (Higgs portal) scenario. We discuss the application of these lifetime constraints in long-lived particle searches that use the LHCb VELO or the ATLAS/CMS inner detectors.


title: Sleptons without Hadrons

abstract: Multilepton searches for electroweakino and slepton pair production at hadron colliders remain some of the best means to test weak-scale supersymmetry. Searches at the CERN Large Hadron Collider, however, are limited by large diboson and top quark pair backgrounds, despite the application of traditional, central jet vetoes. In this context, we report the impact of introducing dynamic jet vetoes in searches for colorless superpartners. As a representative scenario, we consider the Drell-Yan production of a pair of right-handed smuons decaying into a dimuon system accompanied with missing transverse energy. As an exploratory step, we consider several global and local measures of the leptonic and hadronic activity to construct the veto. In most all cases, we find that employing a dynamic jet veto improves the sensitivity, independently of the integrated luminosity. The inclusion of non-perturbative multiple particle interactions and next-to-leading order jet merging does not alter this picture. Directions for further improvements are discussed.