title: Laplace transformed MP2 for three dimensional periodic materials using stochastic orbitals in the plane wave basis and correlated sampling

abstract: We present an implementation and analysis of a stochastic high performance algorithm to calculate the correlation energy of three dimensional periodic systems in second-order M{\o}ller-Plesset perturbation theory (MP2). In particular we measure the scaling behavior of the sample variance and probe whether this stochastic approach is competitive if accuracies well below 1 meV per valence orbital are required, as it is necessary for calculations of adsorption, binding, or surface energies. The algorithm is based on the Laplace transformed MP2 (LTMP2) formulation in the plane wave basis. The time-dependent Hartree-Fock orbitals, appearing in the LTMP2 formulation, are stochastically rotated in the occupied and unoccupied Hilbert space. This avoids a full summation over all combinations of occupied and unoccupied orbitals, as inspired by the work of D. Neuhauser, E. Rabani, and R. Baer in J. Chem. Theory Comput. 9, 24 (2013). Additionally, correlated sampling is introduced, accelerating the statistical convergence significantly.


title: Classification of binary self-dual [76, 38, 14] codes with an automorphism of order 9

abstract: Using the method for constructing binary self-dual codes with an automorphism of order square of a prime number we have classified all binary self-dual codes with length 76 having minimum distance $d=14$ and automorphism of order 9. Up to equivalence, there are six self-dual $[76, 38, 14]$ codes with an automorphism of type $9$-$(8,0,4)$. All codes obtained have new values of the parameter in their weight enumerator thus more than doubling the number of known values.


title: Nonradial solutions of nonlinear scalar field equations

abstract: We prove new results concerning the nonlinear scalar field equation \begin{equation*} \left\{ \begin{array}{ll} -\Delta u = g(u)&\quad \hbox{in }\mathbb{R}^N,\; N\geq 3, u\in H^1(\mathbb{R}^N)& \end{array} \right. \end{equation*} with a nonlinearity $g$ satisfying the general assumptions due to Berestycki and Lions. In particular, we find at least one nonradial solution for any $N\geq 4$ minimizing the energy functional on the Pohozaev constraint in a subspace of $H^1(\mathbb{R}^N)$ consisting of nonradial functions. If in addition $N\neq 5$, then there are infinitely many nonradial solutions. These solutions are sign-changing. The results give a positive answer to a question posed by Berestycki and Lions in [5,6]. Moreover, we build a critical point theory on a topological manifold, which enables us to solve the above equation as well as to treat new elliptic problems.


title: CSWA: Aggregation-Free Spatial-Temporal Community Sensing

abstract: In this paper, we present a novel community sensing paradigm -- {C}ommunity {S}ensing {W}ithout {A}ggregation}. CSWA is designed to obtain the environment information (e.g., air pollution or temperature) in each subarea of the target area, without aggregating sensor and location data collected by community members. CSWA operates on top of a secured peer-to-peer network over the community members and proposes a novel \emph{Decentralized Spatial-Temporal Compressive Sensing} framework based on \emph{Parallelized Stochastic Gradient Descent}. Through learning the \emph{low-rank structure} via distributed optimization, CSWA approximates the value of the sensor data in each subarea (both covered and uncovered) for each sensing cycle using the sensor data locally stored in each member's mobile device. Simulation experiments based on real-world datasets demonstrate that CSWA exhibits low approximation error (i.e., less than $0.2 ^\circ$C in city-wide temperature sensing task and $10$ units of PM2.5 index in urban air pollution sensing) and performs comparably to (sometimes better than) state-of-the-art algorithms based on the data aggregation and centralized computation.


title: Fighting fish and two-stack sortable permutations

abstract: In 2017, Duchi, Guerrini, Rinaldi and Schaeffer proposed a new family of combinatorial objects called "fighting fish", which are counted by the same formula as more classical objects, such as two-stack sortable permutations and non-separable planar maps. In this article, we explore the bijective aspect of fighting fish by establishing a bijection to two-stack sortable permutations, using a new recursive decomposition of these permutations. With our bijection, we give combinatorial explanations of several results on fighting fish proved previously with generating functions. Using the decomposition of two-stack sortable permutations, we also prove the algebraicity of their generating function, extending a result of Bousquet-M\'elou (1998).


title: The interplay between activity and filament flexibility determines the emergent properties of active nematics

