title: Monolayer MoS$_2$ Strained to 1.3\% with a Microelectromechanical System

abstract: We report on a modified transfer technique for atomically thin materials integrated onto microelectromechanical systems (MEMS) for studying strain physics and creating strain-based devices. Our method tolerates the non-planar structures and fragility of MEMS, while still providing precise positioning and crack free transfer of flakes. Further, our method used the transfer polymer to anchor the 2D crystal to the MEMS, which reduces the fabrication time, increases the yield, and allowed us to exploit the strong mechanical coupling between 2D crystal and polymer to strain the atomically thin system. We successfully strained single atomic layers of molybdenum disulfide (MoS$_2$) with MEMS devices for the first time and achieved greater than 1.3\% strain, marking a major milestone for incorporating 2D materials with MEMS We used the established strain response of MoS$_2$ Raman and Photoluminescence spectra to deduce the strain in our crystals and provide a consistency check. We found good comparison between our experiment and literature.


title: Loynes construction for the extended bipartite matching

abstract: We propose an explicit construction of the stationary state of Extended Bipartite Matching (EBM) models, as defined in (Busic et. al., 2013). We use a Loynes-type backwards scheme similar in flavor to that in (Moyal et al., 2017), allowing to show the existence and uniqueness of a bi-infinite perfect matching under various conditions, for a large class of matching policies and of bipartite matching structures. The key algebraic element of our construction is the sub-additivity of a suitable stochastic recursive representation of the model, satisfied under most usual matching policies. By doing so, we also derive stability conditions for the system under general stationary ergodic assumptions, subsuming the classical markovian settings.


title: Back to the Future: The Case for Reversible Computing

abstract: There is one, and only one way, consistent with fundamental physics, that the efficiency of general digital computation can continue increasing indefinitely, and that is to apply the principles of reversible computing. We need to begin intensive development work on this technology soon if we want to maintain advances in computing and the attendant economic growth. NOTE: This paper is an extended author's preprint of the feature article titled "Throwing Computing Into Reverse" (print) or "The Future of Computing Depends on Making it Reversible" (online), published by IEEE Spectrum in Aug.-Sep. 2017. This preprint is based on the original draft manuscript that the author submitted to Spectrum, prior to IEEE edits and feedback from external readers.


title: Photonic quantum information processing: a review

abstract: Photonic quantum technologies represent a promising platform for several applications, ranging from long-distance communications to the simulation of complex phenomena. Indeed, the advantages offered by single photons do make them the candidate of choice for carrying quantum information in a broad variety of areas with a versatile approach. Furthermore, recent technological advances are now enabling first concrete applications of photonic quantum information processing. The goal of this manuscript is to provide the reader with a comprehensive review of the state of the art in this active field, with a due balance between theoretical, experimental and technological results. When more convenient, we will present significant achievements in tables or in schematic figures, in order to convey a global perspective of the several horizons that fall under the name of photonic quantum information.


title: Facebook (A)Live? Are live social broadcasts really broadcasts?

abstract: The era of live-broadcast is back but with two major changes. First, unlike traditional TV broadcasts, content is now streamed over the Internet enabling it to reach a wider audience. Second, due to various user-generated content platforms it has become possible for anyone to get involved, streaming their own content to the world. This emerging trend of going live usually happens via social platforms, where users perform live social broadcasts predominantly from their mobile devices, allowing their friends (and the general public) to engage with the stream in real-time. With the growing popularity of such platforms, the burden on the current Internet infrastructure is therefore expected to multiply. With this in mind, we explore one such prominent platform - Facebook Live. We gather 3TB of data, representing one month of global activity and explore the characteristics of live social broadcast. From this, we derive simple yet effective principles which can decrease the network burden. We then dissect global and hyper-local properties of the video while on-air, by capturing the geography of the broadcasters or the users who produce the video and the viewers or the users who interact with it. Finally, we study the social engagement while the video is live and distinguish the key aspects when the same video goes on-demand. A common theme throughout the paper is that, despite its name, many attributes of Facebook Live deviate from both the concepts of live and broadcast.


title: Lozenge tilings of hexagons with central holes and dents

