title: On the linkability of Zcash transactions

abstract: Zcash is a fork of Bitcoin with optional anonymity features. While transparent transactions are fully linkable, shielded transactions use zero-knowledge proofs to obscure the parties and amounts of the transactions. First, we observe various metrics regarding the usage of shielded addresses. Moreover, we show that most coins sent to shielded addresses are later sent back to transparent addresses. We then search for round-trip transactions, where the same, or nearly the same number of coins are sent from a transparent address, to a shielded address, and back again to a transparent address. We argue that such behavior exhibits high linkability, especially when they occur nearby temporally. Using this heuristic our analysis matched 31.5% of all coins sent to shielded addresses.


title: A Unified Study of Continuous and Discontinuous Galerkin Methods

abstract: A unified study is presented in this paper for the design and analysis of different finite element methods (FEMs), including conforming and nonconforming FEMs, mixed FEMs, hybrid FEMs,discontinuous Galerkin (DG) methods, hybrid discontinuous Galerkin (HDG) methods and weak Galerkin (WG) methods. Both HDG and WG are shown to admit inf-sup conditions that hold uniformly with respect to both mesh and penalization parameters. In addition, by taking the limit of the stabilization parameters, a WG method is shown to converge to a mixed method whereas an HDG method is shown to converge to a primal method. Furthermore, a special class of DG methods, known as the mixed DG methods, is presented to fill a gap revealed in the unified framework.


title: Bosonic Integer Quantum Hall States without Landau Levels on Square Lattice

abstract: We study an interacting two-component hard-core bosons on square lattice for which, in the presence of staggered magnetic flux, the ground state is a bosonic integer quantum Hall (BIQH) state. Using a coupled-wire bosonization approach, we analytically show this model exhibits a BIQH state at total charge half filling associated with a symmetry-protected topological phase under $U(1)$ charge conservation. These theoretical expectations are verified, using the infinite density matrix renormalization group method, by providing numerical evidences for: (i) a quantized Hall conductance $\sigma_{xy}=\pm2$, and (ii) two counter-propagating gapless edge modes. Our model is a bosonic cousin of the fermionic Haldane model and serves as an additional case of analogy between bosonic and fermionic quantum Hall states.


title: An Encoder-Decoder Model for ICD-10 Coding of Death Certificates

abstract: Information extraction from textual documents such as hospital records and healthrelated user discussions has become a topic of intense interest. The task of medical concept coding is to map a variable length text to medical concepts and corresponding classification codes in some external system or ontology. In this work, we utilize recurrent neural networks to automatically assign ICD-10 codes to fragments of death certificates written in English. We develop end-to-end neural architectures directly tailored to the task, including basic encoder-decoder architecture for statistical translation. In order to incorporate prior knowledge, we concatenate cosine similarities vector among the text and dictionary entry to the encoded state. Being applied to a standard benchmark from CLEF eHealth 2017 challenge, our model achieved F-measure of 85.01% on a full test set with significant improvement as compared to the average score of 62.2% for all official participants approaches.


title: TRIFECTA: Security, Energy-Efficiency, and Communication Capacity Comparison for Wireless IoT Devices

abstract: The widespread proliferation of sensor nodes in the era of Internet of Things (IoT) coupled with increasing sensor fidelity and data acquisition modality is expected to generate 3+ Exabytes of data per day by 2018. Since most of these IoT devices will be wirelessly connected at the last few feet, wireless communication is an integral part of the future IoT scenario. The ever-shrinking size of unit computation (Moore's Law) and continued improvements in efficient communication (Shannon's Law) is expected to harness the true potential of the IoT revolution and produce dramatic societal impact. However, reducing size of IoT nodes and lack of significant improvement in energy-storage density leads to reducing energy-availability. Moreover, smaller size and energy means less resources available for securing IoT nodes, making the energy-sparse low-cost leaf nodes of the network as prime targets for attackers. In this paper, we survey six prominent wireless technologies with respect to the three dimensions - security, energy efficiency, and communication capacity. We point out the state-of-the-art, open issues, and the road ahead for promising research directions.


title: Exact solutions, finite time singularities and non-singular universe models from a variety of $\Lambda(t)$ cosmologies

abstract: Cosmological models with time dependent $\Lambda$ (read as $\Lambda (t)$) have been investigated widely in the literature. Models that solve background dynamics analytically, are of special interest. Additionally, the allowance of past or future singularities at finite cosmic time in a specific model signals for a generic test on its viabilities with the current observations. Following these, in this work we consider a variety of $\Lambda (t)$ models focusing on their evolutions and singular behaviour. We found that a series of models in this class can be exactly solved when the background universe is described by a spatially flat Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) line element. The solutions in terms of the scale factor of the FLRW universe offer different universe models, such as power law expansion, oscillating, and the singularity free universe as well. However, we also noticed that a large number of models in this series permits past or future cosmological singularities at finite cosmic time. At last we close the work with a note that the avoidance of future singularities is possible for certain models under some specific restrictions.


title: Large Deviation Function for a Driven Underdamped Particle in a Periodic Potential

abstract: Employing large deviation theory, we explore current fluctuations of underdamped Brownian motion for the paradigmatic example of a single particle in a one dimensional periodic potential. Two different approaches to the large deviation function of the particle current are presented. First, we derive an explicit expression for the large deviation functional of the empirical phase space density, which replaces the level 2.5 functional used for overdamped dynamics. Using this approach, we obtain several bounds on the large deviation function of the particle current. We compare these to bounds for overdamped dynamics that have recently been derived motivated by the thermodynamic uncertainty relation. Second, we provide a method to calculate the large deviation function via the cumulant generating function. We use this method to assess the tightness of the bounds in a numerical case study for a cosine potential.


title: Iterative Deep Learning for Network Topology Extraction

abstract: This paper tackles the task of estimating the topology of filamentary networks such as retinal vessels and road networks. Building on top of a global model that performs a dense semantical classification of the pixels of the image, we design a Convolutional Neural Network (CNN) that predicts the local connectivity between the central pixel of an input patch and its border points. By iterating this local connectivity we sweep the whole image and infer the global topology of the filamentary network, inspired by a human delineating a complex network with the tip of their finger. We perform an extensive and comprehensive qualitative and quantitative evaluation on two tasks: retinal veins and arteries topology extraction and road network estimation. In both cases, represented by two publicly available datasets (DRIVE and Massachusetts Roads), we show superior performance to very strong baselines.


title: Experimental investigation of enviroment--induced entanglement using an all--optical setup

abstract: We investigate the generation of entanglement between two non interacting qubits coupled to a common reservoir. An experimental setup was conceived to encode one qubit on the polarization of an optical beam and another qubit on its transverse mode. The action of the reservoir is implemented as conditional operations on these two qubits, controlled by the longitudinal path as an ancillary degree of freedom. An entanglement witness and the two-qubit concurrence are easily evaluated from direct intensity measurements showing an excellent agreement with the theoretical prediction.


title: Star Cluster Formation in Cosmological Simulations. II. Effects of Star Formation Efficiency and Stellar Feedback

abstract: The implementation of star formation and stellar feedback in cosmological simulations plays a critical role in shaping galaxy properties. In the first paper of the series, we presented a new method to model star formation as a collection of star clusters. In this paper, we improve the algorithm by eliminating accretion gaps, boosting momentum feedback, and introducing a subgrid initial bound fraction, $f_i$, that distinguishes cluster mass from stellar particle mass. We perform a suite of simulations with different star formation efficiency per freefall time $\epsilon_{\rm ff}$ and supernova momentum feedback intensity $f_{\rm boost}$. We find that the star formation history of a Milky Way-sized galaxy is sensitive to $f_{\rm boost}$, which allows us to constrain its value, $f_{\rm boost}\approx5$, in the current simulation setup. Changing $\epsilon_{\rm ff}$ from a few percent to 200\% has little effect on global galaxy properties. However, on smaller scales, the properties of star clusters are very sensitive to $\epsilon_{\rm ff}$. We find that $f_i$ increases with $\epsilon_{\rm ff}$ and cluster mass. Through the dependence on $f_i$, the shape of the cluster initial mass function varies strongly with $\epsilon_{\rm ff}$. The fraction of clustered star formation and maximum cluster mass increase with the star formation rate surface density, with the normalization of both relations dependent on $\epsilon_{\rm ff}$. The cluster formation timescale systematically decreases with increasing $\epsilon_{\rm ff}$. Local variations in the gas accretion history lead to a 0.25~dex scatter for the integral cluster formation efficiency. Joint constraints from all the observables prefer the runs that produce a median integral efficiency of 16\%.


title: Observation of fractional quantum Hall effect in an InAs quantum well

abstract: The two-dimensional electron system in an InAs quantum well has emerged as a prime candidate for hosting exotic quasi-particles with non-Abelian statistics such as Majorana fermions and parafermions. To attain its full promise, however, the electron system has to be clean enough to exhibit electron-electron interaction phenomena. Here we report the observation of fractional quantum Hall effect in a very low disorder InAs quantum well with a well-width of 24 nm, containing a two-dimensional electron system with a density $n=7.8 \times 10^{11}$ cm$^{-2}$ and low-temperature mobility $1.8 \times 10^6$ cm$^2$/Vs. At a temperature of $\simeq35$ mK and $B\simeq24$ T, we observe a deep minimum in the longitudinal resistance, accompanied by a nearly quantized Hall plateau at Landau level filling factor $\nu=4/3$.


title: Stable pairs, flat connections and Gopakumar-Vafa invariants

abstract: Using the interpretation of certain generalised Donaldson-Thomas invariants, including stable pairs curve counts, as the monodromy of a flat connection on a formal principal bundle, we show that the conjectural Gopakumar-Vafa contributions of all genera to the Gromov-Witten partition function appear in the asymptotics of the corresponding flat sections. The Fourier-Laplace integrals used to produce flat sections lead naturally to the GW/DT change of variable -q = e^{i u}.


title: The JKind Model Checker

abstract: JKind is an open-source industrial model checker developed by Rockwell Collins and the University of Minnesota. JKind uses multiple parallel engines to prove or falsify safety properties of infinite state models. It is portable, easy to install, performance competitive with other state-of-the-art model checkers, and has features designed to improve the results presented to users: inductive validity cores for proofs and counterexample smoothing for test-case generation. It serves as the back-end for various industrial applications.


title: Exact controllability for string with attached masses

