title: Energy storage in structural composites by introducing CNT fiber/polymer electrolyte interleaves

abstract: This work presents a method to produce structural composites capable of energy storage. They are produced by integrating thin sandwich structures of CNT fiber veils and an ionic liquid-based polymer electrolyte between carbon fiber plies, followed by infusion and curing of an epoxy resin. The resulting structure behaves simultaneously as an electric double-layer capacitor and a structural composite, with flexural modulus of 60 GPa and flexural strength of 153 MPa, combined with 88 mF/g of specific capacitance and the highest power (30 W/kg) and energy (37.5 mWh/kg) densities reported so far for structural supercapacitors. In-situ electrochemical measurements during 4-point bending show that electrochemical performance is retained up to fracture, with minor changes in equivalent series resistance for interleaves under compressive stress. En route to improving interlaminar properties we produce grid-shaped interleaves that enable mechanical interconnection of plies by the stiff epoxy. Synchrotron 3D X-ray tomography analysis of the resulting hierarchical structure confirms the formation of interlaminar epoxy joints. The manuscript discusses encapsulation role of epoxy, demonstrated by charge-discharge measurements of composites immersed in water, a deleterious agent for ionic liquids. Finally, we show different architectures free of current collector and electrical insulators, in which both CNT fiber and CF act as active electrodes.


title: Large Scale Clustering with Variational EM for Gaussian Mixture Models