abstract: Active nematics are microscopically driven liquid crystals that exhibit dynamical steady states characterized by the creation and annihilation of topological defects. Motivated by experimental realizations of such systems made of biopolymer filaments and motor proteins, we describe a large-scale simulation study of a particle-based computational model that explicitly incorporates the semiflexibility of the biopolymers. We find that energy injected into the system at the particle scale preferentially excites bend deformations, renormalizing the filament bend modulus to smaller values. The emergent characteristics of the active nematic depend on activity and flexibility only through this activity-renormalized bend modulus, demonstrating that material parameters such as the Frank `constants' must explicitly depend on activity in a continuum hydrodynamic description of an active nematic. Further, we present a systematic way to estimate these material parameters from observations of deformation fields and defect shapes in experimental or simulation data.

categories: cond-mat.soft cond-mat.stat-mech physics.bio-ph physics.comp-ph


title: BBQ-Networks: Efficient Exploration in Deep Reinforcement Learning for Task-Oriented Dialogue Systems

abstract: We present a new algorithm that significantly improves the efficiency of exploration for deep Q-learning agents in dialogue systems. Our agents explore via Thompson sampling, drawing Monte Carlo samples from a Bayes-by-Backprop neural network. Our algorithm learns much faster than common exploration strategies such as \epsilon-greedy, Boltzmann, bootstrapping, and intrinsic-reward-based ones. Additionally, we show that spiking the replay buffer with experiences from just a few successful episodes can make Q-learning feasible when it might otherwise fail.


title: Effective alignments as building blocks of flavour models

abstract: Flavour models typically rely on flavons - scalars that break the family symmetry by acquiring vacuum expectation values in specific directions. We develop the idea of effective alignments, i.e.\ cases where the contractions of multiple flavons give rise to directions that are hard or impossible to obtain directly by breaking the family symmetry. Focusing on the example where the symmetry is $S_4$, we list the effective alignments that can be obtained from flavons vacuum expectation values that arise naturally from $S_4$. Using those effective alignments as building blocks, it is possible to construct flavour models, for example by using the effective alignments in constrained sequential dominance models. We illustrate how to obtain several of the mixing schemes in the literature, and explicitly construct renormalizable models for three viable cases, two of which lead to trimaximal mixing scenarios.


title: Variational Bi-LSTMs

abstract: Recurrent neural networks like long short-term memory (LSTM) are important architectures for sequential prediction tasks. LSTMs (and RNNs in general) model sequences along the forward time direction. Bidirectional LSTMs (Bi-LSTMs) on the other hand model sequences along both forward and backward directions and are generally known to perform better at such tasks because they capture a richer representation of the data. In the training of Bi-LSTMs, the forward and backward paths are learned independently. We propose a variant of the Bi-LSTM architecture, which we call Variational Bi-LSTM, that creates a channel between the two paths (during training, but which may be omitted during inference); thus optimizing the two paths jointly. We arrive at this joint objective for our model by minimizing a variational lower bound of the joint likelihood of the data sequence. Our model acts as a regularizer and encourages the two networks to inform each other in making their respective predictions using distinct information. We perform ablation studies to better understand the different components of our model and evaluate the method on various benchmarks, showing state-of-the-art performance.


title: Phase-sensitive evidence for the sign-reversal s+- symmetry of the order parameter in an iron-pnictide superconductor using Nb/Ba1-xNaxFe2As2 Josephson junctions

abstract: Josephson current between two superconductors provides a phase sensitive tool for probing their pairing symmetries. Here we fabricate and study experimentally high-quality Josephson junctions between a conventional s-wave superconductor Nb and a multi-band iron-pnictide Ba$_{1-x}$Na$_x$Fe$_2$As$_2$. Junctions exhibit a large enough critical current density to preclude the d-wave symmetry of the order parameter in the pnictide. However, the $I_cR_n$ product is very small $\simeq 3~\mu$V, which is not consistent with the sign-preserving $s_{++}$ symmetry either. We argue that the small $I_cR_n$ value along with its unusual temperature dependence provide evidence for the $s_{\pm}$ symmetry of the order parameter in Ba$_{1-x}$Na$_x$Fe$_2$As$_2$. We conclude that it is the phase sensitivity of our junctions that leads to an almost complete (bellow a sub-percent) cancellation of opposite supercurrents from the sign-reversal $s_{\pm}$ bands in the pnictide.


title: Classical branches and entanglement structure in the wavefunction of cosmological fluctuations

abstract: The emergence of preferred classical variables within a many-body wavefunction is encoded in its entanglement structure in the form of redundant classical information shared between many spatially local subsystems. We show how such structure can be generated via cosmological dynamics from the vacuum state of a massless field, causing the wavefunction to branch into classical field configurations on large scales. An accelerating epoch first excites the vacuum into a superposition of classical fields as well as a highly sensitive bath of squeezed super-horizon modes. During a subsequent decelerating epoch, gravitational interactions allow these modes to collect information about longer-wavelength fluctuations. This information disperses spatially, creating long-range redundant correlations in the wavefunction. The resulting classical observables, preferred basis for decoherence, and system/environment decomposition emerge from the translation invariant wavefunction and Hamiltonian, rather than being fixed by hand. We discuss the role of squeezing, the cosmological horizon scale, and phase information in the wavefunction, as well as aspects of wavefunction branching for continuous variables and in field theories.


title: Electromagnetically induced transparency in an isotopically purified Nd$^{3+}$:YLiF$_4$ crystal

abstract: We report the first observation of electromagnetically induced transparency (EIT) in an isotopically purified Nd$^{3+}$:YLiF$_4$ crystal. This crystal demonstrates inhomogeneous broadening of optical transitions of about 35 MHz. EIT is observed in a symmetrical $\Lambda$-like system formed by two hyperfine sublevels of the ground state corresponding to a zero first order Zeeman (ZEFOZ) transition and a hyperfine sublevel of the excited state, which is not coupled to other ground-state sublevels. It is found that transmission of the probe field as a function of the two-photon detuning demonstrates a comb-like structure that can be attributed to superhyperfine coupling between Nd$^{3+}$ ions and fluorine nuclei. The observed structure can be resolved only in the vicinity of the ZEFOZ point where the homogeneous linewidth of the spin transition is sufficiently small. The results pave the way for implementing solid-state quantum memories based on off-resonant Raman interaction without spectral tailoring of optical transitions.


title: GUT monopoles, the Witten effect and QCD axion phenomenology

abstract: Taken individually, magnetic monopoles and axions are both well-motivated aspects of physics beyond the Standard Model. We demonstrate that by virtue of the Witten effect, their interplay is furthermore nontrivial; monopoles break the corresponding axial symmetries explicitly and result in an axion potential lacking the usual instanton-derived exponential suppression factor. As axions are distinguished by their extremely weak interactions, any regime where their effects are enhanced may be of interest. As an application of these findings, we demonstrate that a phenomenologically acceptable population of grand-unification monopoles in the early Universe can efficiently suppress both the QCD axion dark matter abundance and CMB isocurvature contribution.


title: Signal for a light singlet scalar at the LHC

abstract: In the general Higgs portal like models, the extra neutral scalar, $S$, can mix with the Standard Model (SM) Higgs boson, $H$. We perform an exploratory study focusing on the direct search for such a light singlet $S$ at the Large Hadron Collider (LHC). After careful study of the SM background, we find the process $pp\rightarrow t\bar{t}S$ followed by $S\rightarrow b\bar{b}$ can be used to investigate $S$ with mass in the $20<M_S<100$ GeV range, which has not been well explored at the LHC. The signal significance becomes meaningful with a luminosity around a few $\mbox{ab}^{-1}$. Also, we study the prospects of finding the light scalar at the future 100 TeV $pp$ collider, the $Z$ and Higgs factories. With similar luminosity, the current LEP limits on the mixing between $S$ and $H$ can be improved by at least one or two orders of magnitude.


title: Evidence of incoherent carriers associated with resonant impurity levels and their influence on superconductivity in the anomalous superconductor Pb$_{1-x}$Tl$_x$Te

abstract: We present a combined experimental and theoretical study of the evolution of the Fermi surface of the anomalous superconductor Pb$_{1-x}$Tl$_x$Te as a function of thallium concentration, drawing on a combination of magnetotransport measurements (Shubnikov de Haas oscillations and Hall coefficient), Angle Resolved Photoemission Spectroscopy (ARPES), and density functional theory (DFT) calculations of the electronic structure. Our results indicate that for Tl concentrations beyond a critical value the Fermi energy coincides with resonant impurity states in Pb$_{1-x}$Tl$_x$Te, and we rule out the presence of an additional valence band maximum at the Fermi energy. Through comparison to non-superconducting Pb$_{1-x}$Na$_x$Te we argue that the presence of these states at the Fermi energy provides the pairing interaction and thus also the anomalously high temperature superconductivity in this material.


title: Exact Limits of Inference in Coalescent Models

abstract: Recovery of population size history from molecular sequence data is an important problem in population genetics. Inference commonly relies on a coalescent model linking the population size history to genealogies. The high computational cost of estimating parameters from these models usually compels researchers to select a subset of the available data or to rely on non-sufficient summary statistics for statistical inference. We consider the problem of recovering the true population size history from two possible alternatives on the basis of coalescent time data. We give exact expressions for the probability of selecting the correct alternative in a variety of biologically interesting cases as a function of the separation between the alternative size histories, the number of individuals, loci, and the sampling times. The results are applied to human population history. This work has significant implications for optimal design when the inferential goal is to test scientific hypotheses about population size trajectories in coalescent models with and without recombination.


title: Extremes of multifractional Brownian motion

abstract: Let $B_{H}(t), t\geq [0,T], T\in(0,\infty)$ be the standard Multifractional Brownian Motion(mBm), in this contribution we are concerned with the exact asymptotics of \begin{eqnarray*} \mathbb{P}\left\{\sup_{t\in[0,T]}B_{H}(t)>u\right\} \end{eqnarray*} as $u\rightarrow\infty$. Mainly depended on the structures of $H(t)$, the results under several important cases are investigated.


title: Markov Decision Processes with Continuous Side Information

abstract: We consider a reinforcement learning (RL) setting in which the agent interacts with a sequence of episodic MDPs. At the start of each episode the agent has access to some side-information or context that determines the dynamics of the MDP for that episode. Our setting is motivated by applications in healthcare where baseline measurements of a patient at the start of a treatment episode form the context that may provide information about how the patient might respond to treatment decisions. We propose algorithms for learning in such Contextual Markov Decision Processes (CMDPs) under an assumption that the unobserved MDP parameters vary smoothly with the observed context. We also give lower and upper PAC bounds under the smoothness assumption. Because our lower bound has an exponential dependence on the dimension, we consider a tractable linear setting where the context is used to create linear combinations of a finite set of MDPs. For the linear setting, we give a PAC learning algorithm based on KWIK learning techniques.


title: Counter-propagating Fractional Hall states in mirror-symmetric Dirac semi-metals

abstract: The Landau bands of mirror symmetric 2D Dirac semi-metals (for example odd-layers of ABA-graphene) can be identified by their parity with respect to mirror symmetry. This symmetry facilitates a new class of counter-propagating Hall states at opposite but equal electron and hole filling factors $|\nu_{\pm}|=1/m$ ({\it m} odd). Here, we propose a Laughlin-like correlated liquid wavefunction, at the charge neutrality point, that exhibits fractionally charged quasi-particle/hole pair excitation of opposite parity. Using a bosonized one-dimensional edge state theory, we show that the longitudinal conductance of this state, $\sigma_{xx} = 2e^2/(m h)$, is robust to short-ranged inter-mode interactions.


title: A new method of measuring center-of-mass velocities of radially pulsating stars from high-resolution spectroscopy

abstract: We present a radial velocity analysis of 20 solar neighborhood RR Lyrae and 3 Population II Cepheids variables. We obtained high-resolution, moderate-to-high signal-to-noise ratio spectra for most stars and obtained spectra were covering different pulsation phases for each star. To estimate the gamma (center-of-mass) velocities of the program stars, we use two independent methods. The first, `classic' method is based on RR Lyrae radial velocity curve templates. The second method is based on the analysis of absorption line profile asymmetry to determine both the pulsational and the gamma velocities. This second method is based on the Least Squares Deconvolution (LSD) technique applied to analyze the line asymmetry that occurs in the spectra. We obtain measurements of the pulsation component of the radial velocity with an accuracy of $\pm$ 3.5 km s$^{-1}$. The gamma velocity was determined with an accuracy $\pm$ 10 km s$^{-1}$, even for those stars having a small number of spectra. The main advantage of this method is the possibility to get the estimation of gamma velocity even from one spectroscopic observation with uncertain pulsation phase. A detailed investigation of the LSD profile asymmetry shows that the projection factor $p$ varies as a function of the pulsation phase -- this is a key parameter which converts observed spectral line radial velocity variations into photospheric pulsation velocities. As a byproduct of our study, we present 41 densely-spaced synthetic grids of LSD profile bisectors that are based on atmospheric models of RR Lyr covering all pulsation phases.


title: Single and multiple recurrence along non-polynomial sequences

abstract: We establish new recurrence and multiple recurrence results for a rather large family $\mathcal{F}$ of non-polynomial functions which includes tempered functions defined in [11], as well as functions from a Hardy field with the property that for some $\ell\in \mathbb{N}\cup\{0\}$, $\lim_{x\to\infty }f^{(\ell)}(x)=\pm\infty$ and $\lim_{x\to\infty }f^{(\ell+1)}(x)=0$. Among other things, we show that for any $f\in\mathcal{F}$, any invertible probability measure preserving system $(X,\mathcal{B},\mu,T)$, any $A\in\mathcal{B}$ with $\mu(A)>0$, and any $\epsilon>0$, the sets of returns $$ R_{\epsilon, A}= \big\{n\in\mathbb{N}:\mu(A\cap T^{-\lfloor f(n)\rfloor}A)>\mu^2(A)-\epsilon\big\} $$ and $$ R^{(k)}_{A}= \big\{ n\in\mathbb{N}: \mu\big(A\cap T^{\lfloor f(n)\rfloor}A\cap T^{\lfloor f(n+1)\rfloor}A\cap\cdots\cap T^{\lfloor f(n+k)\rfloor}A\big)>0\big\} $$ possess somewhat unexpected properties of largeness; in particular, they are thick, i.e., contain arbitrarily long intervals.


title: Perturbations against a Q-ball: Charge, energy, and additivity property

abstract: In the present paper, perturbations against a Q-ball solution are considered. It is shown that if we calculate the U(1) charge and the energy of the modes, which are solutions to linearized equations of motion, up to the second order in perturbations, we will get incorrect results. In particular, for the time-dependent modes we will obtain nonzero terms, which explicitly depend on time, indicating the nonconservation over time of the charge and the energy. It is shown that, as expected, this problem can be resolved by considering nonlinear equations of motion for the perturbations, providing second-order corrections to the solutions of linearized equations of motion. It turns out that contributions of these corrections to the charge and the energy can be taken into account without solving explicitly the nonlinear equations of motion for the perturbations. It is also shown that the use of such nonlinear equations not only recovers the conservation over time of the charge and the energy but also results in the additivity of the charge and the energy of different modes forming the perturbation.


title: Android Malware Detection using Markov Chain Model of Application Behaviors in Requesting System Services

abstract: Widespread growth in Android malwares stimulates security researchers to propose different methods for analyzing and detecting malicious behaviors in applications. Nevertheless, current solutions are ill-suited to extract the fine-grained behavior of Android applications accurately and efficiently. In this paper, we propose ServiceMonitor, a lightweight host-based detection system that dynamically detects malicious applications directly on mobile devices. ServiceMonitor reconstructs the fine-grained behavior of applications based on a novel systematic system service use analysis technique. Using proposed system service use perspective enables us to build a statistical Markov chain model to represent what and how system services are used to access system resources. Afterwards, we consider built Markov chain in the form of a feature vector and use it to classify the application behavior into either malicious or benign using Random Forests classification algorithm. ServiceMonitor outperforms current host-based solutions with evaluating it against 4034 malwares and 10024 benign applications and obtaining 96\% of accuracy rate and negligible overhead and performance penalty.