abstract: We consider the problem of boundary control for a vibrating string with $N$ interior point masses. We assume the control is at the left end, and the string is fixed at the right end. Singularities in waves are "smoothed" out to one order as they cross a point mass. We characterize the reachable set for a $L^2$ control. The control problem is reduced to a moment problem, which is then solved using the theory of exponential divided differences in tandem with unique shape and velocity controllability results.


title: The random gas of hard spheres

abstract: The inconsistency between the time-reversible Liouville equation and time-irreversible Boltzmann equation has been pointed out long ago by Loschmidt. To avoid Loschmidt's objection, here we propose a new dynamical system to model the motion of atoms of gas, with their interactions triggered by a random point process. Despite being random, this model can approximate the collision dynamics of rigid spheres via adjustable parameters. We compute the exact statistical steady state of the system, and determine the form of its marginal distributions for a large number of spheres. We find that the Kullback-Leibler entropy (a generalization of the conventional Boltzmann entropy) of the full system of random gas spheres is a nonincreasing function of time. Unlike the conventional hard sphere model, the proposed random gas model results in a variant of the Enskog equation, which is known to be a more accurate model of dense gas than the Boltzmann equation. We examine the hydrodynamic limit of the derived Enskog equation for spheres of constant mass density, and find that the corresponding Enskog-Euler and Enskog-Navier-Stokes equations acquire additional effects in both the advective and viscous terms. In the dilute gas approximation, the Enskog equation simplifies to the Boltzmann equation, while the Enskog-Euler and Enskog-Navier-Stokes equations become the conventional Euler and Navier-Stokes equations.


title: Almost Quantum Correlations are Inconsistent with Specker's Principle