abstract: Ciucu showed that the number of lozenge tilings of a hexagon in which a chain of equilateral triangles of alternating orientations, called a `\emph{fern}', has been removed in the center is given by a simple product formula (Adv. Math. 2017). In this paper, we present a multi-parameter generalization of this work by giving an explicit tiling enumeration for a hexagon with three ferns removed, besides the middle fern located in the center as in Ciucu's region, we remove two additional ferns from two sides of the hexagon. Our result also implies a counterpart of MacMahon's classical formula of boxed plane partitions, corresponding the \emph{exterior} of the union of three disjoint concave polygons obtained by turning 120 degrees after drawing each side.


title: Topology Learning of Radial Dynamical Systems with Latent Nodes

abstract: In this article, we present a method to reconstruct the topology of a partially observed radial network of linear dynamical systems with bi-directional interactions. Our approach exploits the structure of the inverse power spectral density matrix and recovers edges involving nodes up to four hops away in the underlying topology. We then present an algorithm with provable guarantees, which eliminates the spurious links obtained and also identifies the location of the unobserved nodes in the inferred topology. The algorithm recovers the exact topology of the network by using only time-series of the states at the observed nodes. The effectiveness of the method developed is demonstrated by applying it on a typical distribution system of the electric grid.


title: Demonstration of a novel method for measuring mass-loss rates for massive stars

abstract: The rate at which massive stars eject mass in stellar winds significantly influences their evolutionary path. Cosmic rates of nucleosynthesis, explosive stellar phenomena, and compact object genesis depend on this poorly known facet of stellar evolution. We employ an unexploited observational technique for measuring the mass-loss rates of O- and early-B stars. Our approach, which has no adjustable parameters, uses the principle of pressure equilibrium between the stellar wind and the ambient interstellar medium for a high-velocity star generating an infrared bowshock nebula. Results for twenty bowshock-generating stars show good agreement with two sets of theoretical predictions for O5--O9.5 main-sequence stars, yielding $\dot M=$1.3$\times$10$^{-6}$ to 2$\times$10$^{-9}$ solar masses per year. Although $\dot M$ values derived for this sample are smaller than theoretical expectations by a factor of about two, this discrepancy is greatly reduced compared to canonical mass-loss methods. Bowshock-derived mass-loss rates are factors of ten smaller than H$\alpha$-based measurements (uncorrected for clumping) for similar stellar types and are nearly an order of magnitude larger than P$^{4+}$ and some other UV absorption-line-based diagnostics. Ambient interstellar densities of at least several cm$^{-3}$ appear to be required for formation of a prominent infrared bowshock nebula. $\dot M$ measurements for early-B stars are not yet compelling owing to the small number in our sample and the lack of clear theoretical predictions in the regime of lower stellar luminosities. These results may constitute a partial resolution of the extant "weak-wind problem" for late-O stars. The technique shows promise for determining mass-loss rates in the weak-wind regime.


title: Complexity Growth Rate in Lovelock Gravity

abstract: Using the "Complexity = Action" framework we compute the late time growth of complexity for charged black holes in Lovelock gravity. Our calculation is facilitated by the fact that the null boundaries of the Wheeler-DeWitt patch do not contribute at late times and essential contributions coming from the joints are now understood arXiv:1803.00172. The late time growth rate reduces to a difference of internal energies associated with the inner and outer horizons, and in the limit where the mass is much larger than the charge, we reproduce the celebrated result of $2M/\pi$ with corrections proportional to the highest Lovelock coupling in even (boundary) dimensions. We find in some cases a minimum mass below which complexity remains effectively constant, even if the black hole contains a non-degenerate horizon.


title: Resource Polymorphism

abstract: We present a resource-management model for ML-style programming languages, designed to be compatible with the OCaml philosophy and runtime model. This is a proposal to extend the OCaml language with destructors, move semantics, and resource polymorphism, to improve its safety, efficiency, interoperability, and expressiveness. It builds on the ownership-and-borrowing models of systems programming languages (Cyclone, C++11, Rust) and on linear types in functional programming (Linear Lisp, Clean, Alms). It continues a synthesis of resources from systems programming and resources in linear logic initiated by Baker. It is a combination of many known and some new ideas. On the novel side, it highlights the good mathematical structure of Stroustrup's "Resource acquisition is initialisation" (RAII) idiom for resource management based on destructors, a notion sometimes confused with finalizers, and builds on it a notion of resource polymorphism, inspired by polarisation in proof theory, that mixes C++'s RAII and a tracing garbage collector (GC). The proposal targets a new spot in the design space, with an automatic and predictable resource-management model, at the same time based on lightweight and expressive language abstractions. It is backwards-compatible: current code is expected to run with the same performance, the new abstractions fully combine with the current ones, and it supports a resource-polymorphic extension of libraries. It does so with only a few additions to the runtime, and it integrates with the current GC implementation. It is also compatible with the upcoming multicore extension, and suggests that the Rust model for eliminating data-races applies. Interesting questions arise for a safe and practical type system, many of which have already been thoroughly investigated in the languages and prototypes Cyclone, Rust, and Alms.


title: Non-equilibrium scaling behaviour in driven soft biological assemblies

abstract: Measuring and quantifying non-equilibrium dynamics in active biological systems is a major challenge, because of their intrinsic stochastic nature and the limited number of variables accessible in any real experiment. We investigate what non-equilibrium information can be extracted from non-invasive measurements using a stochastic model of soft elastic networks with a heterogeneous distribution of activities, representing enzymatic force generation. In particular, we use this model to study how the non-equilibrium activity, detected by tracking two probes in the network, scales as a function of the distance between the probes. We quantify the non-equilibrium dynamics through the cycling frequencies, a simple measure of circulating currents in the phase space of the probes. We find that these cycling frequencies exhibit power-law scaling behavior with the distance between probes. In addition, we show that this scaling behavior governs the entropy production rate that can be recovered from the two traced probes. Our results provide insight in to how internal enzymatic driving generates non-equilibrium dynamics on different scales in soft biological assemblies.


title: Optimal Threshold-Based Control Policies for Persistent Monitoring on Graphs

abstract: We consider the optimal multi-agent persistent monitoring problem defined by a team of cooperating agents visiting a set of nodes (targets) on a graph with the objective of minimizing a measure of overall node state uncertainty. The solution to this problem involves agent trajectories defined both by the sequence of nodes to be visited by each agent and the amount of time spent at each node. Since such optimal trajectories are generally intractable, we propose a class of distributed threshold-based parametric controllers through which agent transitions from one node to the next are controlled by threshold parameters on the node uncertainty states. The resulting behavior of the agent-target system can be described by a hybrid dynamic system. This enables the use of Infinitesimal Perturbation Analysis (IPA) to determine on line (locally) optimal threshold parameters through gradient descent methods and thus obtain optimal controllers within this family of threshold-based policies. We further show that in a single-agent case the IPA gradient is monotonic, which implies a simple structure whereby an agent visiting a node should reduce the uncertainty state to zero before moving to the next node. Simulation examples are included to illustrate our results and compare them to optimal solutions derived through dynamic programming when this is possible.


title: Projective Hessian and Sasakian manifolds

abstract: The Hessian geometry is the real analogue of the K\"ahler one. Sasakian geometry is an odd-dimensional counterpart of the K\"ahler geometry. In the paper, we study the connection between projective Hessian and Sasakian manifolds analogous to the one between Hessian and K\"ahler manifolds. In particular, we construct a Sasakian structure on $TM\times \mathbb{R}$ from a projective Hessian structure on $M$. Especially, we are interested in the case of invariant structure on Lie groups. We define semi-Sasakian Lie groups as a generalization of Sasakian Lie groups. Then we construct a semi-Sasakian structure on a group $G\ltimes \mathbb{R}^{n+1}$ for a projective Hessian Lie group $G$. Further, we describe examples of homogeneous Hessian Lie groups and corresponding semi-Sasakian Lie groups. The big class of projective Hessian Lie groups can be constructed by homogeneous regular domains in $\mathbb{R}^n$. The groups $\text{SO}(2)$ and $\text{SU}(2)$ belong to another kind of examples. Using them, we construct semi-Sasakian structures on the group of the Euclidean motions of the real plane and the group of isometries of the complex plane.


title: Long-branch attraction in species tree estimation: inconsistency of partitioned likelihood and topology-based summary methods

abstract: With advances in sequencing technologies, there are now massive amounts of genomic data from across all life, leading to the possibility that a robust Tree of Life can be constructed. However, "gene tree heterogeneity", which is when different genomic regions can evolve differently, is a common phenomenon in multi-locus datasets, and reduces the accuracy of standard methods for species tree estimation that do not take this heterogeneity into account. New methods have been developed for species tree estimation that specifically address gene tree heterogeneity, and that have been proven to converge to the true species tree when the number of loci and number of sites per locus both increase (i.e., the methods are said to be "statistically consistent"). Yet, little is known about the biologically realistic condition where the number of sites per locus is bounded. We show that when the sequence length of each locus is bounded (by any arbitrarily chosen value), the most common approaches to species tree estimation that take heterogeneity into account (i.e., traditional fully partitioned concatenated maximum likelihood and newer approaches, called summary methods, that estimate the species tree by combining gene trees) are not statistically consistent, even when the heterogeneity is extremely constrained. The main challenge is the presence of conditions such as long branch attraction that create biased tree estimation when the number of sites is restricted. Hence, our study uncovers a fundamental challenge to species tree estimation using both traditional and new methods.


title: Re-entrant bimodality in spheroidal chiral swimmers in shear flow

abstract: We use a continuum model to report on the behavior of a dilute suspension of chiral swimmers subject to externally imposed shear in a planar channel. Swimmer orientation in response to the imposed shear can be characterized by two distinct phases of behavior, corresponding to unimodal or bimodal distribution functions for swimmer orientation along the channel. These phases indicate the occurrence (or not) of a population splitting phenomenon changing the swimming direction of a macroscopic fraction of active particles to the exact opposite of that dictated by the imposed flow. We present a detailed quantitative analysis elucidating the complexities added to the population splitting behavior of swimmers when they are chiral. In particular, the transition from unimodal to bimodal and vice versa are shown to display a re-entrant behavior across the parameter space spanned by varying the chiral angular speed. We also present the notable effects of particle aspect ratio and self-propulsion speed on system phase behavior and discuss potential implications of our results in applications such as swimmer separation/sorting.

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


title: Realistic atomistic structure of amorphous silicon from machine-learning-driven molecular dynamics

abstract: Amorphous silicon (a-Si) is a widely studied non-crystalline material, and yet the subtle details of its atomistic structure are still unclear. Here, we show that accurate structural models of a-Si can be obtained by harnessing the power of machine-learning algorithms to create interatomic potentials. Our best a-Si network is obtained by cooling from the melt in molecular-dynamics simulations, at a rate of 10$^{11}$ K/s (that is, on the 10 ns timescale). This structure shows a defect concentration of below 2% and agrees with experiments regarding excess energies, diffraction data, as well as $^{29}$Si solid-state NMR chemical shifts. We show that this level of quality is impossible to achieve with faster quench simulations. We then generate a 4,096-atom system which correctly reproduces the magnitude of the first sharp diffraction peak (FSDP) in the structure factor, achieving the closest agreement with experiments to date. Our study demonstrates the broader impact of machine-learning interatomic potentials for elucidating accurate structures and properties of amorphous functional materials.


title: Neutron star tidal deformabilities constrained by nuclear theory and experiment

abstract: We confront observational data from gravitational wave event GW170817 with microscopic modeling of the cold neutron star equation of state. We develop and employ a Bayesian statistical framework that enables us to implement constraints on the equation of state from laboratory measurements of nuclei and state-of-the-art chiral effective field theory methods. The energy density functionals constructed from the posterior probability distributions are then used to compute consistently the neutron star equation of state from the outer crust to the inner core, assuming a composition consisting of protons, neutrons, electrons, and muons. In contrast to previous studies, we find that the 95% credibility range of predicted neutron star tidal deformabilities ($136 < \Lambda < 519$) for a 1.4 solar-mass neutron star is already consistent with the upper bound deduced from observations of the GW170817 event. However, we find that lower bounds on the neutron star tidal deformability will very strongly constrain microscopic models of the dense matter equation of state. We also demonstrate a strong correlation between the neutron star tidal deformability and the pressure of beta-equilibrated matter at twice saturation density.


title: Severely Constraining Dark Matter Interpretations of the 21-cm Anomaly

abstract: The EDGES Collaboration has recently reported the detection of a stronger-than-expected absorption feature in the global 21-cm spectrum, centered at a frequency corresponding to a redshift of z ~ 17. This observation has been interpreted as evidence that the gas was cooled during this era as a result of scattering with dark matter. In this study, we explore this possibility, applying constraints from the cosmic microwave background, light element abundances, Supernova 1987A, and a variety of laboratory experiments. After taking these constraints into account, we find that the vast majority of the parameter space capable of generating the observed 21-cm signal is ruled out. The only range of models that remains viable is that in which a small fraction, ~ 0.3-2%, of the dark matter consists of particles with a mass of ~ 10-80 MeV and which couple to the photon through a small electric charge, epsilon ~ 10^{-6}-10^{-4}. Furthermore, in order to avoid being overproduced in the early universe, such models must be supplemented with an additional depletion mechanism, such as annihilations through a L_{\mu}-L_{\tau} gauge boson or annihilations to a pair of rapidly decaying hidden sector scalars.

categories: hep-ph astro-ph.CO astro-ph.GA astro-ph.HE astro-ph.IM


title: Rotomagnetic coupling in fine grained multiferroic BiFeO3: Theory and experiment

abstract: Using Landau-Ginzburg-Devonshire (LGD) theory for BiFeO3 dense fine grained ceramics with quasi spherical grains and nanosized inter grain spaces enriched by elastic defects, we calculated a surprisingly strong size-induced increase of the AFM temperature caused by the joint action of rotomagnetic and magnetostrictive coupling. Notably that all parameters included in the LGD functional have been extracted from experiments, not assumed. Complementary we performed experiments for dense BiFeO3 ceramics, which revealed that the shift of antiferromagnetic transition to 690 K instead of 645 K for a single crystal. To explain theoretically the result, we consider the possibility to control antiferromagnetic state of multiferroic BiFeO3 via biquadratic antiferrodistortive rotomagnetic, rotoelectric, magnetostrictive and magnetoelectric couplings. According to our calculations the highest is the rotostriction contribution, the magnetostrictive and electrostriction contributions appeared smaller.


title: Spitzer Infrared Spectrograph Observations of the Galactic Center: Quantifying the Extreme Ultraviolet/so ft X-ray fluxes

abstract: It has long been shown that the extreme ultraviolet spectrum of the ionizing stars of H II regions can be estimated by comparing the observed line emission to detailed models. In the Galactic Center (GC), however, previous observations have shown that the ionizing spectral energy distribution (SED) of the local photon field is strange, producing both very low excitation ionized gas (indicative of ionization by late O stars) and also widespread diffuse emission from atoms too highly ionized to be found in normal H II regions. This paper describes the analysis of all the GC spectra taken by Spitzer's Infrared Spectrograph and downloaded from the Spitzer Heritage Archive. In it, H II region densities and abundances are described, and serendipitously discovered candidate planetary nebulae, compact shocks, and candidate young stellar objects are tabulated. Models were computed with Cloudy, using SEDs from Starburst99 plus additional X-rays, and compared to the observed mid-infrared forbidden and recombination lines. The ages inferred from the model fits do not agree with recent proposed star formation sequences (star formation in the GC occurring along streams of gas with density enhancements caused by close encounters with the black hole, Sgr A*), with Sgr~B1, Sgr~C, and the Arches Cluster being all about the same age, around 4.5 Myr old, with similar X-ray requirements. The fits for the Quintuplet Cluster appear to give a younger age, but that could be caused by higher-energy photons from shocks from stellar winds or from a supernova.


title: Flexible and Efficient Algorithms for Abelian Matching in Strings

abstract: The abelian pattern matching problem consists in finding all substrings of a text which are permutations of a given pattern. This problem finds application in many areas and can be solved in linear time by a naive sliding window approach. In this short communication we present a new class of algorithms based on a new efficient fingerprint computation approach, called Heap-Counting, which turns out to be fast, flexible and easy to be implemented. It can be proved that our solutions have a linear worst case time complexity and, in addition, we present an extensive experimental evaluation which shows that our newly presented algorithms are among the most efficient and flexible solutions in practice for the abelian matching problem in strings.


title: OntoWind: An Improved and Extended Wind Energy Ontology

abstract: Ontologies are critical sources of semantic information for many application domains. Hence, there are ontologies proposed and utilized for domains such as medicine, chemical engineering, and electrical energy. In this paper, we present an improved and extended version of a wind energy ontology previously proposed. First, the ontology is restructured to increase its understandability and coverage. Secondly, it is enriched with new concepts, crisp/fuzzy attributes, and instances to increase its usability in semantic applications regarding wind energy. The ultimate ontology is utilized within a Web-based semantic portal application for wind energy, in order to showcase its contribution in a genuine application. Hence, the current study is a significant to wind and thereby renewable energy informatics, with the presented publicly-available wind energy ontology and the implemented proof-of-concept system.


title: The size of the giant component in random hypergraphs: a short proof

abstract: We consider connected components in $k$-uniform hypergraphs for the following notion of connectedness: given integers $k\ge 2$ and $1\le j \le k-1$, two $j$-sets (of vertices) lie in the same $j$-component if there is a sequence of edges from one to the other such that consecutive edges intersect in at least $j$ vertices. We prove that certain collections of $j$-sets constructed during a breadth-first search process on $j$-components in a random $k$-uniform hypergraph are reasonably regularly distributed with high probability. We use this property to provide a short proof of the asymptotic size of the giant $j$-component shortly after it appears.


title: Kernel and cokernel in the category of augmented involutive stereotype algebras

abstract: We prove several properties of kernels and cokernels in the category of augmented involutive stereotype algebras: 1) the morphisms of the augmented involutive stereotype algebras have kernels and cokernels, 2) the cokernel is preserved under the passage to the group stereotype algebras, and 3) the notion of cokernel allows to prove that the continuous envelope $\operatorname{Env} {\mathcal C}^\star(G)$ of the group algebra ${\mathcal C}^\star(G)$ is an involutive Hopf algebra in the category of stereotype spaces $({\tt Ste},\odot)$, if $G$ has the form $Z\cdot K$, where $Z$ is a commutative locally compact group, and $K$ a compact group. The last result plays an important role in the generalization of the Pontryagin duality for arbitrary Moore groups.


title: Accelerated Methods for Deep Reinforcement Learning

abstract: Deep reinforcement learning (RL) has achieved many recent successes, yet experiment turn-around time remains a key bottleneck in research and in practice. We investigate how to optimize existing deep RL algorithms for modern computers, specifically for a combination of CPUs and GPUs. We confirm that both policy gradient and Q-value learning algorithms can be adapted to learn using many parallel simulator instances. We further find it possible to train using batch sizes considerably larger than are standard, without negatively affecting sample complexity or final performance. We leverage these facts to build a unified framework for parallelization that dramatically hastens experiments in both classes of algorithm. All neural network computations use GPUs, accelerating both data collection and training. Our results include using an entire DGX-1 to learn successful strategies in Atari games in mere minutes, using both synchronous and asynchronous algorithms.


title: On continuous duality for Moore groups

abstract: In this paper we correct the errors of Yu.~N.~Kuznetsova's paper on the continuous duality for Moore groups.


title: New Experimental Limits on Exotic Spin-Spin-Velocity-Dependent Interactions By Using SmCo$_5$ Spin Sources

abstract: We report the latest results of searching for possible new macro-scale spin-spin-velocity-dependent forces (SSVDFs) based on specially designed iron-shielded SmCo$_5$ (ISSC) spin sources and a spin exchange relaxation free (SERF) co-magnetometer. The ISSCs have high net electron spin densities of about $1.7\times 10^{21}$ cm$^{-3}$, which mean high detecting sensitivity; and low magnetic field leakage of about $\sim$mG level due to iron shielding, which means low detecting noise. With help from the ISSCs, the high sensitivity SERF co-magnetometer, and the similarity analysis method, new constraints on SSVDFs with forms of $V_{6+7}$, $V_8$, $V_{15}$, and $V_{16}$ have been obtained, which represent the tightest limits in force range of 5 cm -- 1 km to the best of our knowledge.


title: Addendum to Pontryagin's maximum principle for dynamic systems on time scales

abstract: This note is an addendum to [1,2], pointing out the differences between these papers and raising open questions.


title: Sever: A Robust Meta-Algorithm for Stochastic Optimization

abstract: In high dimensions, most machine learning methods are brittle to even a small fraction of structured outliers. To address this, we introduce a new meta-algorithm that can take in a base learner such as least squares or stochastic gradient descent, and harden the learner to be resistant to outliers. Our method, Sever, possesses strong theoretical guarantees yet is also highly scalable -- beyond running the base learner itself, it only requires computing the top singular vector of a certain $n \times d$ matrix. We apply Sever on a drug design dataset and a spam classification dataset, and find that in both cases it has substantially greater robustness than several baselines. On the spam dataset, with $1\%$ corruptions, we achieved $7.4\%$ test error, compared to $13.4\%-20.5\%$ for the baselines, and $3\%$ error on the uncorrupted dataset. Similarly, on the drug design dataset, with $10\%$ corruptions, we achieved $1.42$ mean-squared error test error, compared to $1.51$-$2.33$ for the baselines, and $1.23$ error on the uncorrupted dataset.


title: Shedding Light on the Dark Corners of the Internet: A Survey of Tor Research

abstract: Anonymity services have seen high growth rates with increased usage in the past few years. Among various services, Tor is one of the most popular peer-to-peer anonymizing service. In this survey paper, we summarize, analyze, classify and quantify 26 years of research on the Tor network. Our research shows that `security' and `anonymity' are the most frequent keywords associated with Tor research studies. Quantitative analysis shows that the majority of research studies on Tor focus on `deanonymization' the design of a breaching strategy. The second most frequent topic is analysis of path selection algorithms to select more resilient paths. Analysis shows that the majority of experimental studies derived their results by deploying private testbeds while others performed simulations by developing custom simulators. No consistent parameters have been used for Tor performance analysis. The majority of authors performed throughput and latency analysis.


title: Stochastic nonlinear Schr\"odinger equations on tori

abstract: We consider the stochastic nonlinear Schr\"odinger equations (SNLS) posed on $d$-dimensional tori with either additive or multiplicative stochastic forcing. In particular, for the one-dimensional cubic SNLS, we prove global well-posedness in $L^2(\mathbb{T})$. As for other power-type nonlinearities, namely (i) (super)quintic when $d = 1$ and (ii) (super)cubic when $d \geq 2$, we prove local well-posedness in all scaling-subcritical Sobolev spaces and global well-posedness in the energy space for the defocusing, energy-subcritical problems.


title: Path Planning and Navigation Inside Off-World Lava Tubes and Caves

abstract: Detailed surface images of the Moon and Mars reveal hundreds of cave-like openings. These cave-like openings are theorized to be remnants of lava-tubes and their interior maybe in pristine conditions. These locations may have well preserved geological records of the Moon and Mars, including evidence of past water flow and habitability. Exploration of these caves using wheeled rovers remains a daunting challenge. These caves are likely to have entrances with caved-in ceilings much like the lava-tubes of Arizona and New Mexico. Thus, the entrances are nearly impossible to traverse even for experienced human hikers. Our approach is to utilize the SphereX robot, a 3 kg, 30 cm diameter robot with computer hardware and sensors of a smartphone attached to rocket thrusters. Each SphereX robot can hop, roll or fly short distances in low gravity, airless or low-pressure environments. Several SphereX robots maybe deployed to minimize single-point failure and exploit cooperative behaviors to traverse the cave. There are some important challenges for navigation and path planning in these cave environments. Localization systems such as GPS are not available nor are they easy to install due to the signal blockage from the rocks. These caves are too dark and too large for conventional sensor such as cameras and miniature laser sensors to perform detailed mapping and navigation. In this paper, we identify new techniques to map these caves by performing localized, cooperative mapping and navigation.


title: Correlated emission and spin-down variability in radio pulsars

abstract: The recent revelation that there are correlated period derivative and pulse shape changes in pulsars has dramatically changed our understanding of timing noise as well as the relationship between the radio emission and the properties of the magnetosphere as a whole. Using Gaussian processes we are able to model timing and emission variability using a regression technique that imposes no functional form on the data. We revisit the pulsars first studied by Lyne et al. (2010). We not only confirm the emission and rotational transitions revealed therein, but reveal further transitions and periodicities in 8 years of extended monitoring. We also show that in many of these objects the pulse profile transitions between two well-defined shapes, coincident with changes to the period derivative. With a view to the SKA and other telescopes capable of higher cadence we also study the detection limitations of period derivative changes.


title: Imitate or innovate: competition of strategy updating attitudes in spatial social dilemma games

abstract: Evolution is based on the assumption that competing players update their strategies to increase their individual payoffs. However, while the applied updating method can be different, most of previous works proposed uniform models where players use identical way to revise their strategies. In this work we explore how imitation-based or learning attitude and innovation-based or myopic best response attitude compete for space in a complex model where both attitudes are available. In the absence of additional cost the best response trait practically dominates the whole snow-drift game parameter space which is in agreement with the average payoff difference of basic models. When additional cost is involved then the imitation attitude can gradually invade the whole parameter space but this transition happens in a highly nontrivial way. However, the role of competing attitudes is reversed in the stag-hunt parameter space where imitation is more successful in general. Interestingly, a four-state solution can be observed for the latter game which is a consequence of an emerging cyclic dominance between possible states. These phenomena can be understood by analyzing the microscopic invasion processes, which reveals the unequal propagation velocities of strategies and attitudes.


title: Winds as the origin of radio emission in $z=2.5$ radio-quiet extremely red quasars

abstract: Most active galactic nuclei (AGNs) are radio-quiet, and the origin of their radio emission is not well-understood. One hypothesis is that this radio emission is a by-product of quasar-driven winds. In this paper, we present the radio properties of 108 extremely red quasars (ERQs) at $z=2-4$. ERQs are among the most luminous quasars ($L_{bol} \sim 10^{47-48}$ erg/s) in the Universe, with signatures of extreme ($\gg 1000$ km/s) outflows in their [OIII]$\lambda$5007 \AA\ emission, making them the best subjects to seek the connection between radio and outflow activity. All ERQs but one are unresolved in the radio on $\sim 10$ kpc scales, and the median radio luminosity of ERQs is $\nu L_\nu [{\rm 6\,GHz}] = 10^{41.0}$ erg/s, in the radio-quiet regime, but one to two orders of magnitude higher than that of other quasar samples. The radio spectra are steep, with a mean spectral index $\langle \alpha \rangle = -1.0$. In addition, ERQs neatly follow the extrapolation of the low-redshift correlation between radio luminosity and the velocity dispersion of [OIII]-emitting ionized gas. Uncollimated winds, with a power of one per cent of the bolometric luminosity, can account for all these observations. Such winds would interact with and shock the gas around the quasar and in the host galaxy, resulting in acceleration of relativistic particles and the consequent synchrotron emission observed in the radio. Our observations support the picture in which ERQs are signposts of extremely powerful episodes of quasar feedback, and quasar-driven winds as a contributor of the radio emission in the intermediate regime of radio luminosity $\nu L_\nu = 10^{39}-10^{42}$ erg/s.


title: De Broglie relations, Gravitational time dilation and weak equivalence principle

abstract: Interplays between quantum physics and gravity have long inspired exciting studies, which reveal subtle connections between quantum laws and the general notion of curved spacetime. One important example is the uniqueness of free-falling motions in both quantum and gravitational physics. In this work, we re-investigate the free-falling motions of quantum test wave packets that are distributed over weakly curved spacetime backgrounds. Except for the de Broglie relations, no assumption of priori given Hamiltonians or least actions satisfied by the quantum system is made. We find that the mean motions of quantum test wave packets can be deduced naturally from the de Broglie relations with a generalized treatment of gravitational time dilations in quantum waves. Such mean motions of quantum test masses in gravitational field are independent of their masses and compositions, and restore exactly the free-falling or geodesic motions of classical test masses in curved spacetime. This suggests a novel perspective that weak equivalence principle, which states the universality of free-fall, may be deeply rooted in quantum physics and be a phenomenon that has emerged from the quantum world.


title: A unified and improved Chebotarev density theorem

abstract: We establish an unconditional effective Chebotarev density theorem that improves uniformly over the well-known result of Lagarias and Odlyzko. As a consequence, we give a new asymptotic form of the Chebotarev density theorem that can count much smaller primes with arbitrary log-power savings, even in the case where a Landau-Siegel zero is present. Our main theorem interpolates the strongest unconditional upper bound for the least prime ideal with a given Artin symbol as well as the Chebotarev analogue of the Brun-Titchmarsh theorem proved by the authors. We also present a new application of our main result that exhibits considerable gains over earlier versions of the Chebotarev density theorem. If $f$ is a positive definite primitive binary quadratic form then we count lattice points $(u,v) \in \mathbb{Z}^2$ such that $f(u,v)$ is prime and $u, v$ have no prime factors $\leq z$ with uniformity in $z$ and the discriminant of $f$.


title: A Detection of the Environmental Dependence of the Sizes and Stellar Haloes of Massive Central Galaxies

abstract: We use ~100 square deg of deep (>28.5 mag arcsec$^{-2}$ in i-band), high-quality (median 0.6 arcsec seeing) imaging data from the Hyper Suprime-Cam (HSC) survey to reveal the halo mass dependence of the surface mass density profiles and outer stellar envelopes of massive galaxies. The i-band images from the HSC survey reach ~4 magnitudes deeper than Sloan Digital Sky Survey and enable us to directly trace stellar mass distributions to 100 kpc without requiring stacking. We conclusively show that, at fixed stellar mass, the stellar profiles of massive galaxies depend on the masses of their dark matter haloes. On average, massive central galaxies with $\log M_{\star, 100\ \mathrm{kpc}}>11.6$ in more massive haloes at 0.3 < z < 0.5 have shallower inner stellar mass density profiles (within ~10-20 kpc) and more prominent outer envelopes. These differences translate into a halo mass dependence of the mass-size relation. Central galaxies in haloes with $\log M_{\rm{Halo}}>14.0$ are ~20% larger in $R_{\mathrm{50}}$ at fixed stellar mass. Such dependence is also reflected in the relationship between the stellar mass within 10 and 100 kpc. Comparing to the mass--size relation, the $\log M_{\star, 100\ \rm{kpc}}$-$\log M_{\star, 10\ \rm{kpc}}$ relation avoids the ambiguity in the definition of size, and can be straightforwardly compared with simulations. Our results demonstrate that, with deep images from HSC, we can quantify the connection between halo mass and the outer stellar halo, which may provide new constraints on the formation and assembly of massive central galaxies.


title: A novel quantitative indicator of the left ventricular contraction based on volume changes of the left ventricular myocardial segments

abstract: Ejection fraction (EF) is commonly measured by echocardiography, by dividing the volume ejected by the heart (stroke volume) by the volume of the filled heart (end-diastolic volume). Utilizing volume changes of left myocardial segments per a cardiac cycle, physical laws and mathematical equations specific echocardiographic data, this paper serves to generalize EF by a novel parameter over the time that it can make available, more detailed, valuable and practical information to fully describe the left ventricular (LV) contractility function.Patients who underwent clinically-directed standard transthoracic echocardiography using 2D conventional echocardiography machines armed to measuring strain components, were asked to estimate displacements and longitudinal, radial and circumferential strains for each LV echocardiographic segments per a cardiac cycle. Volume fractional changes of the LV echocardiographic segments are expanded based on their strain components over the time. Ejected blood volume fraction induced by a left myocardial sample, is computed within a cardiac cycle. Total fraction of the ejected blood volume in the left ventricular cavity was obtained by integrating over the times and LV myocardial segments. EF is an especial value of this total fraction at the end systolic time. The common measurement of EF is only based on LV cavity volumes at the end diastolic and systolic phases. These findings lead to determine detailed aspects of the left ventricular contraction. This generalized parameter has important implications to give the real value of EF in the sever Mitral valve regurgitations.


title: Mori-Zwanzig reduced models for uncertainty quantification

abstract: In many time-dependent problems of practical interest the parameters and/or initial conditions entering the equations describing the evolution of the various quantities exhibit uncertainty. One way to address the problem of how this uncertainty impacts the solution is to expand the solution using polynomial chaos expansions and obtain a system of differential equations for the evolution of the expansion coefficients. We present an application of the Mori-Zwanzig (MZ) formalism to the problem of constructing reduced models of such systems of differential equations. In particular, we construct reduced models for a subset of the polynomial chaos expansion coefficients that are needed for a full description of the uncertainty caused by uncertain parameters or initial conditions. Even though the MZ formalism is exact, its straightforward application to the problem of constructing reduced models for estimating uncertainty involves the computation of memory terms whose cost can become prohibitively expensive. For those cases, we present a Markovian reformulation of the MZ formalism which can lead to approximations that can alleviate some of the computational expense while retaining an accuracy advantage over reduced models that discard the memory altogether. Our results support the conclusion that successful reduced models need to include memory effects.


title: Approaches to the Optimal Nonlinear Analysis of Microcalorimeter Pulses

abstract: We consider how to analyze microcalorimeter pulses for quantities that are nonlinear in the data, while preserving the signal-to-noise advantages of lin- ear optimal filtering. We successfully apply our chosen approach to compute the electrothermal feedback energy deficit (the "Joule energy") of a pulse, which has been proposed as a linear estimator of the deposited photon energy.


title: Orbital Period Ratios and Fibonacci Numbers in Solar Planetary and Satellite Systems and in Exoplanetary Systems

abstract: It is shown that orbital period ratios of successive secondaries in the Solar planetary and giant satellite systems and in exoplanetary systems are preferentially closer to irreducible fractions formed with Fibonacci numbers between 1 and 8 than to other fractions, in a ratio of approximately 60% to 40%. Furthermore, if sets of minor planets are chosen with gradually smaller inclinations and eccentricities, the proximity to Fibonacci fractions of their period ratios with Jupiter or Mars$'$ period tends to increase. Finally, a simple model explains why the resonances with ratios of orbital periods $P_{1}$ and $P_{2}$ of successive secondaries being equal to ratio of small integers $p$ and $(p+q)$, $P_{1}/P_{2}=p/(p+q)$, are stronger and more commonly observed.


title: CSO CO(2-1) and Spitzer IRAC Observations of a Bipolar Outflow in High-mass Star-forming Region IRAS 22506+5944

abstract: We present Caltech Submillimeter Observatory CO (2-1) and Spitzer IRAC observations toward IRAS 22506+5944, which is a 10$^4$ $L_{\odot}$ massive star-forming region. The CO (2-1) maps show an east-west bipolar molecular outflow originating from the 3 mm dust continuum peak. The Spitzer IRAC color-composite image reveals a pair of bow-shaped tips which are prominent in excess 4.5 $\mu$m emission and are located to the leading fronts of the bipolar outflow, providing compelling evidence for the existence of bow-shocks as the driving agents of the molecular outflow. By comparing our CO (2-1) observations with previously published CO (1-0) data, we find that the CO (2-1)/(1-0) line ratio increases from low (~5 km s$^{-1}$) to moderate (~8 - 12 km s$^{-1}$) velocities, and then decreases at higher velocities. This is qualitatively consistent with the scenario that the molecular outflow is driven by multiple bow-shocks. We also revisit the position-velocity diagram of the CO (1-0) data, and find two spur structures along the outflow axis, which are further evidence for the presence of multiple jet bow-shocks. Finally, power-law fittings to the mass spectrum of the outflow gives power law indexes more consistent with the jet bow-shock model than the wide-angle wind model.


title: Charting the AdS Islands of Stability with Multi-oscillators?

abstract: We propose the existence of an infinite-parameter family of solutions in AdS that oscillate on any number of non-commensurate frequencies. Some of these solutions appear stable when perturbed, and we suggest that they can be used to map out the AdS "islands of stability". By numerically constructing two-frequency solutions and exploring their parameter space, we find that both collapse and non-collapse are generic scenarios near AdS. Unlike other approaches, our results are valid on any timescale and do not rely on perturbation theory.


title: Radio transients investigation with VLBI

abstract: The technique of Very Long Baseline Interferometry (VLBI) can provide accurate localization and unique physical information about radio transients. However, it is still underutilized due to the inherent difficulties of VLBI data analysis and practical difficulties of organizing observations on short notice. We present a brief overview of the currently available VLBI arrays and observing strategies used to study long- and short-duration radio transients.


title: Kuiper Belt Analogues in Nearby M-type Planet-host Systems

abstract: We present the results of a Herschel survey of 21 late-type stars that host planets discovered by the radial velocity technique. The aims were to discover new disks in these systems and to search for any correlation between planet presence and disk properties. In addition to the known disk around GJ 581, we report the discovery of two new disks, in the GJ 433 and GJ 649 systems. Our sample therefore yields a disk detection rate of 14%, higher than the detection rate of 1.2% among our control sample of DEBRIS M-type stars with 98% confidence. Further analysis however shows that the disk sensitivity in the control sample is about a factor of two lower in fractional luminosity than for our survey, lowering the significance of any correlation between planet presence and disk brightness below 98%. In terms of their specific architectures, the disk around GJ 433 lies at a radius somewhere between 1 and 30au. The disk around GJ 649 lies somewhere between 6 and 30au, but is marginally resolved and appears more consistent with an edge-on inclination. In both cases the disks probably lie well beyond where the known planets reside (0.06-1.1au), but the lack of radial velocity sensitivity at larger separations allows for unseen Saturn-mass planets to orbit out to $\sim$5au, and more massive planets beyond 5au. The layout of these M-type systems appears similar to Sun-like star + disk systems with low-mass planets.


title: Half-Metallic Superconducting Triplet Spin MultiValves

abstract: We study spin switching effects in finite-size superconducting multivalve structures. We examine $\rm F_1F_2SF_3$ and $\rm F_1F_2SF_3F_4$ hybrids where a singlet superconductor ($\rm S$) layer is sandwiched among ferromagnet ($\rm F$) layers with differing thicknesses and magnetization orientations. Our results reveal a considerable number of experimentally viable spin valve configurations that lead to on-off switching of the superconducting state. For $\rm S$ widths on the order of the superconducting coherence length $\xi_0$, non-collinear magnetization orientations in adjacent $\rm F$ layers with multiple spin-axes leads to a rich variety of triplet spin-valve effects. Motivated by recent experiments, we focus on samples where magnetizations in the $\rm F_1$ and $\rm F_4$ layers exist in a fully spin-polarized half metallic phase, and calculate the superconducting transition temperature, spatially and energy resolved density of states, and the spin-singlet and spin-triplet superconducting correlations. Our findings demonstrate that superconductivity in these devices can be completely switched on or off over a wide range of magnetization misalignment angles due to the generation of equal-spin and opposite-spin triplet pairings.


title: Fundamental limits to quantum channel discrimination

abstract: What is the ultimate performance for discriminating two arbitrary quantum channels acting on a finite-dimensional Hilbert space? Here we address this basic question by deriving a general and fundamental lower bound. More precisely, we investigate the symmetric discrimination of two arbitrary qudit channels by means of the most general protocols based on adaptive (feedback-assisted) quantum operations. In this general scenario, we first show how port-based teleportation can be used to simplify these adaptive protocols into a much simpler non-adaptive form, designing a new type of teleportation stretching. Then, we prove that the minimum error probability affecting the channel discrimination cannot beat a bound determined by the Choi matrices of the channels, establishing a general, yet computable formula for quantum hypothesis testing. As a consequence of this bound, we derive ultimate limits and no-go theorems for adaptive quantum illumination and single-photon quantum optical resolution. Finally, we show how the methodology can also be applied to other tasks, such as quantum metrology, quantum communication and secret key generation.


title: A new solar neutrino channel for grand-unification monopole searches

abstract: We identify a previously untapped discovery channel for grand-unification monopoles, arising from their ability to catalyse the direct decay of protons into monoenergetic 459 MeV antineutrinos within the Sun. Previous analyses omit this possibility as it necessarily involves an electroweak suppression factor, and instead search for the unsuppressed 20-50 MeV neutrinos produced via two-stage proton decays. By accounting for the relative difference in interaction cross section and experimental background at typical neutrino detection experiments, we demonstrate that this new channel in fact possesses greater discovery potential. As a case in point, using 5326 live days of Super-Kamiokande exposure we find that $2\;\sigma$ ($3\;\sigma$) deviations in the 20-50 MeV channel are amplified to $3\;\sigma$ ($4.6\;\sigma$) deviations in the 459 MeV case. Exploiting correlations between these two channels may also offer even greater statistical power.


title: Quantifying quantum invasiveness

abstract: We propose a resource theory of the quantum invasiveness of general quantum operations, i.e., those defined by quantum channels in Leggett-Garg scenarios. We are then able to compare the resource-theoretic framework of quantum invasiveness to the resource theory of coherence. We also show that the Fisher information is a quantifier of quantum invasiveness. This result allows us to establish a direct connection between the concept of quantum invasiveness and quantum metrology, by exploring the utility of the definition of quantum invasiveness in the context of metrological protocols.


title: Primordial Black Holes from String Inflation

abstract: We present a single-field string inflationary model which allows for the generation of primordial black holes in the low mass region where they can account for a significant fraction of the dark matter abundance. The potential is typical of type IIB fibre inflation models and features a plateau at CMB scales and a near inflection point at large distance scales where the power spectrum is enhanced due to a period of ultra slow-roll. The tunability of the underlying parameters is guaranteed by scanning through the string landscape and their stability against quantum corrections is ensured by an effective shift symmetry.


title: Beyond many-body localized states in a spin-disordered Hubbard model with pseudo-spin symmetry

abstract: A prime characterization of many-body localized (MBL) systems is the entanglement of their eigenstates; in contrast to the typical ergodic phase whose eigenstates are volume law, MBL eigenstates obey an area law. In this work, we show that a spin-disordered Hubbard model has both a large number of area-law eigenstates as well as a large number of eigenstates whose entanglement scales logarithmically with system size (log-law). This model, then, is a microscopic Hamiltonian which is neither ergodic nor many-body localized. We establish these results through a combination of analytic arguments based on the eta-pairing operators combined with a numerical analysis of eigenstates. In addition, we describe and simulate a dynamic time evolution approach starting from product states through which one can separately probe the area law and log-law eigenstates in this system.


title: The emergent algebraic structure of RNNs and embeddings in NLP

abstract: We examine the algebraic and geometric properties of a uni-directional GRU and word embeddings trained end-to-end on a text classification task. A hyperparameter search over word embedding dimension, GRU hidden dimension, and a linear combination of the GRU outputs is performed. We conclude that words naturally embed themselves in a Lie group and that RNNs form a nonlinear representation of the group. Appealing to these results, we propose a novel class of recurrent-like neural networks and a word embedding scheme.


title: Interstellar object 'Oumuamua as an extinct fragment of an ejected cometary planetesimal

abstract: 'Oumuamua was discovered passing through our Solar System on a hyperbolic orbit. It presents an apparent contradiction, with colors similar to those of volatile-rich Solar System bodies but with no visible outgassing or activity during its close approach to the Sun. Here we show that this contradiction can be explained by the dynamics of planetesimal ejection by giant planets. We propose that 'Oumuamua is an extinct fragment of a comet-like planetesimal born in a planet-forming disk that also formed Neptune- to Jupiter-mass giant planets. On its pathway to ejection 'Oumuamua's parent body underwent a close encounter with a giant planet and was tidally disrupted into small pieces, similar to comet Shoemaker-Levy 9's disruption after passing close to Jupiter. We use dynamical simulations to show that 0.1-1% of cometary planetesimals undergo disruptive encounters prior to ejection. Rocky asteroidal planetesimals are unlikely to disrupt due to their higher densities. After disruption, the bulk of fragments undergo enough close passages to their host stars to lose their surface volatiles and become extinct. Planetesimal fragments such as 'Oumuamua contain little of the mass in the population of interstellar objects but dominate by number. Our model makes predictions that will be tested in the coming decade by LSST.


title: Affine and bilinear systems on Lie groups

abstract: In this paper we study affine and bilinear systems on Lie groups. We show that there is an intrinsic connection between the solutions of both systems. Such relation allows us to obtain some preliminary controllability results of affne systems on compact and solvable Lie groups. We also show that the controllability property of bilinear systems is very restricted and may only be achieved if the state space G is an Euclidean space.


title: Bisecting measures with hyperplane arrangements

abstract: We show that when $n$ is a power of two, any $nD$ measures in $\mathbb{R}^n$ can be bisected by an arrangement of $D$ hyperplanes.


title: Flip procedure in geometric approximation of multiple-component shapes -- Application to multiple-inclusion detection

abstract: We are interested in geometric approximation by parameterization of two-dimensional multiple-component shapes, in particular when the number of components is a priori unknown. Starting a standard method based on successive shape deformations with a one-component initial shape in order to approximate a multiple-component target shape usually leads the deformation flow to make the boundary evolve until it surrounds all the components of the target shape. This classical phenomenon tends to create double points on the boundary of the approximated shape. In order to improve the approximation of multiple-component shapes (without any knowledge on the number of components in advance), we use in this paper a piecewise B\'ezier parameterization and we consider two procedures called intersecting control polygons detection and flip procedure. The first one allows to prevent potential collisions between two parts of the boundary of the approximated shape, and the second one permits to change its topology by dividing a one-component shape into a two-component shape. For an experimental purpose, we include these two processes in a basic geometrical shape optimization algorithm and test it on the classical inverse obstacle problem. This new approach allows to obtain a numerical approximation of the unknown inclusion, detecting both the topology (i.e. the number of connected components) and the shape of the obstacle. Several numerical simulations are performed.


title: Rydberg atoms based creation of N particle GHZ state using STIRAP

abstract: Schemes for creation of N particle entangled Greenberger-Horne-Zeilinger (GHZ) states are important for understanding multi-particle non-classical correlations. Here, a theoretical scheme for creation of a multi-particle GHZ state implemented on a target ensemble of N, $\Lambda$ three-level Rydberg atoms and a single Rydberg atom as a control using Stimulated Raman Adiabatic Passage (STIRAP) is presented. We work in the Rydberg blockade regime for the ensemble atoms induced due to excitation of the control atom to a high lying Rydberg level. It is shown that using STIRAP, atoms from one ground state of the ensemble can be adiabatically transferred to the other ground state, depending on the state of the control atom with high fidelity. Measurement of the control atom in a specific basis after this conditional transfer facilitates one-step creation of a N particle GHZ state. A thorough analysis of adiabatic conditions for this scheme and the influence of radiative decay from the excited Rydberg levels is presented. We show that this scheme is immune to the decay rate of the excited level in ensemble atoms and provides a robust way of creating GHZ states.


title: Black Hole Echology: The Observer's Manual

abstract: While recent detections of gravitational waves from the mergers of binary black holes match well with the predictions of General Relativity (GR), they cannot directly confirm the existence of event horizons. Exotic compact objects (ECOs) are motivated by quantum models of black holes, and can have exotic structure (or a "wall") just outside the (would-be) horizon. ECOs produce similar ringdown waveforms to the GR black holes, but they are followed by delayed "echoes". By solving linearized Einstein equations we can model these echoes and provide analytic templates that can be used to compare to observations. For concreteness, we consider GW150914 event, detected by the LIGO/Virgo collaboration, and study the model dependence of its echo properties. We find that echoes are reasonably approximated by complex gaussians, with amplitudes that decay as a power law in time, while their width in time (frequency) grows (shrinks) over subsequent echoes. We also show that trapped modes between a perfectly reflecting wall and angular momentum barrier in Kerr metric can exhibit superradiant instability over long times, as expected.


title: BSMPT - Beyond the Standard Model Phase Transitions -A Tool for the Electroweak Phase Transition in Extended Higgs Sectors

abstract: We provide the C++ tool BSMPT for calculating the strength of the electroweak phase transition in extended Higgs sectors. This relies on the loop-corrected effective potential at finite temperature including daisy resummation of the bosonic masses. The program allows to compute the vacuum expectation value (VEV) $v$ of the potential as a function of the temperature, and in particular the critical VEV $v_c$ at the temperature $T_c$ where the phase transition takes place. In addition, the loop-corrected trilinear Higgs self-couplings are provided. We apply an 'on-shell' renormalization scheme in the sense that the loop-corrected masses and mixing angles are required to be equal to their tree-level input values. This allows for efficient scans in the parameter space of the models. The models implemented so far are the CP-conserving and the CP-violating 2-Higgs-Doublet Models (2HDM) and the Next-to-Minimal 2HDM (N2HDM). The program structure is such that the user can easily implement further models. Our tool can be used for the investigation of electroweak baryogenesis in models with extended Higgs sectors and the related Higgs self-couplings. The combination with parameter scans in the respective models allows to study the impact on collider phenomenology and to make a link between collider phenomenology and cosmology. The program package can be downloaded at: https://github.com/phbasler/BSMPT.


title: Multiple Gaps in the Disk of the Class I Protostar GY 91

abstract: We present the highest spatial resolution ALMA observations to date of the Class I protostar GY 91 in the $\rho$ Ophiuchus L1688 molecular cloud complex. Our 870 $\mu$m and 3 mm dust continuum maps show that the GY 91 disk has a radius of $\sim$80 AU, and an inclination of $\sim$40$^{\circ}$, but most interestingly that the disk has three dark lanes located at 10 AU, 40 AU, and 70 AU. We model these features assuming they are gaps in the disk surface density profile and find that their widths are 7 AU, 30 AU, and 10 AU. These gaps bear a striking resemblance to the gaps seen in the HL Tau disk, suggesting that there may be Saturn-mass planets hiding in the disk. To constrain the relative ages of GY 91 and HL Tau, we also model the disk and envelope of HL Tau and find that they are of similar ages, although GY 91 may be younger. Although snow lines and magnetic dead zones can also produce dark lanes, if planets are indeed carving these gaps then Saturn-mass planets must form within the first $\sim$0.5 Myr of the lifetime of protoplanetary disks.


title: The Randomized Kaczmarz Method with Mismatched Adjoint

abstract: This paper investigates the randomized version of the Kaczmarz method to solve linear systems in the case where the adjoint of the system matrix is not exact---a situation we refer to as "mismatched adjoint". We show that the method may still converge both in the over- and underdetermined consistent case under appropriate conditions, and we calculate the expected asymptotic rate of linear convergence. Moreover, we analyze the inconsistent case and obtain results for the method with mismatched adjoint as for the standard method. Finally, we derive a method to compute optimized probabilities for the choice of the rows and illustrate our findings with numerical example.


title: Good Distance Lattices from High Dimensional Expanders

abstract: We show a new framework for constructing good distance lattices from high dimensional expanders. For error-correcting codes, which have a similar flavor as lattices, there is a known framework that yields good codes from expanders. However, prior to our work, there has been no framework that yields good distance lattices directly from expanders. Interestingly, we need the notion of high dimensional expansion (and not only one dimensional expansion) for obtaining large distance lattices which are dense. Our construction is obtained by proving the existence of bounded degree high dimensional cosystolic expanders over any ring, and in particular over $\mathbb{Z}$. Previous bounded degree cosystolic expanders were known only over $\mathbb{F}_2$. The proof of the cosystolic expansion over any ring is composed of two main steps, each of an independent interest: We show that coboundary expansion over any ring of the links of a bounded degree complex implies that the complex is a cosystolic expander over any ring. We then prove that all the links of Ramanujan complexes (which are called spherical buildings) are coboundary expanders over any ring. We follow the strategy of [LMM16] for proving that the spherical building is a coboundary expander over any ring. Besides of generalizing their proof from $\mathbb{F}_2$ to any ring, we present it in a detailed way, which might serve readers with less background who wish to get into the field.


title: Circular and magnetoinduced photocurrents in Weyl semimetals

abstract: We develop a theory of the direct interband and indirect intraband photogalvanic effects in Weyl semimetals belonging to the gyrotropic classes with improper symmetry operations. At zero magnetic field, an excitation of such a material with circularly polarized light leads to a photocurrent whose direction depends on the light helicity. We show that in the semimetals of the C$_{2v}$ symmetry, an allowance for the tilt term in the effective Hamiltonian is enough to prevent cancellation of the photocurrent contributions from the Weyl cones of opposite chiralities. In the case of the C$_{4v}$ symmetry, in addition to the tilt it is necessary to include terms of the second- or third-order in the electron quasi-momentum. For indirect intraband transitions, the helicity-dependent photocurrent generated within each Weyl node takes on a universal value determined by the fundamental constants, the light frequency and electric field. We have complementarily investigated the magneto-gyrotropic photogalvanic effect, i.e. an appearance of a photocurrent under unpolarized excitation in a magnetic field. In quantized magnetic fields, the photocurrent is caused by optical transitions between the one-dimensional magnetic subbands. A value of the photocurrent is particularly high if one of the photocarriers is excited to the chiral subband with the energy below the cyclotron energy.


title: Ion Traps at the Sun: Implications for Elemental Fractionation

abstract: Why the tenuous solar outer atmosphere, or corona, is much hotter than the underlying layers remains one of the greatest challenge for solar modeling. Detailed diagnostics of the coronal thermal structure come from extreme ultraviolet (EUV) emission. The EUV emission is produced by heavy ions in various ionization states and depends on the amount of these ions and on plasma temperature and density. Any nonuniformity of the elemental distribution in space or variability in time affects thermal diagnostics of the corona. Here we theoretically predict ionized chemical element concentrations in some areas of the solar atmosphere, where the electric current is directed upward. We then detect these areas observationally, by comparing the electric current density with the EUV brightness in an active region. We found a significant excess in EUV brightness in the areas with positive current density rather than negative. Therefore, we report the observational discovery of substantial concentrations of heavy ions in current-carrying magnetic flux tubes, which might have important implications for the elemental fractionation in the solar corona known as the first ionization potential (FIP) effect. We call such areas of heavy ion concentration the "ion traps." These traps hold enhanced ion levels until they are disrupted by a flare whether large or small.


title: Value Alignment, Fair Play, and the Rights of Service Robots

abstract: Ethics and safety research in artificial intelligence is increasingly framed in terms of "alignment" with human values and interests. I argue that Turing's call for "fair play for machines" is an early and often overlooked contribution to the alignment literature. Turing's appeal to fair play suggests a need to correct human behavior to accommodate our machines, a surprising inversion of how value alignment is treated today. Reflections on "fair play" motivate a novel interpretation of Turing's notorious "imitation game" as a condition not of intelligence but instead of value alignment: a machine demonstrates a minimal degree of alignment (with the norms of conversation, for instance) when it can go undetected when interrogated by a human. I carefully distinguish this interpretation from the Moral Turing Test, which is not motivated by a principle of fair play, but instead depends on imitation of human moral behavior. Finally, I consider how the framework of fair play can be used to situate the debate over robot rights within the alignment literature. I argue that extending rights to service robots operating in public spaces is "fair" in precisely the sense that it encourages an alignment of interests between humans and machines.


title: A remark on two notions of order of contact

abstract: We recall two measurements of the order of contact of an ideal in the ring of germs of holomorphic functions at a point and we provide a class of examples in which they differ.


title: A family of singular integral operators which control the Cauchy transform

abstract: We study the behaviour of singular integral operators $T_{k_t}$ of convolution type on $\mathbb{C}$ associated with the parametric kernels $$ k_t(z):=\frac{(\Re z)^{3}}{|z|^{4}}+t\cdot \frac{\Re z}{|z|^{2}}, \quad t\in \mathbb{R},\qquad k_\infty(z):=\frac{\Re z}{|z|^{2}}\equiv \Re \frac{1}{z},\quad z\in \mathbb{C}\setminus\{0\}. $$ It is shown that for any positive locally finite Borel measure with linear growth the corresponding $L^2$-norm of $T_{k_0}$ controls the $L^2$-norm of $T_{k_\infty}$ and thus of the Cauchy transform. As a corollary, we prove that the $L^2(\mathcal{H}^1\lfloor E)$-boundedness of $T_{k_t}$ with a fixed $t\in (-t_0,0)$, where $t_0>0$ is an absolute constant, implies that $E$ is rectifiable. This is so in spite of the fact that the usual curvature method fails to be applicable in this case. Moreover, as a corollary of our techniques, we provide an alternative and simpler proof of the bi-Lipschitz invariance of the $L^2$-boundedness of the Cauchy transform, which is the key ingredient for the bi-Lipschitz invariance of analytic capacity.


title: Satisficing in Time-Sensitive Bandit Learning

abstract: Much of the recent literature on bandit learning focuses on algorithms that aim to converge on an optimal action. One shortcoming is that this orientation does not account for time sensitivity, which can play a crucial role when learning an optimal action requires much more information than near-optimal ones. Indeed, popular approaches such as upper-confidence-bound methods and Thompson sampling can fare poorly in such situations. We consider instead learning a satisficing action, which is near-optimal while requiring less information, and propose satisficing Thompson sampling, an algorithm that serves this purpose. We establish a general bound on expected discounted regret and study the application of satisficing Thompson sampling to linear and infinite-armed bandits, demonstrating arbitrarily large benefits over Thompson sampling. We also discuss the relation between the notion of satisficing and the theory of rate distortion, which offers guidance on the selection of satisficing actions.


title: A Passive Re-Directing Van Atta Type Reflector

abstract: The Van Atta retro-reflector can be envisioned as a self-phasing antenna array; a phase gradient, naturally induced by an external incident wave, is inverted so the device re-radiates back towards the incident direction. This Letter demonstrates how to re-direct the re-radiated beam through passive alteration of this phase gradient. For this purpose, crosspropagating isolation is required between the incident and reradiated signal paths. To this end, polarization duplexing can be used to achieve this isolation with a passive and reciprocal structure. To provide further demonstration, a 4-element redirective array is designed and fabricated to offset its re-radiated beam by -15{\deg} from incidence. Performance is then verified using fullwave electromagnetic simulation and experimental radiation pattern measurements. It is found that the -15{\deg} angular offset is effectively retained over a wide range of incidence angles.


title: Gerrymandering and Compactness: Implementation Flexibility and Abuse

abstract: Political districts may be drawn to favor one group or political party over another, or gerrymandered. A number of measurements have been suggested as ways to detect and prevent such behavior. These measures give concrete axes along which districts and districting plans can be compared. However, measurement values are affected by both noise and the compounding effects of seemingly innocuous implementation decisions. Such issues will arise for any measure. As a case study demonstrating the effect, we show that commonly-used measures of geometric compactness for district boundaries are affected by several factors irrelevant to fairness or compliance with civil rights law. We further show that an adversary could manipulate measurements to affect the assessment of a given plan. This instability complicates using these measurements as legislative or judicial standards to counteract unfair redistricting practices. This paper accompanies the release of packages in C++, Python, and R that correctly, efficiently, and reproducibly calculate a variety of compactness scores.


title: Two warm, low-density sub-Jovian planets orbiting bright stars in K2 campaigns 13 and 14

abstract: We report the discovery of two planets transiting the bright stars HD 89345 (EPIC 248777106, $V=9.376$, $K=7.721$) in K2 Campaign 14 and HD 286123 (EPIC 247098361, $V=9.822$, $K=8.434$) in K2 Campaign 13. Both stars are G-type stars, one of which is at or near the end of its main sequence lifetime, and the other that is just over halfway through its main sequence lifetime. HD 89345 hosts a warm sub-Saturn (0.66 $R_J$, 0.11 $M_J$, $T_\mathrm{eq}=1100$ K) in an 11.81-day orbit. The planet is similar in size to WASP-107b, which falls in the transition region between ice giants and gas giants. HD 286123 hosts a Jupiter-sized, low-mass planet (1.06 $R_J$, 0.39 $M_J$, $T_\mathrm{eq}=1000$ K) in an 11.17-day, mildly eccentric orbit, with $e=0.255\pm0.035$. Given that they orbit relatively evolved main-sequence stars and have orbital periods longer than 10 days, these planets are interesting candidates for studies of gas planet evolution, migration, and (potentially) re-inflation. Both planets have spent their entire lifetimes near the proposed stellar irradiation threshold at which giant planets become inflated, and neither shows any sign of radius inflation. They probe the regime where inflation begins to become noticeable and are valuable in constraining planet inflation models. In addition, the brightness of the host stars, combined with large atmospheric scale heights of the planets, makes these two systems favorable targets for transit spectroscopy to study their atmospheres and perhaps provide insight into the physical mechanisms that lead to inflated hot Jupiters.


title: A Runaway Yellow Supergiant Star in the Small Magellanic Cloud

abstract: We recently discovered a yellow supergiant (YSG) in the Small Magellanic Cloud (SMC) with a heliocentric radial velocity of ~300 km/s which is much larger than expected for a star in its location in the SMC. This is the first runaway YSG ever discovered and only the second evolved runaway star discovered in a different galaxy than the Milky Way. We classify the star as G5-8I, and use de-reddened broad-band colors with model atmospheres to determine an effective temperature of 4700+/-250K, consistent with what is expected from its spectral type. The star's luminosity is then L/Lo ~ 4.2+/-0.1, consistent with it being a ~30Myr 9Mo star according to the Geneva evolution models. The star is currently located in the outer portion of the SMC's body, but if the star's transverse peculiar velocity is similar to its peculiar radial velocity, in 10Myr the star would have moved 1.6 degrees across the disk of the SMC, and could easily have been born in one of the SMC's star-forming regions. Based on its large radial velocity, we suggest it originated in a binary system where the primary exploded as a supernovae thus flinging the runaway star out into space. Such stars may provide an important mechanism for the dispersal of heavier elements in galaxies given the large percentage of massive stars that are runaways. In the future we hope to look into additional evolved runaway stars that were discovered as part of our other past surveys.


title: Responses of the chiral-magnetic-effect-sensitive sine observable to resonance backgrounds in heavy-ion collisions

abstract: A new sine observable, $R_{\Psi_2}(\Delta S)$, has been proposed to measure the chiral magnetic effect (CME) in heavy-ion collisions; $\Delta S = \left \langle \sin \varphi_+ \right \rangle - \left \langle \sin \varphi_- \right \rangle$, where $\varphi_\pm$ are azimuthal angles of positively and negatively charged particles relative to the reaction plane and averages are event-wise, and $R_{\Psi_2}(\Delta S)$ is a normalized event probability distribution. Preliminary STAR data reveal concave $R_{\Psi_2}(\Delta S)$ distributions in 200 GeV Au+Au collisions. Studies with a multiphase transport (AMPT) and anomalous-viscous Fluid Dynamics (AVFD) models show concave $R_{\Psi_2}(\Delta S)$ distributions for CME signals and convex ones for typical resonance backgrounds. A recent hydrodynamic study, however, indicates concave shapes for backgrounds as well. To better understand these results, we report a systematic study of the elliptic flow ($v_{2}$) and transverse momentum ($p_{T}$) dependences of resonance backgrounds with toy-model simulations and central limit theorem (CLT) calculations. It is found that the concavity or convexity of $R_{\Psi_2}(\Delta S)$ depends sensitively on the resonance $v_2$ (which yields different numbers of decay $\pi^+\pi^-$ pairs in the in-plane and out-of-plane directions) and $p_T$ (which affects the opening angle of the decay $\pi^+\pi^-$ pair). Qualitatively, low $p_{T}$ resonances decay into large opening-angle pairs and result in more `back-to-back' pairs out-of-plane, mimicking a CME signal, or a concave $R_{\Psi_2}(\Delta S)$. Supplemental studies of $R_{\Psi_3}(\Delta S)$ in terms of the triangular flow ($v_3$), where only backgrounds exist but any CME would average to zero, are also presented.