title: ParaNMT-50M: Pushing the Limits of Paraphrastic Sentence Embeddings with Millions of Machine Translations

abstract: We describe PARANMT-50M, a dataset of more than 50 million English-English sentential paraphrase pairs. We generated the pairs automatically by using neural machine translation to translate the non-English side of a large parallel corpus, following Wieting et al. (2017). Our hope is that ParaNMT-50M can be a valuable resource for paraphrase generation and can provide a rich source of semantic knowledge to improve downstream natural language understanding tasks. To show its utility, we use ParaNMT-50M to train paraphrastic sentence embeddings that outperform all supervised systems on every SemEval semantic textual similarity competition, in addition to showing how it can be used for paraphrase generation.


title: Characterization of Lifshitz transitions in topological nodal line semimetals

abstract: We introduce a two-band model of three-dimensional nodal line semimetals, the Fermi surface of which at half-filling may form various one-dimensional configurations of different topology. We study the symmetries and "drumhead" surface states of the model, and find that the transitions between different configurations, namely, the Lifshitz transitions, can be identified solely by the number of gap-closing points on some high-symmetry planes in the Brillouin zone. A global phase diagram of this model is also obtained accordingly. We then investigate the effect of some extra terms analogous to a two-dimensional Rashba-type spin-orbit coupling. The introduced extra terms open a gap for the nodal line semimetals and can be useful in engineering different topological insulating phases. We demonstrate that the behavior of surface Dirac cones in the resulting insulating system has a clear correspondence with the different configurations of the original nodal lines in the absence of the gap terms.


title: Chipmunk: A Systolically Scalable 0.9 mm${}^2$, 3.08 Gop/s/mW @ 1.2 mW Accelerator for Near-Sensor Recurrent Neural Network Inference

abstract: Recurrent neural networks (RNNs) are state-of-the-art in voice awareness/understanding and speech recognition. On-device computation of RNNs on low-power mobile and wearable devices would be key to applications such as zero-latency voice-based human-machine interfaces. Here we present Chipmunk, a small (<1 mm${}^2$) hardware accelerator for Long-Short Term Memory RNNs in UMC 65 nm technology capable to operate at a measured peak efficiency up to 3.08 Gop/s/mW at 1.24 mW peak power. To implement big RNN models without incurring in huge memory transfer overhead, multiple Chipmunk engines can cooperate to form a single systolic array. In this way, the Chipmunk architecture in a 75 tiles configuration can achieve real-time phoneme extraction on a demanding RNN topology proposed by Graves et al., consuming less than 13 mW of average power.


title: Is 1I/2017 U1 really of interstellar origin ?

abstract: Here we investigate the possibility that the asteroid 1I/2017 U1 actually is a Solar System object, currently expelled from the Solar System by the recent encouter with a Solar System planet. We show that this possibility is extremely unlikeley and that 1I/2017 U1 really is an interstellar object.


title: The ALMA-PILS survey: Formaldehyde deuteration in warm gas on small scales toward IRAS 16293-2422 B

abstract: [abridged] The enhanced degrees of deuterium fractionation observed in envelopes around protostars demonstrate the importance of chemistry at low temperatures, relevant in pre- and protostellar cores. formaldehyde is an important species in the formation of methanol and more complex molecules. Here, we present the first study of formaldehyde deuteration on small scales around the prototypical low-mass protostar IRAS 16293-2422 using high spatial and spectral resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations. Numerous isotopologues of formaldehyde are detected, among them H$_2$C$^{17}$O, and D$_2^{13}$CO for the first time in the ISM. The large range of upper energy levels covered by the HDCO lines help constrain the excitation temperature to 106$\pm$13 K. Using the derived column densities, formaldehyde shows a deuterium fractionation of HDCO/H$_2$CO=6.5$\pm$1%, D$_2$CO/HDCO=12.8$^{+3.3}_{-4.1}$%, and D$_2$CO/H$_2$CO=0.6(4)$\pm$0.1%. The isotopic ratios derived are $^{16}$O/$^{18}$O=805, $^{18}$O/$^{17}$O=3.2 and $^{12}$C/$^{13}$C=56. The HDCO/H$_2$CO ratio is lower than found in previous studies, highlighting the uncertainties involved in interpreting single dish observations of the inner warm regions. The D$_2$CO/HDCO ratio is only slightly larger than the HDCO/H$_2$CO ratio. This is consistent with formaldehyde forming in the ice as soon as CO has frozen onto the grains, with most of the deuteration happening towards the end of the prestellar core phase. A comparison with available time-dependent chemical models indicates that the source is in the early Class 0 stage.


title: Computing mean logarithmic mass from muon counts in air shower experiments

abstract: I discuss the conversion of muon counts in air showers, which are observable by experiments, into mean logarithmic mass, an important variable to express the mass composition of cosmic rays. Stochastic fluctuations in the shower development and statistical fluctuations from muon sampling can subtly bias the conversion. A central theme is that the mean of the logarithm of the muon number is not identical to the logarithm of the mean. It is discussed how that affects the conversion in practice. Simple analytical formulas to quantify and correct such biases are presented, which are applicable to any kind of experiment.


title: The Neural Network Pushdown Automaton: Model, Stack and Learning Simulations

abstract: In order for neural networks to learn complex languages or grammars, they must have sufficient computational power or resources to recognize or generate such languages. Though many approaches have been discussed, one ob- vious approach to enhancing the processing power of a recurrent neural network is to couple it with an external stack memory - in effect creating a neural network pushdown automata (NNPDA). This paper discusses in detail this NNPDA - its construction, how it can be trained and how useful symbolic information can be extracted from the trained network. In order to couple the external stack to the neural network, an optimization method is developed which uses an error function that connects the learning of the state automaton of the neural network to the learning of the operation of the external stack. To minimize the error function using gradient descent learning, an analog stack is designed such that the action and storage of information in the stack are continuous. One interpretation of a continuous stack is the probabilistic storage of and action on data. After training on sample strings of an unknown source grammar, a quantization procedure extracts from the analog stack and neural network a discrete pushdown automata (PDA). Simulations show that in learning deterministic context-free grammars - the balanced parenthesis language, 1*n0*n, and the deterministic Palindrome - the extracted PDA is correct in the sense that it can correctly recognize unseen strings of arbitrary length. In addition, the extracted PDAs can be shown to be identical or equivalent to the PDAs of the source grammars which were used to generate the training strings.


title: Planet-Planet Occultations in TRAPPIST-1 and Other Exoplanet Systems

abstract: We explore the occurrence and detectability of planet-planet occultations (PPOs) in exoplanet systems. These are events during which a planet occults the disk of another planet in the same system, imparting a small photometric signal as its thermal or reflected light is blocked. We focus on the planets in TRAPPIST-1, whose orbital planes we show are aligned to within 0.3 degrees at 90% confidence. We present a photodynamical model for predicting and computing PPOs in TRAPPIST-1 and other systems for various assumptions of the planets' atmospheric states. When marginalizing over the uncertainties on all orbital parameters, we find that the rate of PPOs in TRAPPIST-1 is about 1.4 per day. We investigate the prospects for detection of these events with the James Webb Space Telescope, finding that ~10-20 occultations per year of b and c should be above the noise level at 12-15 microns. Joint modeling of several of these PPOs could lead to a robust detection. Alternatively, observations with the proposed Origins Space Telescope should be able to detect individual PPOs at high signal-to-noise. We show how PPOs can be used to break transit timing variation degeneracies, imposing strong constraints on the eccentricities and masses of the planets, as well as to constrain the longitudes of nodes and thus the complete three-dimensional structure of the system. We further show how modeling of these events can be used to reveal a planet's day/night temperature contrast and construct crude surface maps. We make our photodynamical code available on github.


title: Possibilities for Probabilities

abstract: In ordinary situations involving a small part of the universe, Born's rule seems to work well for calculating probabilities of observations in quantum theory. However, there are a number of reasons for believing that it is not adequate for many cosmological purposes. Here a number of possible generalizations of Born's rule are discussed, explaining why they are consistent with the present statistical support for Born's rule in ordinary situations but can help solve various cosmological problems.


title: Accidental Peccei-Quinn Symmetry from Discrete Flavour Symmetry and Pati-Salam

abstract: We show how an accidental $U(1)$ Peccei-Quinn (PQ) symmetry can arise from a discrete $A_4$ family symmetry combined with a discrete flavour symmetry $ \mathbb{Z}_3 \times \mathbb{Z}_5^2 $, in a realistic Pati-Salam unified theory of flavour. Imposing only these discrete flavour symmetries, the axion solution to the strong $ CP $ problem is protected from PQ-breaking operators to the required degree. A QCD axion arises from a linear combination of $ A_4 $ triplet flavons, which are also responsible for fermion flavour structures due to their vacuum alignments. We find that the requirement of an accidental PQ symmetry arising from a discrete flavour symmetry constrains the form of the Yukawa matrices, providing a link between flavour and the strong $ CP $ problem. Our model predicts specific flavour-violating couplings of the flavourful axion and thus puts a strong limit on the axion scale from kaon decays.


title: Detection strategies for the first supernovae with JWST

abstract: Pair-instability supernovae (PISNe) are very luminous explosions of massive, low metallicity stars. They can potentially be observed out to high redshifts due to their high explosion energies, thus providing a probe of the Universe prior to reionization. The near-infrared camera, NIRCam, on board the James Webb Space Telescope is ideally suited for detecting their redshifted ultraviolet emission. We calculate the photometric signature of high-redshift PISNe and derive the optimal detection strategy for identifying their prompt emission and possible afterglow. We differentiate between PISNe and other sources that could have a similar photometric signature, such as active galactic nuclei or high-redshift galaxies. We demonstrate that the optimal strategy, which maximizes the visibility time of the PISN lightcurve per invested exposure time, consists of the two wide-band filters F200W and F356W with an exposure time of 600s. For such exposures, we expect one PISN at $z \lesssim 7.5$ per at least 50,000 different field of view, which can be accomplished with parallel observations and an extensive archival search. The PISN afterglow, caused by nebular emission and reverberation, is very faint and requires unfeasibly long exposure times to be uniquely identified. However, this afterglow would be visible for several hundred years, about two orders of magnitude longer than the prompt emission, rendering PISNe promising targets for future, even more powerful telescopes.


