title: Composing Diverse Policies for Temporally Extended Tasks

abstract: Robot control policies for temporally extended and sequenced tasks are often characterized by discontinuous switches between different local dynamics. These change-points are often exploited in hierarchical motion planning to build approximate models and to facilitate the design of local, region-specific controllers. However, it becomes combinatorially challenging to implement such a pipeline for complex temporally extended tasks, especially when the sub-controllers work on different information streams, time scales and action spaces. In this paper, we introduce a method that can compose diverse policies comprising motion planning trajectories, dynamic motion primitives and neural network controllers. We introduce a global goal scoring estimator that uses local, per-motion primitive dynamics models and corresponding activation state-space sets to sequence diverse policies in a locally optimal fashion. We use expert demonstrations to convert what is typically viewed as a gradient-based learning process into a planning process without explicitly specifying pre- and post-conditions. We first illustrate the proposed framework using an MDP benchmark to showcase robustness to action and model dynamics mismatch, and then with a particularly complex physical gear assembly task, solved on a PR2 robot. We show that the proposed approach successfully discovers the optimal sequence of controllers and solves both tasks efficiently.


title: Origin of the Anomalous Electronic Shot Noise in Atomic-Scale Junctions

abstract: Fluctuations pose fundamental limitations in making sensitive measurements, yet at the same time, noise unravels properties that are inaccessible at the level of the averaged signal. In electronic devices, shot noise arises from the discrete nature of charge carriers and it increases linearly with the applied voltage according to the celebrated Schottky formula. Nonetheless, measurements of shot noise in atomic-scale junctions at high voltage reveal significant nonlinear (anomalous) behavior, which varies from sample to sample, and has no specific trend. Here, we provide a viable, unifying explanation for these diverse observations based on the theory of quantum coherent transport. Our formula for the anomalous shot noise relies on---and allows us to resolve---two key characteristics of a conducting junction: The structure of the transmission function at the vicinity of the Fermi energy and the asymmetry of the bias voltage drop at the contacts. We tested our theory on high voltage shot noise measurements on Au atomic scale junctions and demonstrated a quantitative agreement, recovering both the enhancement and suppression of shot noise as observed in different junctions. The good theory-experiment correspondence supports our modelling and emphasizes that the asymmetry of the bias drop on the contacts is a key factor in nanoscale electronic transport, which may substantially impact electronic signals even in incomplex structures.


title: Approximations for the boundary crossing probabilities of moving sums of random variables

abstract: In this paper we study approximations for the boundary crossing probabilities of moving sums of i.i.d. normal r.v. We approximate a discrete time problem with a continuous time problem allowing us to apply established theory for stationary Gaussian processes. By then subsequently correcting approximations for discrete time, we show that the developed approximations are very accurate even for small window length. Also, they have high accuracy when the original r.v. are not exactly normal and when the weights in the moving window are not all equal. We then provide accurate and simple approximations for ARL, the average run length until crossing the boundary.


title: Mathematical description of a Frozen Wave beam after passing through a pair of convex lenses with different focal distance

abstract: In this paper, we shall provide an analytical solution describing a Frozen Wave beam after passing through a pair of convex lenses with different focal distances.


title: The closure-complement-frontier problem in saturated polytopological spaces

abstract: Let $X$ be a space equipped with $n$ topologies $\tau_1,...,\tau_n$ which are pairwise comparable and saturated, and for each $1\leq i\leq n$ let $k_i$ and $f_i$ be the associated topological closure and frontier operators, respectively. Inspired by the closure-complement theorem of Kuratowski, we prove that the monoid of set operators $\mathcal{KF}_n$ generated by $\{k_i,f_i:1\leq i\leq n\}\cup\{c\}$ (where $c$ denotes the set complement operator) has cardinality no more than $2p(n)$ where $p(n)=\frac{5}{24}n^4+\frac{37}{12}n^3+\frac{79}{24}n^2+\frac{101}{12}n+2$. The bound is sharp in the following sense: for each $n$ there exists a saturated polytopological space $(X,\tau_1,...,\tau_n)$ and a subset $A\subseteq X$ such that repeated application of the operators $k_i, f_i, c$ to $A$ will yield exactly $2p(n)$ distinct sets. In particular, following the tradition for Kuratowski-type problems, we exhibit an explicit initial set in $\mathbb{R}$, equipped with the usual and Sorgenfrey topologies, which yields $2p(2)=120$ distinct sets under the action of the monoid $\mathcal{KF}_2$.


title: Topological theory of Lieb-Schultz-Mattis theorems in quantum spin systems

abstract: The Lieb-Schultz-Mattis (LSM) theorem states that a spin system with translation and spin rotation symmetry and half-integer spin per unit cell does not admit a gapped symmetric ground state lacking fractionalized excitations. That is, the ground state must be gapless, spontaneously break a symmetry, or be a gapped spin liquid. Thus, such systems are natural spin-liquid candidates if no ordering is found. In this work, we give a much more general criterion that determines when an LSM-type theorem holds in a spin system. For example, we consider quantum magnets with arbitrary space group symmetry and/or spin-orbit coupling. Our criterion is intimately connected to recent work on the general classification of topological phases with spatial symmetries and also allows for the computation of an "anomaly" associated with the existence of an LSM theorem. Moreover, our framework is also general enough to encompass recent works on "SPT-LSM" theorems where the system admits a gapped symmetric ground state without fractionalized excitations, but such a ground state must still be non-trivial in the sense of symmetry-protected topological (SPT) phases.

categories: cond-mat.str-el math-ph math.MP quant-ph


title: Tidal radii of main sequence stars -- I. Physical tidal radius, semi-analytic model and their implications

abstract: A star is tidally disrupted by a supermassive black hole when their separation is shorter than the "tidal radius". This quantity is often estimated on an order-of-magnitude basis without reference to the star's internal structure. Using MESA models for main sequence stars and fully general relativistic dynamics, we find the physical tidal radius for complete disruption $\cal{R}_t$ for a $10^6M_\odot$ black hole (BH). We find that across a factor $\sim20$ in stellar mass $M_*$, i.e., $0.15M_{\odot}\leq M_*\leq3M_\odot$, $\cal{R}_t\sim27\times$(BH's gravitational radius). When comparing $\cal{R}_t$ with the commonly used order-of-magnitude estimate $r_t$, we find that $\cal{R}_t\sim1.05-1.45r_t$ for $0.15M_\odot\leq M_*\leq0.5M_\odot$, but between $0.5 M_\odot$ and $1 M_\odot$, $\cal{R}_t$ drops to $\sim 0.45r_t$, and it remains at this value up to $10 M_\odot$. The near-constancy of $\cal{R}_t$ implies a weaker dependence of the full disruption rate on $M_*$ than when predicted with $r_t$. The characteristic energy width of the debris $\Delta E$ ranges from $\sim1.2\Delta\cal{E}$ for low-mass stars to $\sim 0.35\Delta\cal{E}$ for higher-mass stars, where $\Delta\cal{E}=GM_{\rm BH}R_*/\cal{R}_t^{2}$. We present analytic fits for the $M_*$ dependence of $\cal{R}_t$ and $\Delta E$; these fits lead to analytic expressions for the time of peak mass fallback rate and the maximal mass fallback rate. Our results also bear on the fraction of events leading to fast or slow circularization, as well as on the character of the tidal event occurring when the remnant of a partial disruption returns to the black hole. Using a semi-analytic model, we show that $\cal{R}_t$ is primarily determined by the star's central density rather than its mean density. For high-mass stars, the full disruption rate is roughly 1/4 the partial disruption rate, while this ratio is close to unity for low-mass stars.


title: What is wrong with steady accretion discs?

abstract: In a standard, steady, thin accretion disc, the radial distribution of the dissipation of the accretion energy is determined simply by energy considerations. Here we draw attention to the fact that while the (quasi-)steady discs in dwarf novae in outburst are in agreement with the expected emission distribution, the steady discs in the nova-like variables are not. We note that essentially the only difference between these two sets of discs is the time for which they have been in the high viscosity, high accretion rate state. In such discs, the major process by which angular momentum is transported outwards is MHD turbulence. We speculate that such turbulence gives rise to corona-like structures (here called magnetically controlled zones, or MCZs) which are also able to provide non-negligible angular momentum transport, the magnitude of which depends on the spatial scale $L$ of the magnetic field structures in such zones. For short-lived, high accretion rate discs (such as those in dwarf novae) we expect $L \sim H$ and the MCZ to have little effect. But, with time (such as in the nova-like variables) an inverse cascade in the MHD turbulence enables $L$, and the net effect of the MCZ, to grow. We present a simple toy model which demonstrates that such ideas can provide an explanation for the difference between the dwarf novae and the nova-like variable discs.


title: Conformal Perturbation Theory for Twisted Fields

abstract: We investigate second order conformal perturbation theory for $\mathbb{Z}_2$ orbifolds of conformal field theories in two dimensions. To evaluate the necessary twisted sector correlation functions and their integrals, we map them from the sphere to its torus double cover. We discuss how this relates crossing symmetry to the modular group, and introduce a regularization scheme on the cover that allows to evaluate the integrals numerically. These methods do not require supersymmetry. As an application, we show that in the torus orbifold of 8 and 16 free bosons, $\mathbb{Z}_2$ twist fields are marginal at first order, but stop being marginal at second order.


title: Destruction of long-range order in non-collinear two-dimensional antiferromagnets by random-bond disorder

abstract: We consider frustrated Heisenberg antiferromagnets, whose clean-limit ground state is characterized by non-collinear long-range order with non-zero vector chirality, and study the effects of quenched bond disorder, i.e., random exchange couplings. A single bond defect is known to induce a dipolar texture in the spin background independent of microscopic details. Using general analytical arguments as well as large-scale simulations for the classical triangular-lattice Heisenberg model, we show that any finite concentration of such defects destroys long-range order for spatial dimension $d\leq 2$, in favor of a glassy state whose correlation length in $d=2$ is exponentially large for small randomness. Our results are relevant for a wide range of layered frustrated magnets.