title: LHC Phenomenology of Dark Matter with a Color-Octet Partner

abstract: Colored dark sectors where the dark matter particle is accompanied by colored partners have recently attracted theoretical and phenomenological interest. We explore the possibility that the dark sector consists of the dark matter particle and a color-octet partner, where the interaction with the Standard Model is governed by an effective operator involving gluons. The resulting interactions resemble the color analogues of electric and magnetic dipole moments. Although many phenomenological features of this kind of model only depend on the group representation of the partner under SU(3)$_c$, we point out that interesting collider signatures such as $R$-hadrons are indeed controlled by the interaction operator between the dark and visible sector. We perform a study of the current constraints and future reach of LHC searches, where the complementarity between different possible signals is highlighted and exploited.


title: Approximation algorithms for two-machine flow-shop scheduling with a conflict graph

abstract: Path cover is a well-known intractable problem that finds a minimum number of vertex disjoint paths in a given graph to cover all the vertices. We show that a variant, where the objective function is not the number of paths but the number of length-$0$ paths (that is, isolated vertices), turns out to be polynomial-time solvable. We further show that another variant, where the objective function is the total number of length-$0$ and length-$1$ paths, is also polynomial-time solvable. Both variants find applications in approximating the two-machine flow-shop scheduling problem in which job processing has constraints that are formulated as a conflict graph. For the unit jobs, we present a $4/3$-approximation algorithm for the scheduling problem with an arbitrary conflict graph, based on the exact algorithm for the variants of the path cover problem. For the arbitrary jobs while the conflict graph is the union of two disjoint cliques, that is, all the jobs can be partitioned into two groups such that the jobs in a group are pairwise conflicting, we present a simple $3/2$-approximation algorithm.