title: Corrections to di-Higgs boson production with light stops and modified Higgs couplings

abstract: The Higgs pair production in gluon fusion is a sensitive probe of beyond-Standard Model (BSM) phenomena and its detection is a major goal for the LHC and higher energy hadron collider experiments. In this work we reanalyze the possible modifications of the Higgs pair production cross section within low energy supersymmetry models. We show that the supersymmetric contributions to the Higgs pair production cross section are strongly correlated with the ones of the single Higgs production in the gluon fusion channel. Motivated by the analysis of ATLAS and CMS Higgs production data, we show that the scalar superpartners' contributions may lead to significant modification of the di-Higgs production rate and invariant mass distribution with respect to the SM predictions. We also analyze the combined effects on the di-Higgs production rate of a modification of the Higgs trilinear and top-quark Yukawa couplings in the presence of light stops. In particular, we show that due to the destructive interference of the triangle and box amplitude contributions to the di-Higgs production cross section, even a small modification of the top-quark Yukawa coupling can lead to a significant increase of the di-Higgs production rate.


title: Improving galaxy morphologies for SDSS with Deep Learning

abstract: We present a morphological catalogue for $\sim$ 670,000 galaxies in the Sloan Digital Sky Survey in two flavours: T-Type, related to the Hubble sequence, and Galaxy Zoo 2 (GZ2 hereafter) classification scheme. By combining accurate existing visual classification catalogues with machine learning, we provide the largest and most accurate morphological catalogue up to date. The classifications are obtained with Deep Learning algorithms using Convolutional Neural Networks (CNNs). We use two visual classification catalogues, GZ2 and Nair & Abraham (2010), for training CNNs with colour images in order to obtain T-Types and a series of GZ2 type questions (disk/features, edge-on galaxies, bar signature, bulge prominence, roundness and mergers). We also provide an additional probability enabling a separation between pure elliptical (E) from S0, where the T-Type model is not so efficient. For the T-Type, our results show smaller offset and scatter than previous models trained with support vector machines. For the GZ2 type questions, our models have large accuracy (> 97\%), precision and recall values (> 90\%) when applied to a test sample with the same characteristics as the one used for training. The catalogue is publicly released with the paper.


title: Approximations for the Lowest Energy Eigenstates in a Double Square Well Potential

abstract: Highly accurate closed-form approximations are given for the ground state and first excited state wavefunctions and energies for a nonrelativistic particle in a one-dimensional double square well potential with a square barrier in between (so that the potential is a sequence of five constant pieces that alternate in value from being above and below the ground state and first excited state energies), under the assumption that the barrier is sufficiently wide that the tunneling across it is very small.


title: From Dirac semimetals to topological phases in three dimensions: a coupled wire construction

abstract: Weyl and Dirac (semi)metals in three dimensions have robust gapless electronic band structures. Their massless single-body energy spectra are protected by symmetries such as lattice translation, (screw) rotation and time reversal. In this manuscript, we discuss many-body interactions in these systems. We focus on strong interactions that preserve symmetries and are outside the single-body mean-field regime. By mapping a Dirac (semi)metal to a model based on a three dimensional array of coupled Dirac wires, we show (1) the Dirac (semi)metal can acquire a many-body excitation energy gap without breaking the relevant symmetries, and (2) interaction can enable an anomalous Weyl (semi)metallic phase that is otherwise forbidden by symmetries in the single-body setting and can only be present holographically on the boundary of a four dimensional weak topological insulator. Both of these topological states support fractional gapped (gapless) bulk (resp. boundary) quasiparticle excitations.


title: Exploring Speech Enhancement with Generative Adversarial Networks for Robust Speech Recognition

abstract: We investigate the effectiveness of generative adversarial networks (GANs) for speech enhancement, in the context of improving noise robustness of automatic speech recognition (ASR) systems. Prior work demonstrates that GANs can effectively suppress additive noise in raw waveform speech signals, improving perceptual quality metrics; however this technique was not justified in the context of ASR. In this work, we conduct a detailed study to measure the effectiveness of GANs in enhancing speech contaminated by both additive and reverberant noise. Motivated by recent advances in image processing, we propose operating GANs on log-Mel filterbank spectra instead of waveforms, which requires less computation and is more robust to reverberant noise. While GAN enhancement improves the performance of a clean-trained ASR system on noisy speech, it falls short of the performance achieved by conventional multi-style training (MTR). By appending the GAN-enhanced features to the noisy inputs and retraining, we achieve a 7% WER improvement relative to the MTR system.


title: Phenomenological consequences of superfluid dark matter with baryon-phonon coupling

abstract: Recently, a new form of dark matter has been suggested to naturally reproduce the empirically successful aspects of Milgrom's law in galaxies. The dark matter particle candidates are axion-like, with masses of order eV and strong self-interactions. They Bose-Einstein condense into a superfluid phase in the central regions of galaxy halos. The superfluid phonon excitations in turn couple to baryons and mediate an additional long-range force. For a suitable choice of the superfluid equation of state, this force can mimic Milgrom's law. In this paper we develop in detail some of the main phenomenological consequences of such a formalism, by revisiting the expected dark matter halo profile in the presence of an extended baryon distribution. In particular, we show how rotation curves of both high and low surface brightness galaxies can be reproduced, with a slightly rising rotation curve at large radii in massive high surface brightness galaxies, thus subtly different from Milgrom's law. We finally point out other expected differences with Milgrom's law, in particular in dwarf spheroidal satellite galaxies, tidal dwarf galaxies, and globular clusters, whose Milgromian or Newtonian behavior depends on the position with respect to the superfluid core of the host galaxy. We also expect ultra-diffuse galaxies within galaxy clusters to have velocities slightly above the baryonic Tully-Fisher relation. Finally, we note that, in this framework, photons and gravitons follow the same geodesics, and that galaxy-galaxy lensing, probing larger distances within galaxy halos than rotation curves, should follow predictions closer to the standard cosmological model than those of Milgrom's law.


title: First and second $K$-groups of an elliptic curve over a global field of positive characteristic

abstract: In this paper, we show that the maximal divisible subgroup of groups $K_1$ and $K_2$ of an elliptic curve $E$ over a function field is uniquely divisible. Further those $K$-groups modulo this uniquely divisible subgroup are explicitly computed. We also calculate the motivic cohomology groups of the minimal regular model of $E$, which is an elliptic surface over a finite field.


title: A Multi-wavelength Analysis of Binary-AGN Candidate PSO J334.2028+01.4075

abstract: We present analysis of the first Chandra observation of PSO J334.2028+01.4075 (PSO J334), targeted as a binary-AGN candidate based on periodic variations of the optical flux. With no prior targeted X-ray coverage for PSO J334, our new 40 ksec Chandra observation allows for the opportunity to differentiate between a single or binary-AGN system, and if a binary, can characterize the mode of accretion. Simulations show that the two expected accretion disk morphologies for binary-AGN systems are (i) a "cavity", where the inner region of the accretion disk is mostly empty and emission is truncated blueward of the wavelength associated with the temperature of the innermost ring, or (ii) "minidisks", where there is substantial accretion from the cirumbinary disk onto one or both of the members of the binary, each with their own shock-heated thin-disk accretion system. We find the X-ray emission to be well-fit with an absorbed power-law, incompatible with the simple cavity scenario. Further, we construct an SED of PSO J334 by combining radio through X-ray observations and find that the SED agrees well with that of a normal AGN, most likely incompatible with the minidisk scenario. Other analyses, such as locating the quasar on IR color-color diagrams and analyzing the quasar mass predicted by the fundamental plane of black hole activity, further highlight the similarity of PSO J334 with respect to normal AGN. On the multi-wavelength fronts we investigated, we find no evidence supporting PSO J334 as a binary-AGN system, though our analysis remains insensitive to some binary configurations.


title: The Gaia-ESO Survey: Churning through the Milky Way

abstract: We attempt to determine the relative fraction of stars that have undergone significant radial migration by studying the orbital properties of metal-rich ([Fe/H]$>0.1$) stars within 2 kpc of the Sun using a sample of more than 3,000 stars selected from iDR4 of the Gaia-ESO Survey. We investigate the kinematic properties, such as velocity dispersion and orbital parameters, of stellar populations near the sun as a function of [Mg/Fe] and [Fe/H], which could show evidence of a major merger in the past history of the Milky Way. This was done using the stellar parameters from the Gaia-ESO Survey along with proper motions from PPMXL to determine distances, kinematics, and orbital properties for these stars to analyze the chemodynamic properties of stellar populations near the Sun. Analyzing the kinematics of the most metal-rich stars ([Fe/H]$>0.1$), we find that more than half have small eccentricities ($e<0.2$) or are on nearly circular orbits. Slightly more than 20\% of the metal-rich stars have perigalacticons $R_p>7$ kpc. We find that the highest [Mg/Fe], metal-poor populations have lower vertical and radial velocity dispersions compared to lower [Mg/Fe] populations of similar metallicity by $\sim10$ km s$^{-1}$. The median eccentricity increases linearly with [Mg/Fe] across all metallicities, while the perigalacticon decreases with increasing [Mg/Fe] for all metallicities. Finally, the most [Mg/Fe]-rich stars are found to have significant asymmetric drift and rotate more than 40 km s$^{-1}$ slower than stars with lower [Mg/Fe] ratios. While our results cannot constrain how far stars have migrated, we propose that migration processes are likely to have played an important role in the evolution of the Milky Way, with metal-rich stars migrating from the inner disk toward to solar neighborhood and past mergers potentially driving enhanced migration of older stellar populations in the disk.


title: Quantum Origami: Transversal Gates for Quantum Computation and Measurement of Topological Order