title: Neural Networks for Full Phase-space Reweighting and Parameter Tuning

abstract: Precise scientific analysis in collider-based particle physics is possible because of complex simulations that connect fundamental theories to observable quantities. The significant computational cost of these programs limits the scope, precision, and accuracy of Standard Model measurements and searches for new phenomena. We therefore introduce Deep neural networks using Classification for Tuning and Reweighting (DCTR), a neural network-based approach to reweight and fit simulations using all kinematic and flavor information -- the full phase space. DCTR can perform tasks that are currently not possible with existing methods, such as estimating non-perturbative fragmentation uncertainties. The core idea behind the new approach is to exploit powerful high-dimensional classifiers to reweight phase space as well as to identify the best parameters for describing data. Numerical examples from $e^+e^-\rightarrow\text{jets}$ demonstrate the fidelity of these methods for simulation parameters that have a big and broad impact on phase space as well as those that have a minimal and/or localized impact. The high fidelity of the full phase-space reweighting enables a new paradigm for simulations, parameter tuning, and model systematic uncertainties across particle physics and possibly beyond.


title: Modular Bosonic Subsystem Codes

abstract: We introduce a framework to decompose a bosonic mode into two virtual subsystems-a logical qubit and a gauge mode. This framework allows the entire toolkit of qubit-based quantum information to be applied in the continuous-variable setting. We give a detailed example based on a modular decomposition of the position basis and apply it in two situations. First, we decompose Gottesman-Kitaev-Preskill grid states and find that the encoded logical state can be damaged due to entanglement with the gauge mode. Second, we identify and disentangle qubit cluster states hidden inside of Gaussian continuous-variable cluster states.


title: Data Analysis for Precision 21 cm Cosmology

abstract: The redshifted 21 cm line is an emerging tool in cosmology, in principle permitting three-dimensional surveys of our Universe that reach unprecedentedly large volumes, previously inaccessible length scales, and hitherto unexplored epochs of our cosmic timeline. Large radio telescopes have been constructed for this purpose, and in recent years there has been considerable progress in transforming 21 cm cosmology from a field of considerable theoretical promise to one of observational reality. Increasingly, practitioners in the field are coming to the realization that the success of observational 21cm cosmology will hinge on software algorithms and analysis pipelines just as much as it does on careful hardware design and telescope construction. This review provides a pedagogical introduction to state-of-the-art ideas in 21 cm data analysis, covering a wide variety of steps in a typical analysis pipeline, from calibration to foreground subtraction to mapmaking to power spectrum estimation to parameter estimation.


title: Collapse and revival of quantum many-body scars via Floquet engineering

abstract: The presence of quantum scars, athermal eigenstates of a many-body Hamiltonian with finite energy density, leads to absence of ergodicity and long-time coherent dynamics in closed quantum systems starting from simple initial states. Such non-ergodic coherent dynamics, where the system does not explore its entire phase space, has been experimentally observed in a chain of ultracold Rydberg atoms. We show, via study of a periodically driven Rydberg chain, that the drive frequency acts as a tuning parameter for several reentrant transitions between ergodic and non-ergodic regimes. The former regime shows rapid thermalization of correlation functions and absence of scars in the spectrum of the system's Floquet Hamiltonian. The latter regime, in contrast, has scars in its Floquet spectrum which control the long-time coherent dynamics of correlation functions. Our results open a new possibility of drive frequency-induced tuning between ergodic and non-ergodic dynamics in experimentally realizable disorder-free quantum many-body systems.


title: A Catalog of Hyper-luminous X-ray Sources and Intermediate-Mass Black Hole Candidates out to High Redshifts

abstract: Hyper-luminous X-ray sources (HLXs; L_X>10^41 erg s^-1) are off-nuclear X-ray sources in galaxies and strong candidates for intermediate-mass black holes (IMBHs). We have constructed a sample of 169 HLX candidates by combining X-ray detections from the Chandra Source Catalog (Version 2) with galaxies from the Sloan Digital Sky Survey and registering individual images for improved relative astrometric accuracy. The spatial resolution of Chandra allows for the sample to extend out to z~0.9. Optical counterparts are detected among one-fourth of the sample, one-third of which are consistent with dwarf galaxy stellar masses. The average intrinsic X-ray spectral slope indicates efficient accretion, potentially driven by galaxy mergers, and the column densities suggest one-third of the sample has significant X-ray absorption. We find that 144 of the HLX candidates have X-ray emission that is significantly in excess of the expected contribution from star formation and hot gas, strongly suggesting that they are produced by accretion onto black holes more massive than stars. After correcting for an average background or foreground contamination rate of 8%, we estimate that at least ~20 of the HLX candidates are consistent with IMBH masses, and this estimate is potentially several times higher assuming more efficient accretion. This catalog currently represents the largest sample of uniformly-selected, off-nuclear IMBH candidates. These sources may represent scenarios in which a low-mass galaxy hosting an IMBH has merged with a more massive galaxy and provide an excellent sample for testing models of low-mass BH formation and merger-driven growth.


title: Dynamical formation of a magnetic polaron in a two-dimensional quantum antiferromagnet

abstract: We numerically study the real-time dynamics of a single hole created in the $t-J$ model on a square lattice. Initially, the hole spreads ballistically with a velocity proportional to the hopping matrix element. At intermediate to long times, the dimensionality as well as the spin background determine the hole dynamics. A hole created in the ground state of a two dimensional quantum antiferromagnet propagates again ballistically at long times but with a velocity proportional to the spin exchange coupling, showing the formation of a magnetic polaron. We provide an intuitive explanation of this dynamics in terms of a parton construction, which leads to a good quantitative agreement with the numerical simulations. In the limit of infinite temperature and no spin exchange couplings, the dynamics can be approximated by a quantum random walk on the Bethe lattice. Adding Ising interactions corresponds to an effective disordered potential, which can dramatically slow down the hole propagation, consistent with subdiffusive dynamics.

categories: cond-mat.quant-gas cond-mat.dis-nn cond-mat.str-el quant-ph


title: Flat band of topological states bound to a mobile impurity

abstract: I consider a particle in the topologically non-trivial Su-Schrieffer-Heeger (SSH) model interacting strongly with a mobile impurity, whose quantum dynamics is described by a topologically trivial Hamiltonian. A particle in the SSH model admits a topological zero-energy edge mode when a hard boundary is placed at a given site of the chain, which may be modelled by a static impurity. By solving the two-body problem analytically I show that, when the impurity is mobile, the topological edge states of the Su-Schrieffer-Heeger model remain fully robust and a flat band of bound states at zero energy is formed as long as the continuum spectrum of the two-body problem remains gapped, without the need for any boundaries in the system. This is guaranteed for a sufficiently heavy impurity. As a consequence of the infinite degeneracy of the zero energy modes, it is possible to spatially localise the particle-impurity bound states, effectively making the impurity immobile. These effects can be readily observed using two-dimensional photonic lattices.


title: Measurements of capacitive coupling within a quadruple quantum dot array

abstract: We present measurements of the capacitive coupling energy and the inter-dot capacitances in a linear quadruple quantum dot array in undoped Si/SiGe. With the device tuned to a regime of strong ($>$1 GHz) intra-double dot tunnel coupling, as is typical for double dot qubits, we measure a capacitive coupling energy of $20.9 \pm 0.3$ GHz. In this regime, we demonstrate a fitting procedure to extract all the parameters in the 4D Hamiltonian for two capacitively coupled charge qubits from a 2D slice through the quadruple dot charge stability diagram. We also investigate the tunability of the capacitive coupling energy, using inter-dot barrier gate voltages to tune the inter- and intra-double dot capacitances, and change the capacitive coupling energy of the double dots over a range of 15-32 GHz. We provide a model for the capacitive coupling energy based on the electrostatics of a network of charge nodes joined by capacitors, which shows how the coupling energy should depend on inter-double dot and intra-double dot capacitances in the network, and find that the expected trends agree well with the measurements of coupling energy.


title: High sensitivity maps of molecular ions in L1544: I. Deuteration of N2H+ and HCO+ and first evidence of N2D+ depletion

abstract: Context. The deuterium fraction in low-mass prestellar cores is a good diagnostic indicator of the initial phases of star formations, and it is also a fundamental quantity to infer the ionisation degree in these objects. Aims. With the analysis of multiple transitions of $\rm N_2H^+$, $\rm N_2D^+$, $\rm HC^{18}O^+$ and $\rm DCO^+$ we are able to determine the molecular column density maps and the deuterium fraction in $\rm N_2H^+$ and $\rm HCO^+$ toward the prototypical prestellar core L1544. This is the preliminary step to derive the ionisation degree in the source. Methods. We use a non-local thermodynamic equilibrium (non-LTE) radiative transfer code, combined with the molecular abundances derived from a chemical model, to infer the excitation conditions of all the observed transitions, which allows us to derive reliable maps of each molecule's column density. The ratio between the column density of a deuterated species and its non-deuterated counterpart gives the searched deuteration level. Results. The non-LTE analysis confirms that, for the analysed molecules, higher-J transitions are characterised by excitation temperatures $\approx 1-2\,$K lower than the lower-J ones. The chemical model that provides the best fit to the observational data predict the depletion of $\rm N_2H^+$ and to a lesser extent of $\rm N_2D^+$ in the innermost region. The peak values for the deuterium fraction that we find are $\mathrm{D/H_{N_2H^+}} = 0.26^{+0.15}_{-0.14}$ and $\mathrm{D/H_{HCO^+}} = 0.035^{+0.015}_{-0.012}$, in good agreement with previous estimates in the source.


title: Pion and Kaon Structure at the Electron-Ion Collider