abstract: Ernst Specker considered a particular feature of quantum theory to be especially fundamental, namely that pairwise joint measurability of sharp measurements implies their global joint measurability (https://vimeo.com/52923835). To date, Specker's principle seemed incapable of singling out quantum theory from the space of all general probabilistic theories. In particular, its well-known consequence for experimental statistics, the principle of consistent exclusivity, does not rule out the set of correlations known as almost quantum, which is strictly larger than the set of quantum correlations. Here we show that, contrary to the popular belief, Specker's principle cannot be satisfied in any theory that yields almost quantum correlations.


title: On Capacity-Achieving Distributions for Complex AWGN Channels Under Nonlinear Power Constraints and their Applications to SWIPT

abstract: The capacity of a complex and discrete-time memoryless additive white Gaussian noise (AWGN) channel under three constraints, namely, input average power, input amplitude and output delivered power is studied. The output delivered power constraint is modelled as the average of linear combination of even moments of the channel input being larger than a threshold. It is shown that the capacity of an AWGN channel under transmit average power and receiver delivered power constraints is the same as the capacity of an AWGN channel under an average power constraint. However, depending on the two constraints, the capacity can be either achieved by a Gaussian distribution or arbitrarily approached by using time-sharing between a Gaussian distribution and On-Off Keying. As an application, a simultaneous wireless information and power transfer (SWIPT) problem is studied, where an experimentally-validated nonlinear model of the harvester is used. It is shown that the delivered power depends on higher order moments of the channel input. Two inner bounds, one based on complex Gaussian inputs and the other based on further restricting the delivered power are obtained for the Rate-Power (RP) region. For Gaussian inputs, the optimal inputs are zero mean and a tradeoff between transmitted information and delivered power is recognized by considering asymmetric power allocations between inphase and quadrature subchannels. Through numerical algorithms, it is observed that input distributions (obtained by restricting the delivered power) attain larger RP region compared to Gaussian input counterparts.


title: Determinacy of Schmidt's Game and Other Intersection Games

abstract: Schmidt's game, and other similar intersection games have played an important role in recent years in applications to number theory, dynamics, and Diophantine approximation theory. These games are real games, that is, games in which the players make moves from a complete separable metric space. The determinacy of these games trivially follows from the axiom of determinacy for real games, $\mathsf{AD}_\mathbb{R}$, which is a much stronger axiom than that asserting all integer games are determined, $\mathsf{AD}$. One of our main results is a general theorem which under the hypothesis $\mathsf{AD}$ implies the determinacy of intersection games which have a property allowing strategies to be simplified. In particular, we show that Schmidt's $(\alpha,\beta,\rho)$ game on $\mathbb{R}$ is determined from $\mathsf{AD}$ alone, but on $\mathbb{R}^n$ for $n \geq 3$ we show that $\mathsf{AD}$ does not imply the determinacy of this game. We also prove several other results specifically related to the determinacy of Schmidt's game. These results highlight the obstacles in obtaining the determinacy of Schmidt's game from $\mathsf{AD}$.


title: Why surplus structure is not superfluous

abstract: The idea that gauge theory has 'surplus' structure poses a puzzle: in one much discussed sense, this structure is redundant; but on the other hand, it is also widely held to play an essential role in the theory. In this paper, we employ category-theoretic tools to illuminate an aspect of this puzzle. We precisify what is meant by 'surplus' structure by means of functorial comparisons with equivalence classes of gauge fields, and then show that such structure is essential for any theory that represents a rich collection of physically relevant fields which are 'local' in nature.


title: The SLUGGS Survey: The Inner Dark Matter Density Slope of the Massive Elliptical Galaxy NGC 1407

abstract: We investigate the dark matter density profile of the massive elliptical galaxy, NGC 1407, by constructing spherically symmetric Jeans models of its field star and globular cluster systems. Two major challenges in such models are the degeneracy between the stellar mass and the dark matter halo profiles, and the degeneracy between the orbital anisotropy of the tracer population and the total mass causing the observed motions. We address the first issue by using new measurements of the mass-to-light ratio profile from stellar population constraints that include a radially varying initial mass function. To deal with the mass--anisotropy degeneracy, we make use of multiple kinematic tracers, including two subpopulations of globular clusters in addition to the galaxy's field stars. We create a hierarchical Bayesian model that addresses several often neglected systematic uncertainties, such as the statistical weight given to various datasets and the adopted distance. After sampling the posterior probability distribution with a Markov Chain Monte Carlo method, we find evidence for a central cusp with a log-slope of $\gamma = 1.0^{+0.3}_{-0.4}$. This is lower than expected for dark matter halos that have undergone adiabatic contraction, supporting inferences from gravitational lensing that some process has suppressed the steepening of halos in massive galaxies. We also confirm radially-biased orbits for the metal-rich globular clusters and tangentially-biased orbits for the metal-poor globular clusters, which remains a puzzling finding for an accretion-dominated halo.


title: Charm Penguin in $B^\pm \to K^\pm K^+ K^-$: partonic and hadronic loops

abstract: Charm penguin diagrams are known to be the main contribution to charmless B decay process with strangeness variation equal to minus one, which is the case of $B^\pm \to K^\pm K^+ K^-$ decay. The large phase space available in this and other B three-body decays allows non trivial final state interactions with all sort of rescattering processes and also access high momentum transfers in the central region of the Dalitz plane. In this work we investigate the charm Penguin contribution to $B^\pm \to K^\pm K^+ K^-$, described by a hadronic triangle loop in nonperturbative regions of the phase space, and by a partonic loop at the quasi perturbative region. These nonresonant amplitudes should have a particular structure in the Dalitz plane and their contributions to the final decay amplitude can be confirmed by a data amplitude analysis in this channel. In particular, the hadronic amplitude has a changing sign in the phase at $D\bar{D}$ threshold which can result in a change of sign for the CP asymmetry.


title: Sub-clustering in decomposable graphs and size-varying junction trees

abstract: This paper proposes a novel representation of decomposable graphs based on semi-latent tree-dependent bipartite graphs. The novel representation has two main benefits. First, it enables a form of sub-clustering within maximal cliques of the graph, adding informational richness to the general use of decomposable graphs that could be harnessed in applications with behavioural type of data. Second, it allows for a new node-driven Markov chain Monte Carlo sampler of decomposable graphs that can easily parallelize and scale. The proposed sampler also benefits from the computational efficiency of junction-tree-based samplers of decomposable graphs.


title: Minimal Input Selection for Robust Control

abstract: This paper studies the problem of selecting a minimum-size set of input nodes to guarantee stability of a networked system in the presence of uncertainties and time delays. Current approaches to input selection in networked dynamical systems focus on nominal systems with known parameter values in the absence of delays. We derive sufficient conditions for existence of a stabilizing controller for an uncertain system that are based on a subset of system modes lying within the controllability subspace induced by the set of inputs. We then formulate the minimum input selection problem and prove that it is equivalent to a discrete optimization problem with bounded submodularity ratio, leading to polynomial-time algorithms with provable optimality bounds. We show that our approach is applicable to different types of uncertainties, including additive and multiplicative uncertainties in the system matrices as well as uncertain time delays. We demonstrate our approach in a numerical case study on the IEEE 39-bus test power system.


title: Andreev levels as a first approach to quantum computing with high Tc superconductors

abstract: A computing platform based on low temperature superconductors (LTS) has already been proven both theoretically and experimentally. However, qubits based on high Tc superconductors (HTS) are not yet well understood. Here we study the Andreev bounds states (ABS) in the later materials in order to show that a formal correspondence exists between the Mathieu levels in a Cooper Pair Box qubit built with LTS and the Andreev levels in HTS junctions.


title: Structure of temporal correlations of a qubit

abstract: In quantum mechanics, spatial correlations arising from measurements at separated particles are well studied. This is not the case, however, for the temporal correlations arising from a single quantum system subjected to a sequence of generalized measurements. We first characterize the polytope of temporal quantum correlations coming from the most general measurements. We then show that if the dimension of the quantum system is bounded, only a subset of the most general correlations can be realized and identify the correlations in the simplest scenario that can not be reached by two-dimensional systems. This leads to a temporal inequality for a dimension test, and we discuss a possible implementation using nitrogen-vacancy centers in diamond.


title: On the Real-time Vehicle Placement Problem

abstract: Motivated by ride-sharing platforms' efforts to reduce their riders' wait times for a vehicle, this paper introduces a novel problem of placing vehicles to fulfill real-time pickup requests in a spatially and temporally changing environment. The real-time nature of this problem makes it fundamentally different from other placement and scheduling problems, as it requires not only real-time placement decisions but also handling real-time request dynamics, which are influenced by human mobility patterns. We use a dataset of ten million ride requests from four major U.S. cities to show that the requests exhibit significant self-similarity. We then propose distributed online learning algorithms for the real-time vehicle placement problem and bound their expected performance under this observed self-similarity.


title: Negative Differential Mobility in Interacting Particle Systems

abstract: Driven particles in presence of crowded environment, obstacles or kinetic constraints often exhibit negative differential mobility (NDM) due to their decreased dynamical activity. We propose a new mechanism for complex many-particle systems where slowing down of certain {\it non-driven} degrees of freedom by the external field can give rise to NDM. This phenomenon, resulting from inter-particle interactions, is illustrated in a pedagogical example of two interacting random walkers, one of which is biased by an external field while the same field only slows down the other keeping it unbiased. We also introduce and solve exactly the steady state of several driven diffusive systems, including a two species exclusion model, asymmetric misanthrope and zero-range processes, to show explicitly that this mechanism indeed leads to NDM.


title: Electron scattering wings on lines in interacting supernovae

abstract: We consider the effect of electron scattering on lines emitted as a result of supernova interaction with a circumstellar medium, assuming that the scattering occurs in ionized gas in the preshock circumstellar medium. The single scattering case gives the broad component in the limit of low optical depth, showing a velocity full width half maximum that is close to the thermal velocities of electrons. The line shape is approximately exponential at low velocities and steepens at higher velocities. At higher optical depths, the line profile remains exponential at low velocities, but wings strengthen with increasing optical depth. In addition to the line width, the ratio of narrow to broad (scattered) line strength is a possible diagnostic of the gas. The results depend on the density profile of the circumstellar gas, especially if the scattering and photon creation occur in different regions. We apply the scattering model to a number of supernovae, including Type IIn and Type Ia-CSM events. The asymmetry to the red found in some cases can be explained by scattering in a fast wind region which is indicated by observations.


title: Learning by Asking Questions

abstract: We introduce an interactive learning framework for the development and testing of intelligent visual systems, called learning-by-asking (LBA). We explore LBA in context of the Visual Question Answering (VQA) task. LBA differs from standard VQA training in that most questions are not observed during training time, and the learner must ask questions it wants answers to. Thus, LBA more closely mimics natural learning and has the potential to be more data-efficient than the traditional VQA setting. We present a model that performs LBA on the CLEVR dataset, and show that it automatically discovers an easy-to-hard curriculum when learning interactively from an oracle. Our LBA generated data consistently matches or outperforms the CLEVR train data and is more sample efficient. We also show that our model asks questions that generalize to state-of-the-art VQA models and to novel test time distributions.


title: DAMPE Electron-Positron Excess in Leptophilic $Z'$ model

abstract: Recently the DArk Matter Particle Explorer (DAMPE) has reported an excess in the electron-positron flux of the cosmic rays which is interpreted as a dark matter particle with the mass about $1.5$ TeV. We come up with a leptophilic $Z'$ scenario including a Dirac fermion dark matter candidate which beside explaining the observed DAMPE excess, is able to pass various experimental/observational constraints including the relic density value from the WMAP/Planck, the invisible Higgs decay bound at the LHC, the LEP bounds in electron-positron scattering, the muon anomalous magnetic moment constraint, Fermi-LAT data, and finally the direct detection experiment limits from the XENON1t/LUX. By computing the electron-positron flux produced from a dark matter with the mass about $1.5$ TeV we show that the model predicts the peak observed by the DAMPE.


title: Bragg-induced oscillations in non-PT complex photonic lattices

abstract: When a monochromatic beam of light propagates through a periodic structure with the incident angle satisfying the Bragg condition, its Fourier spatial spectra oscillates between the resonant modes situated on the edges of the Brillouin zones of the lattice with a nontrivial dynamics. We wish to study here these Bragg-induced oscillations in a specific non-PT periodic structure, that is, a periodic medium not invariant under the combined action of parity and time reversal symmetries. We compare our analytic results based on the expansion of the optical field in Bragg-resonant plane waves with direct numerical integration of the paraxial wave equation using a wide Gaussian beam as initial condition without assuming the shallow potential approximation. In particular, we wish to investigate under what conditions a mode trapping phenomenon can still be observed and to verify closely how the energy exchange between the spectral modes takes place during propagation in this more general class of potentials with arbitrary gain and loss profiles.


title: Clustering Stable Instances of Euclidean k-means

abstract: The Euclidean k-means problem is arguably the most widely-studied clustering problem in machine learning. While the k-means objective is NP-hard in the worst-case, practitioners have enjoyed remarkable success in applying heuristics like Lloyd's algorithm for this problem. To address this disconnect, we study the following question: what properties of real-world instances will enable us to design efficient algorithms and prove guarantees for finding the optimal clustering? We consider a natural notion called additive perturbation stability that we believe captures many practical instances. Stable instances have unique optimal k-means solutions that do not change even when each point is perturbed a little (in Euclidean distance). This captures the property that the k-means optimal solution should be tolerant to measurement errors and uncertainty in the points. We design efficient algorithms that provably recover the optimal clustering for instances that are additive perturbation stable. When the instance has some additional separation, we show an efficient algorithm with provable guarantees that is also robust to outliers. We complement these results by studying the amount of stability in real datasets and demonstrating that our algorithm performs well on these benchmark datasets.


title: Elasto-plastic models of the yielding transition with stress-dependent transition rates

abstract: Elasto-plastic models are among the most successful ways to study the critical properties of the plastic yielding transition of amorphous solids. Typically these models are studied under a condition of constant transition rates from one plastic configuration to another, and in this form they predict the existence of well defined critical exponents that display universality, in the same sense that in standard equilibrium phase transitions. I show however that very naturally the transition rates must not be taken as a constant, but dependent of the local stress excess above the critical value. This modification in the model is seen to affect the values of some of the exponents of the transition, concretely, of the dynamical exponents that are related to the speed at which the system is driven. I argue about the reason for this dependence, claiming that it is due to the quasi-mean field nature of the plastic yielding transition originated in the fact that elastic interactions are long range.


title: A variation on bisecting the binomial coefficients

abstract: In this paper, we present an algorithm which allows us to search for all the bisections for the binomial coefficients $\{\binom{n}{k} \}_{k=0,...,n}$ and include a table with the results for all $n\le 154$. Connections with previous work on this topic is included. We conjecture that the probability of having only trivial solutions is $5/6$. \end{abstract}


title: Explaining the DAMPE data with scalar dark matter and gauged $U(1)_{L_e-L_\mu}$ interaction

abstract: Inspired by the peak structure observed by recent DAMPE experiment in $e^+e^-$ cosmic-ray spectrum, we consider a scalar dark matter (DM) model with gauged $U(1)_{L_e-L_\mu}$ symmetry, which is the most economical anomaly-free theory to potentially explain the peak by DM annihilation in nearby subhalo. We utilize the process $\chi \chi \to Z^\prime Z^\prime \to l \bar{l} l^\prime \bar{l}^\prime$, where $\chi$, $Z^\prime$, $l^{(\prime)}$ denote the scalar DM, the new gauge boson and $l^{(\prime)} =e, \mu$, respectively, to generate the $e^+e^-$ spectrum. By fitting the predicted spectrum to the experimental data, we obtain the favored DM mass range $m_\chi \simeq 3060^{+80}_{-100} \, {\rm GeV}$ and $\Delta m \equiv m_\chi - m_{Z^\prime} \lesssim 14 \, {\rm GeV}$ at $68\%$ Confidence Level (C.L.). Furthermore, we determine the parameter space of the model which can explain the peak and meanwhile satisfy the constraints from DM relic abundance, DM direct detection and the collider bounds. We conclude that the model we consider can account for the peak, although there exists a tension with the constraints from the LEP-II bound on $m_{Z^\prime}$ arising from the cross section measurement of $e^+e^- \to Z^{\prime\ast} \to e^+ e^-$.


title: Exponential Generalised Network Descriptors

abstract: In communication networks theory the concepts of networkness and network surplus have recently been defined. Together with transmission and betweenness centrality, they were based on the assumption of equal communication between vertices. Generalised versions of these four descriptors were presented, taking into account that communication between vertices $u$ and $v$ is decreasing as the distance between them is increasing. Therefore, we weight the quantity of communication by $\lambda^{d(u,v)}$ where $\lambda \in \left\langle0,1 \right\rangle$. Extremal values of these descriptors are analysed.


title: Reannihilation of self-interacting dark matter

abstract: We explore the phenomenology of having a second epoch of dark matter annihilation into dark radiation long after the standard thermal freeze-out. Such a hidden reannihilation process could affect visible sectors only gravitationally. As a concrete realization we consider self-interacting dark matter (SIDM) with a light force mediator coupled to dark radiation. We demonstrate how resonantly Sommerfeld enhanced cross sections emerge to induce the reannihilation epoch. The effect is a temporally local modification of the Hubble expansion rate and we show that the Cosmic Microwave Background (CMB) measurements -- as well as other observations -- have a high sensitivity to observe this phenomenon. Special attention is given to the model region where late kinetic decoupling and strong self-interactions can alleviate several small-scale problems in the cold dark matter paradigm at the same time. Interestingly, we find that reannihilation might here also simultaneously lower the tension between CMB and low-redshift astronomical observations of $H_{0}$ and $\sigma_{8}$. Moreover, we identify reannihilation as a clear signature to discriminate between the phenomenologically otherwise almost identical vector and scalar mediator realizations of SIDM.


title: Evidence for a shallow thin magnetic structure and solar dynamo: the driver of torsional oscillations

abstract: The solar dynamo and the solar Global internal Magnetic Structure (GMS) appear to be a thin ($\sim$2 Mm thick) structure near ($\sim$1 Mm below) the solar surface. Evidence for these properties are found from the amplitude of the torsional oscillations and in their velocity contours relationship to solar magnetogram; the power to the chromosphere; power to the corona and the solar wind; the current in the helio-current-sheet measured at the radius of the orbit of Earth; the calculated size ($\sim$1 Mm) of the expanding polar flux when it enters the photosphere; and from the observation that solar magnetic activity is generated near the surface. A thin stable minimum energy state seems to be covering most of the solar surface just below the photosphere. The magnetic field lines should be parallel to the solar surface and rotate with distance from the surface for 2$\pi$ radians in $\sim$2 Mm. Resistive diffusion helps to push the magnetic fields to the surface and the GMS seems to lose $\pi$ radians every 11 years, causing the observed 180$^\circ$ flipping of the solar magnetic fields including the flipping of the polar flux. Further evidence for this GMS and its loss is that solar prominences are made of thin sheets of magnetized plasma, which are, likely, remnants of the lost thin sheet of the GMS. The loss process is consistent with the butterfly pattern of the sunspots and with the differences observed between solar maximum and solar minimum in the corona. The solar dynamo drives current parallel to the polar flux, which, in turn, sustains the GMS using cross-field current drive. For completeness, the formation of sunspots, CMEs and flares is discussed.


title: Dynamo cycles in global convection simulations of solar-like stars

abstract: Several solar-like stars exhibit cyclic magnetic activity similar to the Sun as found in photospheric and chromospheric emission. We want to understand the physical mechanism involved in rotational dependence of these activity cycle periods. We use three-dimensional magnetohydrodynamical simulations of global convective dynamos models of solar-like stars to investigate the rotational dependency of dynamos. We further apply the test-field method to determine the $\alpha$ effect in these simulations. We find dynamo with clear oscillating mean magnetic fields for moderately and rapidly rotating runs. For slower rotation, the field is constant or exhibit irregular cycles. In the moderately and rapidly rotating regime the cycle periods increase weakly with rotation. This behavior can be well explained with a Parker-Yoshimura dynamo wave traveling equatorward. Even though the $\alpha$ effect becomes stronger for increasing rotation, the shear decreases steeper, causing this weak dependence on rotation. Similar as other numerical studies, we find no indication of activity branches as suggested by Brandenburg et al. (1998). However, our simulation seems to agree more with the transitional branch suggested by Distefano et al. (2017) and Olspert et al. (2017). If the Sun exhibit a dynamo wave similar as we find in our simulations, it would operate deep inside the convection zone.


title: On Out-of-Band Emissions of Quantized Precoding in Massive MU-MIMO-OFDM

abstract: We analyze out-of-band (OOB) emissions in the massive multi-user (MU) multiple-input multiple-output (MIMO) downlink. We focus on systems in which the base station (BS) is equipped with low-resolution digital-to-analog converters (DACs) and orthogonal frequency-division multiplexing (OFDM) is used to communicate to the user equipments (UEs) over frequency-selective channels. We demonstrate that analog filtering in combination with simple frequency-domain digital predistortion (DPD) at the BS enables a significant reduction of OOB emissions, but degrades the signal-to-interference-noise-and-distortion ratio (SINDR) at the UEs and increases the peak-to-average power ratio (PAR) at the BS. We use Bussgang's theorem to characterize the tradeoffs between OOB emissions, SINDR, and PAR, and to study the impact of analog filters and DPD on the error-rate performance of the massive MU-MIMO-OFDM downlink. Our results show that by carefully tuning the parameters of the analog filters, one can achieve a significant reduction in OOB emissions with only a moderate degradation of error-rate performance and PAR.


title: The algebraic geometry of Kazhdan-Lusztig-Stanley polynomials

abstract: Kazhdan-Lusztig-Stanley polynomials are a combinatorial generalization of Kazhdan-Lusztig polynomials of for Coxeter groups that include g-polynomials of polytopes and Kazhdan-Lusztig polynomials of matroids. In the cases of Weyl groups, rational polytopes, and realizable matroids, one can count points over finite fields on flag varieties, toric varieties, or reciprocal planes to obtain cohomological interpretations of these polynomials. We survey these results and unite them under a single geometric framework.


title: Critical $\mathrm{L}^p$-differentiability of $\mathrm{BV}^{\mathbb{A}}$-maps and canceling operators

abstract: We give a generalization of Dorronsoro's Theorem on critical $\mathrm{L}^p$-Taylor expansions for $\mathrm{BV}^k$-maps on $\mathbb{R}^n$, i.e., we characterize homogeneous linear differential operators $\mathbb{A}$ of $k$-th order such that $D^{k-j}u$ has $j$-th order $\mathrm{L}^{n/(n-j)}$-Taylor expansion a.e. for all $u\in\mathrm{BV}^\mathbb{A}_{\text{loc}}$ (here $j=1,\ldots, k$, with an appropriate convention if $j\geq n$). The space $\mathrm{BV}^\mathbb{A}_{\text{loc}}$ consists of those locally integrable maps $u$ such that $\mathbb{A} u$ is a Radon measure on $\mathbb{R}^n$. A new $\mathrm{L}^\infty$-Sobolev inequality is established to cover higher order expansions. Lorentz refinements are also considered. The main results can be seen as pointwise regularity statements for linear elliptic systems with measure-data.


title: An Equivalence of Fully Connected Layer and Convolutional Layer

abstract: This article demonstrates that convolutional operation can be converted to matrix multiplication, which has the same calculation way with fully connected layer. The article is helpful for the beginners of the neural network to understand how fully connected layer and the convolutional layer work in the backend. To be concise and to make the article more readable, we only consider the linear case. It can be extended to the non-linear case easily through plugging in a non-linear encapsulation to the values like this $\sigma(x)$ denoted as $x^{\prime}$.


title: Implementation of Multilayer Perceptron Network with Highly Uniform Passive Memristive Crossbar Circuits

abstract: The progress in the field of neural computation hinges on the use of hardware more efficient than the conventional microprocessors. Recent works have shown that mixed-signal integrated memristive circuits, especially their passive ('0T1R') variety, may increase the neuromorphic network performance dramatically, leaving far behind their digital counterparts. The major obstacle, however, is relatively immature memristor technology so that only limited functionality has been demonstrated to date. Here we experimentally demonstrate operation of one-hidden layer perceptron classifier entirely in the mixed-signal integrated hardware, comprised of two passive 20x20 metal-oxide memristive crossbar arrays, board-integrated with discrete CMOS components. The demonstrated multilayer perceptron network, whose complexity is almost 10x higher as compared to previously reported functional neuromorphic classifiers based on passive memristive circuits, achieves classification fidelity within 3 percent of that obtained in simulations, when using ex-situ training approach. The successful demonstration was facilitated by improvements in fabrication technology of memristors, specifically by lowering variations in their I-V characteristics.


title: Bounds on Neutrino Mass in Viscous Cosmology

abstract: Effective field theory of dark matter fluid on large scales predicts the presence of viscosity of the order of $10^{-6} H_0 M_P^2$. It has been shown that this magnitude of viscosities can resolve the discordance between large scale structure observations and Planck CMB data in the $\sigma_8$-$\Omega_m^0$ and $H_0$-$\Omega_m^0$ parameters space. Massive neutrinos suppresses the matter power spectrum on the small length scales similar to the viscosities. We show that by including the effective viscosity, which arises from summing over non linear perturbations at small length scales, severely constrains the cosmological bound on neutrino masses. Under a joint analysis of Planck CMB and different large scale observation data, we find that upper bound on the sum of the neutrino masses at 2-$\sigma$ level, decreases from $\sum m_\nu \le 0.396\,$eV (normal hierarchy) and $\sum m_\nu \le 0.378 \,$eV (inverted hierarchy) to $\sum m_\nu \le 0.267\,$eV (normal hierarchy) and $\sum m_\nu \le 0.146\,$eV (inverted hierarchy) when the effective viscosities are included.


title: Upper Tail Large Deviations in First Passage Percolation

abstract: For first passage percolation on $\mathbb{Z}^2$ with i.i.d. bounded edge weights, we consider the upper tail large deviation event; i.e., the rare situation where the first passage time between two points at distance $n$, is macroscopically larger than typical. It was shown by Kesten (1986) that the probability of this event decays as $\exp (-\Theta(n^2))$. However the question of existence of the rate function i.e., whether the log-probability normalized by $n^2$ tends to a limit, had remained open. We show that under some additional mild regularity assumption on the passage time distribution, the rate function for upper tail large deviation indeed exists. Our proof can be generalized to work in higher dimensions and for the corresponding problem in last passage percolation as well. The key intuition behind the proof is that a limiting metric structure which is atypical causes the upper tail large deviation event. The formal argument then relies on an approximate version of the above which allows us to dilate the large deviation environment to compare the upper tail probabilities for various values of $n.$

categories: math.PR cond-mat.stat-mech math-ph math.MP


title: IR properties of chiral effects in pionic matter

abstract: Chiral effects exhibit peculiar universality in idealized theoretical limits. However, they are known to be infrared sensitive and get modified in more realistic settings. In this work, we study how the corresponding conductivities vary with the constituent mass. We concentrate on a pionic realization of chiral effects which provides a better control over infrared properties of the theory. The pionic medium is considered at finite vector and axial isospin chemical potentials in the presence of an external magnetic field. This system supports electric and axial isospin currents along the magnetic field which correspond to chiral magnetic and chiral separation effects. We show that these currents are sensitive to the finite mass of the constituents, but the conductivities follow a simple scaling with the corresponding charge densities as one would expect for polarization effects. It is argued that this relation can capture the dependence of chiral effects on other infrared parameters. Finally, we briefly comment on the realization of the 't Hooft matching condition in pionic media at finite densities.


title: Development of Statewide AADT Estimation Model from Short-Term Counts: A Comparative Study for South Carolina

abstract: Annual Average Daily Traffic (AADT) is an important parameter used in traffic engineering analysis. Departments of Transportation (DOTs) continually collect traffic count using both permanent count stations (i.e., Automatic Traffic Recorders or ATRs) and temporary short-term count stations. In South Carolina, 87% of the ATRs are located on interstates and arterial highways. For most secondary highways (i.e., collectors and local roads), AADT is estimated based on short-term counts. This paper develops AADT estimation models for different roadway functional classes with two machine learning techniques: Artificial Neural Network (ANN) and Support Vector Regression (SVR). The models aim to predict AADT from short-term counts. The results are first compared against each other to identify the best model. Then, the results of the best model are compared against a regression method and factor-based method. The comparison reveals the superiority of SVR for AADT estimation for different roadway functional classes over all other methods. Among all developed models for different functional roadway classes, the SVR-based model shows a minimum root mean square error (RMSE) of 0.22 and a mean absolute percentage error (MAPE) of 11.3% for the interstate/expressway functional class. This model also shows a higher R-squared value compared to the traditional factor-based model and regression model. SVR models are validated for each roadway functional class using the 2016 ATR data and selected short-term count data collected by the South Carolina Department of Transportation (SCDOT). The validation results show that the SVR-based AADT estimation models can be used by the SCDOT as a reliable option to predict AADT from the short-term counts.


title: A pedagogical overview on 2D and 3D Toric Codes and the origin of their topological orders

abstract: In this work, we will show how the topological order of the Toric Code appears when the lattice on which it is defined discretizes a three-dimensional torus. In order to do this, we will present a pedagogical review of the traditional two-dimensional Toric Code, with an emphasis on how its quasiparticles are conceived and transported. With that, we want to make clear not only how all these same quasiparticle conception and transportation fit into this three-dimensional model, but to make it clear how topology controls the degeneracy of ground state in this new situation.


title: A Perceptual Measure for Deep Single Image Camera Calibration

abstract: Most current single image camera calibration methods rely on specific image features or user input, and cannot be applied to natural images captured in uncontrolled settings. We propose directly inferring camera calibration parameters from a single image using a deep convolutional neural network. This network is trained using automatically generated samples from a large-scale panorama dataset, and considerably outperforms other methods, including recent deep learning-based approaches, in terms of standard L2 error. However, we argue that in many cases it is more important to consider how humans perceive errors in camera estimation. To this end, we conduct a large-scale human perception study where we ask users to judge the realism of 3D objects composited with and without ground truth camera calibration. Based on this study, we develop a new perceptual measure for camera calibration, and demonstrate that our deep calibration network outperforms other methods on this measure. Finally, we demonstrate the use of our calibration network for a number of applications including virtual object insertion, image retrieval and compositing.


title: Continuous Matrix Product States for Inhomogeneous Quantum Field Theories: a Basis-Spline Approach

abstract: Continuous Matrix Product States (cMPS) are powerful variational ansatz states for ground states of continuous quantum field theories in (1+1) dimension. In this paper we introduce a novel parametrization of the cMPS wave function based on basis-spline functions, which we coin spline-based MPS (spMPS), and develop novel regauging techniques for inhomogeneous cMPS. We extend a recently developed ground-state optimization algorithm for translational invariant cMPS [M. Ganahl, J. Rinc\'on, G.Vidal. Phys.Rev.Lett. 118,220402 (2017)] to the case of inhomogeneous cMPS and, as proof-of-principle, use it to obtain the ground-state of a gas of Lieb-Liniger bosons in a periodic potential. The proposed method provides a first working implementation of a cMPS optimization for non-translational invariant continuous Hamiltonians.


title: SfSNet: Learning Shape, Reflectance and Illuminance of Faces in the Wild

abstract: We present SfSNet, an end-to-end learning framework for producing an accurate decomposition of an unconstrained human face image into shape, reflectance and illuminance. SfSNet is designed to reflect a physical lambertian rendering model. SfSNet learns from a mixture of labeled synthetic and unlabeled real world images. This allows the network to capture low frequency variations from synthetic and high frequency details from real images through the photometric reconstruction loss. SfSNet consists of a new decomposition architecture with residual blocks that learns a complete separation of albedo and normal. This is used along with the original image to predict lighting. SfSNet produces significantly better quantitative and qualitative results than state-of-the-art methods for inverse rendering and independent normal and illumination estimation.


title: Compatibility Family Learning for Item Recommendation and Generation

abstract: Compatibility between items, such as clothes and shoes, is a major factor among customer's purchasing decisions. However, learning "compatibility" is challenging due to (1) broader notions of compatibility than those of similarity, (2) the asymmetric nature of compatibility, and (3) only a small set of compatible and incompatible items are observed. We propose an end-to-end trainable system to embed each item into a latent vector and project a query item into K compatible prototypes in the same space. These prototypes reflect the broad notions of compatibility. We refer to both the embedding and prototypes as "Compatibility Family". In our learned space, we introduce a novel Projected Compatibility Distance (PCD) function which is differentiable and ensures diversity by aiming for at least one prototype to be close to a compatible item, whereas none of the prototypes are close to an incompatible item. We evaluate our system on a toy dataset, two Amazon product datasets, and Polyvore outfit dataset. Our method consistently achieves state-of-the-art performance. Finally, we show that we can visualize the candidate compatible prototypes using a Metric-regularized Conditional Generative Adversarial Network (MrCGAN), where the input is a projected prototype and the output is a generated image of a compatible item. We ask human evaluators to judge the relative compatibility between our generated images and images generated by CGANs conditioned directly on query items. Our generated images are significantly preferred, with roughly twice the number of votes as others.


title: Data-Driven Spatio-Temporal Analysis of Curbside Parking Demand: A Case-Study in Seattle

abstract: Due to rapid expansion of urban areas in recent years, management of curbside parking has become increasingly important. To mitigate congestion, while meeting a city's diverse needs, performance-based pricing schemes have received a significant amount of attention. However, several recent studies suggest location, time-of-day, and awareness of policies are the primary factors that drive parking decisions. In light of this, we provide an extensive data-driven study of the spatio-temporal characteristics of curbside parking. This work advances the understanding of where and when to set pricing policies, as well as where to target information and incentives to drivers looking to park. Harnessing data provided by the Seattle Department of Transportation, we develop a Gaussian mixture model based technique to identify zones with similar spatial parking demand as quantified by spatial autocorrelation. In support of this technique, we introduce a metric based on the repeatability of our Gaussian mixture model to investigate temporal consistency.


title: A Context-aware Recommender System for Hyperlocal News: A Conceptual Framework

abstract: Recommender systems (RSs) have been popular in variety of application domains due to the increased demand for filtering and sorting items and information. Today, there is a numerous approaches and algorithms of data filtering and recommendations. This works presents a conceptual framework for constructing a mobile RS in hyper-local news domain. The mobile RS is designed to deal with specific requirements of news readers, such as spatial- temporal relevance, recency, real-time update and validated news. The implementation of the RS in a distributed file system is also discussed.


title: A Multiple Observer Probability Analysis for Bell Scenarios in Special Relativity

abstract: Here we present a Multiple Observer Probability Analysis (MOPA) for the purpose of clarifying topics in experimental Bell scenarios. Because Bell scenarios are interested in quantum effects between nonlocal measurement devices, we assign an observer to each device: Alice and Bob. Given that the observers are stationary and space-like separated, each observer is privy to different information along their shared equi-temporal lines due to permutations in the order they observe events. Therefore, each observer is inclined to assign different probability distributions to the same set of propositions due to these informational differences. The observers are obligated to update their probability distributions on the basis of locally observed events, and in this sense, factuality is informational locality. In this framework, only local variables or detections may be factual, but nothing prevents an observer from inquiring or making if-then inferences on the counterfactual basis of a nonlocal proposition being true. Indeed the objects pertaining to these nonlocal counterfactual propositions may be far outside an observer's light cone. The MOPA arrives at the conclusion that the CHSH inequality is only nonlocally violated "counterfactually" by each observer whereas local violations of the CHSH may be factual or counterfactual. We believe the MOPA to better gel probability theory (and thus QM) with Special Relativity than does the standard locality conditions imposed in the Bell and CHSH inequalities. The no-signaling condition is reinterpreted, and perhaps further clarified, in the MOPA and statements about counterfactuality and observer dependent QM are made.


title: On scattering for the defocusing nonlinear Schr\"odinger equation on waveguide $ \mathbb{R}^m$ $\times$ $\mathbb{T}$ (when $m=2,3$)

abstract: In the article, we prove the large data scattering for two problems, i.e. the defocusing quintic nonlinear Schr{\"o}dinger equation on $\mathbb{R}^2$ $\times$ $\mathbb{T}$ and the defocusing cubic nonlinear Schr{\"o}dinger equation on $\mathbb{R}^3$ $\times$ $\mathbb{T}$. Both of the two equations are mass supercritical and energy critical. The main ingredients of the proofs contain global Stricharz estimate, profile decomposition and energy induction method. This paper is the second project of our series work (two papers, together with [36]) on large data scattering for the defocusing critical NLS with integer index nonlinearity on low dimensional waveguides. At this point, that type of problems are almost solved except for two remaining resonant system conjectures and the quintic NLS problem on $\mathbb{R}\times \mathbb{T}$.


title: Comment on "Unified view of quantum correlations and quantum coherence"

abstract: We show that contrary to what it is claimed in Phys. Rev. A \textbf{94}, 022329 (2016), in general the local projective measurement that induces maximal coherence loss is \textit{not} the projection onto the local basis that defines the coherence of the system, at least for all quantum-incoherent states.


title: An interior spherical static solution of Einstein-Maxwell equations

abstract: The exact interior static solution of Einstein-Maxwell equations in Bondi's coordinates for the electrically charged fluid ball is obtained as a source of the Reissner-Nordstr\"om solution.


title: Two-weight, weak type norm inequalities for a class of sublinear operators on weighted Morrey and amalgam spaces

abstract: Let $\mathcal T_\alpha~(0\leq\alpha<n)$ be a class of sublinear operators satisfying certain size conditions introduced by Soria and Weiss, and let $[b,\mathcal T_\alpha]~(0\leq\alpha<n)$ be the commutators generated by $\mathrm{BMO}(\mathbb R^n)$ functions and $\mathcal T_\alpha$. This paper is concerned with two-weight, weak type norm estimates for these sublinear operators and their commutators on the weighted Morrey and amalgam spaces. Some boundedness criterions for such operators are given, under the assumptions that weak-type norm inequalities on weighted Lebesgue spaces are satisfied. As applications of our main results, we can obtain the weak-type norm inequalities for several integral operators as well as the corresponding commutators in the framework of weighted Morrey and amalgam spaces.


title: Finite-time and exact Lyapunov dimension of the Henon map

abstract: This work is devoted to further consideration of the Henon map with negative values of the shrinking parameter and the study of transient oscillations, multistability, and possible existence of hidden attractors. The computation of the finite-time Lyapunov exponents by different algorithms is discussed. A new adaptive algorithm for the finite-time Lyapunov dimension computation in studying the dynamics of dimension is used. Analytical estimates of the Lyapunov dimension using the localization of attractors are given. A proof of the conjecture on the Lyapunov dimension of self-excited attractors and derivation of the exact Lyapunov dimension formula are revisited.


title: On the 2-part of the Birch and Swinnerton-Dyer conjecture for quadratic twists of elliptic curves

abstract: In the present paper, we prove, for a large class of elliptic curves defined over $\mathbb{Q}$, the existence of an explicit infinite family of quadratic twists with analytic rank $0$. In addition, we establish the $2$-part of the conjecture of Birch and Swinnerton-Dyer for many of these infinite families of quadratic twists. Recently, Xin Wan has used our results to prove for the first time the full Birch--Swinnerton-Dyer conjecture for some explicit infinite families of elliptic curves defined over $\mathbb{Q}$ without complex multiplication.


title: Layer-wise Learning of Stochastic Neural Networks with Information Bottleneck

abstract: Information Bottleneck (IB) is a generalization of rate-distortion theory that naturally incorporates compression and relevance trade-offs for learning. Though the original IB has been extensively studied, there has not been much understanding of multiple bottlenecks which better fit in the context of neural networks. In this work, we propose Information Multi-Bottlenecks (IMBs) as an extension of IB to multiple bottlenecks which has a direct application to training neural networks by considering layers as multiple bottlenecks and weights as parameterized encoders and decoders. We show that the multiple optimality of IMB is not simultaneously achievable for stochastic encoders. We thus propose a simple compromised scheme of IMB which in turn generalizes maximum likelihood estimate (MLE) principle in the context of stochastic neural networks. We demonstrate the effectiveness of IMB on classification tasks and adversarial robustness in MNIST and CIFAR10.


title: Phase sensitive imaging of 10 GHz vibrations in an AlN microdisk resonator

abstract: We demonstrate a high frequency phase-sensitive heterodyne vibrometer, operating up to 10 GHz. Using this heterodyne vibrometer, the amplitude and phase fields of the fundamental thickness mode, the radial fundamental and the 2nd-order modes of an AlN optomechanical microdisk resonator are mapped with a displacement sensitivity of around 0.36 picometer at 1 Hz resolution bandwidth. The simultaneous amplitude and phase measurement allows precise mode identification and characterization. The recorded modal frequencies and proles are consistent with numerical simulations. This vibrometer will be of great significance for the development of high frequency mechanical devices.


title: A Type D Non-Vacuum Spacetime with Causality Violating Curves, and its Physical Interpretation

abstract: We present a topologically trivial, non-vacuum solution of the Einstein's field equations in four-dimensions, which is regular everywhere. The metric admits circular closed timelike curves, which appear beyond the null curve, and these timelike curves are linearly stable under linear perturbations. Additionally, the spacetime admits null geodesics curve which are not closed, and the metric is of type D in the Petrov classification scheme. The stress-energy tensor anisotropic fluid satisfy the different energy conditions and a generalization of Equation-of-State parameter of perfect fluid $p=\omega\,\rho$. The metric admits a twisting, shearfree, non-exapnding timelike geodesic congruence. Finally, the physical interpretation of this solution, based on the study of the equation of the geodesics deviation, will be presented.


title: A Deeper Look at Experience Replay

abstract: Recently experience replay is widely used in various deep reinforcement learning (RL) algorithms, in this paper we rethink the utility of experience replay. It introduces a new hyper-parameter, the memory buffer size, which needs carefully tuning. However unfortunately the importance of this new hyper-parameter has been underestimated in the community for a long time. In this paper we did a systematic empirical study of experience replay under various function representations. We showcase that a large replay buffer can significantly hurt the performance. Moreover, we propose a simple O(1) method to remedy the negative influence of a large replay buffer. We showcase its utility in both simple grid world and challenging domains like Atari games.


title: The case for preserving our knowledge and data in physics experiments

abstract: This proceeding covers tools and technologies at our disposal for scientific data preservation and shows that this extends the scientific reach of our experiments. It is cost-efficient to warehouse data from completed experiments on the tape archives of our national and international laboratories. These subject-specific data stores also offer the technologies to capture and archive knowledge about experiments in the form of technical notes, electronic logs, websites, etc. Furthermore, it is possible to archive our source code and computing environments. The paper illustrates these challenges with experience from preserving the LEP data for the long term.


title: A coherence study on EEG and EMG signals

abstract: The aim of this study is to investigate bursts- related EEG signals in a focal hand dystonia patient. Despite of considering time domain and frequency domain techniques as mutually exclusive analysis, in this contribution we have taken advantage from both of them: particularly, in the frequency domain, coherence was used to identify the most likely frequency bands of interaction between brain and muscles, then, in the time domain, cross-correlation was exploited to verify the physiological reliability of such a relationship in terms of signal transmission delay from the centre to the periphery. Our preliminary results suggest - in line with recent literature - that activity in the high beta band (around 30 Hz) could represent an electroencephalographic correlate for the pathological electromyographic bursts affecting the focal hand dystonia condition. Even though a future study on a larger sample is needed to statistically support these preliminary findings, this contribution allows to think of new kinds of rehabilitation from focal hand dystonia that could target the actual electroencephalographic correlate of the pathology, i.e. phenotypically expressed by bursts, with the consequence of a relevant functional improvement.


title: Noise characterization for resonantly-enhanced polarimetric vacuum magnetic-birefringence experiments

abstract: In this work we present data characterizing the sensitivity of the Bir\'{e}fringence Magnetique du Vide (BMV) instrument. BMV is an experiment attempting to measure vacuum magnetic birefringence (VMB) via the measurement of an ellipticity induced in a linearly polarized laser field propagating through a birefringent region of vacuum in the presence of an external magnetic field. Correlated measurements of laser noise alongside the measurement in the main detection channel allow us to separate measured sensing noise from the inherent birefringence noise of the apparatus. To this end we model different sources of sensing noise for cavity-enhanced polarimetry experiments, such as BMV. Our goal is to determine the main sources of noise, clarifying the limiting factors of such an apparatus. We find our noise models are compatible with the measured sensitivity of BMV. In this context we compare the phase sensitivity of separate-arm interferometers to that of a polarimetry apparatus for the discussion of current and future VMB measurements.


title: ${SO(10)}$ inspired $Z'$ models at the LHC

abstract: We study and compare various $Z'$ models arising from $SO(10)$, focussing in particular on the Abelian subgroup $U(1)_{R} \times U(1)_{B-L}$, broken at the TeV scale to Standard Model hypercharge $U(1)_{Y}$. The gauge group $U(1)_{R} \times U(1)_{B-L}$, which is equivalent to the $U(1)_{Y}\times U(1)_{\chi}$ in a different basis, is well motivated from $SO(10)$ breaking and allows neutrino mass via the linear seesaw mechanism. Assuming supersymmetry, we consider single step gauge unification to predict the gauge couplings, then consider the detection and characterisation prospects of the resulting $Z'$ at the LHC by studying its possible decay modes into di-leptons as well as into Higgs bosons. The main new result here is to analyse in detail the expected leptonic forward-backward asymmetry at the high luminosity LHC and show that it may be used to discriminate the $U(1)_{R} \times U(1)_{B-L}$ model from the usual $B-L$ model based on $U(1)_{Y}\times U(1)_{B-L}$.


title: The Dispersion of Fast Radio Bursts from a Structured Intergalactic Medium at Redshifts z < 1.5

abstract: We analyze the sources of free electrons that produce the large dispersion measures, DM $\approx 300-1600$ (in units cm$^{-3}$ pc), observed toward fast radio bursts (FRBs). Individual galaxies typically produce DM $\sim 25-60$ cm$^{-3}$ pc from ionized gas in their disk, disk-halo interface, and circumgalactic medium. Toward an FRB source at redshift $z$, a homogeneous IGM containing a fraction $f_{\rm IGM}$ of cosmological baryons will produce DM $= (935~{\rm cm}^{-3}~{\rm pc}) f_{\rm IGM} \, h_{70}^{-1} I(z)$, where $I(z) = (2/3 \Omega_m)[ \{ \Omega_m(1+z)^3 + \Omega_{\Lambda} \}^{1/2} - 1 ]$. A structured IGM of photoionized Ly-alpha absorbers in the cosmic web produces similar dispersion, modeled from the observed distribution, $f_b(N,z)$, of H I (Lya-forest) absorbers in column density and redshift with ionization corrections and scaling relations from cosmological simulations. An analytic formula for DM($z$) applied to observed FRB dispersions suggests that $z_{\rm FRB} \approx 0.2-1.5$ for an IGM containing a significant baryon fraction, $f_{\rm IGM} = 0.6\pm0.1$. Future surveys of the statistical distribution, DM($z)$, of FRBs identified with specific galaxies and redshifts, can be used to calibrate the IGM baryon fraction and distribution of Ly-alpha absorbers. Fluctuations in DM at the level $\pm10$ cm$^{-3}$ pc will arise from filaments and voids in the cosmic web.


title: Studying tidal effects in planetary systems with Posidonius. A N-body simulator written in Rust

abstract: Planetary systems with several planets in compact orbital configurations such as TRAPPIST-1 are surely affected by tidal effects. Its study provides us with important insight about its evolution. We developed a second generation of a N-body code based on the tidal model used in Mercury-T, re-implementing and improving its functionalities using Rust as programming language (including a Python interface for easy use) and the WHFAST integrator. The new open source code ensures memory safety, reproducibility of numerical N-body experiments, it improves the spin integration compared to Mercury-T and allows to take into account a new prescription for the dissipation of tidal inertial waves in the convective envelope of stars. Posidonius is also suitable for binary system simulations with evolving stars.

categories: astro-ph.EP astro-ph.IM astro-ph.SR cs.PL


title: Hessian eigenvalue distribution in a random Gaussian landscape

abstract: The energy landscape of multiverse cosmology is often modeled by a multi-dimensional random Gaussian potential. The physical predictions of such models crucially depend on the eigenvalue distribution of the Hessian matrix at potential minima. In particular, the stability of vacua and the dynamics of slow-roll inflation are sensitive to the magnitude of the smallest eigenvalues. The Hessian eigenvalue distribution has been studied earlier, using the saddle point approximation, in the leading order of $1/N$ expansion, where $N$ is the dimensionality of the landscape. This approximation, however, is insufficient for the small eigenvalue end of the spectrum, where sub-leading terms play a significant role. We extend the saddle point method to account for the sub-leading contributions. We also develop a new approach, where the eigenvalue distribution is found as an equilibrium distribution at the endpoint of a stochastic process (Dyson Brownian motion). The results of the two approaches are consistent in cases where both methods are applicable. We discuss the implications of our results for vacuum stability and slow-roll inflation in the landscape.

categories: hep-th cond-mat.stat-mech math-ph math.MP


title: Spatially resolved kinematics in the central 1 kpc of a compact star-forming galaxy at z=2.3 from ALMA CO observations

abstract: We present high spatial resolution (FWHM$\sim$0.14'') observations of the CO($8-7$) line in GDS-14876, a compact star-forming galaxy at $z=2.3$ with total stellar mass of $\log(M_{\star}/M_{\odot})=10.9$. The spatially resolved velocity map of the inner $r\lesssim1$~kpc reveals a continous velocity gradient consistent with the kinematics of a rotating disk with $v_{\rm rot}(r=1\rm kpc)=163\pm5$ km s$^{-1}$ and $v_{\rm rot}/\sigma\sim2.5$. The gas-to-stellar ratios estimated from CO($8-7$) and the dust continuum emission span a broad range, $f^{\rm CO}_{\rm gas}=M_{\rm gas}/M_{\star}=13-45\%$ and $f^{\rm cont}_{\rm gas}=50-67\%$, but are nonetheless consistent given the uncertainties in the conversion factors. The dynamical modeling yields a dynamical mass of$\log(M_{\rm dyn}/M_{\odot})=10.58^{+0.5}_{-0.2}$ which is lower, but still consistent with the baryonic mass, $\log$(M$_{\rm bar}$= M$_{\star}$ + M$^{\rm CO}_{\rm gas}$/M$_{\odot}$)$=11.0$, if the smallest CO-based gas fraction is assumed. Despite a low, overall gas fraction, the small physical extent of the dense, star-forming gas probed by CO($8-7$), $\sim3\times$ smaller than the stellar size, implies a strong concentration that increases the gas fraction up to $f^{\rm CO, 1\rm kpc}_{\rm gas}\sim 85\%$ in the central 1 kpc. Such a gas-rich center, coupled with a high star-formation rate, SFR$\sim$ 500 M$_{\odot}$ yr$^{-1}$, suggests that GDS-14876 is quickly assembling a dense stellar component (bulge) in a strong nuclear starburst. Assuming its gas reservoir is depleted without replenishment, GDS-14876 will quickly ($t_{\rm depl}\sim27$ Myr) become a compact quiescent galaxy that could retain some fraction of the observed rotational support.


title: M13 multiple stellar populations seen with the eyes of Str\"omgren photometry

abstract: We present a photometric study of M13 multiple stellar populations over a wide field of view, covering approximately 6.5 half-light radii, using archival Isaac Newton Telescope observations to build an accurate multi-band Str\"omgren catalogue. The use of the Str\"omgren index $c_{y}$ permits us to separate the multiple populations of M13 on the basis of their position on the red giant branch. The comparison with medium and high resolution spectroscopic analysis confirms the robustness of our selection criterion. To determine the radial distribution of stars in M13, we complemented our dataset with Hubble Space Telescope observations of the cluster core, to compensate for the effect of incompleteness affecting the most crowded regions. From the analysis of the radial distributions we do not find any significant evidence of spatial segregation. Some residuals may be present in the external regions where we observe only a small number of stars. This finding is compatible with the short dynamical timescale of M13 and represents, to date, one of the few examples of fully spatially mixed multiple populations in a massive globular cluster.


title: Spitzer Secondary Eclipses of Qatar-1b

abstract: Previous secondary eclipse observations of the hot Jupiter Qatar-1b in the Ks band suggest that it may have an unusually high day side temperature, indicative of minimal heat redistribution. There have also been indications that the orbit may be slightly eccentric, possibly forced by another planet in the system. We investigate the day side temperature and orbital eccentricity using secondary eclipse observations with Spitzer. We observed the secondary eclipse with Spitzer/IRAC in subarray mode, in both 3.6 and 4.5 micron wavelengths. We used pixel-level decorrelation to correct for Spitzer's intra-pixel sensitivity variations and thereby obtain accurate eclipse depths and central phases. Our 3.6 micron eclipse depth is 0.149 +/- 0.051% and the 4.5 micron depth is 0.273 +/- 0.049%. Fitting a blackbody planet to our data and two recent Ks band eclipse depths indicates a brightness temperature of 1506 +/- 71K. Comparison to model atmospheres for the planet indicates that its degree of longitudinal heat redistribution is intermediate between fully uniform and day side only. The day side temperature of the planet is unlikely to be as high (1885K) as indicated by the ground-based eclipses in the Ks band, unless the planet's emergent spectrum deviates strongly from model atmosphere predictions. The average central phase for our Spitzer eclipses is 0.4984 +/- 0.0017, yielding e cos(omega) = -0.0028 +/- 0.0027. Our results are consistent with a circular orbit, and we constrain e cos(omega) much more strongly than has been possible with previous observations.


title: Multiple Stellar Populations in Globular Clusters

abstract: Globular Clusters (GCs) exhibit star-to-star variations in specific elements (e.g., He, C, N, O, Na, Al) that bare the hallmark of high temperature H burning. These abundance variations can be observed spectroscopically and also photometrically, with the appropriate choice of filters, due to the changing of spectral features within the band pass. This phenomenon is observed in nearly all of the ancient GCs, although, to date, has not been found in any massive cluster younger than 2~Gyr. Many scenarios have been suggested to explain this phenomenon, with most invoking multiple epochs of star-formation within the cluster, however all have failed to reproduce various key observations, in particular when a global view of the GC population is taken. We review the state of current observations, and outline the successes and failures of each of the main proposed models. The traditional idea of using the stellar ejecta from a 1st generation of stars to form a 2nd generation of stars, while conceptually straight forward, has failed to reproduce an increasing number of observational constraints. We conclude that the puzzle of multiple populations remains unsolved, hence alternative theories are needed.


title: The Near-Infrared CO Absorption Band as a Probe to the Innermost Part of an AGN Obscuring Material

abstract: We performed a systematic analysis of the 4.67 $\mu$m CO ro-vibrational absorption band toward nearby active galactic nuclei (AGNs) and analyzed the absorption profiles of ten nearby galaxies collected from the AKARI and Spitzer spectroscopic observations that show the CO absorption feature by fitting a plane-parallel local thermal equilibrium gas model. We found that CO gas is warm (200--500 K) and has a large column density ($N_\mathrm{H}\gtrsim10^{23}~\mathrm{cm^{-2}}$). The heating of the gas is not explicable by either UV heating or shock heating because these processes cannot represent the large column densities of the warm gas. Instead, X-ray photons from the nuclei, which can produce large columns of warm gas with up to $N_\mathrm{H}\sim10^{24}~\mathrm{cm^{-2}}$, are the most convincing power source. The hydrogen column density estimated from the CO band is smaller than that inferred from X-ray observations. These results indicate that the region probed by the near-infrared CO absorption is in the vicinity of the nuclei and is located outside the X-ray emitting region. Furthermore, the covering factors close to unity required by the observed deep absorption profiles suggest that the probed region is close to the continuum source, which can be designated as the inner rim of the obscuring material around the AGN.


title: Instanton counting in Class $\mathcal{S}_k$

abstract: We compute the instanton partition functions of $\mathcal{N}=1$ SCFTs in class $\mathcal{S}_k$. We obtain this result via orbifolding Dp/D(p-4) brane systems and calculating the partition function of the supersymmetric gauge theory on the worldvolume of $K$ D(p-4) branes. Starting with D5/D1 setups probing a $\mathbb{Z}_\ell\times \mathbb{Z}_k$ orbifold singularity we obtain the $K$ instanton partition functions of 6d $(1,0)$ theories on $\mathbb{R}^4 \times T^2$ in the presence of orbifold defects on $T^2$ via computing the 2d superconformal index of the worldvolume theory on $K$ D1 branes wrapping the $T^2$. We then reduce our results to the 5d and to the 4d instanton partition functions. For $k=1$ we check that we reproduce the known elliptic, trigonometric and rational Nekrasov partition functions. Finally, we show that the instanton partition functions of $SU(N)$ quivers in class $\mathcal{S}_k$ can be obtained from the class $\mathcal{S}$ mother theory partition functions with $SU(kN)$ gauge factors via imposing the `orbifold condition' $a_{\mathcal{A}} \rightarrow a_A e^{2\pi i j/k}$ with $\mathcal{A}=jA$ and $A=1,\dots, N$, $j=1,\dots, k$ on the Coulomb moduli and the mass parameters.


title: Consistent hydrodynamic theory of chiral electrons in Weyl semimetals

abstract: The complete set of Maxwell's and hydrodynamic equations for the chiral electrons in Weyl semimetals is presented. The formulation of the Euler equation takes into account the explicit breaking of the Galilean invariance by the ion lattice. It is shown that the Chern-Simons (or Bardeen-Zumino) contributions should be added to the electric current and charge densities in Maxwell's equations that provide the information on the separation of Weyl nodes in energy and momentum. On the other hand, these topological contributions do not directly affect the Euler equation and the energy conservation relation for the electron fluid. By making use of the proposed consistent hydrodynamic framework, we show that the Chern-Simons contributions strongly modify the dispersion relations of collective modes in Weyl semimetals. This is reflected, in particular, in the existence of distinctive anomalous Hall waves, which are sustained by the local anomalous Hall currents.


title: Spitzer Observations of the North Ecliptic Pole

abstract: We present a photometric catalog for Spitzer Space Telescope warm mission observations of the North Ecliptic Pole (NEP; centered at $\rm R.A.=18^h00^m00^s$, $\rm Decl.=66^d33^m38^s.552$). The observations are conducted with IRAC in 3.6 $\mu$m and 4.5 $\mu$m bands over an area of 7.04 deg$^2$ reaching 1$\sigma$ depths of 1.29 $\mu$Jy and 0.79 $\mu$Jy in the 3.6 $\mu$m and 4.5 $\mu$m bands respectively. The photometric catalog contains 380,858 sources with 3.6 $\mu$m and 4.5 $\mu$m band photometry over the full-depth NEP mosaic. Point source completeness simulations show that the catalog is 80% complete down to 19.7 AB. The accompanying catalog can be utilized in constraining the physical properties of extra-galactic objects, studying the AGN population, measuring the infrared colors of stellar objects, and studying the extra-galactic infrared background light.


title: Machine learning for predictive estimation of qubit dynamics subject to dephasing

abstract: Decoherence remains a major challenge in quantum computing hardware and a variety of physical-layer controls provide opportunities to mitigate the impact of this phenomenon through feedback and feedforward control. In this work, we compare a variety of machine learning algorithms derived from diverse fields for the task of state estimation (retrodiction) and forward prediction of future qubit state evolution for a single qubit subject to classical, non-Markovian dephasing. Our approaches involve the construction of a dynamical model capturing qubit dynamics via autoregressive or Fourier-type protocols using only a historical record of projective measurements. A detailed comparison of achievable prediction horizons, model robustness, and measurement-noise-filtering capabilities for Kalman Filters (KF) and Gaussian Process Regression (GPR) algorithms is provided. We demonstrate superior performance from the autoregressive KF relative to Fourier-based KF approaches and focus on the role of filter optimization in achieving suitable performance. Finally, we examine several realizations of GPR using different kernels and discover that these approaches are generally not suitable for forward prediction. We highlight the underlying failure mechanism in this application and identify ways in which the output of the algorithm may be misidentified numerical artefacts.


title: The VANDELS survey: Dust attenuation in star-forming galaxies at $\mathbf{z=3-4}$

abstract: We present the results of a new study of dust attenuation at redshifts $3 < z < 4$ based on a sample of $236$ star-forming galaxies from the VANDELS spectroscopic survey. Motivated by results from the First Billion Years (FiBY) simulation project, we argue that the intrinsic spectral energy distributions (SEDs) of star-forming galaxies at these redshifts have a self-similar shape across the mass range $8.2 \leq$ log$(M_{\star}/M_{\odot}) \leq 10.6$ probed by our sample. Using FiBY data, we construct a set of intrinsic SED templates which incorporate both detailed star formation and chemical abundance histories, and a variety of stellar population synthesis (SPS) model assumptions. With this set of intrinsic SEDs, we present a novel approach for directly recovering the shape and normalization of the dust attenuation curve. We find, across all of the intrinsic templates considered, that the average attenuation curve for star-forming galaxies at $z\simeq3.5$ is similar in shape to the commonly-adopted Calzetti starburst law, with an average total-to-selective attenuation ratio of $R_{V}=4.18\pm0.29$. We show that the optical attenuation ($A_V$) versus stellar mass ($M_{\star}$) relation predicted using our method is consistent with recent ALMA observations of galaxies at $2<z<3$ in the \emph{Hubble} \emph{Ultra} \emph{Deep} \emph{Field} (HUDF), as well as empirical $A_V - M_{\star}$ relations predicted by a Calzetti-like law. Our results, combined with other literature data, suggest that the $A_V - M_{\star}$ relation does not evolve over the redshift range $0<z<5$, at least for galaxies with log$(M_{\star}/M_{\odot}) \gtrsim 9.5$. Finally, we present tentative evidence which suggests that the attenuation curve may become steeper at log$(M_{\star}/M_{\odot}) \lesssim 9.0$.


title: Probabilistic treatment of the uncertainty from the finite size of weighted Monte Carlo data

abstract: Parameter estimation in HEP experiments often involves Monte-Carlo simulation to model the experimental response function. A typical application are forward-folding likelihood analyses with re-weighting, or time-consuming minimization schemes with a new simulation set for each parameter value. Problematically, the finite size of such Monte Carlo samples carries intrinsic uncertainty that can lead to a substantial bias in parameter estimation if it is neglected and the sample size is small. We introduce a probabilistic treatment of this problem by replacing the usual likelihood functions with novel generalized probability distributions that incorporate the finite statistics via suitable marginalization. These new PDFs are analytic, and can be used to replace the Poisson, multinomial, and sample-based unbinned likelihoods, which covers many use cases in high-energy physics. In the limit of infinite statistics, they reduce to the respective standard probability distributions. In the general case of arbitrary Monte Carlo weights, the expressions involve the fourth Lauricella function $F_D$, for which we find a new finite-sum representation in a certain parameter setting. The result also represents an exact form for Carlson's Dirichlet average $R_n$ with $n>0$, and thereby an efficient way to calculate the probability generating function of the Dirichlet-multinomial distribution, the extended divided difference of a monomial, or arbitrary moments of univariate B-splines. We demonstrate the bias reduction of our approach with a typical toy Monte Carlo problem, estimating the normalization of a peak in a falling energy spectrum, and compare the results with previously published methods from the literature.

categories: physics.data-an astro-ph.IM hep-ex math.ST stat.TH


title: A New Public Release of the GIZMO Code

abstract: We describe a major update to the public GIZMO code. GIZMO has been used in simulations of cosmology; galaxy and star formation and evolution; black hole accretion and feedback; proto-stellar disk dynamics and planet formation; fluid dynamics and plasma physics; dust-gas dynamics; giant impacts and solid-body interactions; collisionless gravitational dynamics; and more. This release of the public code supports: hydrodynamics (using various mesh-free finite-volume Godunov methods or SPH), ideal and non-ideal MHD, anisotropic conduction and viscosity, radiative cooling and chemistry, star and black hole formation and feedback, sink particles, dust-gas (aero)-dynamics (with or without magnetic fields), elastic/plastic dynamics, arbitrary (gas, stellar, degenerate, solid/liquid material) equations of state, passive scalar/turbulent diffusion, large-eddy and shearing boxes, self-gravity with fully-adaptive force softenings, arbitrary cosmological expansion, and on-the-fly group-finding. It is massively-parallel with hybrid MPI+OpenMP scaling verified up to >1 million threads. The code is extensively documented, with test problems and tutorials provided for these different physics modules.

categories: astro-ph.IM astro-ph.CO astro-ph.EP astro-ph.GA physics.flu-dyn


title: Leptoquarks meet $\varepsilon'/\varepsilon$ and rare Kaon processes

abstract: We analyse for the first time the CP violating ratio $\varepsilon'/\varepsilon$ in $K\to \pi\pi$ decays in leptoquark (LQ) models. Assuming a mass gap to the electroweak (EW) scale, the main mechanism for LQs to contribute to $\varepsilon'/\varepsilon$ is EW gauge-mixing of semi-leptonic into non-leptonic operators, which we treat in the Standard Model effective theory (SMEFT). We perform also the one-loop decoupling for scalar LQs, finding that in all models with both left-handed and right-handed LQ couplings box-diagrams generate numerically strongly enhanced EW-penguin operators $Q_{8,8'}$ already at the LQ scale. We then investigate correlations of $\varepsilon'/\varepsilon$ with rare Kaon processes ($K_L\to\pi^0\nu\bar\nu$, $K^+\to\pi^+\nu\bar\nu$, $K_L\to\pi^0\ell\bar\ell$, $K_S\to\mu\bar\mu$, $\Delta M_K$ and $\varepsilon_K$) and find that even imposing only a moderate enhancement of $(\varepsilon'/\varepsilon)_{\rm NP} = 5 \times 10^{-4}$ to explain the current anomaly hinted by the Dual QCD approach and RBC-UKQCD lattice QCD calculations leads to conflicts with experimental upper bounds on rare Kaon processes. They exclude all LQ models with only a single coupling as an explanation of the $\varepsilon'/\varepsilon$ anomaly and put strong-to-serious constraints on parameter spaces of the remaining models. Future results on $K^+\to\pi^+\nu\bar\nu$ from the NA62 collaboration, $K_L\to\pi^0\nu\bar\nu$ from the KOTO experiment and $K_S\to\mu\bar\mu$ from LHCb will even stronger exhibit the difficulty of LQ models in explaining the measured $\varepsilon'/\varepsilon$, in case the $\varepsilon'/\varepsilon$ anomaly will be confirmed by improved lattice QCD calculations. Hopefully also improved measurements of $K_L\to\pi^0\ell\bar\ell$ decays will one day help in this context.


title: Neutral hydrogen in the post-reionization universe

abstract: The evolution of neutral hydrogen (HI) across redshifts is a powerful probe of cosmology, large scale structure in the universe and the intergalactic medium. Using a data-driven halo model to describe the distribution of HI in the post-reionization universe ($z \sim $ 5 to 0), we obtain the best-fitting parameters from a rich sample of observational data: low redshift 21-cm emission line studies, intermediate redshift intensity mapping experiments, and higher redshift Damped Lyman Alpha (DLA) observations. Our model describes the abundance and clustering of neutral hydrogen across redshifts 0 - 5, and is useful for investigating different aspects of galaxy evolution and for comparison with hydrodynamical simulations. The framework can be applied for forecasting future observations with neutral hydrogen, and extended to the case of intensity mapping with molecular and other line transitions at intermediate redshifts.


title: The Bulge Metallicity Distribution from the APOGEE Survey

abstract: The Apache Point Observatory Galactic Evolution Experiment (APOGEE) provides spectroscopic information of regions of the inner Milky Way inaccessible to optical surveys. We present the first large study of the metallicity distribution of the innermost Galactic regions based on homogeneous measurements from the SDSS Data Release 12 for 7545 red giant stars within 4.5 kpc of the Galactic center, with the goal to shed light on the structure and origin of the Galactic bulge. Stellar metallicities are found, through multiple-Gaussian decompositions, to be distributed in several components indicative of the presence of various stellar populations such as the bar, or the thin and the thick disk. A super-solar ([Fe/H]=+0.32) and a solar ([Fe/H]=+0.00) metallicity components, tentatively associated with the thin disk and the Galactic bar, respectively, seem to be the major contributors near the midplane. The solar-metallicity component extends outwards in the midplane but is not observed in the innermost regions. The central regions (within 3 kpc of the Galactic center) reveal, on the other hand, the presence of a significant metal-poor population ([Fe/H] =-0.46), tentatively associated with the thick disk, and which becomes the dominant component far from the midplane ($|Z| > +0.75 kpc). Varying contributions from these different components produce a transition region at +0.5 kpc < |Z| < +1.0 kpc characterized by a significant vertical metallicity gradient.


title: Using time-dependent indirect $CP$ asymmetries to measure $T$ and $CPT$ violation in $B^0$-${\bar B}^0$ mixing

abstract: Quantum field theory, which is the basis for all of particle physics, requires that all processes respect $CPT$ invariance. It is therefore of paramount importance to test the validity of $CPT$ conservation. In this Letter, we show that the time-dependent, indirect $CP$ asymmetries involving $B$ decays to a $CP$ eigenstate contain enough information to measure $T$ and $CPT$ violation in $B^0$-${\bar B}^0$ mixing, in addition to the standard $CP$-violating weak phases. Entangled $B^0{\bar B}^0$ states are not required (so that this analysis can be carried out at LHCb, as well as at the $B$ factories), penguin pollution need not be neglected, and the measurements can be made even if the $B^0$-${\bar B}^0$ width difference vanishes.


title: Future prospects of mass-degenerate Higgs bosons in the $CP$-conserving two-Higgs-doublet model

abstract: The scenario of two mass-degenerate Higgs bosons within the general two-Higgs-doublet model (2HDM) is revisited. We focus on the global picture when two $CP$-even Higgs bosons of $h$ and $H$ are nearly mass-degenerate. A global fit to the signal strength of the 125 GeV Higgs measured at the LHC is performed. Based on the best-fit result of the 2HDM mixing angles $(\alpha,\beta)$, theoretical constraints, charged and $CP$-odd Higgs boson direct search constraints and the electroweak precision constraints are imposed to the 2HDM parameter space. We present the signal predictions of the $(4b\,, 2b\,2\gamma)$ channels for the benchmark models at the LHC 14 TeV runs. We also study the direct Higgs boson pair productions at the LHC, and the Z-associated Higgs boson pair production search at the ILC 500 GeV runs, as well as the indirect probes at the CEPC 250 GeV run. We find that the mass-degenerate Higgs boson scenario in the Type-II 2HDM can be fully probed by these future experimental searches.