abstract: This paper represents a preliminary (pre-reviewing) version of a sublinear variational algorithm for isotropic Gaussian mixture models (GMMs). Further developments of the algorithm for GMMs with diagonal covariance matrices (instead of isotropic clusters) and their corresponding benchmarking results have been published by TPAMI (doi:10.1109/TPAMI.2021.3133763) in the paper "A Variational EM Acceleration for Efficient Clustering at Very Large Scales". We kindly refer the reader to the TPAMI paper instead of this much earlier arXiv version (the TPAMI paper is also open access). Publicly available source code accompanies the paper (see https://github.com/variational-sublinear-clustering). Please note that the TPAMI paper does not contain the benchmark on the 80 Million Tiny Images dataset anymore because we followed the call of the dataset creators to discontinue the use of that dataset. The aim of the project (which resulted in this arXiv version and the later TPAMI paper) is the exploration of the current efficiency and large-scale limits in fitting a parametric model for clustering to data distributions. To reduce computational complexity, we used a clustering objective based on truncated variational EM (which reduces complexity for many clusters) in combination with coreset objectives (which reduce complexity for many data points). We used efficient coreset construction and efficient seeding to translate the theoretical sublinear complexity gains into an efficient algorithm. In applications to standard large-scale benchmarks for clustering, we then observed substantial wall-clock speedups compared to already highly efficient clustering approaches. To demonstrate that the observed efficiency enables applications previously considered unfeasible, we clustered the entire and unscaled 80 Million Tiny Images dataset into up to 32,000 clusters.


title: Deep sequential models for sampling-based planning

abstract: We demonstrate how a sequence model and a sampling-based planner can influence each other to produce efficient plans and how such a model can automatically learn to take advantage of observations of the environment. Sampling-based planners such as RRT generally know nothing of their environments even if they have traversed similar spaces many times. A sequence model, such as an HMM or LSTM, guides the search for good paths. The resulting model, called DeRRT*, observes the state of the planner and the local environment to bias the next move and next planner state. The neural-network-based models avoid manual feature engineering by co-training a convolutional network which processes map features and observations from sensors. We incorporate this sequence model in a manner that combines its likelihood with the existing bias for searching large unexplored Voronoi regions. This leads to more efficient trajectories with fewer rejected samples even in difficult domains such as when escaping bug traps. This model can also be used for dimensionality reduction in multi-agent environments with dynamic obstacles. Instead of planning in a high-dimensional space that includes the configurations of the other agents, we plan in a low-dimensional subspace relying on the sequence model to bias samples using the observed behavior of the other agents. The techniques presented here are general, include both graphical models and deep learning approaches, and can be adapted to a range of planners.


title: Molecular and Ionic Dipole Effects on the Electronic Properties of Si/SiO2 Grafted Alkylamine Monolayers

abstract: In this work, we demonstrate the tunability of electronic properties of Si/SiO2 substrate by molecular and ionic surface modifications. The change in the electronic properties such as the work function (WF) and electron affinity (EA), were experimentally measured by contact potential difference (CPD) technique and theoretically supported by DFT calculations. We attribute these molecular electronic effects mainly to the variations of molecular and surface dipoles of the ionic and neutral species. We have previously showed that for the alkylhalide monolayers, changing the tail group from Cl to I decreased the work function of the substrate. Here we report on the opposite trend of WF changes, i.e. increase of the WF, obtained by using the anions of those halides from Cl$^{-}$ to I$^{-}$. This trend was observed on self-assembled alkylamonium halide (-NH3$^{+}$ X$^{-}$, where X$^{-}$=Cl$^{-}$, Br$^{-}$, I$^{-}$) monolayers modified substrates. The monolayer' formation was supported by Ellipsometry measurements, X-Ray Photoelectron Spectroscopy and Atomic Force Microscopy. Comparison of the theoretical and experimental data suggests that ionic surface dipole depends mainly on the polarizability and the position of the counter halide anion along with the organization and packaging of the layer. The described ionic modification can be easily used for facile tailoring and design of the electronic properties Si/SiO2 substrates for various device applications.


title: On Isoclasses of Maximal Subalgebras Determined by Automorphisms

abstract: Let $k$ be an algebraically-closed field, and let $B = kQ/I$ be a basic, finite-dimensional associative $k$-algebra with $n := \dim_kB < \infty$. Previous work shows that the collection of maximal subalgebras of $B$ carries the structure of a projective variety, denoted by $\operatorname{msa} (Q)$, which only depends on the underlying quiver $Q$ of $B$. The automorphism group $\operatorname{Aut}_k(B)$ acts regularly on $\operatorname{msa} (Q)$. Since $\operatorname{msa} (Q)$ does not depend on the admissible ideal $I$, it is not necessarily easy to tell when two points of $\operatorname{msa} (Q)$ actually correspond to isomorphic subalgebras of $B$. One way to gain insight into this problem is to study $\operatorname{Aut}_k(B)$-orbits of $\operatorname{msa} (Q)$, and attempt to understand how isoclasses of maximal subalgebras decompose as unions of $\operatorname{Aut}_k(B)$-orbits. This paper investigates the problem for $B = kQ$, where $Q$ is a type $\mathbb{A}$ Dynkin quiver. We show that for such $B$, two maximal subalgebras with connected Ext quivers are isomorphic if and only if they lie in the same $\operatorname{Aut}_k(B)$-orbit of $\operatorname{msa} (Q)$.


title: Decay of plasmonic waves in Josephson junction chains

abstract: We study the damping of plasma waves in linear Josephson junction chains as well as in two capacitively coupled chains. In the parameter regime where the ground capacitance can be neglected, the theory of the antisymmetric mode in the double chain can be mapped onto the theory of a single chain. We consider two sources of relaxation: the scattering from quantum phase slips (QPS) and the interaction among plasmons related to the nonlinearity of the Josephson potential. The contribution to the relaxation rate $1/\tau$ from the nonlinearity scales with the fourth power of frequency $\omega$, while the phase-slip contribution behaves as a power law with a non-universal exponent. In the parameter regime where the charging energy related to the junction capacitance is much smaller than the Josephson energy, the amplitude of QPS is strongly suppressed. This makes the relaxation mechanism related to QPS efficient only at very low frequencies. As a result, for chains that are in the infrared limit on the insulating side of the superconductor-insulator transition, the quality factor $\omega\tau$ shows a strongly non-monotonic dependence on frequency, as was observed in a recent experiment.


title: Algebraic varieties are homeomorphic to varieties defined over number fields

abstract: We show that every affine or projective algebraic variety defined over the field of real or complex numbers is homeomorphic to a variety defined over the field of algebraic numbers. We construct such a homeomorphism by choosing a small deformation of the coefficients of the original equations. This method is based on the properties of Zariski equisingular families of varieties. Moreover we construct an algorithm, that, given a system of equations defining a variety $V$, produces a system of equations with algebraic coefficients of a variety homeomorphic to $V$


title: Stability and instability of self-gravitating relativistic matter distributions

abstract: We consider steady state solutions of the massive, asymptotically flat, spherically symmetric Einstein-Vlasov system, i.e., relativistic models of galaxies or globular clusters, and steady state solutions of the Einstein-Euler system, i.e., relativistic models of stars. Such steady states are embedded into one-parameter families parameterized by their central redshift $\kappa>0$. We prove their linear instability when $\kappa$ is sufficiently large, i.e., when they are strongly relativistic, and that the instability is driven by a growing mode. Our work confirms the scenario of dynamic instability proposed in the 1960s by Zel'dovich \& Podurets (for the Einstein-Vlasov system) and by Harrison, Thorne, Wakano, \& Wheeler (for the Einstein-Euler system). Our results are in sharp contrast to the corresponding non-relativistic, Newtonian setting. We carry out a careful analysis of the linearized dynamics around the above steady states and prove an exponential trichotomy result and the corresponding index theorems for the stable/unstable invariant spaces. Finally, in the case of the Einstein-Euler system we prove a rigorous version of the turning point principle which relates the stability of steady states along the one-parameter family to the winding points of the so-called mass-radius curve.


title: Multilevel Adaptive Sparse Grid Quadrature for Monte Carlo models

abstract: Many problems require to approximate an expected value by some kind of Monte Carlo (MC) sampling, e.g. molecular dynamics (MD) or simulation of stochastic reaction models (also termed kinetic Monte Carlo (kMC)). Often, we are furthermore interested in some integral of the MC model's output over the input parameters. We present a Multilevel Adaptive Sparse Grid strategy for the numerical integration of such problems where the integrand is implicitly defined by a Monte Carlo model. In this approach, we exploit different levels of sampling accuracy in the Monte Carlo model to reduce the overall computational costs compared to a single level approach. Unlike existing approaches for Multilevel Numerical Quadrature, our approach is not based on a telescoping sum, but we rather utilize the intrinsic multilevel structure of the sparse grids and the employed locally supported, piecewise linear basis functions. Besides illustrative toy models, we demonstrate the methodology on a realistic kMC model for CO oxidation. We find significant savings compared to the single level approach - often orders of magnitude.


title: Caterpillars in Erd\H{o}s-Hajnal

abstract: Let $T$ be a tree such that all its vertices of degree more than two lie on one path, that is, $T$ is a caterpillar subdivision. We prove that there exists $\epsilon>0$ such that for every graph $G$ with $|V(G)|\ge 2$ not containing $T$ as an induced subgraph, either some vertex has at least $\epsilon|V(G)|$ neighbours, or there are two disjoint sets of vertices $A,B$, both of cardinality at least $\epsilon|V(G)|$, where there is no edge joining $A$ and $B$. A consequence is: for every caterpillar subdivision $T$, there exists $c>0$ such that for every graph $G$ containing neither of $T$ and its complement as an induced subgraph, $G$ has a clique or stable set with at least $|V(G)|^c$ vertices. This extends a theorem of Bousquet, Lagoutte and Thomass\'e [JCTB 2015], who proved the same when $T$ is a path, and a recent theorem of Choromanski, Falik, Liebenau, Patel and Pilipczuk [Electron. J. Combin. 2018], who proved it when $T$ is a "hook".


title: Hypergraph polynomials and the Bernardi process

abstract: Recently O. Bernardi gave a formula for the Tutte polynomial $T(x,y)$ of a graph, based on spanning trees and activities just like the original definition, but using a fixed ribbon structure to order the set of edges in a different way for each tree. The interior polynomial $I$ is a generalization of $T(x,1)$ to hypergraphs. We supply a Bernardi-type description of $I$ using a ribbon structure on the underlying bipartite graph $G$. Our formula works because it is determined by the Ehrhart polynomial of the root polytope of $G$ in the same way as $I$ is. To prove this we interpret the Bernardi process as a way of dissecting the root polytope into simplices, along with a shelling order. We also show that our generalized Bernardi process gives a common extension of bijections (and their inverses) constructed by Baker and Wang between spanning trees and break divisors.


title: Exact spectral solution of two interacting run-and-tumble particles on a ring lattice

abstract: Exact solutions of interacting random walk models, such as 1D lattice gases, offer precise insight into the origin of nonequilibrium phenomena. Here, we study a model of run-and-tumble particles on a ring lattice interacting via hardcore exclusion. We present the exact solution for one and two particles using a generating function technique. For two particles, the eigenvectors and eigenvalues are explicitly expressed using two parameters reminiscent of Bethe roots, whose numerical values are determined by polynomial equations which we derive. The spectrum depends in a complicated way on the ratio of direction reversal rate to lattice jump rate, $\omega$. For both one and two particles, the spectrum consists of separate real bands for large $\omega$, which mix and become complex-valued for small $\omega$. At exceptional values of $\omega$, two or more eigenvalues coalesce such that the Markov matrix is non-diagonalizable. A consequence of this intricate parameter dependence is the appearance of dynamical transitions: non-analytic minima in the longest relaxation times as functions of $\omega$ (for a given lattice size). Exceptional points are theoretically and experimentally relevant in, e.g., open quantum systems and multichannel scattering. We propose that the phenomenon should be a ubiquitous feature of classical nonequilibrium models as well, and of relevance to physical observables in this context.


title: Local disc model in view of Gaia DR1 and RAVE data

abstract: We test the performance of the semi-analytic self-consistent Just-Jahrei{\ss} disc model (JJ model) with the astrometric data from the Tycho-Gaia Astrometric Solution (TGAS) sub-catalogue of the first Gaia data release (Gaia DR1), as well as the radial velocities from the fifth data release of the Radial Velocity Experiment survey (RAVE DR5). We use a sample of 19,746 thin disc stars from the TGAS$\times$RAVE cross-match selected in the local solar cylinder of 300 pc radius and 1 kpc height below the Galactic plane and simulate this sample via the forward modelling technique. First, we convert the predicted vertical density laws of the thin disc populations into a mock sample. Then the obtained mock populations are reddened with a 3D dust map and are subjected to the selection criteria corresponding to the RAVE and TGAS observational limitations as well as to additional cuts applied to the data sample. We calculate the quantities of interest separately at different heights above the Galactic plane taking into account the distance error effects separately in horizontal and vertical directions. We investigate the simulated sample in terms of the vertical number density profiles, Hess diagrams and velocity distribution functions. Basing on a good agreement of our simulations with the data, we conclude that our fiducial disc model confidently reproduces the vertical trends in the thin disc stellar population properties. Thus, it can serve as a starting point for the future extension of the JJ model to other Galactocentric distances.


title: Monoclinic SrIrO$_3$ - A Dirac semimetal produced by non-symmorphic symmetry and spin-orbit coupling

abstract: SrIrO$_3$ crystallizes in a monoclinic structure of distorted hexagonal perovskite at ambient pressure. The transport measurements show that the monoclinic SrIrO$_3$ is a low-carrier density semimetal, as in the orthorhombic perovskite polymorph. The electronic structure calculation indicates a semimetallic band structure with Dirac bands at two high-symmetry points of Brillouin zone only when spin-orbit coupling is incorporated, suggesting that the semimetallic state is produced by the strong spin-orbit coupling. We argue that the Dirac bands are protected by the non-symmorphic symmetry of lattice.


title: Cholinergic switch between two types of slow waves in cerebral cortex

abstract: Sleep slow waves are known to participate in memory consolidation, yet slow waves occurring under anesthesia present no positive effects on memory. Here, we shed light onto this paradox, based on a combination of extracellular recordings in vivo, in vitro, and computational models. We find two types of slow waves, based on analyzing the temporal patterns of successive slow-wave events. The first type is consistently observed in natural slow-wave sleep, while the second is shown to be ubiquitous under anesthesia. Network models of spiking neurons predict that the two slow wave types emerge due to a different gain on inhibitory vs excitatory cells and that different levels of spike-frequency adaptation in excitatory cells can account for dynamical distinctions between the two types. This prediction was tested in vitro by varying adaptation strength using an agonist of acetylcholine receptors, which demonstrated a neuromodulatory switch between the two types of slow waves. Finally, we show that the first type of slow-wave dynamics is more sensitive to external stimuli, which can explain how slow waves in sleep and anesthesia differentially affect memory consolidation, as well as provide a link between slow-wave dynamics and memory diseases.


title: On the maximum number of copies of H in graphs with given size and order

abstract: We study the maximum number $ex(n,e,H)$ of copies of a graph $H$ in graphs with given number of vertices and edges. We show that for any fixed graph $H$, $ex(n,e,H)$ is asymptotically realized by the quasi-clique provided that the edge density is sufficiently large. We also investigate a variant of this problem, when the host graph is bipartite.


title: RGB-D Object Detection and Semantic Segmentation for Autonomous Manipulation in Clutter

abstract: Autonomous robotic manipulation in clutter is challenging. A large variety of objects must be perceived in complex scenes, where they are partially occluded and embedded among many distractors, often in restricted spaces. To tackle these challenges, we developed a deep-learning approach that combines object detection and semantic segmentation. The manipulation scenes are captured with RGB-D cameras, for which we developed a depth fusion method. Employing pretrained features makes learning from small annotated robotic data sets possible. We evaluate our approach on two challenging data sets: one captured for the Amazon Picking Challenge 2016, where our team NimbRo came in second in the Stowing and third in the Picking task, and one captured in disaster-response scenarios. The experiments show that object detection and semantic segmentation complement each other and can be combined to yield reliable object perception.


title: QCD calculations of $B \to \pi, K$ form factors with higher-twist corrections

abstract: We update QCD calculations of $B \to \pi, K$ form factors at large hadronic recoil by including the subleading-power corrections from the higher-twist $B$-meson light-cone distribution amplitudes (LCDAs) up to the twist-six accuracy and the strange-quark mass effects at leading-power in $\Lambda/m_b$ from the twist-two $B$-meson LCDA $\phi_B^{+}(\omega, \mu)$. The higher-twist corrections from both the two-particle and three-particle $B$-meson LCDAs are computed from the light-cone QCD sum rules (LCSR) at tree level. In particular, we construct the local duality model for the twist-five and -six $B$-meson LCDAs, in agreement with the corresponding asymptotic behaviours at small quark and gluon momenta, employing the QCD sum rules in heavy quark effective theory at leading order in $\alpha_s$. The strange quark mass effects in semileptonic $B \to K$ form factors yield the leading-power contribution in the heavy quark expansion, consistent with the power-counting analysis in soft-collinear effective theory, and they are also computed from the LCSR approach due to the appearance of the rapidity singularities. We further explore the phenomenological aspects of the semileptonic $B \to \pi \ell \nu$ decays and the rare exclusive processes $B \to K \nu \nu$, including the determination of the CKM matrix element $|V_{ub}|$, the normalized differential $q^2$ distributions and precision observables defined by the ratios of branching fractions for the above-mentioned two channels in the same intervals of $q^2$.


title: Recent Results from ANITA

abstract: The ANtarctic Impulsive Transient Antenna (ANITA) long-duration balloon payload searches for Askaryan radio emission from ultra-high-energy ($>10^{18}$ eV) neutrinos interacting in Antarctic ice. ANITA is also sensitive to geomagnetic radio emission from extensive air showers (EAS). This talk summarizes recently released results from the third flight of ANITA, which flew during the 2014-2015 Austral summer. The most sensitive search from ANITA-III identified one neutrino candidate with an a priori background estimate of 0.7$^{+0.5}_{-0.3}$. When combined with previous flights, ANITA sets the best limits on diffuse neutrino flux at energies above $\sim10^{19.5}$ eV. Additionally, ANITA-III searches identified nearly 30 EAS candidates. One unusual event appears to correspond to an upward-going air shower, similar to an event from ANITA-I.


title: Variational Discriminator Bottleneck: Improving Imitation Learning, Inverse RL, and GANs by Constraining Information Flow

abstract: Adversarial learning methods have been proposed for a wide range of applications, but the training of adversarial models can be notoriously unstable. Effectively balancing the performance of the generator and discriminator is critical, since a discriminator that achieves very high accuracy will produce relatively uninformative gradients. In this work, we propose a simple and general technique to constrain information flow in the discriminator by means of an information bottleneck. By enforcing a constraint on the mutual information between the observations and the discriminator's internal representation, we can effectively modulate the discriminator's accuracy and maintain useful and informative gradients. We demonstrate that our proposed variational discriminator bottleneck (VDB) leads to significant improvements across three distinct application areas for adversarial learning algorithms. Our primary evaluation studies the applicability of the VDB to imitation learning of dynamic continuous control skills, such as running. We show that our method can learn such skills directly from \emph{raw} video demonstrations, substantially outperforming prior adversarial imitation learning methods. The VDB can also be combined with adversarial inverse reinforcement learning to learn parsimonious reward functions that can be transferred and re-optimized in new settings. Finally, we demonstrate that VDB can train GANs more effectively for image generation, improving upon a number of prior stabilization methods.


title: Nonlinear Electrodynamics and Optics of Graphene

abstract: Graphene is a two-dimensional material with strongly nonlinear electrodynamics and optical properties. We present some of our recent theoretical results on the quantum and non-perturbative quasi-classical theories of nonlinear effects in graphene, influence of substrates on graphene nonlinearities, plasma oscillations in graphene in the nonlinear regime and other effects.


title: Approximating mixed H\"older functions using random samples

abstract: Suppose $f : [0,1]^2 \rightarrow \mathbb{R}$ is a $(c,\alpha)$-mixed H\"older function that we sample at $l$ points $X_1,\ldots,X_l$ chosen uniformly at random from the unit square. Let the location of these points and the function values $f(X_1),\ldots,f(X_l)$ be given. If $l \ge c_1 n \log^2 n$, then we can compute an approximation $\tilde{f}$ such that $$ \|f - \tilde{f} \|_{L^2} = \mathcal{O}(n^{-\alpha} \log^{3/2} n), $$ with probability at least $1 - n^{2 -c_1}$, where the implicit constant only depends on the constants $c > 0$ and $c_1 > 0$.


title: Three plots about the Cremona groups

abstract: The first group of results of this paper concerns the compressibility of finite subgroups of the Cremona groups. The second concerns the embeddability of other groups in the Cremona groups and, conversely, the Cremona groups in other groups. The third concerns the connectedness of the Cremona groups.


title: Set Transformer: A Framework for Attention-based Permutation-Invariant Neural Networks

abstract: Many machine learning tasks such as multiple instance learning, 3D shape recognition, and few-shot image classification are defined on sets of instances. Since solutions to such problems do not depend on the order of elements of the set, models used to address them should be permutation invariant. We present an attention-based neural network module, the Set Transformer, specifically designed to model interactions among elements in the input set. The model consists of an encoder and a decoder, both of which rely on attention mechanisms. In an effort to reduce computational complexity, we introduce an attention scheme inspired by inducing point methods from sparse Gaussian process literature. It reduces the computation time of self-attention from quadratic to linear in the number of elements in the set. We show that our model is theoretically attractive and we evaluate it on a range of tasks, demonstrating the state-of-the-art performance compared to recent methods for set-structured data.


title: How Powerful are Graph Neural Networks?

abstract: Graph Neural Networks (GNNs) are an effective framework for representation learning of graphs. GNNs follow a neighborhood aggregation scheme, where the representation vector of a node is computed by recursively aggregating and transforming representation vectors of its neighboring nodes. Many GNN variants have been proposed and have achieved state-of-the-art results on both node and graph classification tasks. However, despite GNNs revolutionizing graph representation learning, there is limited understanding of their representational properties and limitations. Here, we present a theoretical framework for analyzing the expressive power of GNNs to capture different graph structures. Our results characterize the discriminative power of popular GNN variants, such as Graph Convolutional Networks and GraphSAGE, and show that they cannot learn to distinguish certain simple graph structures. We then develop a simple architecture that is provably the most expressive among the class of GNNs and is as powerful as the Weisfeiler-Lehman graph isomorphism test. We empirically validate our theoretical findings on a number of graph classification benchmarks, and demonstrate that our model achieves state-of-the-art performance.


title: Five Driving Forces of Multi-Access Edge Computing

abstract: The emergence of Multi-Access Edge Computing (MEC) technology aims at extending cloud computing capabilities to the edge of the wireless access networks. MEC provides real-time, high-bandwidth, low-latency access to radio network resources, allowing operators to open their networks to a new ecosystem and value chain. Moreover, it will provide a new insight to the design of future 5th Generation (5G) wireless systems. This paper describes five key technologies, including Network Function Vitalization (NFV), Software Defined Networking (SDN), Network Slicing, Information Centric Networking (ICN) and Internet of Things (IoT), that intensify the widespread of MEC and its adoption. Our goal is to provide the associativity between MEC and these five driving technologies in 5G context while identifying the open challenges, future directions, and tangible integration paths.


title: Singularity, Misspecification, and the Convergence Rate of EM

abstract: A line of recent work has analyzed the behavior of the Expectation-Maximization (EM) algorithm in the well-specified setting, in which the population likelihood is locally strongly concave around its maximizing argument. Examples include suitably separated Gaussian mixture models and mixtures of linear regressions. We consider over-specified settings in which the number of fitted components is larger than the number of components in the true distribution. Such misspecified settings can lead to singularity in the Fisher information matrix, and moreover, the maximum likelihood estimator based on $n$ i.i.d. samples in $d$ dimensions can have a non-standard $\mathcal{O}((d/n)^{\frac{1}{4}})$ rate of convergence. Focusing on the simple setting of two-component mixtures fit to a $d$-dimensional Gaussian distribution, we study the behavior of the EM algorithm both when the mixture weights are different (unbalanced case), and are equal (balanced case). Our analysis reveals a sharp distinction between these two cases: in the former, the EM algorithm converges geometrically to a point at Euclidean distance of $\mathcal{O}((d/n)^{\frac{1}{2}})$ from the true parameter, whereas in the latter case, the convergence rate is exponentially slower, and the fixed point has a much lower $\mathcal{O}((d/n)^{\frac{1}{4}})$ accuracy. Analysis of this singular case requires the introduction of some novel techniques: in particular, we make use of a careful form of localization in the associated empirical process, and develop a recursive argument to progressively sharpen the statistical rate.


title: Adaptive Game-Theoretic Decision Making for Autonomous Vehicle Control at Roundabouts

abstract: In this paper, we propose a decision making algorithm for autonomous vehicle control at a roundabout intersection. The algorithm is based on a game-theoretic model representing the interactions between the ego vehicle and an opponent vehicle, and adapts to an online estimated driver type of the opponent vehicle. Simulation results are reported.


title: A Kinematic Evolution Equation for the Dynamic Contact Angle and some Consequences

abstract: We investigate the moving contact line problem for two-phase incompressible flows with a kinematic approach. The key idea is to derive an evolution equation for the contact angle in terms of the transporting velocity field. It turns out that the resulting equation has a simple structure and expresses the time derivative of the contact angle in terms of the velocity gradient at the solid wall. Together with the additionally imposed boundary conditions for the velocity, it yields a more specific form of the contact angle evolution. Thus, the kinematic evolution equation is a tool to analyze the evolution of the contact angle. Since the transporting velocity field is required only on the moving interface, the kinematic evolution equation also applies when the interface moves with its own velocity independent of the fluid velocity. We apply the developed tool to a class of moving contact line models which employ the Navier slip boundary condition. We derive an explicit form of the contact angle evolution for sufficiently regular solutions, showing that such solutions are unphysical. Within the simplest model, this rigorously shows that the contact angle can only relax to equilibrium if some kind of singularity is present at the contact line. Moreover, we analyze more general models including surface tension gradients at the contact line, slip at the fluid-fluid interface and mass transfer across the fluid-fluid interface.


title: Nonlinear Electrodynamic Properties of Graphene and Other Two-Dimensional Materials

abstract: Graphene, the first truly two-dimensional (one atom thin) material, possesses strongly nonlinear electrodynamic and optical properties. At low (microwave, terahertz) frequencies this results from the unique electronic property of graphene - the "relativistic", linear energy dispersion of its electrons and holes. At high (infrared, visible) frequencies its nonlinear response functions have resonances due to the inter-band transitions between the electron and hole bands. The position of these resonances on the frequency axis depends on the Fermi energy, or the charge carrier density, which can be varied by applying a gate voltage in the graphene field-effect transistor. This opens up the unique opportunities to electrically control the nonlinear graphene response, - the property which was not available in conventional (three-dimensional) nonlinear optical materials. In this paper we give a short overview of the state of the art in this rapidly developing area of physics and present some of our recent results in this field. In particular, we discuss essential differences between the nonlinear response functions and parameters of traditional three-dimensional and new two-dimensional crystals.


abstract: Data integration, or the strategic analysis of multiple sources of data simultaneously, can often lead to discoveries that may be hidden in individualistic analyses of a single data source. We develop a new unsupervised data integration method named Integrated Principal Components Analysis (iPCA), which is a model-based generalization of PCA and serves as a practical tool to find and visualize common patterns that occur in multiple data sets. The key idea driving iPCA is the matrix-variate normal model, whose Kronecker product covariance structure captures both individual patterns within each data set and joint patterns shared by multiple data sets. Building upon this model, we develop several penalized (sparse and non-sparse) covariance estimators for iPCA, and using geodesic convexity, we prove that our non-sparse iPCA estimator converges to the global solution of a non-convex problem. We also demonstrate the practical advantages of iPCA through extensive simulations and a case study application to integrative genomics for Alzheimer's disease. In particular, we show that the joint patterns extracted via iPCA are highly predictive of a patient's cognition and Alzheimer's diagnosis.


title: Monotone Lagrangians in $\mathbb{CP}^n$ of minimal Maslov number $n+1$

abstract: We show that a monotone Lagrangian $L$ in $\mathbb{CP}^n$ of minimal Maslov number $n + 1$ is homeomorphic to a double quotient of a sphere, and thus homotopy equivalent to $\mathbb{RP}^n$. To prove this we use Zapolsky's canonical pearl complex for $L$ with coefficients in $\mathbb{Z}$, and various twisted versions thereof, where the twisting is determined by connected covers of $L$. The main tool is the action of the quantum cohomology of $\mathbb{CP}^n$ on the resulting Floer homologies.


title: Controlling a Quantum System via its Boundary Conditions

abstract: We numerically study a particle in a box with moving walls. In the case where the walls are oscillating sinusoidally with small amplitude, we show that states up to the fourth state can be populated with more than 80 percent population, while higher-lying states can also be selectively excited. This work introduces a way of controlling quantum systems which does not rely on (dipole) selection rules.


title: Automated Discovery of Jet Substructure Analyses

abstract: The study of the substructure of collimated particles from quarks and gluons, or jets, has the promise to reveal the details how color charges interact with the QCD plasma medium created in colliders such as RHIC and the LHC. Traditional jet substructure observables have been constructed using expert knowledge, and are largely transplanted, unmodified, from the high-energy physics, where the goal is primarily the study of boosted hadronic decays. A novel neural network architecture is described that is capable of examining theoretical models, and constructs, on its own, an analysis procedure that is sensitive to the internal model features. This architecture, in combination with symbolic regression, further allows the extraction of closed-form algebraic expressions from the learned result -- enabling the automatically constructed jet substructure analysis to be subsequently understood and reproduced by humans. This system is then tasked to construct an analysis that infers the plasma temperature from observing jets, which is demonstrated using both JEWEL and the Linearized Boltzmann Transport model, and at the presence of a realistic remnant of the plasma, or underlying event, that the measurement has to overcome. In a demonstration how algorithms can produce original research in direct competition to human experts, the resulting jet substructure variables and analyses are capable of determining the initial temperature of the plasma medium from analyzing 1200--2500 jets, a performance not seen in existing, manually designed analyses. Comparison of an incidentally discovered observable with the existing literature further indicates that the system described is capable of examining the model phase spaces to a detail at least comparable to the current field of human experts.


title: Input from the charm threshold for the measurement of $\gamma$

abstract: A brief overview of the inputs from charm threshold that are essential to the determination of the Unitarity Triangle angle $\gamma$ is presented. The focus is on the measurements of four-body final states that have not previously been considered: $D^0\to K_{\rm S}^0 \pi^+ \pi^- \pi^0 $ and $D^0 \to \pi^+ \pi^- \pi^+ \pi^-$.


title: Graphene for Emission of Terahertz Radiation

abstract: Graphene is a two-dimensional crystal consisting of a monatomic layer of carbon atoms. Electrons and holes in graphene behave as quasi-relativistic particles with zero effective mass and large (as compared to semiconductors) Fermi velocity. These unique physical properties can be used for designing detectors, emitters and modulators of terahertz radiation. Here we discuss several ideas of using graphene based structures for emission of terahertz radiation.


title: Multimodal Interactive Learning of Primitive Actions

abstract: We describe an ongoing project in learning to perform primitive actions from demonstrations using an interactive interface. In our previous work, we have used demonstrations captured from humans performing actions as training samples for a neural network-based trajectory model of actions to be performed by a computational agent in novel setups. We found that our original framework had some limitations that we hope to overcome by incorporating communication between the human and the computational agent, using the interaction between them to fine-tune the model learned by the machine. We propose a framework that uses multimodal human-computer interaction to teach action concepts to machines, making use of both live demonstration and communication through natural language, as two distinct teaching modalities, while requiring few training samples.


title: Network Modeling and Pathway Inference from Incomplete Data ("PathInf")

abstract: In this work, we developed a network inference method from incomplete data ("PathInf") , as massive and non-uniformly distributed missing values is a common challenge in practical problems. PathInf is a two-stages inference model. In the first stage, it applies a data summarization model based on maximum likelihood to deal with the massive distributed missing values by transforming the observation-wise items in the data into state matrix. In the second stage, transition pattern (i.e. pathway) among variables is inferred as a graph inference problem solved by greedy algorithm with constraints. The proposed method was validated and compared with the state-of-art Bayesian network method on the simulation data, and shown consistently superior performance. By applying the PathInf on the lymph vascular metastasis data, we obtained the holistic pathways of the lymph node metastasis with novel discoveries on the jumping metastasis among nodes that are physically apart. The discovery indicates the possible presence of sentinel node groups in the lung lymph nodes which have been previously speculated yet never found. The pathway map can also improve the current dissection examination protocol for better individualized treatment planning, for higher diagnostic accuracy and reducing the patients trauma.


title: Universality of Toda equation in ${\cal N}=2$ superconformal field theories

abstract: We show that extremal correlators in all Lagrangian ${\cal N}=2$ superconformal field theories with a simple gauge group are governed by the same universal Toda system of equations, which dictates the structure of extremal correlators to all orders in the perturbation series. A key point is the construction of a convenient orthogonal basis for the chiral ring, by arranging towers of operators in order of increasing dimension, which has the property that the associated two-point functions satisfy decoupled Toda chain equations. We explicitly verify this in all known SCFTs based on $\mathrm{SU}(N)$ gauge groups as well as in superconformal QCD based on orthogonal and symplectic groups. As a by-product, we find a surprising non-renormalization property for the ${\cal N}=2$ $\mathrm{SU}(N)$ SCFT with one hypermultiplet in the rank-2 symmetric representation and one hypermultiplet in the rank-2 antisymmetric representation, where the two-loop terms of a large class of supersymmetric observables identically vanish.


title: Precision measurements of the CKM parameters (mainly $\gamma/\phi_{3}$ measurements)

abstract: The CKM angle $\gamma/\phi_{3}$ is the only one that is accessible with tree level decays in a theoretically clean way such that it provides a precision test of $CP$ violation in the standard model. The Belle II experiment is a substantial upgrade of the Belle detector and will operate at the SuperKEKB asymmetric-energy $e^{+}e^{-}$ collider. The accelerator has already successfully completed the first phase of commissioning in 2016 and first $e^{+}e^{-}$ collisions in Belle II happened during April 2018. The design luminosity of SuperKEKB is 8 $\times$ 10$^{35}$ cm$^{-2}$s$^{-1}$ and the Belle II experiment aims to record 50 ab$^{-1}$ of data, a factor of 50 more than its predecessor (Belle). The key method to measure $\phi_{3}$ is through interference between the $B^- \to D^0 K^-$ and $B^- \to \bar{D^{0}} K^{-}$ decays which occurs if the final state of the charm-meson decay is accessible to both the $D^0$ and $\bar{D^{0}}$ mesons. To achieve the best sensitivity, a large variety of $D$ and $B$ decay modes are required, which is possible at the Belle II experiment as almost any final state can be reconstructed including those with photons. With the ultimate Belle II data sample of 50 ab$^{-1}$, a determination of $\phi_{3}$ with a precision of $1^{\circ}$ or better is foreseen. We explain herein the details of the planned measurements at Belle II.


title: Generalized, energy-conserving numerical simulations of particles in general relativity. II. Test particles in electromagnetic fields and GRMHD

abstract: Direct observations of compact objects, in the form of radiation spectra, gravitational waves from VIRGO/LIGO, and forthcoming direct imaging, are currently one of the primary source of information on the physics of plasmas in extreme astrophysical environments. The modeling of such physical phenomena requires numerical methods that allow for the simulation of microscopic plasma dynamics in presence of both strong gravity and electromagnetic fields. In Bacchini et al. (2018) we presented a detailed study on numerical techniques for the integration of free geodesic motion. Here we extend the study by introducing electromagnetic forces in the simulation of charged particles in curved spacetimes. We extend the Hamiltonian energy-conserving method presented in Bacchini et al. (2018) to include the Lorentz force and we test its performance compared to that of standard explicit Runge-Kutta and implicit midpoint rule schemes against analytic solutions. Then, we show the application of the numerical schemes to the integration of test particle trajectories in general relativistic magnetohydrodynamic (GRMHD) simulations, by modifying the algorithms to handle grid-based electromagnetic fields. We test this approach by simulating ensembles of charged particles in a static GRMHD configuration obtained with the Black Hole Accretion Code (BHAC).


title: Nearly all the sky is covered by Lyman-alpha emission around high redshift galaxies

abstract: Galaxies are surrounded by large reservoirs of gas, mostly hydrogen, fed by inflows from the intergalactic medium and by outflows due to galactic winds. Absorption-line measurements along the sightlines to bright and rare background quasars indicate that this circumgalactic medium pervades far beyond the extent of starlight in galaxies, but very little is known about the spatial distribution of this gas. A new window into circumgalactic environments was recently opened with the discovery of ubiquitous extended Lyman-alpha emission from hydrogen around high-redshift galaxies, facilitated by the extraordinary sensitivity of the MUSE instrument at the ESO Very Large Telescope. Due to the faintness of this emission, such measurements were previously limited to especially favourable systems or to massive statistical averaging. Here we demonstrate that low surface brightness Lyman-alpha emission surrounding faint galaxies at redshifts between 3 and 6 adds up to a projected sky coverage of nearly 100%. The corresponding rate of incidence (the mean number of Lyman-alpha emitters penetrated by any arbitrary line of sight) is well above unity and similar to the incidence rate of high column density absorbers frequently detected in the spectra of distant quasars. This similarity suggests that most circumgalactic atomic hydrogen at these redshifts has now been detected also in emission.


title: A New SU(2) Anomaly

abstract: A familiar anomaly affects SU(2) gauge theory in four dimensions: a theory with an odd number of fermion multiplets in the spin 1/2 representation of the gauge group, and more generally in representations of spin 2r+1/2, is inconsistent. We describe here a more subtle anomaly that can affect SU(2) gauge theory in four dimensions under the condition that fermions transform with half-integer spin under SU(2) and bosons with integer spin. Such a theory, formulated in a way that requires no choice of spin structure, and with an odd number of fermion multiplets in representations of spin 4r+3/2, is inconsistent. The theory is consistent if one picks a spin or spin_c structure. Under Higgsing to U(1), the new SU(2) anomaly reduces to a known anomaly of "all-fermion electrodynamics." Like that theory, an SU(2) theory with an odd number of fermion multiplets in representations of spin 4r+3/2 can provide a boundary state for a five-dimensional gapped theory whose partition function on a closed five-manifold Y is $(-1)^{\int_Y w_2w_3}$. All statements have analogs with SU(2) replaced by Sp(2N). There is also an analog in five dimensions.


title: CHET: Compiler and Runtime for Homomorphic Evaluation of Tensor Programs

abstract: Fully Homomorphic Encryption (FHE) refers to a set of encryption schemes that allow computations to be applied directly on encrypted data without requiring a secret key. This enables novel application scenarios where a client can safely offload storage and computation to a third-party cloud provider without having to trust the software and the hardware vendors with the decryption keys. Recent advances in both FHE schemes and implementations have moved such applications from theoretical possibilities into the realm of practicalities. This paper proposes a compact and well-reasoned interface called the Homomorphic Instruction Set Architecture (HISA) for developing FHE applications. Just as the hardware ISA interface enabled hardware advances to proceed independent of software advances in the compiler and language runtimes, HISA decouples compiler optimizations and runtimes for supporting FHE applications from advancements in the underlying FHE schemes. This paper demonstrates the capabilities of HISA by building an end-to-end software stack for evaluating neural network models on encrypted data. Our stack includes an end-to-end compiler, runtime, and a set of optimizations. Our approach shows generated code, on a set of popular neural network architectures, is faster than hand-optimized implementations.


title: Classification from Positive, Unlabeled and Biased Negative Data

abstract: In binary classification, there are situations where negative (N) data are too diverse to be fully labeled and we often resort to positive-unlabeled (PU) learning in these scenarios. However, collecting a non-representative N set that contains only a small portion of all possible N data can often be much easier in practice. This paper studies a novel classification framework which incorporates such biased N (bN) data in PU learning. We provide a method based on empirical risk minimization to address this PUbN classification problem. Our approach can be regarded as a novel example-weighting algorithm, with the weight of each example computed through a preliminary step that draws inspiration from PU learning. We also derive an estimation error bound for the proposed method. Experimental results demonstrate the effectiveness of our algorithm in not only PUbN learning scenarios but also ordinary PU learning scenarios on several benchmark datasets.


title: Exploring the Performance Boundaries of NB-IoT

abstract: NarrowBand-IoT has just joined the LPWAN community. Unlike most of its competitors, NB-IoT did not emerge from a blank slate. Indeed, it is closely linked to LTE, from which it inherits many of the features that undoubtedly determine its behavior. In this paper, we empirically explore the boundaries of this technology, analyzing from a user's point of view critical characteristics such as energy consumption, reliability and delays. The results show that its performance in terms of energy is comparable and even outperforms, in some cases, an LPWAN reference technology like LoRa, with the added benefit of guaranteeing delivery. However, the high variability observed in both energy expenditure and network delays call into question its suitability for some applications, especially those subject to service-level agreements.


title: N-, P-, As-triphenylene-graphdiyne: Strong and stable 2D semiconductors with outstanding capacities as anodes for Li-ion batteries

abstract: Since the first report of graphdiyne nanomembranes synthesis in 2010, different novel graphdiyne nanosheets have been fabricated. In a latest experimental advance, triphenylene-graphdiyne (TpG), a novel two-dimensional (2D) material was fabricated using the liquid/liquid interfacial method. In this study, we employed extensive first-principles simulations to investigate the mechanical/failure, thermal stability, electronic and optical properties of single-layer TpG. In addition, we predicted and explored the properties of nitrogenated-, phosphorated- and arsenicated-TpG monolayers. Our results reveal that TpG, N-TpG, P-TpG and As-TpG nanosheets can exhibit outstanding thermal stability. These nanomembranes moreover were found to yield linear elasticity with considerable tensile strengths. Notably, it was predicted that monolayer TpG, As-TpG, P-TpG and N-TpG show semiconducting electronic characters with direct band-gaps of 1.94 eV, 0.88 eV, 1.54 eV and 1.91 eV, respectively, along with highly attractive optical properties. We particularly analyzed the application prospect of these novel 2D materials as anodes for Li-ion batteries. Remarkably, P-TpG and N-TpG nanosheets were predicted to yield ultrahigh charge capacities of 1979 mAh/g and 2664 mAh/g, respectively, for Li-ions storage. The acquired results by this work suggest TpG based nanomembranes as highly promising candidates for the design of flexible nanoelectronics and energy storage devices.


title: LEGEND: The Large Enriched Germanium Experiment for Neutrinoless Double-Beta Decay

abstract: The lepton number violating process of neutrinoless double-beta decay could result from the physics beyond the Standard Model needed to generate the neutrino masses. Taking different approaches, the current generation of $^{76}$Ge experiments, the MAJORANA DEMONSTRATOR and GERDA, lead the field in both the ultra-low background and energy resolution achieved. The next generation of neutrinoless double-beta decay experiments requires increased mass and further reduction of backgrounds to maximize discovery potential. Building on the successes of the MAJORANA DEMONSTRATOR and GERDA, the LEGEND collaboration has been formed to pursue a tonne-scale $^{76}$Ge experiment, with discovery potential at a half-life beyond $10^{28}$ years. The collaboration aims to develop a phased neutrinoless double-beta decay experimental program, starting with a 200 kg measurement using the existing GERDA cryostat at LNGS. These proceedings discuss the plans and physics reach of LEGEND, and the combination of R&D efforts and existing resources being employed to expedite physics results.


title: Augmented Mitotic Cell Count using Field Of Interest Proposal

abstract: Histopathological prognostication of neoplasia including most tumor grading systems are based upon a number of criteria. Probably the most important is the number of mitotic figures which are most commonly determined as the mitotic count (MC), i.e. number of mitotic figures within 10 consecutive high power fields. Often the area with the highest mitotic activity is to be selected for the MC. However, since mitotic activity is not known in advance, an arbitrary choice of this region is considered one important cause for high variability in the prognostication and grading. In this work, we present an algorithmic approach that first calculates a mitotic cell map based upon a deep convolutional network. This map is in a second step used to construct a mitotic activity estimate. Lastly, we select the image segment representing the size of ten high power fields with the overall highest mitotic activity as a region proposal for an expert MC determination. We evaluate the approach using a dataset of 32 completely annotated whole slide images, where 22 were used for training of the network and 10 for test. We find a correlation of r=0.936 in mitotic count estimate.


title: Solvability for a drift-diffusion system with Robin boundary conditions

abstract: This paper focuses on a drift-diffusion system subjected to boundedly non dissipative Robin boundary conditions. A general existence result with large initial conditions is established by using suitable L1, L2 and trace estimates. Finally, two examples coming from the corrosion and the self-gravitation model are analyzed.


title: Raster Grid Pathology and the Cure

abstract: Blind ptychography is a phase retrieval method using multiple coded diffraction patterns from different, overlapping parts of the unknown extended object illuminated with an unknown window function. The window function is also known as the probe in the optics literature. As such blind ptychography is an inverse problem of simultaneous recovery of the object and the window function given the intensities of the windowed Fourier transform and has a multi-scale set-up in which the probe has an intermediate scale between the pixel scale and the macro-scale of the extended object. Uniqueness problem for blind ptychography is analyzed rigorously for the raster scan (of a constant step size {\tau}) and its variants, in which another scale comes into play: the overlap between adjacent blocks (the shifted windows). The block phases are shown to form an arithmetic progression and the complete characterization of the raster scan ambiguities is given, including: First, the periodic raster grid pathology of degrees of freedom proportional to {\tau}^2 and, second, a non-periodic, arithmetically progressing phase shift from block to block. Finally irregularly perturbed raster scans are shown to remove all ambiguities other than the inherent ambiguities of the scaling factor and the affine phase ambiguity under the minimum requirement of roughly 50% overlap ratio.


title: Simulations of $ SU(2) $ lattice gauge theory with dynamical reduced staggered fermions

abstract: We simulate $ SU(2) $ lattice gauge theory using dynamical reduced staggered fermions. The latter lead to two rather than four Dirac fermions in the continuum limit. We review the derivation and properties of reduced staggered fermions and show that in the case of fields in the fundamental representation of $SU(2)$ the theory does not exhibit a sign problem and can be simulated using the RHMC algorithm. We present results on lattices up to $16^4$ for a wide range of bare fermion masses. We find a single site condensate appears at strong coupling that spontaneously breaks the one global $U(1)$ symmetry remaining in the reduced fermion action.


title: Black hole remnants may exist if Starobinsky inflation occurred

abstract: Zero temperature black hole solutions to the semiclassical backreaction equations are investigated. Evidence is provided that certain components of the stress-energy tensors for free quantum fields at the horizon only depend on the local geometry near the horizon. This allows the semiclassical backreaction equations to be solved near the horizon. It is found that macroscopic uncharged zero temperature black hole solutions to the equations may exist if the coefficient of one of the higher derivative terms in the gravitational Lagrangian is large enough and of the right sign for Starobinsky inflation to have occurred in the early Universe.


title: Ultra-Wideband Detection of 22 Coherent Radio Bursts on M Dwarfs

abstract: Coherent radio bursts detected from M dwarfs have some analogy with solar radio bursts, but reach orders of magnitude higher luminosities. These events trace particle acceleration, powered by magnetic reconnection, shock fronts (such as formed by coronal mass ejections, CMEs), and magnetospheric currents, in some cases offering the only window into these processes in stellar atmospheres. We conducted a 58-hour, ultra-wideband survey for coherent radio bursts on 5 active M dwarfs. We used the Karl G. Jansky Very Large Array (VLA) to observe simultaneously in three frequency bands covering a subset of 224-482 MHz and 1-6 GHz, achieving the widest fractional bandwidth to date for any observations of stellar radio bursts. We detected 22 bursts across 13 epochs, providing the first large sample of wideband dynamic spectra of stellar coherent radio bursts. The observed bursts have diverse morphology, with durations ranging from seconds to hours, but all share strong (40-100%) circular polarization. No events resemble solar Type II bursts (often associated with CMEs), but we cannot rule out the occurrence of radio-quiet stellar CMEs. The hours-long bursts are all polarized in the sense of the x-mode of the star's large-scale magnetic field, suggesting they are cyclotron maser emission from electrons accelerated in the large-scale field, analogous to auroral processes on ultracool dwarfs. The duty cycle of luminous coherent bursts peaks at 25% at 1-1.4 GHz, declining at lower and higher frequencies, indicating source regions in the low corona. At these frequencies, active M dwarfs should be the most common galactic transient source.


title: Equivariant prequantisation of the super-0-brane in ${\rm AdS}_2\times\mathbb{S}^2$ - a toy model for supergerbe theory on curved spaces

abstract: The paper is another step towards a realisation of the goal, advanced in articles 1706.05682 [hep-th] and 1808.04470 [hep-th], of a systematic supersymmetry-equivariant geometrisation of physically distinguished Green-Schwarz super-(p+2)-cocycles defining classes in the supersymmetry-invariant refinement of the de Rham cohomology of homogeneous spaces of (supersymmetry) Lie supergroups, associated with reductive decompositions of their Lie superalgebras. It deals with a correlated geometrisation of a pair of such super-(p+2)-cocycles on spaces in correspondence under a blow-up transformation dual to the \.In\"on\"u-Wigner contraction that relates the respective (supersymmetry) Lie superalgebras, the latter correspondence being taken as the organising principle of the geometrisation procedure that exploits the link between the Cartan-Eilenberg cohomology of the supersymmetry group and the Chevalley-Eilenberg cohomology of its Lie superalgebra, alongside a cohomological classification of central extensions thereof. A general scheme of a correlated geometrisation compatible with the contraction is laid out and illustrated on the nontrivial example of a pair of consistent super-0-brane backgrounds: the super-Minkowski space $sMink^{3,1|8}$ with the standard N=2 Green-Schwarz super-2-cocycle on it and the super-$AdS_2 \times S^2$ space with Zhou's super-2-cocycle on it, asymptoting to the former in the limit of an infinite common radius of the $AdS_2$ and the $S^2$ in the body $AdS_2\times S^2$ of the supertarget. The geometrisation yields the respective supersymmetry-equivariant super-0-gerbes, i.e., the prequantum bundles of the associated Green-Schwarz super-$\sigma$-models in the Nambu-Goto formulation. Upon passing to the equivalent Hughes-Polchinski formulation of the models, the relevant extended super-0-gerbes are verified to possess a weak $\kappa$-equivariant structure.


title: Measurement of the electron-antineutrino correlation in neutron beta decay: aCORN experiment

abstract: The aCORN experiment uses a novel asymmetry method to measure the electron-antineutrino correlation (a-coefficient) in free neutron decay that does not require precision proton spectroscopy. aCORN completed two physics runs at the NIST Center for Neutron Research. The first run on the NG-6 beam line in 2013--2014 obtained the result a = 0.1090 +/- 0.0030 (stat) +/- 0.0028 (sys), a total uncertainty of 3.8%. The second run on the new NG-C high flux beam line promises an improvement in precision to <2%.


title: Giant Intrinsic Spin Hall Effect in W$_3$Ta and other A15 Superconductors

abstract: The spin Hall effect (SHE) is the conversion of charge current to spin current, and non-magnetic metals with large SHEs are extremely sought after for spintronic applications, but their rarity has stifled widespread use. Here we predict and explain the large intrinsic SHE in $\beta$-W and the A15 family of superconductors: W$_3$Ta, Ta$_3$Sb, and Cr$_3$Ir having spin hall conductivities (SHC) of -2250, -1400, and 1210 $\frac{\hbar}{e}(\Omega cm)^{-1}$, respectively. Combining concepts from topological physics with the dependence of the SHE on the spin Berry curvature (SBC) of the electronic bands, we propose a simple strategy to rapidly search for materials with large intrinsic SHEs based on the following ideas: high symmetry combined with heavy atoms gives rise to multiple Dirac-like crossings in the electronic structure, without sufficient symmetry protection these crossings gap due to spin orbit coupling (SOC), and gapped Dirac crossings create large spin Berry curvature.


title: Dynamic Sparse Graph for Efficient Deep Learning

abstract: We propose to execute deep neural networks (DNNs) with dynamic and sparse graph (DSG) structure for compressive memory and accelerative execution during both training and inference. The great success of DNNs motivates the pursuing of lightweight models for the deployment onto embedded devices. However, most of the previous studies optimize for inference while neglect training or even complicate it. Training is far more intractable, since (i) the neurons dominate the memory cost rather than the weights in inference; (ii) the dynamic activation makes previous sparse acceleration via one-off optimization on fixed weight invalid; (iii) batch normalization (BN) is critical for maintaining accuracy while its activation reorganization damages the sparsity. To address these issues, DSG activates only a small amount of neurons with high selectivity at each iteration via a dimension-reduction search (DRS) and obtains the BN compatibility via a double-mask selection (DMS). Experiments show significant memory saving (1.7-4.5x) and operation reduction (2.3-4.4x) with little accuracy loss on various benchmarks.


title: CMS precision timing physics impact for the HL-LHC upgrade

abstract: As part of the Phase II upgrade program, the Compact Muon Solenoid (CMS) detector will incorporate a new timing layer designed to measure minimum ionizing particles (MIPs) with a time resolution of $\sim$30 ps. Precision timing will mitigate the impact of the challenging levels of pileup expected at the High Luminosity LHC. The time information assigned to each track will enable the use of 4D-vertexing which will render a 5-fold pile-up reduction, thus recovering the current conditions. Precision timing will also enable new time-based isolations and improved $b$-tagging algorithms. All of this translates into a $\sim20\%$ gain in effective luminosity when looking at di-Higgs boson events decaying to a pair of $b$-quarks and two photons. We present the expected improvements in physics performance with precision timing with the upgraded CMS detector.


title: ProxQuant: Quantized Neural Networks via Proximal Operators

abstract: To make deep neural networks feasible in resource-constrained environments (such as mobile devices), it is beneficial to quantize models by using low-precision weights. One common technique for quantizing neural networks is the straight-through gradient method, which enables back-propagation through the quantization mapping. Despite its empirical success, little is understood about why the straight-through gradient method works. Building upon a novel observation that the straight-through gradient method is in fact identical to the well-known Nesterov's dual-averaging algorithm on a quantization constrained optimization problem, we propose a more principled alternative approach, called ProxQuant, that formulates quantized network training as a regularized learning problem instead and optimizes it via the prox-gradient method. ProxQuant does back-propagation on the underlying full-precision vector and applies an efficient prox-operator in between stochastic gradient steps to encourage quantizedness. For quantizing ResNets and LSTMs, ProxQuant outperforms state-of-the-art results on binary quantization and is on par with state-of-the-art on multi-bit quantization. For binary quantization, our analysis shows both theoretically and experimentally that ProxQuant is more stable than the straight-through gradient method (i.e. BinaryConnect), challenging the indispensability of the straight-through gradient method and providing a powerful alternative.


title: Neutron Lifetime Discrepancy as a Sign of a Dark Sector?

abstract: We summarize our recent proposal of explaining the discrepancy between the bottle and beam measurements of the neutron lifetime through the existence of a dark sector, which the neutron can decay to with a branching fraction 1%. We show that viable particle physics models for such neutron dark decays can be constructed and we briefly comment on recent developments in this area.


title: Convergence rates for quantum evolution & entropic continuity bounds in infinite dimensions

abstract: By extending the concept of energy-constrained diamond norms, we obtain continuity bounds on the dynamics of both closed and open quantum systems in infinite-dimensions, which are stronger than previously known bounds. We extensively discuss applications of our theory to quantum speed limits, attenuator and amplifier channels, the quantum Boltzmann equation, and quantum Brownian motion. Next, we obtain explicit log-Lipschitz continuity bounds for entropies of infinite-dimensional quantum systems, and classical capacities of infinite-dimensional quantum channels under energy-constraints. These bounds are determined by the high energy spectrum of the underlying Hamiltonian and can be evaluated using Weyl's law.


title: A unified theory of variational and quantum Monte Carlo methods and beyond

abstract: We present a unified theory of the variational Monte Carlo (VMC) and determinant quantum Monte Carlo (DQMC) methods using a novel density matrix formulation of VMC. We introduce an efficient algorithm for VMC to compute correlation functions and expectation values based on the auxiliary field Hirsch-Hubbard-Stratonovic transformation. We show that this new approach to VMC converges significantly faster than its traditional implementations. Furthermore, we generalize the Trotter-Suzuki decomposition to finite imaginary time steps $\tau \sim O(1)$ and develop a variational quantum Monte Carlo (VQMC) method accordingly, which is more accurate than VMC and can incorporate quantum fluctuations more efficiently. The two extreme limits of the VQMC method, namely infinitesimal and infinite imaginary time steps, correspond to the DQMC and VMC techniques, respectively. We demonstrate that our VQMC allows us to access lower temperatures in comparison with the conventional DQMC before the sign problem comes into play. We finally show that our VQMC can also enhance the accuracy of the projector Monte Carlo methods by providing better and less biased candidates for the trial wave functions, requiring shorter projection times for a given accuracy and alleviating the sign problem further.


title: Perturbation Gadgets: Arbitrary Energy Scales from a Single Strong Interaction

abstract: In this work we propose a many-body Hamiltonian construction which introduces only a single separate energy scale of order $\Theta(1/N^{2+\delta})$, for a small parameter $\delta>0$, and for $N$ terms in the target Hamiltonian. In its low-energy subspace, the construction can approximate any normalized target Hamiltonian $H_\mathrm{t}=\sum_{i=1}^N h_i$ with norm ratios $r=\max_{i,j\in\{1,\ldots,N\}}\|h_i\| / \| h_j \|=O(\exp(\exp(\mathrm{poly} n)))$ to within relative precision $O(N^{-\delta})$. This comes at the expense of increasing the locality by at most one, and adding an at most poly-sized ancilliary system for each coupling; interactions on the ancilliary system are geometrically local, and can be translationally-invariant. As an application, we discuss implications for QMA-hardness of the local Hamiltonian problem, and argue that "almost" translational invariance-defined as arbitrarily small relative variations of the strength of the local terms-is as good as non-translational-invariance in many of the constructions used throughout Hamiltonian complexity theory. We furthermore show that the choice of geared limit of many-body systems, where e.g. width and height of a lattice are taken to infinity in a specific relation, can have different complexity-theoretic implications: even for translationally-invariant models, changing the geared limit can vary the hardness of finding the ground state energy with respect to a given promise gap from computationally trivial, to QMAEXP-, or even BQEXPSPACE-complete.


title: A cavity-QED simulator of slow and fast scrambling

abstract: We study information scrambling, as diagnosed by the out-of-time order correlations (OTOCs), in a system of large spins collectively interacting via spatially inhomogeneous and incommensurate exchange couplings. The model is realisable in a cavity QED system in the dispersive regime. Fast scrambling, signalled by an exponential growth of the OTOCs, is observed when the couplings do not factorise into the product of a pair of local interaction terms, and at the same time the state of the spins points initially coplanar to the equator of the Bloch sphere. When one of these conditions is not realised, OTOCs grow algebraically with an exponent sensitive to the orientation of the spins in the initial state. The impact of initial conditions on the scrambling dynamics is attributed to the presence of a global conserved quantity, which critically slows down the evolution for initial states close to the poles of the Bloch sphere.


title: Domain-Adversarial Multi-Task Framework for Novel Therapeutic Property Prediction of Compounds

abstract: With the rapid development of high-throughput technologies, parallel acquisition of large-scale drug-informatics data provides huge opportunities to improve pharmaceutical research and development. One significant application is the purpose prediction of small molecule compounds, aiming to specify therapeutic properties of extensive purpose-unknown compounds and to repurpose novel therapeutic properties of FDA-approved drugs. Such problem is very challenging since compound attributes contain heterogeneous data with various feature patterns such as drug fingerprint, drug physicochemical property, drug perturbation gene expression. Moreover, there is complex nonlinear dependency among heterogeneous data. In this paper, we propose a novel domain-adversarial multi-task framework for integrating shared knowledge from multiple domains. The framework utilizes the adversarial strategy to effectively learn target representations and models their nonlinear dependency. Experiments on two real-world datasets illustrate that the performance of our approach obtains an obvious improvement over competitive baselines. The novel therapeutic properties of purpose-unknown compounds we predicted are mostly reported or brought to the clinics. Furthermore, our framework can integrate various attributes beyond the three domains examined here and can be applied in the industry for screening the purpose of huge amounts of as yet unidentified compounds. Source codes of this paper are available on Github.


title: Truly Concurrent Process Algebra Is Reversible

abstract: Based on our previous process algebra for concurrency APTC, we prove that it is reversible with a little modifications. The reversible algebra has four parts: Basic Algebra for Reversible True Concurrency (BARTC), Algebra for Parallelism in Reversible True Concurrency (APRTC), recursion and abstraction.


title: Training Machine Learning Models by Regularizing their Explanations

abstract: Neural networks are among the most accurate supervised learning methods in use today. However, their opacity makes them difficult to trust in critical applications, especially when conditions in training may differ from those in practice. Recent efforts to develop explanations for neural networks and machine learning models more generally have produced tools to shed light on the implicit rules behind predictions. These tools can help us identify when models are right for the wrong reasons. However, they do not always scale to explaining predictions for entire datasets, are not always at the right level of abstraction, and most importantly cannot correct the problems they reveal. In this thesis, we explore the possibility of training machine learning models (with a particular focus on neural networks) using explanations themselves. We consider approaches where models are penalized not only for making incorrect predictions but also for providing explanations that are either inconsistent with domain knowledge or overly complex. These methods let us train models which can not only provide more interpretable rationales for their predictions but also generalize better when training data is confounded or meaningfully different from test data (even adversarially so).


title: Truncated Dynamics, Ring Molecules and Mechanical Time Crystals

abstract: In applications of mechanics, including quantum mechanics, we often consider complex systems, where complete solutions of the underlying "fundamental" equations is both impractical and unnecessary to describe appropriate observations accurately. For example, practical chemistry, including even precision first-principles quantum chemistry, is never concerned with the behavior of the subnuclear quarks and gluons. Instead, we often focus on a few key variables, and construct a so-called effective theory for those. Such effective theories can become complicated and non-local, even for fairly simple systems. But in many circumstances, when there is a separation of scales, we can treat the reduced set of variables as a conventional dynamical system in its own right, governed by an energy conserving Lagrangian or Hamiltonian, in a useful approximation. The structure of that emergent description can display qualitatively new features, notably including reduced dimensionality, manifested through unconventional Poisson brackets. Here we discuss the physical meaning and consequences of such truncated dynamics. We propose physically realizable toy models of molecular rings, wherein time crystals emerge at the classical level. We propose that such behavior occurs in the effective theory of highly diamagnetic aromatic ring molecules, and could be widespread.


title: Automatic Skin Lesion Segmentation Using GrabCut in HSV Colour Space

abstract: Skin lesion segmentation is one of the first steps towards automatic Computer-Aided Diagnosis of skin cancer. Vast variety in the appearance of the skin lesion makes this task very challenging. The contribution of this paper is to apply a power foreground extraction technique called GrabCut for automatic skin lesion segmentation with minimal human interaction in HSV color space. Preprocessing was performed for removing the outer black border. Jaccard Index was measured to evaluate the performance of the segmentation method. On average, 0.71 Jaccard Index was achieved on 1000 images from ISIC challenge 2017 Training Dataset.


title: Quantum Gravity and background field formalism

abstract: We analyze the problem of general covariance for quantum gravity theories in the background field formalism with respect to gauge fixing procedure. We prove that the background effective action is not invariant under general coordinate transformations of background metric tensor in non-linear gauges.


title: Compiling Stan to Generative Probabilistic Languages and Extension to Deep Probabilistic Programming

abstract: Stan is a probabilistic programming language that is popular in the statistics community, with a high-level syntax for expressing probabilistic models. Stan differs by nature from generative probabilistic programming languages like Church, Anglican, or Pyro. This paper presents a comprehensive compilation scheme to compile any Stan model to a generative language and proves its correctness. We use our compilation scheme to build two new backends for the Stanc3 compiler targeting Pyro and NumPyro. Experimental results show that the NumPyro backend yields a 2.3x speedup compared to Stan in geometric mean over 26 benchmarks. Building on Pyro we extend Stan with support for explicit variational inference guides and deep probabilistic models. That way, users familiar with Stan get access to new features without having to learn a fundamentally new language.


title: Experimental demonstration of the violations of Mermin's and Svetlichny's inequalities for W- and GHZ-class of states

abstract: Violation of Mermin's and Svetlichny's inequalities can rule out the predictions of local hidden variable theory and can confirm the existence of true nonlocal correlation for n-particle pure quantum systems. Here we demonstrate the experimental violation of the above inequalities for W- and GHZ-class of states. We use IBM's five-qubit quantum computer for experimental implementation of these states and illustration of inequalities' violations. Our results clearly show the violations of both Mermin's and Svetlichny's inequalities for W and GHZ states respectively. Being a superconducting qubit-based quantum computer, the platform used here opens up the opportunity to explore multipartite inequalities which is beyond the reach of other existing technologies.


title: Solving 3SAT By Reduction To Testing For Odd Hole

abstract: An algorithm is given for finding the solutions to 3SAT problems. The algorithm uses Bienstock's reduction from 3SAT to existence of induced odd cycle of length greater than three, passing through a prescribed node in the constructed graph. The algorithm proceeds to find what will be called the hole complexes of the graph. The set of the boundary nodes of the hole complex containing the prescribed node is then searched for the subsets of 8 nodes corresponding to the 3SAT's literals. If a complete set of literals is contained in the boundary then the 3SAT is solvable.


title: Graph Isomorphism by Conversion to Chordal (6, 3) Graphs

abstract: Babel has shown that for an extended class of chordal (6, 3) graphs the coarsest regular simplicial partition is equivalent to the graph's automorphism partition. We give a reversible transformation for any graph to one of these graph by using Booth's reduction of a graph to a chordal graph and elimination of Babel's forbidden subgraphs for these graphs by adding edges to them.


title: Privacy and Utility Tradeoff in Approximate Differential Privacy

abstract: We characterize the minimum noise amplitude and power for noise-adding mechanisms in $(\epsilon, \delta)$-differential privacy for single real-valued query function. We derive new lower bounds using the duality of linear programming, and new upper bounds by proposing a new class of $(\epsilon,\delta)$-differentially private mechanisms, the \emph{truncated Laplacian} mechanisms. We show that the multiplicative gap of the lower bounds and upper bounds goes to zero in various high privacy regimes, proving the tightness of the lower and upper bounds and thus establishing the optimality of the truncated Laplacian mechanism. In particular, our results close the previous constant multiplicative gap in the discrete setting. Numeric experiments show the improvement of the truncated Laplacian mechanism over the optimal Gaussian mechanism in all privacy regimes.


title: VISION - Vienna survey in Orion. III. Young stellar objects in Orion A

abstract: We have extended and refined the existing young stellar object (YSO) catalogs for the Orion A molecular cloud, the closest massive star-forming region to Earth. This updated catalog is driven by the large spatial coverage (18.3 deg$^2$, $\sim$950 pc$^2$), seeing limited resolution ($\sim$0.7$"$), and sensitivity ($K_s<19$ mag) of the ESO-VISTA near-infrared survey of the Orion A cloud (VISION). Combined with archival mid- to far-infrared data, the VISTA data allow for a refined and more robust source selection. We estimate that among previously known protostars and pre-main-sequence stars with disks, source contamination levels (false positives) are at least $\sim$6.4% and $\sim$2.3%, respectively, mostly due to background galaxies and nebulosities. We identify 274 new YSO candidates using VISTA/Spitzer based selections within previously analyzed regions, and VISTA/WISE based selections to add sources in the surroundings, beyond previously analyzed regions. The WISE selection method recovers about 59% of the known YSOs in Orion A's low-mass star-forming part L1641, which shows what can be achieved by the all-sky WISE survey in combination with deep near-infrared data in regions without the influence of massive stars. The new catalog contains 2980 YSOs, which were classified based on the de-reddened mid-infrared spectral index into 188 protostars, 185 flat-spectrum sources, and 2607 pre-main-sequence stars with circumstellar disks. We find a statistically significant difference in the spatial distribution of the three evolutionary classes with respect to regions of high dust column-density, confirming that flat-spectrum sources are at a younger evolutionary phase compared to Class IIs, and are not a sub-sample seen at particular viewing angles.


title: The Universal Geometry of Heterotic Vacua

abstract: We consider a family of perturbative heterotic string backgrounds. These are complex threefolds X with c_1 = 0, each with a gauge field solving the Hermitian Yang-Mill's equations and compatible B and H fields that satisfy the anomaly cancellation conditions. Our perspective is to consider a geometry in which these backgrounds are fibred over a parameter space. If the manifold X has coordinates x, and parameters are denoted by y, then it is natural to consider coordinate transformations x \to \tilde{x}(x,y) and y \to \tilde{y}(y). Similarly, gauge transformations of the gauge field and B field also depend on both x and y. In the process of defining deformations of the background fields that are suitably covariant under these transformations, it turns out to be natural to extend the gauge field A to a gauge field \IA on the extended (x,y)-space. Similarly, the B, H, and other fields are also extended. The total space of the fibration of the heterotic structures is the Universal Geometry of the title. The extension of gauge fields has been studied in relation to Donaldson theory and monopole moduli spaces. String vacua furnish a richer application of these ideas. One advantage of this point of view is that previously disparate results are unified into a simple tensor formulation. In a previous paper, by three of the present authors, the metric on the moduli space of heterotic theories was derived, correct through order \alpha', and it was shown how this was related to a simple Kahler potential. With the present formalism, we are able to rederive the results of this previously long and involved calculation, in less than a page.


title: Baryogenesis and Dark Matter from $B$ Mesons

abstract: We present a new mechanism of Baryogenesis and dark matter production in which both the dark matter relic abundance and the baryon asymmetry arise from neutral $B$ meson oscillations and subsequent decays. This set-up is testable at hadron colliders and $B$-factories. In the early Universe, decays of a long lived particle produce $B$ mesons and anti-mesons out of thermal equilibrium. These mesons/anti-mesons then undergo CP violating oscillations before quickly decaying into visible and dark sector particles. Dark matter will be charged under Baryon number so that the visible sector baryon asymmetry is produced without violating the total baryon number of the Universe. The produced baryon asymmetry will be directly related to the leptonic charge asymmetry in neutral $B$ decays; an experimental observable. Dark matter is stabilized by an unbroken discrete symmetry, and proton decay is simply evaded by kinematics. We will illustrate this mechanism with a model that is unconstrained by di-nucleon decay, does not require a high reheat temperature, and would have unique experimental signals -- a positive leptonic asymmetry in $B$ meson decays, a new decay of $B$ mesons into a baryon and missing energy, and a new decay of $b$-flavored baryons into mesons and missing energy. These three observables are testable at current and upcoming collider experiments, allowing for a distinct probe of this mechanism.


title: Revealing a Highly-Dynamic Cluster Core in Abell 1664 with Chandra

abstract: We present new, deep (245 ks) Chandra observations of the galaxy cluster Abell 1664 ($z = 0.1283$). These images reveal rich structure, including elongation and accompanying compressions of the X-ray isophotes in the NE-SW direction, suggesting that the hot gas is sloshing in the gravitational potential. This sloshing has resulted in cold fronts, at distances of 55, 115 and 320 kpc from the cluster center. Our results indicate that the core of A1664 is highly disturbed, as the global metallicity and cooling time flatten at small radii, implying mixing on large scales. The central AGN appears to have recently undergone a mechanical outburst, as evidenced by our detection of cavities. These cavities are the X-ray manifestations of radio bubbles inflated by the AGN, and may explain the motion of cold molecular CO clouds previously observed with ALMA. The estimated mechanical power of the AGN, using the minimum energy required to inflate the cavities as a proxy, is $P_{\rm cav} = (1.1 \pm 1.0) \times 10^{44} $ erg s$^{-1}$, which may be enough to drive the molecular gas flows, and offset the cooling luminosity of the ICM, at $L_{\rm cool} = (1.90 \pm0.01)\times 10^{44}$ erg s$^{-1}$. This mechanical power is orders of magnitude higher than the measured upper limit on the X-ray luminosity of the central AGN, suggesting that its black hole may be extremely massive and/or radiatively inefficient. We map temperature variations on the same spatial scale as the molecular gas, and find that the most rapidly cooling gas is mostly coincident with the molecular gas reservoir centered on the BCG's systemic velocity observed with ALMA and may be fueling cold accretion onto the central black hole.


title: Evolution of the cold gas properties of simulated post-starburst galaxies

abstract: Post-starburst galaxies are typically considered to be a transition population, en route to the red sequence after a recent quenching event. Despite this, recent observations have shown that these objects typically have large reservoirs of cold molecular gas. In this paper we study the star-forming gas properties of a large sample of post-starburst galaxies selected from the cosmological, hydrodynamical EAGLE simulations. These objects resemble observed high-mass post-starburst galaxies both spectroscopically and in terms of their space density, stellar mass distribution and sizes. We find that the vast majority of simulated post-starburst galaxies have significant gas reservoirs, with star-forming gas masses of ~10$^9$ M$_{\odot}$, in good agreement with those seen in observational samples. The simulation reproduces the observed time evolution of the gas fraction of the post-starburst galaxy population, with the average galaxy losing ~90 per cent of its star-forming interstellar medium in only ~600 Myr. A variety of gas consumption/loss processes are responsible for this rapid evolution, including mergers and environmental effects, while active galactic nuclei play only a secondary role. The fast evolution in the gas fraction of post-starburst galaxies is accompanied by a clear decrease in the efficiency of star formation, due to a decrease in the dense gas fraction. We predict that forthcoming ALMA observations of the gas reservoirs of low-redshift post-starburst galaxies will show that the molecular gas is typically compact and has disturbed kinematics, reflecting the disruptive nature of many of the evolutionary pathways that build up the post-starburst galaxy population.


title: Bardeen-Petterson Alignment, Jets and Magnetic Truncation in GRMHD Simulations of Tilted Thin Accretion Discs

abstract: Prevalent around luminous accreting black holes, thin discs are challenging to resolve in numerical simulations. When the disc and black hole angular momentum vectors are misaligned, the challenge becomes extreme, requiring adaptive meshes to follow the disc proper as it moves through the computational grid. With our new high-performance general relativistic magnetohydrodynamic (GRMHD) code H-AMR we have simulated the thinnest accretion disc to date, of aspect ratio H/R~0.03, around a rapidly spinning (a=0.9375) black hole, using a cooling function. Initially tilted at 10 degrees, the disc warps inside r~5 r_g into alignment with the black hole, where r_g is the gravitational radius. This is the first demonstration of Bardeen-Petterson alignment in MHD with viscosity self-consistently generated by magnetized turbulence. The disc develops a low-density high-viscosity (alpha_eff ~ 1.0) magnetic-pressure--dominated inner region at r<25 r_g that rapidly empties itself into the black hole. This inner region may in reality, due to thermal decoupling of ions and electrons, evaporate into a radiatively inefficient accretion flow if, as we propose, the cooling time exceeds the accretion time set by the order unity effective viscosity. We furthermore find the unexpected result that even our very thin disc can sustain large-scale vertical magnetic flux on the black hole, which launches powerful relativistic jets that carry 20-50% of the accretion power along the angular momentum vector of the outer tilted disc, providing a potential explanation for the origin of jets in radio-loud quasars.


title: Pressure-Induced Metal-Insulator Transition in Twisted Bi-layer Graphene

abstract: Recent experiments [arXiv: 1808.07865] on twisted bilayer graphene (TBLG) show that under hydrostatic pressure, an insulating state at quarter-filling of the moir\'e superlattice (i.e., one charge per supercell) emerges, in sharp contrast with the previous ambient pressure measurements of Cao et al. where the quarter--filling state (QFS) is a metal [Nature 556, 43 & 80 (2018)]. In fact, the insulating state at the other commensurate fillings of two and three charges per supercell is also enhanced under applied pressure. Based on realistic computations of the band structure for TBLG which show that the bandwidth first shrinks and then expands with increasing hydrostatic pressure, we compute the ratio of the potential to the kinetic energy, $r_s$. We find an experimentally relevant window of pressure for which $r_s$ crosses the threshold for a triangular Wigner crystal, thereby corroborating our previous work [Nano Lett. (2018)] that the insulating states in TBLG are due to Wigner rather than Mott physics. A key prediction of this work is that the window for the onset of the hierarchy of Wigner states that obtains at commensurate fillings is dome-shaped as a function of the applied pressure, which can be probed experimentally. Theoretically, we find a peak for crystallization around $1.5$ GPa relative to the experimental optimal pressure of $1.33$ GPa for the observation of the insulating states. Consequently, TBLG provides a new platform for the exploration of Wigner physics and its relationship with superconductivity.


title: Broadband Intensity Tomography: Spectral Tagging of the Cosmic UV Background

abstract: Cosmic photons can be efficiently collected by broadband intensity mapping but information on their emission redshift and frequency is largely lost. We introduce a technique to statistically recover these otherwise collapsed dimensions by exploiting information in spatial fluctuations and apply it to the Galaxy Evolution Explorer (GALEX) All Sky and Medium Imaging Surveys. By spatially cross-correlating photons in the GALEX far-UV (1500A) and near-UV (2300A) bands with a million spectroscopic objects in the Sloan Digital Sky Survey as a function of redshift, we robustly detect the redshift-dependent intensity of the UV background (UVB) modulated by its clustering bias up to $z\sim2$. These measurements clearly reveal the imprints of UVB spectral features redshifting through the filters. Using a simple parameterization, we simultaneously fit a UVB emissivity and clustering bias factor to these observations and constrain the main spectral features of the UVB spectrum: (i) the Lyman break, (ii) the non-ionizing UV continuum, which agrees with the Haardt & Madau model but does not rely on any assumption regarding the nature of the sources, and (iii) the Ly$\alpha$ emission, whose luminosity density is consistent with estimates of the combined galaxy and AGN contributions at $z\sim1$. Because the technique probes the total background including low surface brightness emission, we place constraints on the amount of UV light originating from the diffuse intergalactic medium (IGM). Finally, the clustering bias of UV photons is found to be chromatic and evolving. Our frequency- and redshift-dependent UVB measurement delivers a summary statistic of the universe's net radiation output from stars, black holes, and the IGM combined.


title: Strong cosmic censorship: The nonlinear story

abstract: A satisfactory formulation of the laws of physics entails that the future evolution of a physical system should be determined from appropriate initial conditions. The existence of Cauchy horizons in solutions of the Einstein field equations is therefore problematic, and expected to be an unstable artifact of General Relativity. This is asserted by the Strong Cosmic Censorship Conjecture, which was recently put into question by an analysis of the linearized equations in the exterior of charged black holes in an expanding universe. Here, we numerically evolve the nonlinear Einstein-Maxwell-scalar field equations with a positive cosmological constant, under spherical symmetry, and provide strong evidence that mass inflation indeed does not occur in the near extremal regime. This shows that nonlinear effects might not suffice to save the Strong Cosmic Censorship Conjecture.