abstract: Understanding the origin and dynamics of hadron structure and in turn that of atomic nuclei is a central goal of nuclear physics. This challenge entails the questions of how does the roughly 1 GeV mass-scale that characterizes atomic nuclei appear; why does it have the observed value; and, enigmatically, why are the composite Nambu-Goldstone (NG) bosons in quantum chromodynamics (QCD) abnormally light in comparison? In this perspective, we provide an analysis of the mass budget of the pion and proton in QCD; discuss the special role of the kaon, which lies near the boundary between dominance of strong and Higgs mass-generation mechanisms; and explain the need for a coherent effort in QCD phenomenology and continuum calculations, in exa-scale computing as provided by lattice QCD, and in experiments to make progress in understanding the origins of hadron masses and the distribution of that mass within them. We compare the unique capabilities foreseen at the electron-ion collider (EIC) with those at the hadron-electron ring accelerator (HERA), the only previous electron-proton collider; and describe five key experimental measurements, enabled by the EIC and aimed at delivering fundamental insights that will generate concrete answers to the questions of how mass and structure arise in the pion and kaon, the Standard Model's NG modes, whose surprisingly low mass is critical to the evolution of our Universe.


title: Physics of relativistic collisionless shocks: II Dynamics of the background plasma

abstract: In this second paper of a series, we discuss the dynamics of a plasma entering the precursor of an unmagnetized, relativistic collisionless pair shock. We discuss how this background plasma is decelerated and heated through its interaction with a microturbulence that results from the growth of a current filamentation instability (CFI) in the shock precursor. We make use, in particular, of the reference frame $\mathcal R_{\rm w}$ in which the turbulence is mostly magnetic. This frame moves at relativistic velocities towards the shock front at rest, decelerating gradually from the far to the near precursor. In a first part, we construct a fluid model to derive the deceleration law of the background plasma expected from the scattering of suprathermal particles off the microturbulence. This law leads to the relationship $\gamma_{\rm p}\,\sim\,\xi_{\rm b}^{-1/2}$ between the background plasma Lorentz factor $\gamma_{\rm p}$ and the normalized pressure of the beam $\xi_{\rm b}$; it is found to match nicely the spatial profiles observed in large-scale 2D3V particle-in-cell simulations. In a second part, we model the dynamics of the background plasma at the kinetic level, incorporating the inertial effects associated with the deceleration of $\mathcal R_{\rm w}$ into a Vlasov-Fokker-Planck equation for pitch-angle diffusion. We show how the effective gravity in $\mathcal R_{\rm w}$ drives the background plasma particles through friction on the microturbulence, leading to efficient plasma heating. Finally, we compare a Monte Carlo simulation of our model with dedicated PIC simulations and conclude that it can satisfactorily reproduce both the heating and the deceleration of the background plasma in the shock precursor, thereby providing a successful 1D description of the shock transition at the microscopic level.


title: Modified swarm-based metaheuristics enhance Gradient Descent initialization performance: Application for EEG spatial filtering

abstract: Gradient Descent (GD) approximators often fail in the solution space with multiple scales of convexities, i.e., in subspace learning and neural network scenarios. To handle that, one solution is to run GD multiple times from different randomized initial states and select the best solution over all experiments. However, this idea is proved impractical in plenty of cases. Even Swarm-based optimizers like Particle Swarm Optimization (PSO) or Imperialistic Competitive Algorithm (ICA), as commonly used GD initializers, have failed to find optimal solutions in some applications. In this paper, Swarm-based optimizers like ICA and PSO are modified by a new optimization framework to improve GD optimization performance. This improvement is for applications with high number of convex localities in multiple scales. Performance of the proposed method is analyzed in a nonlinear subspace filtering objective function over EEG data. The proposed metaheuristic outperforms commonly used baseline optimizers as GD initializers in both the EEG classification accuracy and EEG loss function fitness. The optimizers have been also compared to each other in some of CEC 2014 benchmark functions, where again our method outperforms other algorithms.


title: If the Weak Were Strong and the Strong Were Weak

abstract: We explore the phase structure of the Standard Model as the relative strengths of the SU(2) weak force and SU(3) strong force are varied. With a single generation of fermions, the structure of chiral symmetry breaking suggests that there is no phase transition as we interpolate between the SU(3)-confining phase and the SU(2)-confining phase. Remarkably, the massless left-handed neutrino, familiar in our world, morphs smoothly into a massless right-handed down-quark. With multiple generations, a similar metamorphosis occurs, but now proceeding via a phase transition. In the second half of the paper we introduce a two-parameter extension of the Standard Model, a chiral gauge theory with gauge group U(1) x Sp(r) x SU(N). We again explore the phase structure of the theory as the relative strengths of the Sp(r) and SU(N) gauge couplings vary.


title: Random Walk with Memory on Complex Networks

abstract: We study random walk on complex networks with transition probabilities which depend on the current and previously visited nodes. By using an absorbing Markov chain we derive an exact expression for the mean first passage time between pairs of nodes, for a random walk with a memory of one step. We have analyzed one particular model of random walk, where the transition probabilities depend on the number of paths to the second neighbors. The numerical experiments on paradigmatic complex networks verify the validity of the theoretical expressions, and also indicate that the flattening of the stationary occupation probability accompanies a nearly optimal random search.


title: Renormalized Circuit Complexity

abstract: We propose a modification to Nielsen's circuit complexity for Hamiltonian simulation using the Suzuki-Trotter (ST) method, which provides a network like structure for the quantum circuit. This leads to an optimized gate counting linear in the geodesic distance and spatial volume, unlike in the original proposal. The optimized ST iteration order is correlated with the error tolerance and plays the role of an anti-de Sitter (AdS) radial coordinate. The density of gates is shown to be monotonic with the tolerance and a holographic interpretation using path-integral optimization is given.


title: The one clean qubit model without entanglement is classically simulable

abstract: Entanglement has been shown to be necessary for pure state quantum computation to have an advantage over classical computation. However, it remains open whether entanglement is necessary for quantum computers that use mixed states to also have an advantage. The one clean qubit model is a form of quantum computer in which the input is the maximally mixed state plus one pure qubit. Previous work has shown that there is a limited amount of entanglement present in these computations, despite the fact that they can efficiently solve some problems that are seemingly hard to solve classically. This casts doubt on the notion that entanglement is necessary for quantum speedups. In this work we show that entanglement is indeed crucial for efficient computation in this model, because without it the one clean qubit model is efficiently classically simulable.


title: Dynamical Distance Learning for Semi-Supervised and Unsupervised Skill Discovery

abstract: Reinforcement learning requires manual specification of a reward function to learn a task. While in principle this reward function only needs to specify the task goal, in practice reinforcement learning can be very time-consuming or even infeasible unless the reward function is shaped so as to provide a smooth gradient towards a successful outcome. This shaping is difficult to specify by hand, particularly when the task is learned from raw observations, such as images. In this paper, we study how we can automatically learn dynamical distances: a measure of the expected number of time steps to reach a given goal state from any other state. These dynamical distances can be used to provide well-shaped reward functions for reaching new goals, making it possible to learn complex tasks efficiently. We show that dynamical distances can be used in a semi-supervised regime, where unsupervised interaction with the environment is used to learn the dynamical distances, while a small amount of preference supervision is used to determine the task goal, without any manually engineered reward function or goal examples. We evaluate our method both on a real-world robot and in simulation. We show that our method can learn to turn a valve with a real-world 9-DoF hand, using raw image observations and just ten preference labels, without any other supervision. Videos of the learned skills can be found on the project website: https://sites.google.com/view/dynamical-distance-learning.


title: Who is Afraid of Big Bad Minima? Analysis of Gradient-Flow in a Spiked Matrix-Tensor Model

abstract: Gradient-based algorithms are effective for many machine learning tasks, but despite ample recent effort and some progress, it often remains unclear why they work in practice in optimising high-dimensional non-convex functions and why they find good minima instead of being trapped in spurious ones. Here we present a quantitative theory explaining this behaviour in a spiked matrix-tensor model. Our framework is based on the Kac-Rice analysis of stationary points and a closed-form analysis of gradient-flow originating from statistical physics. We show that there is a well defined region of parameters where the gradient-flow algorithm finds a good global minimum despite the presence of exponentially many spurious local minima. We show that this is achieved by surfing on saddles that have strong negative direction towards the global minima, a phenomenon that is connected to a BBP-type threshold in the Hessian describing the critical points of the landscapes.

categories: cs.LG cond-mat.dis-nn math.ST stat.ML stat.TH


title: Matching coefficients in NRQCD to two-loop accuracy

abstract: We consider the Lagrange density of non-relativistic Quantum Chromodynamics expanded up to order $1/m^2$, where $m$ is the heavy quark mass, and compute several matching coefficients up to two-loop order. Our results are building blocks for next-to-next-to-next-to-leading logarithmic and next-to-next-to-next-to-next-to-leading order corrections to the threshold production of top quark pairs and the decay of heavy quarkonia. We describe the techniques used for the calculation and provide analytic results for a general covariant gauge.


title: TED-On: A Total Error Framework for Digital Traces of Human Behavior on Online Platforms

abstract: Peoples' activities and opinions recorded as digital traces online, especially on social media and other web-based platforms, offer increasingly informative pictures of the public. They promise to allow inferences about populations beyond the users of the platforms on which the traces are recorded, representing real potential for the Social Sciences and a complement to survey-based research. But the use of digital traces brings its own complexities and new error sources to the research enterprise. Recently, researchers have begun to discuss the errors that can occur when digital traces are used to learn about humans and social phenomena. This article synthesizes this discussion and proposes a systematic way to categorize potential errors, inspired by the Total Survey Error (TSE) Framework developed for survey methodology. We introduce a conceptual framework to diagnose, understand, and document errors that may occur in studies based on such digital traces. While there are clear parallels to the well-known error sources in the TSE framework, the new "Total Error Framework for Digital Traces of Human Behavior on Online Platforms" (TED-On) identifies several types of error that are specific to the use of digital traces. By providing a standard vocabulary to describe these errors, the proposed framework is intended to advance communication and research concerning the use of digital traces in scientific social research.


title: A trusted-node-free eight-user metropolitan quantum communication network