title: Voltage driven, local, and efficient excitation of nitrogen-vacancy centers in diamond

abstract: Magnetic sensing technology has found widespread application in industries as diverse as transportation, medicine, and resource exploration. Such use cases often require highly sensitive instruments to measure the extremely small magnetic fields involved, relying on difficult to integrate Superconducting Quantum Interference Device (SQUID) and Spin-Exchange Relaxation Free (SERF) magnetometers. A potential alternative, nitrogen vacancy (NV) centers in diamond, has shown great potential as a high sensitivity and high resolution magnetic sensor capable of operating in an unshielded, room-temperature environment. Transitioning NV center based sensors into practical devices, however, is impeded by the need for high power RF excitation to manipulate them. Here we report an advance that combines two different physical phenomena to enable a highly efficient excitation of the NV centers: magnetoelastic drive of ferromagnetic resonance (FMR) and NV-magnon coupling. Our work demonstrates a new pathway to combine acoustics and magnonics that enables highly energy efficient and local excitation of NV centers without the need for any external RF excitation, and thus could lead to completely integrated, on-chip, atomic sensors.


title: Search for third-generation scalar leptoquarks decaying to a top quark and a $\tau$ lepton at $\sqrt{s}=$ 13 TeV

abstract: A search for pair production of heavy scalar leptoquarks (LQs), each decaying into a top quark and a $\tau$ lepton, is presented. The search considers final states with an electron or a muon, one or two $\tau$ leptons that decayed to hadrons, and additional jets. The data were collected in 2016 in proton-proton collisions at $\sqrt{s}=$ 13 TeV with the CMS detector at the LHC, and correspond to an integrated luminosity of 35.9 fb$^{-1}$. No evidence for pair production of LQs is found. Assuming a branching fraction of unity for the decay LQ$ \to $t$\tau$, upper limits on the production cross section are set as a function of LQ mass, excluding masses below 900 GeV at 95% confidence level. These results provide the most stringent limits to date on the production of scalar LQs that decay to a top quark and a $\tau$ lepton.