abstract: In topology, a torus remains invariant under certain non-trivial transformations known as modular transformations. In the context of topologically ordered quantum states of matter, these transformations encode the braiding statistics and fusion rules of emergent anyonic excitations and thus serve as a diagnostic of topological order. Moreover, modular transformations of higher genus surfaces, e.g. a torus with multiple handles, can enhance the computational power of a topological state, in many cases providing a universal fault-tolerant set of gates for quantum computation. However, due to the intrusive nature of modular transformations, which abstractly involve global operations and manifold surgery, physical implementations of them in local systems have remained elusive. Here, we show that by folding manifolds, modular transformations can be applied in a single shot by independent local unitaries, providing a novel class of transversal logic gates for fault-tolerant quantum computation. Specifically, we demonstrate that multi-layer topological states with appropriate boundary conditions and twist defects allow modular transformations to be effectively implemented by a finite sequence of local SWAP gates between the layers. We further provide methods to directly measure the modular matrices, and thus the fractional statistics of anyonic excitations, providing a novel way to directly measure topological order.

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


title: Thomas Precession for Dressed Particles

abstract: We consider a particle dressed with boundary gravitons in three-dimensional Minkowski space. The existence of BMS transformations implies that the particle's wavefunction picks up a Berry phase when subjected to changes of reference frames that trace a closed path in the asymptotic symmetry group. We evaluate this phase and show that, for BMS superrotations, it provides a gravitational generalization of Thomas precession. In principle, such phases are observable signatures of asymptotic symmetries.

categories: hep-th cond-mat.str-el gr-qc math.RT quant-ph


title: Spaces of Types in Positive Model Theory

abstract: We introduce a notion of the space of types in positive model theory based on Stone duality for distributive lattices. We show that this space closely mirrors the Stone space of types in the full first-order model theory with negation (Tarskian model theory). We use this to generalise some classical results on countable models from the Tarskian setting to positive model theory.


title: Bottom and Charm Mass determinations from global fits to $Q\bar{Q}$ bound states at N$^3$LO

abstract: The bottomonium spectrum up to $n = 3$ is studied within Non-Relativistic Quantum Chromodynamics up to N$^3$LO. We consider finite charm quark mass effects both in the QCD potential and the $\overline{\mathrm{MS}}$-pole mass relation up to third order in the $\Upsilon$-scheme counting. The $u = 1/2$ renormalon of the static potential is canceled by expressing the bottom quark pole mass in terms of the MSR mass. A careful investigation of scale variation reveals that, while $n = 1, 2$ states are well behaved within perturbation theory, $n = 3$ bound states are no longer reliable. We carry out our analysis in the $n_\ell = 3$ and $n_\ell = 4$ schemes and conclude that, as long as finite $m_c$ effects are smoothly incorporated in the MSR mass definition, the difference between the two schemes is rather small. Performing a fit to $b\bar{b}$ bound states we find $\overline{m}_b(\overline{m}_b) = 4.216\pm 0.039$ GeV. We extend our analysis to the lowest lying charmonium states finding $\overline{m}_c(\overline{m}_c)=1.273 \pm 0.054$ GeV. Finally, we perform simultaneous fits for $\overline{m}_b$ and $\alpha_s$ finding $\alpha_s^{(n_f=5)}(m_Z)=0.1178\pm 0.0051$. Additionally, using a modified version of the MSR mass with lighter massive quarks we are able to predict the uncalculated $\mathcal{O}(\alpha_s^4)$ virtual massive quark corrections to the relation between the $\overline{\mathrm{MS}}$ and pole masses.


title: An excess of low-mass X-ray binaries in the outer halo of NGC 4472