abstract: Quantum communication is rapidly gaining popularity due to its high security and technological maturity. However, most implementations are limited to just two communicating parties (users). Quantum communication networks aim to connect a multitude of users. Here we present a fully connected quantum communication network on a city wide scale without active switching or trusted nodes. We demonstrate simultaneous and secure connections between all 28 pairings of 8 users. Our novel network topology is easily scalable to many users, allows traffic management features and minimises the infrastructure as well as the user hardware needed.


title: How many different clonotypes do immune repertoires contain?

abstract: Immune repertoires rely on diversity of T-cell and B-cell receptors to protect us against foreign threats. The ability to recognize a wide variety of pathogens is linked to the number of different clonotypes expressed by an individual. Out of the estimated $\sim 10^{12}$ different B and T cells in humans, how many of them express distinct receptors? We review current and past estimates for these numbers. We point out a fundamental limitation of current methods, which ignore the tail of small clones in the distribution of clone sizes. We show that this tail strongly affects the total number of clones, but it is impractical to access experimentally. We propose that combining statistical models with mechanistic models of lymphocyte clonal dynamics offers possible new strategies for estimating the number of clones.


title: Dark matter in Inert Doublet Model with one scalar singlet and $U(1)_X$ gauge symmetry

abstract: We study Dark Matter (DM) abundance in the framework of the extension of the Standard Model (SM) with an additional $U(1)_X$ gauge symmetry. One complex singlet is included to break the $U(1)_X$ gauge symmetry, meanwhile one of the doublets is considered inert to introduce a DM candidate. The stability of the DM candidate is analyzed with a continuous $U(1)_X$ gauge symmetry as well as discrete $Z_2$ symmetry. We find allowed regions for the free model parameters which are in agreement with the most up-to-date experimental results reported by CMS and ATLAS collaborations, the upper limit on WIMP-nucleon cross section imposed by XENON1T collaboration and the upper limit on the production cross-section of a $Z^{\prime}$ gauge boson times the branching ratio of the $Z^{\prime}$ boson decaying into $\ell^-\ell^+$. We also obtain allowed regions for the DM candidate mass from the relic density reported by the PLANCK collaboration including light, intermediate and heavy masses; depending mainly on two parameters of the scalar potential, $\lambda_{2x}$ and $\lambda_{345}=\lambda_3+\lambda_4+2\lambda_5$. We find that trough $pp\rightarrow \chi\chi \gamma$ production, it may only be possible for a future hadron-hadron Circular Collider (FCC-hh) to be able to detect a DM candidate within the range of masses 10-60 GeV.


title: A note on the classification of Gamma factors

abstract: One of the earliest invariants introduced in the study of finite von Neumann algebras is the property Gamma of Murray and von Neumann. In this note we prove that it is not possible to classify separable $\rm{II}_1$ factors satisfying the property Gamma up to isomorphism by a Borel measurable assignment of countable structures as invariants. We also show that the same holds true for the full $\rm{II}_1$ factors.


title: Analyses of protein cores reveal fundamental differences between solution and crystal structures

abstract: There have been several studies suggesting that protein structures solved by NMR spectroscopy and x-ray crystallography show significant differences. To understand the origin of these differences, we assembled a database of high-quality protein structures solved by both methods. We also find significant differences between NMR and crystal structures---in the root-mean-square deviations of the C$_{\alpha}$ atomic positions, identities of core amino acids, backbone and sidechain dihedral angles, and packing fraction of core residues. In contrast to prior studies, we identify the physical basis for these differences by modelling protein cores as jammed packings of amino-acid-shaped particles. We find that we can tune the jammed packing fraction by varying the degree of thermalization used to generate the packings. For an athermal protocol, we find that the average jammed packing fraction is identical to that observed in the cores of protein structures solved by x-ray crystallography. In contrast, highly thermalized packing-generation protocols yield jammed packing fractions that are even higher than those observed in NMR structures. These results indicate that thermalized systems can pack more densely than athermal systems, which suggests a physical basis for the structural differences between protein structures solved by NMR and x-ray crystallography.


title: Proton impact on ground and excited states of atomic hydrogen

abstract: The processes of electron excitation, capture, and ionization were investigated in proton collisions with atomic hydrogen in the initial $n=1$ and $n=2$ states at impact energies from 1 to 300 keV. The theoretical analysis is based on the close-coupling two-center basis generator method in the semiclassical approximation. Calculated cross sections are compared with previous results which include data obtained from classical-trajectory Monte Carlo, convergent close-coupling, and other two-center atomic orbital expansion approaches. There is an overall good agreement in the capture and excitation cross sections while there are some discrepancies in the ionization results at certain impact energies. These discrepancies in the present results can be partially understood through the use of a $1/n^{3}$ scaling model.


title: Doubly-Adaptive Artificial Compression Methods for Incompressible Flow

abstract: This report presents adaptive artificial compression methods in which the time-step and artificial compression parameter $\varepsilon $ are independently adapted. The resulting algorithms are supported by analysis and numerical tests. The first and second-order methods are embedded. As a result, the computational, cognitive and space complexities of the adaptive $% \varepsilon ,k$ algorithms are negligibly greater than that of the simplest, first-order, constant $\varepsilon ,$ constant $k$ artificial compression method.


title: SLATE: A Super-Lightweight Annotation Tool for Experts

abstract: Many annotation tools have been developed, covering a wide variety of tasks and providing features like user management, pre-processing, and automatic labeling. However, all of these tools use Graphical User Interfaces, and often require substantial effort to install and configure. This paper presents a new annotation tool that is designed to fill the niche of a lightweight interface for users with a terminal-based workflow. Slate supports annotation at different scales (spans of characters, tokens, and lines, or a document) and of different types (free text, labels, and links), with easily customisable keybindings, and unicode support. In a user study comparing with other tools it was consistently the easiest to install and use. Slate fills a need not met by existing systems, and has already been used to annotate two corpora, one of which involved over 250 hours of annotation effort.


title: A Systematic Approach to Detect Hierarchical Healthcare Cost Drivers and Interpretable Change Patterns

abstract: There is strong interest among payers to identify emerging healthcare cost drivers to support early intervention. However, many challenges arise in analyzing large, high dimensional, and noisy healthcare data. In this paper, we propose a systematic approach that utilizes hierarchical and multi-resolution search strategies using enhanced statistical process control (SPC) algorithms to surface high impact cost drivers. Our approach aims to provide interpretable, detailed, and actionable insights of detected change patterns attributing to multiple demographic and clinical factors. We also proposed an algorithm to identify comparable treatment offsets at the population level and quantify the cost impact on their utilization changes.


title: Non-Fermi-liquid behaviors associated with a magnetic quantum-critical point in Sr(Co{1-x}Ni{x})2As2 single crystals

abstract: Electron-doped Sr(Co{1-x}Ni{x})2As2 single crystals with compositions x = 0 to 0.9 were grown out of self-flux and SrNi2As2 single crystals out of Bi flux. The crystals were characterized using single-crystal x-ray diffraction (XRD), magnetic susceptibility chi(H,T), isothermal magnetization M(H,T), heat capacity Cp(H,T), and electrical resistivity ho(H,T) measurements versus applied magnetic field H and temperature T. The chi(T) data show that the crystals exhibit an antiferromagnetic (AFM) ground state almost immediately upon Ni doping on the Co site. Ab-initio electronic-structure calculations for x = 0 and x = 0.15 indicate that a flat band with a peak in the density of states just above the Fermi energy is responsible for this initial magnetic-ordering behavior on Ni doping. The Curie-Weiss-like T dependence of \chi in the paramagnetic (PM) state indicates dominant ferromagnetic (FM) interactions. The small ordered moments ~0.1 muB per transition metal atom and the values of the Rhodes-Wohlfarth ratio indicate that the magnetism is itinerant. The Cp(T) at low T exhibits Fermi-liquid behavior for 0 < x < 0.15 whereas an evolution to a logarithmic non-Fermi-liquid (NFL) behavior is found for x = 0.2 to 0.3. The logarithmic dependence is suppressed in an applied magnetic field. The low- T rho(H = 0,T) data show a T^2 dependence for 0 < x < 0.20 and a power-law dependence with n < 2 for x = 0.20 and 0.30. These low-T NFL behaviors observed in the Cp and rho measurements are most evident near the quantum-critical concentration x ~ 0.3 at which a T = 0 composition-induced transition from the AFM phase to the PM phase occurs.


title: S-parts of values of univariate polynomials