title: WNGrad: Learn the Learning Rate in Gradient Descent

abstract: Adjusting the learning rate schedule in stochastic gradient methods is an important unresolved problem which requires tuning in practice. If certain parameters of the loss function such as smoothness or strong convexity constants are known, theoretical learning rate schedules can be applied. However, in practice, such parameters are not known, and the loss function of interest is not convex in any case. The recently proposed batch normalization reparametrization is widely adopted in most neural network architectures today because, among other advantages, it is robust to the choice of Lipschitz constant of the gradient in loss function, allowing one to set a large learning rate without worry. Inspired by batch normalization, we propose a general nonlinear update rule for the learning rate in batch and stochastic gradient descent so that the learning rate can be initialized at a high value, and is subsequently decreased according to gradient observations along the way. The proposed method is shown to achieve robustness to the relationship between the learning rate and the Lipschitz constant, and near-optimal convergence rates in both the batch and stochastic settings ($O(1/T)$ for smooth loss in the batch setting, and $O(1/\sqrt{T})$ for convex loss in the stochastic setting). We also show through numerical evidence that such robustness of the proposed method extends to highly nonconvex and possibly non-smooth loss function in deep learning problems.Our analysis establishes some first theoretical understanding into the observed robustness for batch normalization and weight normalization.

categories: stat.ML cs.AI cs.LG cs.NA math.NA math.OC