abstract: We present new Chandra observations of the outer halo of the giant elliptical galaxy NGC 4472 (M49) in the Virgo Cluster. The data extend to 130 kpc (28'), and have a combined exposure time of 150 ks. After identifying optical counterparts using the Next Generation Virgo Cluster Survey to remove background active galactic nuclei and globular cluster (GC) sources, and correcting for completeness, we find that the number of field low-mass X-ray binaries (LMXBs) per unit stellar V-band light increases significantly with galactocentric radius. Because the flux limit of the complete sample corresponds to the Eddington limit for neutron stars in NGC 4472, many of the ~90 field LMXBs in this sample could host black holes. The excess of field LMXBs at large galactocentric radii may be partially caused by natal kicks on black holes and neutron stars in binary systems in the inner part of the galaxy. Furthermore, since the metallicity in the halo of NGC 4472 strongly decreases towards larger galactocentric radii, the number of field LMXBs per unit stellar mass is anti-correlated with metallicity, opposite to what is observed in GCs. Another way to explain the spatial distribution of field LMXBs is therefore a reversed metallicity effect, although we have not identified a mechanism to explain this in terms of stellar and binary evolution.


title: The Impact of LHC Jet Data on the MMHT PDF Fit at NNLO

abstract: We investigate the impact of the high precision ATLAS and CMS 7 TeV measurements of inclusive jet production on the MMHT global PDF analysis at next-to-next-to-leading order (NNLO). This is made possible by the recent completion of the long-term project to calculate the NNLO corrections to the hard cross section. We find that a good description of the ATLAS data is not possible with the default treatment of experimental systematic errors, and develop a simplified solution that retains the dominant physical information of the data. We then investigate the fit quality and the impact on the gluon PDF central value and uncertainty when the ATLAS and CMS data are included in a MMHT fit. We consider both common choices for the factorization and renormalization scale, namely the inclusive jet transverse momentum, $p_\perp$, and the leading jet $p_\perp$, as well as the different jet radii for which the ATLAS and CMS data are made available. We find that the impact of these data on the gluon is relatively insensitive to these inputs, in particular the scale choice, while the inclusion of NNLO corrections tends to improve the data description somewhat, and gives a qualitatively similar though not identical impact on the gluon in comparison to NLO.


title: Hemispheric Preference and Cyclic Variation of Solar Filament Chirality from 2000 to 2016

abstract: It is well known that solar filaments are features in the solar atmosphere which show a hemispheric preference in their chirality. The hemispheric preference is such that the dextral chirality dominates in the northern hemisphere while the sinistral chirality dominates in the southern. Determining the strength and cyclic variation of the degree of this hemispheric preference, however, is challenging and tedious and thus needs to be automated. In this paper, we follow Dr. Pietro Bernasconi's algorithm (Bernasconi et al. 2005) to detect filament chirality with two parallel channels of application. The algorithm is applied to H-alpha images by "Advanced Automated Filament Detection and Characterization Code" (AAFDCC) (Bernasconi et al. 2005) or the algorithm is explained to a human and the human determines the chirality of the solar filament. We have conducted this exercise on data during the month of August from years 2000 to 2016 and we found that 83 % of our visually determined filaments follow the hemispheric chirality preference, while 58 % of automatically determined filaments follow it. Our visually compiled results have over 90 % agreement with those of Pevtsov et al. (2003), yet the visually determined chiralities of filaments disagree with automated determinations significantly. We find that the hemispheric preference remained the same between solar cycles 23 and 24 but the preference is very difficult to determine during the solar minimum of 2008-2010 primarily due to the absence of filaments.


title: Angle-dependent electron spin resonance of YbRh$_2$Si$_2$ measured with planar microwave resonators and in-situ rotation

abstract: We present a new experimental approach to investigate the magnetic properties of the anisotropic heavy-fermion system YbRh$_2$Si$_2$ as a function of crystallographic orientation. Angle-dependent electron spin resonance (ESR) measurements are performed at a low temperature of 1.6 K and at an ESR frequency of 4.4 GHz utilizing a superconducting planar microwave resonator in a $^4$He-cryostat in combination with in-situ sample rotation. The obtained ESR g-factor of YbRh$_2$Si$_2$ as a function of the crystallographic angle is consistent with results of previous measurements using conventional ESR spectrometers at higher frequencies and fields. Perspectives to implement this experimental approach into a dilution refrigerator and to reach the magnetically ordered phase of YbRh$_2$Si$_2$ are discussed.


title: Scaled Boundary Parametrizations in Isogeometric Analysis

abstract: This paper deals with a special class of parametrizations for Isogeometric Analysis (IGA). The so-called scaled boundary parametrizations are easy to construct and particularly attractive if only a boundary description of the computational domain is available. The idea goes back to the Scaled Boundary Finite Element Method (SB-FEM), which has recently been extended to IGA. We take here a different viewpoint and study these parametrizations as bivariate or trivariate B-spline functions that are directly suitable for standard Galerkin-based IGA. Our main results are first a general framework for this class of parametrizations, including aspects such as smoothness and regularity as well as generalizations to domains that are not star-shaped. Second, using the Poisson equation as example, we explain the relation between standard Galerkin-based IGA and the Scaled Boundary IGA by means of the Laplace-Beltrami operator. Further results concern the separation of integrals in both approaches and an analysis of the singularity in the scaling center. Among the computational examples we present a planar rotor geometry that stems from a screw compressor machine and compare different parametrization strategies.


title: Interstellar communication. III. Optimal frequency to maximize data rate

abstract: The optimal frequency for interstellar communication, using "Earth 2017" technology, was derived in papers I and II of this series (arXiv:1706.03795, arXiv:1706.05570). The framework included models for the loss of photons from diffraction (free space), interstellar extinction, and atmospheric transmission. A major limit of current technology is the focusing of wavelengths $\lambda<300\,$nm (UV). When this technological constraint is dropped, a physical bound is found at $\lambda\approx1\,$nm ($E\approx\,$keV) for distances out to kpc. While shorter wavelengths may produce tighter beams and thus higher data rates, the physical limit comes from surface roughness of focusing devices at the atomic level. This limit can be surpassed by beam-forming with electromagnetic fields, e.g. using a free electron laser, but such methods are not energetically competitive. Current lasers are not yet cost efficient at nm wavelengths, with a gap of two orders of magnitude, but future technological progress may converge on the physical optimum. We recommend expanding SETI efforts towards targeted (at us) monochromatic (or narrow band) X-ray emission at 0.5-2 keV energies.


title: Random gradient extrapolation for distributed and stochastic optimization

abstract: In this paper, we consider a class of finite-sum convex optimization problems defined over a distributed multiagent network with $m$ agents connected to a central server. In particular, the objective function consists of the average of $m$ ($\ge 1$) smooth components associated with each network agent together with a strongly convex term. Our major contribution is to develop a new randomized incremental gradient algorithm, namely random gradient extrapolation method (RGEM), which does not require any exact gradient evaluation even for the initial point, but can achieve the optimal ${\cal O}(\log(1/\epsilon))$ complexity bound in terms of the total number of gradient evaluations of component functions to solve the finite-sum problems. Furthermore, we demonstrate that for stochastic finite-sum optimization problems, RGEM maintains the optimal ${\cal O}(1/\epsilon)$ complexity (up to a certain logarithmic factor) in terms of the number of stochastic gradient computations, but attains an ${\cal O}(\log(1/\epsilon))$ complexity in terms of communication rounds (each round involves only one agent). It is worth noting that the former bound is independent of the number of agents $m$, while the latter one only linearly depends on $m$ or even $\sqrt m$ for ill-conditioned problems. To the best of our knowledge, this is the first time that these complexity bounds have been obtained for distributed and stochastic optimization problems. Moreover, our algorithms were developed based on a novel dual perspective of Nesterov's accelerated gradient method.


title: The University of Michigan Centimeter-Band All Stokes Blazar Monitoring Program: Single-Dish Polarimetry as a Probe of Parsec-Scale Magnetic Fields

abstract: The University of Michigan 26-m paraboloid was dedicated to obtaining linear polarization and total flux density observations of blazars from the mid-1960s until June 2012 providing an unprecedented record tracking centimeter-band variability over decades at 14.5, 8.0, and 4.8 GHz for both targeted objects and members of flux-limited samples. In the mid-1970s through the mid-1980s, and during the last decade of the program, observations were additionally obtained of circular polarization for a small sample of radio-bright (S>5Jy), active sources. Key program results include evidence supporting class-dependent differences in the magnetic field geometry of BL Lac and QSO jets, identification of linear polarization changes temporally associated with flux outbursts supporting a shock-in-jet scenario, and determination of the spectral evolution of the Stokes V amplitude and polarity for testing proposed models. Recent radiative transfer modeling during large flares supports a jet scenario with a kinetically-dominated, relativistic flow at parsec scales with embedded turbulent magnetic fields and dynamically-weak ordered components which may be helical; the circular polarization observations are consistent with linear-to-circular mode conversion within this turbulent jet environment.


title: Online Allocation with Traffic Spikes: Mixing Adversarial and Stochastic Models

abstract: Motivated by Internet advertising applications, online allocation problems have been studied extensively in various adversarial and stochastic models. While the adversarial arrival models are too pessimistic, many of the stochastic (such as i.i.d or random-order) arrival models do not realistically capture uncertainty in predictions. A significant cause for such uncertainty is the presence of unpredictable traffic spikes, often due to breaking news or similar events. To address this issue, a simultaneous approximation framework has been proposed to develop algorithms that work well both in the adversarial and stochastic models; however, this framework does not enable algorithms that make good use of partially accurate forecasts when making online decisions. In this paper, we propose a robust online stochastic model that captures the nature of traffic spikes in online advertising. In our model, in addition to the stochastic input for which we have good forecasting, an unknown number of impressions arrive that are adversarially chosen. We design algorithms that combine a stochastic algorithm with an online algorithm that adaptively reacts to inaccurate predictions. We provide provable bounds for our new algorithms in this framework. We accompany our positive results with a set of hardness results showing that our algorithms are not far from optimal in this framework. As a byproduct of our results, we also present improved online algorithms for a slight variant of the simultaneous approximation framework.


title: Constraints on core-collapse supernova progenitors from explosion site integral field spectroscopy

abstract: Observationally, supernovae (SNe) are divided into subclasses pertaining to their distinct characteristics. This diversity reflects the diversity in the progenitor stars. It is not entirely clear how different evolutionary paths leading massive stars to become a SN are governed by fundamental parameters such as progenitor initial mass and metallicity. This paper places constraints on progenitor initial mass and metallicity in distinct core-collapse SN subclasses, through a study of the parent stellar populations at the explosion sites. Integral field spectroscopy (IFS) of 83 nearby SN explosion sites with a median distance of 18 Mpc has been collected and analysed, enabling detection and spectral extraction of the parent stellar population of SN progenitors. From the parent stellar population spectrum, the initial mass and metallicity of the coeval progenitor are derived by means of comparison to simple stellar population models and strong-line methods. Additionally, near-infrared IFS was employed to characterise the star formation history at the explosion sites. No significant metallicity differences are observed among distinct SN types. The typical progenitor mass is found to be highest for SN Ic, followed by type Ib, then types IIb and II. SN IIn is the least associated with young stellar populations and thus massive progenitors. However, statistically significant differences in progenitor initial mass are observed only when comparing SNe IIn with other subclasses. Stripped-envelope SN progenitors with initial mass estimate lower than 25~$M_\odot$ are found; these are thought to be the result of binary progenitors. Confirming previous studies, these results support the notion that core-collapse SN progenitors cannot arise from single-star channel only, and both single and binary channels are at play in the production of core-collapse SNe. [ABRIDGED]

categories: astro-ph.SR astro-ph.CO astro-ph.GA astro-ph.HE


title: Fast Predictive Simple Geodesic Regression

abstract: Deformable image registration and regression are important tasks in medical image analysis. However, they are computationally expensive, especially when analyzing large-scale datasets that contain thousands of images. Hence, cluster computing is typically used, making the approaches dependent on such computational infrastructure. Even larger computational resources are required as study sizes increase. This limits the use of deformable image registration and regression for clinical applications and as component algorithms for other image analysis approaches. We therefore propose using a fast predictive approach to perform image registrations. In particular, we employ these fast registration predictions to approximate a simplified geodesic regression model to capture longitudinal brain changes. The resulting method is orders of magnitude faster than the standard optimization-based regression model and hence facilitates large-scale analysis on a single graphics processing unit (GPU). We evaluate our results on 3D brain magnetic resonance images (MRI) from the ADNI datasets.


title: Predicting vehicular travel times by modeling heterogeneous influences between arterial roads

abstract: Predicting travel times of vehicles in urban settings is a useful and tangible quantity of interest in the context of intelligent transportation systems. We address the problem of travel time prediction in arterial roads using data sampled from probe vehicles. There is only a limited literature on methods using data input from probe vehicles. The spatio-temporal dependencies captured by existing data driven approaches are either too detailed or very simplistic. We strike a balance of the existing data driven approaches to account for varying degrees of influence a given road may experience from its neighbors, while controlling the number of parameters to be learnt. Specifically, we use a NoisyOR conditional probability distribution (CPD) in conjunction with a dynamic bayesian network (DBN) to model state transitions of various roads. We propose an efficient algorithm to learn model parameters. We propose an algorithm for predicting travel times on trips of arbitrary durations. Using synthetic and real world data traces we demonstrate the superior performance of the proposed method under different traffic conditions.


title: Towards a non-perturbative calculation of Weak Hamiltonian Wilson coefficients

abstract: We propose a method to compute the Wilson coefficients of the weak effective Hamiltonian to all orders in the strong coupling constant using Lattice QCD simulations. We perform our calculations adopting an unphysically light weak boson mass of around $2~\mathrm{GeV}$. We demonstrate that systematic errors for the Wilson coefficients $C_1$ and $C_2$, related to the current-current four-quark operators, can be controlled and present a path towards precise determinations in subsequent works.


title: PackNet: Adding Multiple Tasks to a Single Network by Iterative Pruning

abstract: This paper presents a method for adding multiple tasks to a single deep neural network while avoiding catastrophic forgetting. Inspired by network pruning techniques, we exploit redundancies in large deep networks to free up parameters that can then be employed to learn new tasks. By performing iterative pruning and network re-training, we are able to sequentially "pack" multiple tasks into a single network while ensuring minimal drop in performance and minimal storage overhead. Unlike prior work that uses proxy losses to maintain accuracy on older tasks, we always optimize for the task at hand. We perform extensive experiments on a variety of network architectures and large-scale datasets, and observe much better robustness against catastrophic forgetting than prior work. In particular, we are able to add three fine-grained classification tasks to a single ImageNet-trained VGG-16 network and achieve accuracies close to those of separately trained networks for each task. Code available at https://github.com/arunmallya/packnet


title: Black phosphorus: A new bandgap tuning knob

abstract: Modern electronics rely on devices whose functionality can be adjusted by the end-user with an external 'knob'. A new tuning knob to modify the band gap of black phosphorus has been experimentally demonstrated.


title: Hierarchical On-Surface Synthesis of Deterministic Graphene Nanoribbon Heterojunctions

abstract: Bottom-up graphene nanoribbon (GNR) heterojunctions are nanoscale strips of graphene whose electronic structure abruptly changes across a covalently bonded interface. Their rational design offers opportunities for profound technological advancements enabled by their extraordinary structural and electronic properties. Thus far the most critical aspect of their synthesis, the control over sequence and position of heterojunctions along the length of a ribbon, has been plagued by randomness in monomer sequences emerging from step-growth copolymerization of distinct monomers. All bottom-up GNR heterojunction structures created so far have exhibited random sequences of heterojunctions and, while useful for fundamental scientific studies, are difficult to incorporate into functional nanodevices as a result. Here we describe a new hierarchical fabrication strategy that allows deterministic growth of bottom-up GNRs that preferentially exhibit a single heterojunction interface rather than a random statistical sequence of junctions along the ribbon. Such heterojunctions provide a viable new platform that could be directly used in functional GNR-based device applications at the molecular scale. Our hierarchical GNR fabrication strategy is based on differences in the dissociation energies of C-Br and C-I bonds that allow control over the growth sequence of the block-copolymers from which GNRs are formed, and consequently yields a significantly higher proportion of single-junction GNR heterostructures. Scanning tunnelling spectroscopy and density functional theory calculations confirm that hierarchically-grown heterojunctions between chevron GNR (cGNR) and binaphthyl-cGNR segments exhibit straddling Type I band alignment in structures that are only one atomic layer thick and 3 nm in width.


title: Latent Constraints: Learning to Generate Conditionally from Unconditional Generative Models

abstract: Deep generative neural networks have proven effective at both conditional and unconditional modeling of complex data distributions. Conditional generation enables interactive control, but creating new controls often requires expensive retraining. In this paper, we develop a method to condition generation without retraining the model. By post-hoc learning latent constraints, value functions that identify regions in latent space that generate outputs with desired attributes, we can conditionally sample from these regions with gradient-based optimization or amortized actor functions. Combining attribute constraints with a universal "realism" constraint, which enforces similarity to the data distribution, we generate realistic conditional images from an unconditional variational autoencoder. Further, using gradient-based optimization, we demonstrate identity-preserving transformations that make the minimal adjustment in latent space to modify the attributes of an image. Finally, with discrete sequences of musical notes, we demonstrate zero-shot conditional generation, learning latent constraints in the absence of labeled data or a differentiable reward function. Code with dedicated cloud instance has been made publicly available (https://goo.gl/STGMGx).


title: Asymptotic charges from soft scalars in even dimensions

abstract: We study asymptotic charges associated to a spin-zero analogue of Weinberg's soft photon and graviton theorems in even dimensions. Simple spacetime expressions for the charges are given, but unlike gravity or electrodynamics, the symmetry interpretation for the charges remains elusive. This work is a higher dimensional extension of the four dimensional case studied in arXiv:1703.07885 .


title: Solutions of the D-dimensional Schrodinger equation with the hyperbolic Poschl Teller potential plus modified ring shaped term

abstract: In this paper, we solve the D-dimensional Schr\"odinger equation with hyperbolic Poschl-Teller potential plus a generalized ring-shaped potential. After the separation of variable in the hyperspherical coordinate. We used Nikiforov-Uvarov (NU) method to solve the resulting radial equation and obtain explicitly the energy level and the corresponding wave function in closed form. The solutions to the angular part are solved using the NU approach as well.


title: End-to-end Training for Whole Image Breast Cancer Diagnosis using An All Convolutional Design

abstract: We develop an end-to-end training algorithm for whole-image breast cancer diagnosis based on mammograms. It requires lesion annotations only at the first stage of training. After that, a whole image classifier can be trained using only image level labels. This greatly reduced the reliance on lesion annotations. Our approach is implemented using an all convolutional design that is simple yet provides superior performance in comparison with the previous methods. On DDSM, our best single-model achieves a per-image AUC score of 0.88 and three-model averaging increases the score to 0.91. On INbreast, our best single-model achieves a per-image AUC score of 0.96. Using DDSM as benchmark, our models compare favorably with the current state-of-the-art. We also demonstrate that a whole image model trained on DDSM can be easily transferred to INbreast without using its lesion annotations and using only a small amount of training data. Code availability: https://github.com/lishen/end2end-all-conv


title: Reconstructing general plane quartics from their inflection lines

abstract: Let $C$ be a general plane quartic and let ${\rm Fl}(C)$ denote the configuration of inflection lines of $C$. We show that if $D$ is any plane quartic with the same configuration of inflection lines ${\rm Fl}(C)$, then the quartics $C$ and $D$ coincide.


title: Brittle fracture of polymer transient networks

abstract: We study the fracture of reversible double transient networks, constituted of water suspensions of entangled surfactant wormlike micelles reversibly linked by various amounts of telechelic polymers. We provide a state diagram that delineates the regime of fracture without necking of the filament from the regime where no fracture or break-up has been observed. We show that filaments fracture when stretched at a rate larger than the inverse of the slowest relaxation time of the networks. We quantitatively demonstrate that dissipation processes are not relevant in our experimental conditions and that, depending on the density of nodes in the networks, fracture occurs in the linear viscoelastic regime or in a non-linear regime. In addition, analysis of the crack opening profiles indicates deviations from a parabolic shape close to the crack tip for weakly connected networks. We demonstrate a direct correlation between the amplitude of the deviation from the parabolic shape and the amount of non linear viscoelasticity.


title: On the values of unipotent characters in bad characteristic

abstract: Let $G(q)$ be a Chevalley group over a finite field $F_q$. By Lusztig's and Shoji's work, the problem of computing the values of the unipotent characters of $G(q)$ is solved, in principle, by the theory of character sheaves; one issue in this solution is the determination of certain scalars relating two types of class functions on $G(q)$. We show that this issue can be reduced to the case where $q$ is a prime, which opens the way to use computer algebra methods. Here, and in a sequel to this article, we use this approach to solve a number of cases in groups of exceptional type which seemed hitherto out of reach.


title: Large-scale Analysis of Opioid Poisoning Related Hospital Visits in New York State

abstract: Opioid related deaths are increasing dramatically in recent years, and opioid epidemic is worsening in the United States. Combating opioid epidemic becomes a high priority for both the U.S. government and local governments such as New York State. Analyzing patient level opioid related hospital visits provides a data driven approach to discover both spatial and temporal patterns and identity potential causes of opioid related deaths, which provides essential knowledge for governments on decision making. In this paper, we analyzed opioid poisoning related hospital visits using New York State SPARCS data, which provides diagnoses of patients in hospital visits. We identified all patients with primary diagnosis as opioid poisoning from 2010-2014 for our main studies, and from 2003-2014 for temporal trend studies. We performed demographical based studies, and summarized the historical trends of opioid poisoning. We used frequent item mining to find co-occurrences of diagnoses for possible causes of poisoning or effects from poisoning. We provided zip code level spatial analysis to detect local spatial clusters, and studied potential correlations between opioid poisoning and demographic and social-economic factors.


title: Detecting Egregious Conversations between Customers and Virtual Agents

abstract: Virtual agents are becoming a prominent channel of interaction in customer service. Not all customer interactions are smooth, however, and some can become almost comically bad. In such instances, a human agent might need to step in and salvage the conversation. Detecting bad conversations is important since disappointing customer service may threaten customer loyalty and impact revenue. In this paper, we outline an approach to detecting such egregious conversations, using behavioral cues from the user, patterns in agent responses, and user-agent interaction. Using logs of two commercial systems, we show that using these features improves the detection F1-score by around 20% over using textual features alone. In addition, we show that those features are common across two quite different domains and, arguably, universal.


title: Disentangling the Galactic Halo with APOGEE: I. Chemical and Kinematical Investigation of Distinct Metal-Poor Populations

abstract: We find two chemically distinct populations separated relatively cleanly in the [Fe/H] - [Mg/Fe] plane, but also distinguished in other chemical planes, among metal-poor stars (primarily with metallicities [Fe/H] $< -0.9$) observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and analyzed for Data Release 13 (DR13) of the Sloan Digital Sky Survey. These two stellar populations show the most significant differences in their [X/Fe] ratios for the $\alpha$-elements, C+N, Al, and Ni. In addition to these populations having differing chemistry, the low metallicity high-Mg population (which we denote the HMg population) exhibits a significant net Galactic rotation, whereas the low-Mg population (or LMg population) has halo-like kinematics with little to no net rotation. Based on its properties, the origin of the LMg population is likely as an accreted population of stars. The HMg population shows chemistry (and to an extent kinematics) similar to the thick disk, and is likely associated with $\it in$ $\it situ$ formation. The distinction between the LMg and HMg populations mimics the differences between the populations of low- and high-$\alpha$ halo stars found in previous studies, suggesting that these are samples of the same two populations.


title: Integral Formulas for Higher Order Conformally Invariant Fermionic Operators

abstract: In this paper, we establish higher order Borel-Pompeiu formulas for conformally invariant fermionic operators in higher spin theory, which is the theory of functions on m-dimensional Euclidean space taking values in arbitrary irreducible representations of the Spin group. As applications, we provide higher order Cauchy's integral formulas for those fermionic operators. This paper continues the work of building up basic integral formulas for conformally invariant differential operators in higher spin theory.


title: High-harmonic generation in solids with and without topological edge states

abstract: High-harmonic generation (HHG) in the two topological phases of a finite, one-dimensional, periodic structure is investigated using a self-consistent time-dependent density functional theory (TDDFT) approach. For harmonic photon energies smaller than the band gap, the harmonic yield is found to differ up to fourteen orders of magnitude for the two topological phases. This giant topological effect is explained by the degree of destructive interference in the harmonic emission of all valence-band (and edge-state) electrons, which strongly depends on whether topological edge states are present or not. The combination of strong-field laser physics with topological condensed matter opens up new possibilities to electronically control strong-field-based light or particle sources or---vice versa---to steer by all optical means topological electronics.


title: Identifying the community structure of the international food-trade multi network

abstract: Achieving international food security requires improved understanding of how international trade networks connect countries around the world through the import-export flows of food commodities. The properties of food trade networks are still poorly documented, especially from a multi-network perspective. In particular, nothing is known about the community structure of food networks, which is key to understanding how major disruptions or 'shocks' would impact the global food system. Here we find that the individual layers of this network have densely connected trading groups, a consistent characteristic over the period 2001 to 2011. We also fit econometric models to identify social, economic and geographic factors explaining the probability that any two countries are co-present in the same community. Our estimates indicate that the probability of country pairs belonging to the same food trade community depends more on geopolitical and economic factors -- such as geographical proximity and trade agreements co-membership -- than on country economic size and/or income. This is in sharp contrast with what we know about bilateral-trade determinants and suggests that food country communities behave in ways that can be very different from their non-food counterparts.


title: Laser pulse probe of chirality of Cooper pairs

abstract: The internal chirality of Cooper pairs is shown to modify strongly the response of a superconductor to the local heating by a laser beam. The suppression of the chiral order parameter inside the hot spot appears to induce the supercurrents flowing around the spot region. The chirality affects also the sequential stage of thermal quench developing according to the Kibble-Zurek scenario: besides the generation of vortex-antivortex pairs the quench facilitates the formation of superconducting domains with different chirality. These fingerprints of the chiral superconducting state can be probed by any experimental techniques sensitive to the local magnetic field. The supercurrents encircling the hot spot originate from the inhomogeneity of the state with the broken time reversal symmetry and their detection would provide a convenient alternative to the search of the spontaneous edge currents sensitive to the boundary properties. Thus, the suggested setup can help to resolve the long-standing problem of unambiguous detection of type of pairing in Sr2RuO4 considered as a good candidate for chiral superconductivity.


title: A Parameter Estimation Method Using Linear Response Statistics: Numerical Scheme

abstract: This paper presents a numerical method to implement the parameter estimation method using response statistics that was recently formulated by the authors. The proposed approach formulates the parameter estimation problem of It\^o drift diffusions as a nonlinear least-squares problem. To avoid solving the model repeatedly when using an iterative scheme in solving the resulting least-squares problems, a polynomial surrogate model is employed on appropriate response statistics with smooth dependence on the parameters. The existence of minimizers of the approximate polynomial least-squares problems that converge to the solution of the true least square problem is established under appropriate regularity assumption of the essential statistics as functions of parameters. Numerical implementation of the proposed method is conducted on two prototypical examples that belong to classes of models with wide range of applications, including the Langevin dynamics and the stochastically forced gradient flows. Several important practical issues, such as the selection of the appropriate response operator to ensure the identifiability of the parameters and the reduction of the parameter space, are discussed. From the numerical experiments, it is found that the proposed approach is superior compared to the conventional approach that uses equilibrium statistics to determine the parameters.


title: WebRelate: Integrating Web Data with Spreadsheets using Examples

abstract: Data integration between web sources and relational data is a key challenge faced by data scientists and spreadsheet users. There are two main challenges in programmatically joining web data with relational data. First, most websites do not expose a direct interface to obtain tabular data, so the user needs to formulate a logic to get to different webpages for each input row in the relational table. Second, after reaching the desired webpage, the user needs to write complex scripts to extract the relevant data, which is often conditioned on the input data. Since many data scientists and end-users come from diverse backgrounds, writing such complex regular-expression based logical scripts to perform data integration tasks is unfortunately often beyond their programming expertise. We present WebRelate, a system that allows users to join semi-structured web data with relational data in spreadsheets using input-output examples. WebRelate decomposes the web data integration task into two sub-tasks of i) URL learning and ii) input-dependent web extraction. The first sub-task generates the URLs for the webpages containing the desired data for all rows in the relational table. WebRelate achieves this by learning a string transformation program using a few example URLs. The second sub-task uses examples of desired data to be extracted from the corresponding webpages and learns a program to extract the data for the other rows. We design expressive domain-specific languages for URL generation and web data extraction, and present efficient synthesis algorithms for learning programs in these DSLs from few input-output examples. We evaluate WebRelate on 88 real-world web data integration tasks taken from online help forums and Excel product team, and show that WebRelate can learn the desired programs within few seconds using only 1 example for the majority of the tasks.