abstract: Let $S=\{p_1,\dots,p_s\}$ be a finite non-empty set of distinct prime numbers, let $f\in \mathbb{Z}[X]$ be a polynomial of degree $n\ge 1$, and let $S'\subseteq S$ be the subset of all $p\in S$ such that $f$ has a root in $\mathbb{Z}_p$. For any non-zero integer $y$, write $y=p_1^{k_1}\dots p_s^{k_s}y_0$, where $k_1,\dots,k_s$ are non-negative integers and $y_0$ is an integer coprime to $p_1,\dots,p_s$. We define the $f$-normalized $S$-part of $y$ by $[y]_{f,S}:=p_1^{k_1 r_{p_1,S}(f)}\dots p_s^{k_s r_{p_s,S}(f)}$, with $r_{p,S}(f)=1$ if $p\in S\setminus S'$ and $r_{p,S}(f)=R_{S'}(f)/R_{p}(f)$ if $p\in S'$, where $R_p(f)$ denotes the largest multiplicity of a root of $f$ in $\mathbb{Z}_p$ and $R_{S'}(f):=\max_{p\in S'} R_p(f)$. For positive real numbers $\varepsilon, B$ with $\varepsilon<R_{S'}(f)/n$, we consider the number $\widetilde{N}(f,S,\varepsilon,B)$ of integers $x$ such that $|x|\le B$ and $0<|f(x)|^{\varepsilon}\le [f(x)]_{f,S}$. We prove that if $s':=\#S'\ge 1$, then $\widetilde{N}(f,S,\varepsilon,B)\asymp_{f,S,\varepsilon} B^{1-(n\varepsilon)/R_{S'}(f)}(\log B)^{s'-1}$ as $B\to \infty$. Moreover, if $f$ has no multiple roots in $\mathbb{Z}_p$ for any $p\in S'$ and $s':=\#S'\ge 2$, then there exists a constant $C(f,S,\varepsilon)>0$ such that $\widetilde{N}(f,S,\varepsilon,B)\sim C(f,S,\varepsilon)\,B^{1-n\varepsilon}(\log B)^{s'-1}$ as $B\to \infty$.


title: Experimental test of exchange fluctuation relations in an open quantum system

abstract: Elucidating the energy transfer between a quantum system and a reservoir is a central issue in quantum non-equilibrium thermodynamics, which could provide novel tools to engineer quantum-enhanced heat engines. The lack of information on the reservoir inherently limits the practical insight that can be gained on the exchange process of open quantum systems. Here, we investigate the energy transfer for an open quantum system in the framework of quantum fluctuation relations. As a novel toolbox, we employ a nitrogen-vacancy center spin qubit in diamond, subject to repeated quantum projective measurements and a tunable dissipation channel. In the presence of energy fluctuations originated by dissipation and quantum projective measurements, the experimental results, supplemented by numerical simulations, show the validity of the energy exchange fluctuation relation, where the energy scale factor encodes missing reservoir information in the system out-of-equilibrium steady state properties. This result is complemented by a theoretical argument showing that, also for an open three-level quantum system, the existence of an out-of-equilibrium steady state dictates a unique time-independent value of the energy scale factor for which the fluctuation relation is verified. Our findings pave the way to the investigation of energy exchange mechanisms in arbitrary open quantum systems.


title: Spontaneous symmetry breaking and Nambu-Goldstone modes in open classical and quantum systems

abstract: We discuss spontaneous symmetry breaking of open classical and quantum systems. When a continuous symmetry is spontaneously broken in an open system, a gapless excitation mode appears corresponding to the Nambu-Goldstone mode. Unlike isolated systems, the gapless mode is not always a propagation mode, but it is a diffusion one. Using the Ward-Takahashi identity and the effective action formalism, we establish the Nambu-Goldstone theorem in open systems, and derive the low-energy coefficients that determine the dispersion relation of Nambu-Goldstone modes. Using these coefficients, we classify the Nambu-Goldstone modes into four types: type-A propagation, type-A diffusion, type-B propagation, and type-B diffusion modes.

categories: hep-th cond-mat.mes-hall cond-mat.quant-gas hep-ph


title: Misplaced Confidence in Observed Power

abstract: A recently published randomized controlled trial in JAMA investigated the impact of the selective serotonin reuptake inhibitor, escitalopram, on the risk of major adverse events (MACE). The authors estimated a hazard ratio (HR) of 0.69 (95% CI: 0.49, 0.96; $p$ = 0.03) and then attempted to calculate how much statistical power their study (test) had attained, and used this measure to assess how reliable their results were. Here, we discuss why this approach, along with other post-hoc power analyses, are highly misleading.


title: Joint Learning of Named Entity Recognition and Entity Linking

abstract: Named entity recognition (NER) and entity linking (EL) are two fundamentally related tasks, since in order to perform EL, first the mentions to entities have to be detected. However, most entity linking approaches disregard the mention detection part, assuming that the correct mentions have been previously detected. In this paper, we perform joint learning of NER and EL to leverage their relatedness and obtain a more robust and generalisable system. For that, we introduce a model inspired by the Stack-LSTM approach (Dyer et al., 2015). We observe that, in fact, doing multi-task learning of NER and EL improves the performance in both tasks when comparing with models trained with individual objectives. Furthermore, we achieve results competitive with the state-of-the-art in both NER and EL.


title: Four Small Planets Buried in K2 Systems: What Can We Learn for TESS?

abstract: The Kepler, K2, and Transiting Exoplanet Survey Satellite (TESS) missions have provided a wealth of confirmed exoplanets, benefiting from a huge effort from the planet-hunting and follow-up community. With careful systematics mitigation, these missions provide precise photometric time series, which enable detection of transiting exoplanet signals. However, exoplanet hunting can be confounded by several factors, including instrumental noise, search biases, and host star variability. In this Letter, we discuss strategies to overcome these challenges using newly emerging techniques and tools. We demonstrate the power of new, fast open-source community tools (e.g., lightkurve, starry, celerite, exoplanet), and discuss four high signal-to-noise ratio (S/N) exoplanets that showcase specific challenges present in planet detection: K2-43c, K2-168c, K2-198c, and K2-198d. These planets have been undetected in several large K2 planet searches, despite having transit signals with S/N > 10. Two of the planets discussed here are new discoveries. In this work we confirm all four as true planets. Alongside these planet systems, we discuss three key challenges in finding small transiting exoplanets. The aim of this Letter is to help new researchers understand where planet detection efficiency gains can be made, and to encourage the continued use of K2 archive data. The considerations presented in this Letter are equally applicable to Kepler, K2, and TESS, and the tools discussed here are available for the community to apply to improve exoplanet discovery and fitting.


title: Bayesian Variable Selection for Gaussian copula regression models

abstract: We develop a novel Bayesian method to select important predictors in regression models with multiple responses of diverse types. A sparse Gaussian copula regression model is used to account for the multivariate dependencies between any combination of discrete and/or continuous responses and their association with a set of predictors. We utilize the parameter expansion for data augmentation strategy to construct a Markov chain Monte Carlo algorithm for the estimation of the parameters and the latent variables of the model. Based on a centered parametrization of the Gaussian latent variables, we design a fixed-dimensional proposal distribution to update jointly the latent binary vectors of important predictors and the corresponding non-zero regression coefficients. For Gaussian responses and for outcomes that can be modeled as a dependent version of a Gaussian response, this proposal leads to a Metropolis-Hastings step that allows an efficient exploration of the predictors' model space. The proposed strategy is tested on simulated data and applied to real data sets in which the responses consist of low-intensity counts, binary, ordinal and continuous variables.


title: Finding First and Most-Beautiful Queens by Integer Programming

abstract: The n-queens puzzle is a well-known combinatorial problem that requires to place n queens on an n x n chessboard so that no two queens can attack each other. Since the 19th century, this problem was studied by many mathematicians and computer scientists. While finding any solution to the n-queens puzzle is rather straightforward, it is very challenging to find the lexicographically first (or smallest) feasible solution. Solutions for this type are known in the literature for n <= 55, while for some larger chessboards only partial solutions are known. The present paper was motivated by the question of whether Integer Linear Programming (ILP) can be used to compute solutions for some open instances. We describe alternative ILP-based solution approaches, and show that they are indeed able to compute (sometimes in unexpectedly-short computing times) many new lexicographically optimal solutions for n ranging from 56 to 115. One of the proposed algorithms is a pure cutting plane method based on a combinatorial variant of classical Gomory cuts. We also address an intriguing "lexicographic bottleneck" (or min-max) variant of the problem that requires finding a most beautiful (in a well defined sense) placement, and report its solution for n up to 176.


title: The abelian complexity of infinite words and the Frobenius problem

abstract: We study the following problem, first introduced by Dekking. Consider an infinite word x over an alphabet {0,1,...,k-1} and a semigroup homomorphism S:{0,1,...,k-1}* -> N. Let L_x denote the set of factors of x. What conditions on S and the abelian complexity of x guarantee that S(L_x) contains all but finitely many elements of N? We examine this question for some specific infinite words x having different abelian complexity functions.


title: Tuning spinterface properties in Iron/Fullerene thin films

abstract: In ferromagnetic (FM) metal/organic semiconductor (OSC) heterostructures charge transfer can occur which leads to induction of magnetism in the non-magnetic OSC. This phenomenon has been described by the change in the density of states in the OSC which leads to a finite magnetic moment at the OSC interface and it is called the "spinterface". One of the main motivation in this field of organic spintronics is how to control the magnetic moment in the spinterface. In this regard, there are several open questions such as (i) which combination of FM and OSC can lead to more moment at the spinterface? (ii) Is the thickness of OSC also important? (iii) How does the spinterface moment vary with the FM thickness? (iv) Does the crystalline quality of the FM matters? (v) What is the effect of spinterface on magnetization reversal, domain structure and anisotropy? In this context, we have tried to answer the last three issues in this paper by studying Fe/C$_{60}$ bilayers of variable Fe thickness deposited on Si substrates. We find that both the induced moment and thickness of the spinterface vary proportionally with the Fe thickness. Such behavior is explained in terms of the growth quality of the Fe layer on the native oxide of the Si (100) substrate. The magnetization reversal, domain structure and anisotropy of these bilayer samples were studied and compared with their respective reference samples without having the C$_{60}$ layer. It is observed that the formation of spinterface leads to reduction in uniaxial anisotropy in Fe/C$_{60}$ on Si (100) in comparison to their reference samples.


title: Basis Function Method for Numerical Loop Quantum Cosmology: The Schwarzschild Black Hole Interior

abstract: Loop quantum cosmology is a symmetry-reduced application of loop quantum gravity that has led to the resolution of classical singularities such as the big bang, and those at the center of black holes. This can be seen through numerical simulations involving the quantum Hamiltonian constraint that is a partial difference equation. The equation allows one to study the evolution of sharply-peaked Gaussian wave packets that generically exhibit a quantum "bounce" or a non-singular passage through the classical singularity, thus offering complete singularity resolution. In addition, von-Neumann stability analysis of the difference equation -- treated as a stencil for a numerical solution that steps through the triad variables -- yields useful constraints on the model and the allowed space of states. In this paper, we develop a new method for the numerical solution of loop quantum cosmology models using a set of basis functions that offer a number of advantages over computing a solution by stepping through the triad variables. We use the Corichi and Singh model for the Schwarzschild interior as the main case study in this effort. The main advantage of this new method is computational efficiency and the ease of parallelization. In addition, we also discuss how the stability analysis appears in the context of this new approach.


title: Multi-Objective Optimization for Energy-and Spectral-Efficiency Tradeoff in In-band Full-Duplex (IBFD) Communication

abstract: The problem of joint power and sub-channel allocation to maximize energy efficiency (EE) and spectral efficiency (SE) simultaneously in in-band full-duplex (IBFD) orthogonal frequency-division multiple access (OFDMA) network is addressed considering users' QoS in both uplink and downlink. The resulting optimization problem is a non-convex mixed-integer non-linear program (MINLP) which is generally difficult to solve. In order to strike a balance between the EE and SE, we restate this problem as a multi-objective optimization problem (MOOP) which aims at maximizing system's throughput and minimizing system's power consumption, simultaneously. To this end, the \epsilon constraint method is adopted to transform the MOOP into single-objective optimization problem (SOOP). The underlying problem is solved via an efficient solution based on the majorization minimization (MM) approach. Furthermore, in order to handle binary subchannel allocation variable constraints, a penalty function is introduced. Simulation results unveil interesting tradeoffs between EE and SE.


title: Responsibility Analysis by Abstract Interpretation

abstract: Given a behavior of interest in the program, statically determining the corresponding responsible entity is a task of critical importance, especially in program security. Classical static analysis techniques (e.g. dependency analysis, taint analysis, slicing, etc.) assist programmers in narrowing down the scope of responsibility, but none of them can explicitly identify the responsible entity. Meanwhile, the causality analysis is generally not pertinent for analyzing programs, and the structural equations model (SEM) of actual causality misses some information inherent in programs, making its analysis on programs imprecise. In this paper, a novel definition of responsibility based on the abstraction of event trace semantics is proposed, which can be applied in program security and other scientific fields. Briefly speaking, an entity ER is responsible for behavior B, if and only if ER is free to choose its input value, and such a choice is the first one that ensures the occurrence of B in the forthcoming execution. Compared to current analysis methods, the responsibility analysis is more precise. In addition, our definition of responsibility takes into account the cognizance of the observer, which, to the best of our knowledge, is a new innovative idea in program analysis.


title: Message passing on networks with loops

abstract: In this paper we offer a solution to a long-standing problem in the study of networks. Message passing is a fundamental technique for calculations on networks and graphs. The first versions of the method appeared in the 1930s and over the decades it has been applied to a wide range of foundational problems in mathematics, physics, computer science, statistics, and machine learning, including Bayesian inference, spin models, coloring, satisfiability, graph partitioning, network epidemiology, and the calculation of matrix eigenvalues. Despite its wide use, however, it has long been recognized that the method has a fundamental flaw: it only works on networks that are free of short loops. Loops introduce correlations that cause the method to give inaccurate answers at best, and to fail completely in the worst cases. Unfortunately, almost all real-world networks contain many short loops, which limits the usefulness of the message passing approach. In this paper we demonstrate how to rectify this shortcoming and create message passing methods that work on any network. We give two example applications, one to the percolation properties of networks and the other to the calculation of the spectra of sparse matrices.


title: Low power (mW) nonlinearities of polarization maintaining fibers

abstract: Polarization maintaining (PM) fibers are meant to maintain linear polarization along a preferred axis. A PM fiber can be seen as the fiber version of a very high order waveplate, designed with different refractive indices along two orthogonal axes. It is shown that monitoring the polarization of initially circularly polarized light sent through a PM fiber, leads not only to new sensing methods, but also to power control, saturable absorption, and optical path stabilization. Even at peak power levels not exceeding a few mW, nonlinear transmission is detected, with time constants in the microsecond range.


title: Wave Excitation Force Estimation of Wave Energy Floats Using Extended Kalman Filters

abstract: In many advanced control strategies the wave excitation force is key to determining the control input. However, it is often difficult to measure the excitation force on a Wave Energy Converter (WEC). The use of Kalman filters to estimate the wave excitation force based on readily available measurement data can potentially fill the gap between the development of WEC control strategies and the data that is available. Two different estimation methods using an nonlinear Extended Kalman Filter are tested on experimental wave tank data for a heaving semi-submerged float. The first method relies on directly including the excitation force as a state in the first order dynamics---which allows the "random walk" of the Kalman filter to identify an estimate of the excitation force. The second method of estimation involves modeling the wave excitation force as a harmonic oscillator comprised of sinusoidal components. Both methods are evaluated for a variety of incident waves and additional sensitivity analyses are performed to investigate the susceptibility of these estimation methods to changes in the model, measurement noise, and sampling rate.


title: Cohomology and deformations of dendriform coalgebras

abstract: Dendriform coalgebras are the dual notion of dendriform algebras and are splitting of associative coalgebras. In this paper, we define a cohomology theory for dendriform coalgebras based on some combinatorial maps. We show that the cohomology with self coefficients governs the formal deformation of the structure. We also relate this cohomology with the cohomology of dendriform algebras, coHochschild (Cartier) cohomology of associative coalgebras and cohomology of Rota-Baxter coalgebras which we introduce in this paper. Finally, using those combinatorial maps, we introduce homotopy analogue of dendriform coalgebras and study some of their properties.


title: Impact of DCSB and dynamical diquark correlations on proton GPDs

abstract: We calculate the leading-twist, helicity-independent generalized parton distributions (GPDs) of the proton, at finite skewness, in the Nambu--Jona-Lasinio (NJL) model of quantum chomodynamics (QCD). The NJL model reproduces low-energy characteristics of QCD, including dynamical chiral symmetry breaking (DCSB). The proton bound-state amplitude is solved for using the Faddeev equation in a quark-diquark approximation, including both dynamical scalar and axial vector diquarks. GPDs are calculated using a dressed non-local correlator, consistent with DCSB, which is obtained by solving a Bethe-Salpeter equation. The model and approximations used observe Lorentz covariance, and as a consequence the GPDs obey polynomiality sum rules. Extractions of electromagnetic and gravitational form factors are performed. We find a D-term of $-1.08$ when the non-local correlator is properly dressed, and $0.85$ when the bare correlator is used instead, suggesting that within this framework proton stability requires the constituent quarks to be dressed consistently with DCSB. We also find that the anomalous gravitomagnetic vanishes, as required by Poincar\'{e} symmetry.


title: On the Monniaux Problem in Abstract Interpretation

abstract: The Monniaux Problem in abstract interpretation asks, roughly speaking, whether the following question is decidable: given a program $P$, a safety (\emph{e.g.}, non-reachability) specification $\varphi$, and an abstract domain of invariants $\mathcal{D}$, does there exist an inductive invariant $I$ in $\mathcal{D}$ guaranteeing that program $P$ meets its specification $\varphi$. The Monniaux Problem is of course parameterised by the classes of programs and invariant domains that one considers. In this paper, we show that the Monniaux Problem is undecidable for unguarded affine programs and semilinear invariants (unions of polyhedra). Moreover, we show that decidability is recovered in the important special case of simple linear loops.


title: GBDT and explicit solutions for the matrix coupled dispersionless equations (local and nonlocal cases)

abstract: We introduce matrix coupled (local and nonlocal) dispersionless equations, construct wide classes of explicit multipole solutions, give explicit expressions for the corresponding Darboux and wave matrix valued functions and consider their asymptotics in some interesting cases. We consider the scalar cases of coupled, complex coupled and nonlocal dispersionless equations as well.


title: WriterForcing: Generating more interesting story endings

abstract: We study the problem of generating interesting endings for stories. Neural generative models have shown promising results for various text generation problems. Sequence to Sequence (Seq2Seq) models are typically trained to generate a single output sequence for a given input sequence. However, in the context of a story, multiple endings are possible. Seq2Seq models tend to ignore the context and generate generic and dull responses. Very few works have studied generating diverse and interesting story endings for a given story context. In this paper, we propose models which generate more diverse and interesting outputs by 1) training models to focus attention on important keyphrases of the story, and 2) promoting generation of non-generic words. We show that the combination of the two leads to more diverse and interesting endings.