title: Torsional instability in the single-chain limit of a transition metal trichalcogenide

abstract: We report the synthesis of the quasi-one-dimensional transition metal trichalcogenide NbSe3 in the few chain limit, including the realization of isolated single chains. The chains are encapsulated in protective boron nitride or carbon nanotube sheaths to prevent oxidation and to facilitate characterization. Transmission electron microscopy reveals static and dynamic structural torsional waves not found in bulk NbSe3 crystals. Electronic structure calculations indicate that charge transfer drives the torsional wave instability. Very little covalent bonding is found between the chains and the nanotube sheath, leading to relatively unhindered longitudinal and torsional dynamics for the encapsulated chains.


title: Phase transitions for a model with uncountable spin space on the Cayley tree: the general case

abstract: In this paper we complete the analysis of a statistical mechanics model on Cayley trees of any degree, started in [EsHaRo12,EsRo10,BoEsRo13,JaKuBo14,Bo17]. The potential is of nearest-neighbor type and the local state space is compact but uncountable. Based on the system parameters we prove existence of a critical value $\theta_{\rm c}$ such that for $\theta\le \theta_{\rm c}$ there is a unique translation-invariant splitting Gibbs measure. For $\theta_{\rm c}<\theta$ there is a phase transition with exactly three translation-invariant splitting Gibbs measures. The proof rests on an analysis of fixed points of an associated non-linear Hammerstein integral operator for the boundary laws.