title: Quantile Markov Decision Process

abstract: The goal of a traditional Markov decision process (MDP) is to maximize expected cumulativereward over a defined horizon (possibly infinite). In many applications, however, a decision maker may beinterested in optimizing a specific quantile of the cumulative reward instead of its expectation. In this paperwe consider the problem of optimizing the quantiles of the cumulative rewards of a Markov decision process(MDP), which we refer to as a quantile Markov decision process (QMDP). We provide analytical resultscharacterizing the optimal QMDP value function and present a dynamic programming-based algorithm tosolve for the optimal policy. The algorithm also extends to the MDP problem with a conditional value-at-risk(CVaR) objective. We illustrate the practical relevance of our model by evaluating it on an HIV treatmentinitiation problem, where patients aim to balance the potential benefits and risks of the treatment.


title: CMU LiveMedQA at TREC 2017 LiveQA: A Consumer Health Question Answering System

abstract: In this paper, we present LiveMedQA, a question answering system that is optimized for consumer health question. On top of the general QA system pipeline, we introduce several new features that aim to exploit domain-specific knowledge and entity structures for better performance. This includes a question type/focus analyzer based on deep text classification model, a tree-based knowledge graph for answer generation and a complementary structure-aware searcher for answer retrieval. LiveMedQA system is evaluated in the TREC 2017 LiveQA medical subtask, where it received an average score of 0.356 on a 3 point scale. Evaluation results revealed 3 substantial drawbacks in current LiveMedQA system, based on which we provide a detailed discussion and propose a few solutions that constitute the main focus of our subsequent work.