title: A geometric approach to the transport of discontinuous densities

abstract: Different observations of a relation between inputs ("sources") and outputs ("targets") are often reported in terms of histograms (discretizations of the source and the target densities). Transporting these densities to each other provides insight regarding the underlying relation. In (forward) uncertainty quantification, one typically studies how the distribution of inputs to a system affects the distribution of the system responses. Here, we focus on the identification of the system (the transport map) itself, once the input and output distributions are determined, and suggest a modification of current practice by including data from what we call "an observation process". We hypothesize that there exists a smooth manifold underlying the relation; the sources and the targets are then partial observations (possibly projections) of this manifold. Knowledge of such a manifold implies knowledge of the relation, and thus of "the right" transport between source and target observations. When the source-target observations are not bijective (when the manifold is not the graph of a function over both observation spaces, either because folds over them give rise to density singularities, or because it marginalizes over several observables), recovery of the manifold is obscured. Using ideas from attractor reconstruction in dynamical systems, we demonstrate how additional information in the form of short histories of an observation process can help us recover the underlying manifold. The types of additional information employed and the relation to optimal transport based solely on density observations is illustrated and discussed, along with limitations in the recovery of the true underlying relation.


title: Solenoidal scaling laws for compressible mixing

abstract: Mixing of passive scalars in compressible turbulence does not obey the same classical Reynolds number scaling as its incompressible counterpart. We first show from a large database of direct numerical simulations that even the solenoidal part of the velocity field fails to follow the classical incompressible scaling when the forcing includes a substantial dilatational component. Though the dilatational effects on the flow remain significant, our main results are that both the solenoidal energy spectrum and the passive scalar spectrum scale assume incompressible forms, and that the scalar gradient aligns with the most compressive eigenvalue of the solenoidal part, provided that only the solenoidal components are used for scaling in a consistent manner. Minor modifications to this statement are also pointed out.


title: Resonant Assisted Annihilation

abstract: Assisted annihilation is a novel mechanism to generate viable sub-GeV thermal dark matter, where a pair of stable dark matter annihilates with an assister to Standard Model states. Typically such $3 \to 2$ annihilation topologies are flux suppressed compared to $2 \to 2$ processes. In this paper, we explore the possibility of a resonant $3 \to 2$ assisted annihilation dominantly driving the freeze-out of dark matter. We demonstrate that in a simple multipartite scalar extension of the Standard Model this can be realized in certain regions of parameter space to provide viable dark matter relic density, in agreement with observation. We demonstrate that for photophilic assisters parts of the parameter space are already constrained by indirect detection experiments and the measurements of CMB anisotropies while substantial regions remain beyond the present limit.