title: Looking for Galactic Diffuse Dark Matter in INO-MagICAL Detector

abstract: The Weakly Interacting Massive Particle (WIMP) is a popular particle physics candidate for the dark matter (DM). It can annihilate and/or decay to neutrino and antineutrino pair. The proposed 50 kt Magnetized Iron CALorimeter (MagICAL) detector at the India-based Neutrino Observatory (INO) can observe these pairs over the conventional atmospheric neutrino and antineutrino fluxes. If we do not see any excess of events in ten years, then INO-Magical can place competitive limits on self-annihilation cross-section ($\langle\sigma v\rangle$) and decay lifetime ($\tau$) of dark matter at 90\% C.L.: $\langle\sigma v\rangle\leq 1.87\,\times\,10^{-24}$ cm$^3$ s$^{-1}$ and $\tau\geq 4.8\,\times\,10^{24}$ s for $m_\chi$ = 10 GeV assuming the NFW as DM density profile.


title: Computing Bottleneck Distance for Multi-parameter Interval Decomposable Persistence Modules

abstract: Computation of the interleaving distance between persistence modules is a central task in topological data analysis. For $1$-parameter persistence modules, thanks to the isometry theorem, this can be done by computing the bottleneck distance with known efficient algorithms. The question is open for most $n$-parameter persistence modules, $n>1$, because of the well recognized complications of the indecomposables. Here, we consider a reasonably complicated class called {\em $n$-parameter interval decomposable} modules whose indecomposables may have a description of non-constant complexity. We present a polynomial time algorithm to compute the bottleneck distance for these modules from indecomposables, which bounds the interleaving distance from above, and give another algorithm to compute a new distance called {\em dimension distance} that bounds it from below. An earlier version of this paper considered only the $2$-parameter interval decomposable modules~\cite{DeyCheng18}.