title: Rigid Graph Compression: Motif-based rigidity analysis for disordered fiber networks

abstract: Using particle-scale models to accurately describe property enhancements and phase transitions in macroscopic behavior is a major engineering challenge in composite materials science. To address some of these challenges, we use the graph theoretic property of rigidity to model me- chanical reinforcement in composites with stiff rod-like particles. We develop an efficient algorithmic approach called rigid graph compression (RGC) to describe the transition from floppy to rigid in disordered fiber networks ('rod-hinge systems'), which form the reinforcing phase in many composite systems. To establish RGC on a firm theoretical foundation, we adapt rigidity matroid theory to identify primitive topological network motifs that serve as rules for composing interacting rigid par- ticles into larger rigid components. This approach is computationally efficient and stable, because RGC requires only topological information about rod interactions (encoded by a sparse unweighted network) rather than geometrical details such as rod locations or pairwise distances (as required in rigidity matroid theory). We conduct numerical experiments on simulated two-dimensional rod-hinge systems to demonstrate that RGC closely approximates the rigidity percolation threshold for such systems, through comparison with the pebble game algorithm (which is exact in two dimensions). Importantly, whereas the pebble game is derived from Laman's condition and is only valid in two dimensions, the RGC approach naturally extends to higher dimensions.


title: Maintaining The Humanity of Our Models

abstract: Artificial intelligence and machine learning have been major research interests in computer science for the better part of the last few decades. However, all too recently, both AI and ML have rapidly grown to be media frenzies, pressuring companies and researchers to claim they use these technologies. As ML continues to percolate into daily life, we, as computer scientists and machine learning researchers, are responsible for ensuring we clearly convey the extent of our work and the humanity of our models. Regularizing ML for mass adoption requires a rigorous standard for model interpretability, a deep consideration for human bias in data, and a transparent understanding of a model's societal effects.


title: Aggregated Wasserstein Metric and State Registration for Hidden Markov Models

abstract: We propose a framework, named Aggregated Wasserstein, for computing a dissimilarity measure or distance between two Hidden Markov Models with state conditional distributions being Gaussian. For such HMMs, the marginal distribution at any time position follows a Gaussian mixture distribution, a fact exploited to softly match, aka register, the states in two HMMs. We refer to such HMMs as Gaussian mixture model-HMM (GMM-HMM). The registration of states is inspired by the intrinsic relationship of optimal transport and the Wasserstein metric between distributions. Specifically, the components of the marginal GMMs are matched by solving an optimal transport problem where the cost between components is the Wasserstein metric for Gaussian distributions. The solution of the optimization problem is a fast approximation to the Wasserstein metric between two GMMs. The new Aggregated Wasserstein distance is a semi-metric and can be computed without generating Monte Carlo samples. It is invariant to relabeling or permutation of states. The distance is defined meaningfully even for two HMMs that are estimated from data of different dimensionality, a situation that can arise due to missing variables. This distance quantifies the dissimilarity of GMM-HMMs by measuring both the difference between the two marginal GMMs and that between the two transition matrices. Our new distance is tested on tasks of retrieval, classification, and t-SNE visualization of time series. Experiments on both synthetic and real data have demonstrated its advantages in terms of accuracy as well as efficiency in comparison with existing distances based on the Kullback-Leibler divergence.


title: On the proximity of large primes

abstract: By a sphere-packing argument, we show that there are infinitely many pairs of primes that are close to each other for some metrics on the integers. In particular, for any numeration basis $q$, we show that there are infinitely many pairs of primes the base $q$ expansion of which differ in at most two digits. Likewise, for any fixed integer $t,$ there are infinitely many pairs of primes, the first $t$ digits of which are the same. In another direction, we show that, there is a constant $c$ depending on $q$ such that for infinitely many integers $m$ there are at least $c\log \log m$ primes which differ from $m$ by at most one base $q$ digit.


title: Spin-orbit coupled bosons in one dimension: emergent gauge field and Lifshitz transition

abstract: In the presence of strong spin-independent interactions and spin-orbit coupling, we show that the spinor Bose liquid confined to one spatial dimension undergoes an interaction- or density-tuned quantum phase transition similar to one theoretically proposed for itinerant magnetic solid-state systems. The order parameter describes broken $Z_2$ inversion symmetry, with the ordered phase accompanied by non-vanishing momentum which is generated by fluctuations of an emergent dynamical gauge field at the phase transition. This quantum phase transition has dynamical critical exponent $z \simeq 2$, typical of a Lifshitz transition, but is described by a nontrivial interacting fixed point. From direct numerical simulation of the microscopic model, we extract previously unknown critical exponents for this fixed point. Our model describes a realistic situation of 1D ultracold atoms with Raman-induced spin-orbit coupling, establishing this system as a platform for studying exotic critical behavior of the Hertz-Millis type.


title: Finer Grained Entity Typing with TypeNet

abstract: We consider the challenging problem of entity typing over an extremely fine grained set of types, wherein a single mention or entity can have many simultaneous and often hierarchically-structured types. Despite the importance of the problem, there is a relative lack of resources in the form of fine-grained, deep type hierarchies aligned to existing knowledge bases. In response, we introduce TypeNet, a dataset of entity types consisting of over 1941 types organized in a hierarchy, obtained by manually annotating a mapping from 1081 Freebase types to WordNet. We also experiment with several models comparable to state-of-the-art systems and explore techniques to incorporate a structure loss on the hierarchy with the standard mention typing loss, as a first step towards future research on this dataset.


title: A rank 18 Waring decomposition of $sM_{\langle 3\rangle}$ with 432 symmetries

abstract: The recent discovery that the exponent of matrix multiplication is determined by the rank of the symmetrized matrix multiplication tensor has invigorated interest in better understanding symmetrized matrix multiplication. I present an explicit rank 18 Waring decomposition of $sM_{\langle 3\rangle}$ and describe its symmetry group.