title: Manipulation of Multimode Squeezing in a Coupled Waveguide Array

abstract: We present a scheme for generating and manipulating three-mode squeezed states with genuine tripartite entanglement by injecting single-mode squeezed light into an array of coupled optical waveguides. We explore the possibility to selectively generate single-mode squeezing or multimode squeezing at the output of an elliptical waveguides array, determined solely by the input light polarization. We study the effect of losses in the waveguides array and show that quantum correlations and squeezing are preserved for realistic parameters. Our results show that arrays of optical waveguides are suitable platforms for generating multimode quantum light, which could lead to novel applications in quantum metrology.


title: A general approach to the assessment of uncertainty in volumes by using the multi-Gaussian model

abstract: The goal of this research is to derive an approach to assess uncertainty in an arbitrary volume conditioned by sampling data, without using geostatistical simulation. We have accomplished this goal by deriving an numerical tool suitable for any probabilistic distribution of the sample data. For this, we have worked with an extension of the traditional multi-Gaussian model, allowing us to obtain a formulation that makes explicit the dependence of the uncertainty in the arbitrary volume from the grades within the volume, the spatial correlation of the data and the conditioning values. A Kriging of the Gaussian values is the only requirement to obtain not only conditional local means and variances but also the complete local distributions at any support, in an easy and straightforward way.


title: Neutral hydrogen gas within and around NGC 1316

abstract: We present MeerKAT observations of neutral hydrogen gas (HI) in the nearby merger remnant NGC 1316 (Fornax A), the brightest member of a galaxy group which is falling into the Fornax cluster. We find HI on a variety of scales, from the galaxy centre to its large-scale environment. For the first time we detect HI at large radii (70 - 150 kpc in projection), mostly distributed on two long tails associated with the galaxy. Gas in the tails dominates the HI mass of NGC 1316: 7e+8 Msun -- 14 times more than in previous observations. The total HI mass is comparable to the amount of neutral gas found inside the stellar body, mostly in molecular form. The HI tails are associated with faint optical tidal features thought to be the remnant of a galaxy merger occurred a few billion years ago. They demonstrate that the merger was gas-rich. During the merger, tidal forces pulled some gas and stars out to large radii, where we now detect them in the form of optical tails and, thanks to our new data, HI tails; while torques caused the remaining gas to flow towards the centre of the remnant, where it was converted into molecular gas and fuelled the starburst revealed by the galaxy's stellar populations. Several of the observed properties of NGC 1316 can be reproduced by a ~ 10:1 merger between a dominant, gas-poor early-type galaxy and a smaller, gas-rich spiral occurred 1 - 3 Gyr ago, likely followed by subsequent accretion of satellite galaxies.


title: Q-type asteroids: Possibility of non-fresh weathered surfaces

abstract: Itokawa particles, which are the recovered samples from the S-complex asteroid 25143 Itokawa by the Hayabusa spacecraft, demonstrate that S-complex asteroids are parent bodies of ordinary chondrite meteorites. Furthermore, they clarify that the space weathering age of the Itokawa surface is of the order of several thousand years. Traditionally, Q-type asteroids have been considered fresh-surfaced. However, as the space weathering timescale is approximately three orders of magnitude lesser than the conventionally considered age, the previously proposed formation mechanisms of Q-type asteroids cannot sufficiently explain the surface refreshening. In this study, we propose a new hypothesis on the surface state of Q-type asteroids: Q-type asteroids have a non-fresh weathered surface with a paucity of fine particles. For verifying this hypothesis, laboratory experiments on the space weathering of ordinary chondrites are performed. Based on the results of these experiments, we found that large (more than 100 {\mu}m) ordinary chondritic particles with space weathering exhibit spectra consistent with Q-type asteroids.


title: Implementation of a Hamming-Distance-Like Genomic Quantum Classifier Using Inner Products on IBMQX4 and IBMQX16

abstract: Motivated by the problem of classifying individuals with a disease versus controls using functional genomic attributes as input, we encode the input as a string of 1s (presence) or 0s (absence) of the genomic attribute across the genome. Blocks of physical regions in the subdivided genome serve as the feature dimensions, which takes full advantage of remaining in the computational basis of a quantum computer. Given that a natural distance between two binary strings is the Hamming distance and that this distance shares properties with an inner product between qubits, we developed two Hamming-distance-like classifiers which apply two different kinds of inner products ("active" and "symmetric") to directly and efficiently classify a test input into either of two training classes. To account for multiple disease and control samples, each train class can be composed of an arbitrary number of bit strings (i.e., number of samples) that can be compressed into one input string per class. Thus, our circuits require the same number of qubits regardless of the number of training samples. These algorithms, which implement a training bisection decision plane, were simulated and implemented on IBMQX4 and IBMQX16. The latter allowed for encoding of $64$ training features across the genome for $2$ (disease and control) training classes.


title: Discrete Object Generation with Reversible Inductive Construction

abstract: The success of generative modeling in continuous domains has led to a surge of interest in generating discrete data such as molecules, source code, and graphs. However, construction histories for these discrete objects are typically not unique and so generative models must reason about intractably large spaces in order to learn. Additionally, structured discrete domains are often characterized by strict constraints on what constitutes a valid object and generative models must respect these requirements in order to produce useful novel samples. Here, we present a generative model for discrete objects employing a Markov chain where transitions are restricted to a set of local operations that preserve validity. Building off of generative interpretations of denoising autoencoders, the Markov chain alternates between producing 1) a sequence of corrupted objects that are valid but not from the data distribution, and 2) a learned reconstruction distribution that attempts to fix the corruptions while also preserving validity. This approach constrains the generative model to only produce valid objects, requires the learner to only discover local modifications to the objects, and avoids marginalization over an unknown and potentially large space of construction histories. We evaluate the proposed approach on two highly structured discrete domains, molecules and Laman graphs, and find that it compares favorably to alternative methods at capturing distributional statistics for a host of semantically relevant metrics.


title: Transit timing variations in the WASP-4 planetary system

abstract: Transits in the planetary system WASP-4 were recently found to occur 80s earlier than expected in observations from the TESS satellite. We present 22 new times of mid-transit that confirm the existence of transit timing variations, and are well fitted by a quadratic ephemeris with period decay dP/dt = -9.2 +/- 1.1 ms/yr. We rule out instrumental issues, stellar activity and the Applegate mechanism as possible causes. The light-time effect is also not favoured due to the non-detection of changes in the systemic velocity. Orbital decay and apsidal precession are plausible but unproven. WASP-4b is only the third hot Jupiter known to show transit timing variations to high confidence. We discuss a variety of observations of this and other planetary systems that would be useful in improving our understanding of WASP-4 in particular and orbital decay in general.


title: Energy-momentum and angular-momentum of a gyratonic pp-waves spacetime

abstract: Gyratonic plane fronted gravitational waves are exact solutions of Einstein's field equations, which correspond to gravitational waves that carry momentum and angular-momentum. Using the definitions of the Hamiltonian formulation of the Teleparallel Equivalent of General Relativity, we explicitly evaluate the general expressions of the energy-momentum and angular-momentum of these space-times. In order to better understand the additional properties of these gravitational waves, we consider the motion of particles in this space-time and obtain an interesting relation between the angular-momentum of the particles and that of the gravitational waves.


title: Astro2020 Project White Paper: The Cosmic Accelerometer

abstract: We propose an experiment, the Cosmic Accelerometer, designed to yield velocity precision of $\leq 1$ cm/s with measurement stability over years to decades. The first-phase Cosmic Accelerometer, which is at the scale of the Astro2020 Small programs, will be ideal for precision radial velocity measurements of terrestrial exoplanets in the Habitable Zone of Sun-like stars. At the same time, this experiment will serve as the technical pathfinder and facility core for a second-phase larger facility at the Medium scale, which can provide a significant detection of cosmological redshift drift on a 6-year timescale. This larger facility will naturally provide further detection/study of Earth twin planet systems as part of its external calibration process. This experiment is fundamentally enabled by a novel low-cost telescope technology called PolyOculus, which harnesses recent advances in commercial off the shelf equipment (telescopes, CCD cameras, and control computers) combined with a novel optical architecture to produce telescope collecting areas equivalent to standard telescopes with large mirror diameters. Combining a PolyOculus array with an actively-stabilized high-precision radial velocity spectrograph provides a unique facility with novel calibration features to achieve the performance requirements for the Cosmic Accelerometer.


title: Weak Adversarial Networks for High-dimensional Partial Differential Equations

abstract: Solving general high-dimensional partial differential equations (PDE) is a long-standing challenge in numerical mathematics. In this paper, we propose a novel approach to solve high-dimensional linear and nonlinear PDEs defined on arbitrary domains by leveraging their weak formulations. We convert the problem of finding the weak solution of PDEs into an operator norm minimization problem induced from the weak formulation. The weak solution and the test function in the weak formulation are then parameterized as the primal and adversarial networks respectively, which are alternately updated to approximate the optimal network parameter setting. Our approach, termed as the weak adversarial network (WAN), is fast, stable, and completely mesh-free, which is particularly suitable for high-dimensional PDEs defined on irregular domains where the classical numerical methods based on finite differences and finite elements suffer the issues of slow computation, instability and the curse of dimensionality. We apply our method to a variety of test problems with high-dimensional PDEs to demonstrate its promising performance.


title: Astro2020 Project White Paper: PolyOculus -- Low-cost Spectroscopy for the Community