title: Speech Enhancement Based on Non-stationary Noise-driven Geometric Spectral Subtraction and Phase Spectrum Compensation

abstract: In this paper, a speech enhancement method based on noise compensation performed on short time magnitude as well phase spectra is presented. Unlike the conventional geometric approach (GA) to spectral subtraction (SS), here the noise estimate to be subtracted from the noisy speech spectrum is proposed to be determined by exploiting the low frequency regions of current frame of noisy speech rather than depending only on the initial silence frames. This approach gives the capability of tracking non-stationary noise thus resulting in a non-stationary noise-driven geometric approach of spectral subtraction for speech enhancement. The noise compensated magnitude spectrum from the GA step is then recombined with unchanged phase of noisy speech spectrum and used in phase compensation to obtain an enhanced complex spectrum, which is used to produce an enhanced speech frame. Extensive simulations are carried out using speech files available in the NOIZEUS database shows that the proposed method consistently outperforms some of the recent methods of speech enhancement when employed on the noisy speeches corrupted by street or babble noise at different levels of SNR in terms of objective measures, spectrogram analysis and formal subjective listening tests.


title: Higgs Vacuum Decay from Particle Collisions?

abstract: We examine the effect of large extra dimensions on black hole seeded vacuum decay using the Randall-Sundrum model as a prototype for warped extra dimensions. We model the braneworld black hole by a tidal solution, and solve the Higgs equations of motion for the instanton on the brane. Remarkably, the action of the static instanton can be shown to be the difference in the bulk areas of the seed and remnant black holes, and we estimate these areas assuming the black holes are small compared to the bulk AdS radius. Comparing to the Hawking evaporation rate shows that small black hole seeds preferentially catalyse vacuum decay, thus extending our previous results to higher dimensional braneworld scenarios. The parameter ranges do not allow for Standard Model Higgs decay from collider black holes, but they can be relevant for cosmic ray collisions.


title: The nested structural organization of the worldwide trade multi-layer network

abstract: Nestedness has traditionally been used to detect assembly patterns in meta-communities and networks of interacting species. Attempts have also been made to uncover nested structures in international trade, typically represented as bipartite networks in which connections can be established between countries (exporters or importers) and industries. A bipartite representation of trade, however, inevitably neglects transactions between industries. To fully capture the organization of the global value chain, we draw on the World Input-Output Database and construct a multi-layer network in which the nodes are the countries, the layers are the industries, and links can be established from sellers to buyers within and across industries. We define the buyers' and sellers' participation matrices in which the rows are the countries and the columns are all possible pairs of industries, and then compute nestedness based on buyers' and sellers' involvement in transactions between and within industries. Drawing on appropriate null models that preserve the countries' or layers' degree distributions in the original multi-layer network, we uncover variations of country- and transaction-based nestedness over time, and identify the countries and industries that most contributed to nestedness. We discuss the implications of our findings for the study of the international production network and other real-world systems.


title: Ultrafast lattice deformations studied by means of time-resolved electron and x-ray diffraction

abstract: Ultrafast lattice deformation of tens to hundreds of nanometer thick metallic crystals, after femtosecond laser excitation, was measured directly using 8.04 keV subpicosecond x-ray and 59 keV femtosecond electron pulses. Coherent phonons were generated in both single crystal and polycrystalline films. Lattice compression was observed within the first few picoseconds after laser irradiation in single crystal aluminum, which was attributed to the generation of a blast force and the propagation of elastic waves. The different time scale of lattice heating for tens and hundreds nanometer thick films are clearly distinguished by electron and x-ray pulse diffraction. The electron and lattice heating due to ultrafast deposition of photon energy was numerically simulated using the two-temperature model (TTM) and the results agreed with experimental observations. The ultrafast heating described by TTM was also discussed from an electrical circuit perspective, which may provide new insights on the possible connection between thermal and electrical processes. This study demonstrates that the combination of two complimentary ultrafast time-resolved methods, ultrafast x-ray and electron diffraction will provide a panoramic picture of the transient atomic motions and structure in crystals.


title: Entanglement in a dephasing model and many-body localization

abstract: We study entanglement dynamics in a diagonal dephasing model in which the strength of interaction decays exponentially with distance -- the so-called l-bit model of many-body localization. We calculate the exact expression for entanglement growth with time, finding in addition to a logarithmic growth, a sublogarithmic correction. Provided the l-bit picture correctly describes the many-body localized phase this implies that the entanglement in such systems does not grow (just) as a logarithm of time, as believed so far.