abstract: As astronomy moves into the era of large-scale time-domain surveys, we are seeing a flood of new transient and variable sources which will reach biblical proportions with the advent of LSST. A key strategic challenge for astronomy in this era is the lack of suitable spectroscopic followup facilities. In response to this need, we have developed the PolyOculus approach for producing large-area-equivalent telescopes by using fiber optics to link modules of multiple semi-autonomous, small, inexpensive, commercial-off-the-shelf telescopes. Crucially, this scalable design has construction costs which are $>10x$ lower than equivalent traditional large-area telescopes. In addition, PolyOculus is inherently highly automated and well-suited for remote operations. Development of this technology will enable the expansion of major research efforts in the LSST era to a host of smaller universities and colleges, including primarily-undergraduate institutions, for budgets consistent with their educational expenditures on similar facilities. We propose to develop and deploy a 1.6-m prototype demonstrator at the Mt. Laguna Observatory in California, followed by a full-scale 5-meter-class PolyOculus facility for linkage to existing and upcoming time-domain surveys.


title: Unitary operators with decomposable corners

abstract: We study pairs $(U,L_0)$, where $U$ is a unitary operator in $H$ and $L_0\subset H$ is a closed subspace, such that $$ P_{L_0}U|_{L_0}:L_0\to L_0 $$ has a singular value decomposition. Abstract characterizations of this condition are given, as well as relations to the geometry of projections and pairs of projections. Several concrete examples are examined.


title: The Purity Conjecture in Type $C$

abstract: A collection $\mathcal{C}$ of $k$-element subsets of $\{1,2,\ldots,m\}$ is weakly separated if for each $I, J \in \mathcal{C}$, when the integers $1,2,\ldots,m$ are arranged around in a circle, there is a chord separating $I \backslash J$ from $J \backslash I$. Oh, Postnikov and Speyer constructed a correspondence between weakly separated collections which are maximal by inclusion and reduced plabic graphs, a class of networks defined by Postnikov which give coordinate charts on the Grassmannian of $k$-planes in m-space. As a corollary, they proved Scott's Purity Conjecture, which states that a weakly separated collection is maximal by inclusion if and only if it is maximal by size. In this note, we describe maximal weakly separated collections corresponding to symmetric plabic graphs, which give coordinate charts on the Lagrangian Grassmannian, and prove a symmetric version of the Purity Conjecture.


title: An AI-based, Multi-stage detection system of banking botnets

abstract: Banking Trojans, botnets are primary drivers of financially-motivated cybercrime. In this paper, we first analyzed how an APT-based banking botnet works step by step through the whole lifecycle. Specifically, we present a multi-stage system that detects malicious banking botnet activities which potentially target the organizations. The system leverages Cyber Data Lake as well as multiple artificial intelligence techniques at different stages. The evaluation results using public datasets showed that Deep Learning based detections were highly successful compared with baseline models.


title: The interplay between tissue healing and bone remodeling around immediately loaded tooth replacement implants

abstract: Long-term bone healing/adaptation after a dental implant treatment starts with diffusion of mesenchymal stem cells to the fracture callus and their subsequent differentiation. The healing phase is followed by the bone-remodeling phase. In this work, a mechano-regulatory cellular differentiation model was used to simulate tissue healing around an immediately loaded dental implant. All tissue types were modeled as poroelastic in the healing phase. Material properties of the healing region were updated after each loading cycle for 30 cycles (days). The tissue distribution in the healed state was then used as the initial condition for the remodeling phase during which regions healed into bone adapt their internal density with respect to a homeostatic remodeling stimulus. The short- and long-term effects of micro-motion on bone healing and remodeling were studied. Development of soft tissue was observed both in the coronal region due to high fluid velocity, and on the vertical sides of the healing-callus due to high shear stress. In cases with small implant micromotion, tissue between the implant threads differentiated into bone during the healing phase, but resorbed during remodeling. In cases with large implant micromotion, higher percentage of the healing region differentiated into soft tissue resulting in less volume available for bone remodeling. But, the remaining bone region developed higher density bone tissue. It was concluded that an optimal range of controlled micromotion could be designed for a given implant in order to achieve the desired functional properties.


title: Fog Function: Serverless Fog Computing for Data Intensive IoT Services

abstract: Fog computing can support IoT services with fast response time and low bandwidth usage by moving computation from the cloud to edge devices. However, existing fog computing frameworks have limited flexibility to support dynamic service composition with a data-oriented approach. Function-as-a-Service (FaaS) is a promising programming model for fog computing to enhance flexibility, but the current event- or topic-based design of function triggering and the separation of data management and function execution result in inefficiency for data-intensive IoT services. To achieve both flexibility and efficiency, we propose a data-centric programming model called Fog Function and also introduce its underlying orchestration mechanism that leverages three types of contexts: data context, system context, and usage context. Moreover, we showcase a concrete use case for smart parking where Fog Function allows service developers to easily model their service logic with reduced learning efforts compared to a static service topology. Our performance evaluation results show that the Fog Function can be scaled to hundreds of fog nodes. Fog Function can improve system efficiency by saving 95% of the internal data traffic over cloud function and it can reduce service latency by 30% over edge function.


title: Statistical Clear Sky Fitting Algorithm

abstract: We present an algorithm that estimates a clear sky performance signal from the measured power of a PV system. The algorithm uses only observed power output, and assumes no knowledge of weather, irradiance data, or system configuration metadata. This is a novel approach to understanding the clear sky behavior of an installed PV system, that does not rely on traditional atmospheric and geometric modeling techniques.


title: Performance and near-wake characterization of a tidal current turbine in elevated levels of free stream turbulence

abstract: Tidal turbines are deployed in sites which have elevated levels of free stream turbulence (FST). Accounting for elevated FST on their operation become vital from a design standpoint. Detailed experimental measurements of the dynamic near-wake of a tidal turbine model in elevated FST environments is presented; an active grid turbulence generator developed by our group was used to seed in the elevated FST and evaluate the influence of turbulence intensity (Ti) and inflow integral length scale (L) on the near-wake of the turbine. Three inflow conditions are tested: a quasi-laminar flow with Ti ~ 2.2% and two elevated Ti (~12-14%) cases, one with L ~ 0.4D (D is the turbine diameter) and the other where L~ D. Elevated Ti cases was found to increase the standard deviation of rotor torque by 4.5 times the value in quasi-laminar flow. Energy recovery was also found to be accelerated; at X/D=4, the percentage of inflow energy recovered was 37% and was twice the corresponding value in quasi-laminar flow. Elevated FST was observed to disrupt the rotational character of the wake; the drop in swirl number ranged between 12% at X/D=0.5 to 71% at X/D=4. Elevated Ti also resulted in L that were considerably larger (> 2 times) than the quasi-laminar flow case. An increase in inflow integral length scale (from 0.4D to D) was observed to result in enhanced wake Ti, wake structures and anisotropy; however, no noticeable influence was found on the rate of wake recovery.


title: Entanglement, squeezing and state-swap in cholesteric liquid crystals

abstract: Phenomena such as state-swap, quadrature squeezing, entanglement and violation of entanglement inequalities are frequently reported to occur in quantum systems only. It is shown here that these effects can also occur in cholesteric liquid crystals in the presence of an applied magnetic field. The nature of these effects as well as the conditions for their observation are analysed.


title: Acoustic Multifunctional Logic Gates and Amplifier based on Passive Parity-Time Symmetry

abstract: Acoustic analogue computation and signal processing are of great significance, however, it's challenging to realize the acoustic computing devices because of their limitations of single function and complex structure. In this paper, an acoustic multifunctional device, which can gate or amplify acoustic waves without resorting to altering the frequency and structure using a passive acoustic parity-time (PT)-symmetric metamaterial, is realized theoretically and experimentally. The metamaterial is constructed by five lossless-loss periodically distributed media which are modulated to achieve the passive PT symmetry. At the coherent perfect absorber (CPA)-emitter point in the broken PT-symmetric phase, the logic gates (AND, OR, XOR and NOT) and small signal amplifier are realized in a single system by adjusting the phase and amplitude differences between two incoming beams, respectively. This work provides a new route for the connection between the PT symmetry and the acoustic metamaterial, which has great potential applications in acoustic modulation and acoustic multifunctional device.


title: Public Plug-in Electric Vehicles + Grid Data: Is a New Cyberattack Vector Viable?

abstract: High-wattage demand-side appliances such as Plug-in Electric Vehicles (PEVs) are proliferating. As a result, information on the charging patterns of PEVs is becoming accessible via smartphone applications, which aggregate real-time availability and historical usage of public PEV charging stations. Moreover, information on the power grid infrastructure and operations has become increasingly available in technical documents and real-time dashboards of the utilities, affiliates, and the power grid operators. The research question that this study explores is: Can one combine high-wattage demand-side appliances with public information to launch cyberattacks on the power grid? To answer this question and report a proof of concept demonstration, the study scrapes data from public sources for Manhattan, NY using the electric vehicle charging station smartphone application and the power grid data circulated by the US Energy Information Administration, New York Independent System Operator, and the local utility in New York City. It then designs a novel data-driven cyberattack strategy using state-feedback based partial eigenvalue relocation, which targets frequency stability of the power grid. The study establishes that while such an attack is not possible at the current penetration level of PEVs, it will be practical once the number of PEVs increases.


title: The Simons Observatory: Astro2020 Decadal Project Whitepaper

abstract: The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.


title: Patterns of primes in the Sato-Tate conjecture

abstract: Fix a non-CM elliptic curve $E/\mathbb{Q}$, and let $a_E(p) = p + 1 - \#E(\mathbb{F}_p)$ denote the trace of Frobenius at $p$. The Sato-Tate conjecture gives the limiting distribution $\mu_{ST}$ of $a_E(p)/(2\sqrt{p})$ within $[-1, 1]$. We establish bounded gaps for primes in the context of this distribution. More precisely, given an interval $I\subseteq [-1, 1]$, let $p_{I,n}$ denote the $n$th prime such that $a_E(p)/(2\sqrt{p})\in I$. We show $\liminf_{n\to\infty}(p_{I,n+m}-p_{I,n}) < \infty$ for all $m\ge 1$ for "most" intervals, and in particular, for all $I$ with $\mu_{ST}(I)\ge 0.36$. Furthermore, we prove a common generalization of our bounded gap result with the Green-Tao theorem. To obtain these results, we demonstrate a Bombieri-Vinogradov type theorem for Sato-Tate primes.