title: Cosmic ray heating in cool core clusters - II. Self-regulation cycle and non-thermal emission

abstract: Self-regulated feedback by active galactic nuclei (AGNs) appears to be critical in balancing radiative cooling of the low-entropy gas at the centres of galaxy clusters and in regulating star formation in central galaxies. In a companion paper, we found steady-state solutions of the hydrodynamic equations that are coupled to the CR energy equation for a large cluster sample. In those solutions, radiative cooling in the central region is balanced by streaming CRs through the generation and dissipation of resonantly generated Alfv{\'e}n waves and by thermal conduction at large radii. Here we demonstrate that the predicted non-thermal emission resulting from hadronic CR interactions in the intra-cluster medium exceeds observational radio (and gamma-ray) data in a subsample of clusters that host radio mini halos (RMHs). In contrast, the predicted non-thermal emission is well below observational data in cooling galaxy clusters without RMHs. These are characterised by exceptionally large AGN radio fluxes, indicating high CR yields and associated CR heating rates. We suggest a self-regulation cycle of AGN feedback in which non-RMH clusters are heated by streaming CRs homogeneously throughout the central cooling region. We predict {\em radio micro halos} surrounding the AGNs of these CR-heated clusters in which the primary emission may predominate the hadronically generated emission. Once the CR population has streamed sufficiently far and lost enough energy, the cooling rate increases, which explains the increased star formation rates in clusters hosting RMHs. Those could be powered hadronically by CRs that have previously heated the cluster core.


title: Automated Visual Fin Identification of Individual Great White Sharks

abstract: This paper discusses the automated visual identification of individual great white sharks from dorsal fin imagery. We propose a computer vision photo ID system and report recognition results over a database of thousands of unconstrained fin images. To the best of our knowledge this line of work establishes the first fully automated contour-based visual ID system in the field of animal biometrics. The approach put forward appreciates shark fins as textureless, flexible and partially occluded objects with an individually characteristic shape. In order to recover animal identities from an image we first introduce an open contour stroke model, which extends multi-scale region segmentation to achieve robust fin detection. Secondly, we show that combinatorial, scale-space selective fingerprinting can successfully encode fin individuality. We then measure the species-specific distribution of visual individuality along the fin contour via an embedding into a global `fin space'. Exploiting this domain, we finally propose a non-linear model for individual animal recognition and combine all approaches into a fine-grained multi-instance framework. We provide a system evaluation, compare results to prior work, and report performance and properties in detail.


title: The Interior of Saturn

abstract: We review our current understanding of the interior structure and thermal evolution of Saturn, with a focus on recent results in the Cassini era. There has been important progress in understanding physical inputs, including equations of state of planetary materials and their mixtures, physical parameters like the gravity field and rotation rate, and constraints on Saturnian free oscillations. At the same time, new methods of calculation, including work on the gravity field of rotating fluid bodies, and the role of interior composition gradients, should help to better constrain the state of Saturn's interior, now and earlier in its history. However, a better appreciation of modeling uncertainties and degeneracies, along with a greater exploration of modeling phase space, still leave great uncertainties in our understanding of Saturn's interior. Further analysis of Cassini data sets, as well as precise gravity field measurements from the Cassini Grand Finale orbits, will further revolutionize our understanding of Saturn's interior over the next few years.


title: Exact distributions of cover times for $N$ independent random walkers in one dimension

abstract: We study the probability density function (PDF) of the cover time $t_c$ of a finite interval of size $L$, by $N$ independent one-dimensional Brownian motions, each with diffusion constant $D$. The cover time $t_c$ is the minimum time needed such that each point of the entire interval is visited by at least one of the $N$ walkers. We derive exact results for the full PDF of $t_c$ for arbitrary $N \geq 1$, for both reflecting and periodic boundary conditions. The PDFs depend explicitly on $N$ and on the boundary conditions. In the limit of large $N$, we show that $t_c$ approaches its average value $\langle t_c \rangle \approx L^2/(16\, D \, \ln N)$, with fluctuations vanishing as $1/(\ln N)^2$. We also compute the centered and scaled limiting distributions for large $N$ for both boundary conditions and show that they are given by nontrivial $N$-independent scaling functions.


title: Probing nonstandard neutrino cosmology with terrestrial neutrino experiments

abstract: Neutrino masses and the number of light neutrino species can be tested in a variety of laboratory experiments and also can be constrained by particle astrophysics and precision cosmology. A conflict between these various results could be an indication of new physics in the neutrino sector. In this paper we explore the possibility for reconciliation of otherwise discrepant results in a simple model containing a light scalar field which produces Mass Varying Neutrinos (MaVaNs). We extend previous work on MaVaNs to consider issues of neutrino clumping, the effects of additional contributions to neutrino mass, and reconciliation of eV mass sterile neutrinos with cosmology.


title: Temporal Map Labeling: A New Unified Framework with Experiments

abstract: The increased availability of interactive maps on the Internet and on personal mobile devices has created new challenges in computational cartography and, in particular, for label placement in maps. Operations like rotation, zoom, and translation dynamically change the map over time and make a consistent adaptation of the map labeling necessary. In this paper, we consider map labeling for the case that a map undergoes a sequence of operations over a specified time span. We unify and generalize several preceding models for dynamic map labeling into one versatile and flexible model. In contrast to previous research, we completely abstract from the particular operations (e.g., zoom, rotation, etc.) and express the labeling problem as a set of time intervals representing the labels' presences, activities, and conflicts. The model's strength is manifested in its simplicity and broad range of applications. In particular, it supports label selection both for map features with fixed position as well as for moving entities (e.g., for tracking vehicles in logistics or air traffic control). Through extensive experiments on OpenStreetMap data, we evaluate our model using algorithms of varying complexity as a case study for navigation systems. Our experiments show that even simple (and thus, fast) algorithms achieve near-optimal solutions in our model with respect to an intuitive objective function.


title: Multiwavelength monitoring and X-ray brightening of Be X-ray binary PSR J2032+4127/MT91 213 on its approach to periastron

abstract: The radio and gamma-ray pulsar PSR J2032+4127 was recently found to be in a decades-long orbit with the Be star MT91 213, with the pulsar moving rapidly towards periastron. This binary shares many similar characteristics with the previously unique binary system PSR B1259-63/LS 2883. Here, we describe radio, X-ray, and optical monitoring of PSR J2032+4127/MT91 213. Our extended orbital phase coverage in radio, supplemented with Fermi gamma-ray data, allows us to update and refine the orbital period to 45-50 yr and time of periastron passage to 2017 November. We analyse archival and recent Chandra and Swift observations and show that PSR J2032+4127/MT91 213 is now brighter in X-rays by a factor of ~70 since 2002 and ~20 since 2010. While the pulsar is still far from periastron, this increase in X-rays is possibly due to collisions between pulsar and Be star winds. Optical observations of the Halpha emission line of the Be star suggest that the size of its circumstellar disc may be varying by ~2 over timescales as short as 1-2 months. Multiwavelength monitoring of PSR J2032+4127/MT91 213 will continue through periastron passage, and the system should present an interesting test case and comparison to PSR B1259-63/LS 2883.


title: Many-Body Theory of Trion Absorption Features in Two-Dimensional Semiconductors

abstract: Recent optical studies of monolayer transition metal dechalcogenides have demonstrated that their excitonic absorption feature splits into two widely separated peaks at finite carrier densities. The additional peak is usually attributed to the presence of trions, bound states of two electrons and a hole or an electron and two holes. Here we argue that in the density range over which the trion peak is well resolved, it cannot be interpreted in terms of weakly coupled three-body systems, and that the appropriate picture is instead one in which excitons are dressed by interactions with a Fermi sea of excess carriers. This coupling splits the exciton spectrum into a lower energy attractive exciton-polaron branch, normally identified as a trion branch, and a higher energy repulsive exciton-polaron branch, normally identified as an exciton branch. We have calculated the frequency and doping dependence of the optical conductivity and found that: (i) the splitting varies linearly with the Fermi energy of the excess quasiparticles; (ii) the trion peak is dominant at high carrier densities; (iii) and the trion peak width is considerably smaller than that of the excitonic peak. Our results are in good agreement with recent experiments.


title: A Time-Dependent Wave-Thermoelastic Solid Interaction

abstract: This paper presents a combined field and boundary integral equation method for solving the time-dependent scattering problem of a thermoelastic body immersed in a compressible, inviscid and homogeneous fluid. The approach here is a generalization of the coupling procedure employed by the authors for the treatment of the time-dependent fluid-structure interaction problem. Using an integral representation of the solution in the infinite exterior domain occupied by the fluid, the problem is reduced to one defined only over the finite region occupied by the solid, with nonlocal boundary conditions. The nonlocal boundary problem is analyzed with Lubich's approach for time-dependent boundary integral equations. Existence and uniqueness results are established in terms of time-domain data with the aid of Laplace-domain techniques. Galerkin semi-discretization approximations are derived and error estimates are obtained. A full discretization based on the Convolution Quadrature method is also outlined. Some numerical experiments in 2D are also included in order to demonstrate the accuracy and efficiency of the procedure.


title: Max-affine estimators for convex stochastic programming

abstract: In this paper, we consider two sequential decision making problems with a convexity structure, namely an energy storage optimization task and a multi-product assembly example. We formulate these problems in the stochastic programming framework and discuss an approximate dynamic programming technique for their solutions. As the cost-to-go functions are convex in these cases, we use max-affine estimates for their approximations. To train such a max-affine estimate, we provide a new convex regression algorithm, and evaluate it empirically for these planning scenarios.


title: Radiative Transfer Modeling of Radio-band Linear Polarization Observations as a Probe of the Physical Conditions in the Jets of Gamma-ray Flaring Blazars

abstract: Since the mid-1980s the shock-in-jet model has been the preferred paradigm to explain radio-band flaring in blazar jets. We describe our radiative transfer model incorporating relativistically-propagating shocks, and illustrate how the 4.8, 8, and 14.5 GHz linear polarization and total flux density data from the University of Michigan monitoring program, in combination with the model, constrain jet flow conditions and shock attributes. Results from strong Fermi-era flares in 4 blazars with widely-ranging properties are presented. Additionally, to investigate jet evolution on decadal time scales we analyze 3 outbursts in OT 081 spanning nearly 3 decades and find intrinsic changes attributable to flow changes at a common spatial location, or, alternatively, to a change in the jet segment viewed. The model's success in reproducing these data supports a scenario in which relativistic shocks compress a plasma with an embedded passive, initially-turbulent magnetic field, with additional ordered magnetic field components, one of which may be helical.


title: Tailoring Thermal Conductivity of Single-stranded Carbon-chain Polymers through Atomic Mass Modification

abstract: Tailoring the thermal conductivity of polymers is central to enlarge their applications in the thermal management of flexible integrated circuits. Progress has been made over the past decade by fabricating materials with various nanostructures, but a clear relationship between various functional groups and thermal properties of polymers remains to be established. Here, we numerically study the thermal conductivity of single-stranded carbon-chain polymers with multiple substituents of hydrogen atoms through atomic mass modification. We find that their thermal conductivity can be tuned by atomic mass modifications as revealed through molecular dynamics simulations. The simulation results suggest that heavy homogeneous substituents do not assist heat transport and trace amounts of heavy substituents can in fact hinder heat transport substantially. Our analysis indicates that carbon chain has the biggest contribution (over 80%) to the thermal conduction in single-stranded carbon-chain polymers. We further demonstrate that atomic mass modifications influence the phonon bands of bonding carbon atoms, and the discrepancies of phonon bands between carbon atoms are responsible for the remarkable drops in thermal conductivity and large thermal resistances in carbon chains. Our study provides fundamental insight into how to tailor the thermal conductivity of polymers through variable substituents.


title: Blockage induced condensation controlled by a local reaction

abstract: We consider the set-up of stationary Zero Range models and discuss the onset of condensation induced by a local blockage on the lattice. We show that the introduction of a local feedback on the hopping rates allows to control the particle fraction in the condensed phase. This phenomenon results in a current vs. blockage parameter curve characterized by two non--analyticity points.


title: Unsupervised learning of transcriptional regulatory networks via latent tree graphical models

abstract: Gene expression is a readily-observed quantification of transcriptional activity and cellular state that enables the recovery of the relationships between regulators and their target genes. Reconstructing transcriptional regulatory networks from gene expression data is a problem that has attracted much attention, but previous work often makes the simplifying (but unrealistic) assumption that regulator activity is represented by mRNA levels. We use a latent tree graphical model to analyze gene expression without relying on transcription factor expression as a proxy for regulator activity. The latent tree model is a type of Markov random field that includes both observed gene variables and latent (hidden) variables, which factorize on a Markov tree. Through efficient unsupervised learning approaches, we determine which groups of genes are co-regulated by hidden regulators and the activity levels of those regulators. Post-processing annotates many of these discovered latent variables as specific transcription factors or groups of transcription factors. Other latent variables do not necessarily represent physical regulators but instead reveal hidden structure in the gene expression such as shared biological function. We apply the latent tree graphical model to a yeast stress response dataset. In addition to novel predictions, such as condition-specific binding of the transcription factor Msn4, our model recovers many known aspects of the yeast regulatory network. These include groups of co-regulated genes, condition-specific regulator activity, and combinatorial regulation among transcription factors. The latent tree graphical model is a general approach for analyzing gene expression data that requires no prior knowledge of which possible regulators exist, regulator activity, or where transcription factors physically bind.


title: Performing Hong-Ou-Mandel-type Numerical Experiments with Repulsive Condensates: The case of Dark and Dark-bright Solitons

abstract: The Hong-Ou-Mandel experiment leads indistinguishable photons simultaneously reaching a 50:50 beam splitter to emerge on the same port through two-photon interference. Motivated by this phenomenon, we consider numerical experiments of the same flavor for classical, {\it wave} objects in the setting of repulsive condensates. We examine dark solitons interacting with a repulsive barrier, a case in which we find no significant asymmetries in the emerging waves after the collision, presumably due to their topological nature. We also consider case examples of two-component systems, where the dark solitons trap a bright structure in the second-component (dark-bright solitary waves). For these, pronounced asymmetries upon collision {\it are} possible for the non-topological bright component. We also show an example of a similar phenomenology for ring dark-bright structures in two dimensions.


title: Hydrocarbon emission rings in protoplanetary disks induced by dust evolution

abstract: We report observations of resolved C2H emission rings within the gas-rich protoplanetary disks of TW Hya and DM Tau using the Atacama Large Millimeter Array (ALMA). In each case the emission ring is found to arise at the edge of the observable disk of mm-sized grains (pebbles) traced by (sub)mm-wave continuum emission. In addition, we detect a C3H2 emission ring with an identical spatial distribution to C2H in the TW Hya disk. This suggests that these are hydrocarbon rings (i.e. not limited to C2H). Using a detailed thermo-chemical model we show that reproducing the emission from C2H requires a strong UV field and C/O > 1 in the upper disk atmosphere and outer disk, beyond the edge of the pebble disk. This naturally arises in a disk where the ice-coated dust mass is spatially stratified due to the combined effects of coagulation, gravitational settling and drift. This stratification causes the disk surface and outer disk to have a greater permeability to UV photons. Furthermore the concentration of ices that transport key volatile carriers of oxygen and carbon in the midplane, along with photochemical erosion of CO, leads to an elemental C/O ratio that exceeds unity in the UV-dominated disk. Thus the motions of the grains, and not the gas, lead to a rich hydrocarbon chemistry in disk surface layers and in the outer disk midplane.


title: Assessing the Impact of Retreat Mechanisms in a Simple Antarctic Ice Sheet Model Using Bayesian Calibration

abstract: The response of the Antarctic ice sheet (AIS) to changing climate forcings is an important driver of sea-level changes. Anthropogenic climate change may drive a sizeable AIS tipping point response with subsequent increases in coastal flooding risks. Many studies analyzing flood risks use simple models to project the future responses of AIS and its sea-level contributions. These analyses have provided important new insights, but they are often silent on the effects of potentially important processes such as Marine Ice Sheet Instability (MISI) or Marine Ice Cliff Instability (MICI). These approximations can be well justified and result in more parsimonious and transparent model structures. This raises the question of how this approximation impacts hindcasts and projections. Here, we calibrate a previously published and relatively simple AIS model, which neglects the effects of MICI and regional characteristics, using a combination of observational constraints and a Bayesian inversion method. Specifically, we approximate the effects of missing MICI by comparing our results to those from expert assessments with more realistic models and quantify the bias during the last interglacial when MICI may have been triggered. Our results suggest that the model can approximate the process of MISI and reproduce the projected median melt from some previous expert assessments in the year 2100. Yet, our mean hindcast is roughly 3/4 of the observed data during the last interglacial period and our mean projection is roughly 1/6 and 1/10 of the mean from a model accounting for MICI in the year 2100. These results suggest that missing MICI and/or regional characteristics can lead to a low-bias during warming period AIS melting and hence a potential low-bias in projected sea levels and flood risks.


title: Profiting from correlations: Adjusted estimators for categorical data

abstract: To take sample biases and skewness in the observations into account, practitioners frequently weight their observations according to some marginal distribution. The present paper demonstrates that such weighting can indeed improve the estimation. Studying contingency tables, estimators for marginal distributions are proposed under the assumption that another marginal is known. It is shown that the weighted estimators have a strictly smaller asymptotic variance whenever the two marginals are correlated. The finite sample performance is illustrated in a simulation study. As an application to traffic accident data the method allows for correcting a well-known bias in the observed injury severity distribution.


title: Pattern Recognition In Non-Kolmogorovian Structures

abstract: We present a generalization of the problem of pattern recognition to arbitrary probabilistic models. This version deals with the problem of recognizing an individual pattern among a family of different species or classes of objects which obey probabilistic laws which do not comply with Kolmogorov's axioms. We show that such a scenario accommodates many important examples, and in particular, we provide a rigorous definition of the classical and the quantum pattern recognition problems, respectively. Our framework allows for the introduction of non-trivial correlations (as entanglement or discord) between the different species involved, opening the door to a new way of harnessing these physical resources for solving pattern recognition problems. Finally, we present some examples and discuss the computational complexity of the quantum pattern recognition problem, showing that the most important quantum computation algorithms can be described as non-Kolmogorovian pattern recognition problems.


title: Markov Random Field Model-Based Salt and Pepper Noise Removal

abstract: Problem of impulse noise reduction is a very well studied problem in image processing community and many different approaches have been proposed to tackle this problem. In the current work, the problem of fixed value impulse noise (salt and pepper) removal from images is investigated by use of a Markov Random Field (MRF) models with smoothness priors. After the formulation of the problem as an inpainting problem, graph cuts with $\alpha$-expansion moves are considered for minimization of the energy functional. As for comparisons, several other minimization techniques that are widely used for MRF models' optimization are considered and the results are compared using Peak-Signal-to-Noise-Ratio (PSNR) and Structural Similarity Index (SSIM) as metrics. The investigations show the superiority of graph cuts with $\alpha$-expansion moves over the other techniques both in terms of PSNR and also computational times.


title: Finding Linear-Recurrent Solutions to Hofstadter-Like Recurrences Using Symbolic Computation

abstract: The Hofstadter Q-sequence, with its simple definition, has defied all attempts at analyzing its behavior. Defined by a simple nested recurrence and an initial condition, the sequence looks approximately linear, though with a lot of noise. But, nobody even knows whether the sequence is infinite. In the years since Hofstadter published his sequence, various people have found variants with predictable behavior. Oftentimes, the resulting sequence looks complicated but provably grows linearly. Other times, the sequences are eventually linear recurrent. Proofs describing the behaviors of both types of sequence are inductive. In the first case, the inductive hypotheses are fairly ad-hoc, but the proofs in the second case are highly automatable. This suggests that a search for more sequences like these may be fruitful. In this paper, we develop a step-by-step symbolic algorithm to search for these sequences. Using this algorithm, we determine that such sequences come in infinite families that are themselves plentiful. In fact,there are hundreds of easy to describe families based on the Hofstadter Q-recurrence alone.


title: The Asymptotics of Representations for Cyclic Opers

abstract: Given a Riemann surface $X = (\Sigma, J)$ we find an expression for the dominant term for the asymptotics of the holonomy of opers over that Riemann surface corresponding to rays in the Hitchin base of the form $(0,0,\cdots,t\omega_n)$. Moreover, we find an associated equivariant map from the universal cover $(\tilde{\Sigma},\tilde{J})$ to the symmetric space SL$_n(\mathbb{C}) / \mbox{SU}(n)$ and show that limits of these maps tend to a sub-building in the asymptotic cone. That sub-building is explicitly constructed from the local data of $\omega_n$.


title: Controlling the Multiport Nature of Bragg Diffraction in Atom Interferometry

abstract: Bragg diffraction has been used in atom interferometers because it allows signal enhancement through multiphoton momentum transfer and suppression of systematics by not changing the internal state of atoms. Its multi-port nature, however, can lead to parasitic interferometers, allows for intensity-dependent phase shifts in the primary interferometers, and distorts the ellipses used for phase extraction. We study and suppress these unwanted effects. Specifically, phase extraction by ellipse fitting and the resulting systematic phase shifts are calculated by Monte Carlo simulations. Phase shifts arising from the thermal motion of the atoms are controlled by spatial selection of atoms and an appropriate choice of Bragg intensity. In these simulations, we found that Gaussian Bragg pulse shapes yield the smallest systematic shifts. Parasitic interferometers are suppressed by a "magic" Bragg pulse duration. The sensitivity of the apparatus was improved by the addition of AC Stark shift compensation, which permits direct experimental study of sub-part-per-billion (ppb) systematics. This upgrade allows for a 310\,$\hbar k$ momentum transfer, giving an unprecedented 6.6\,Mrad measured in a Ramsey-Bord\'e interferometer.


title: Entwicklung eines Reputationssystems f\"ur cyber-physikalische Systeme am Beispiel des inHMotion Forschungsprojektes

abstract: Through the continued progress and advances of networking technology and also through efforts like Industry 4.0, SmartGrid, SmartCities, systems, that have been self-sustaining for the most part are now forming networked Cyber-physical Systems (CPS). Networking does not only have undoubtely benefits but also increase the often overlooked attack surface of these systems. A basic problem in information technology systems is building and sustaining trust. This work looks at the usage of reputation systems in the field of CPS and whether they could be used to assert the trustworthyness of communication partners. This is done in context of a research project inHMotion of the University of applied Sciences Munich, which is in the field of Intelligent Transportation Systems (ITS). A risk based analysis is performed for this project. Based on solicited threats mitigations are proposed. A possible solution might be the usage of a reputation system, which is designed and a prototype implementation for a simulation ist done. The concept of a meta reputation system, which combines different reputation systems and verification systems can achieve promising results in a simple simulation. In a concluding discussion criticism and possible improvements are given.


title: An approximate version of a conjecture of Aharoni and Berger

abstract: Aharoni and Berger conjectured that in every proper edge-colouring of a bipartite multigraph by $n$ colours with at least $n+1$ edges of each colour there is a rainbow matching using every colour. This conjecture generalizes a longstanding problem of Brualdi and Stein about transversals in Latin squares. Here an approximate version of the Aharoni-Berger Conjecture is proved---it is shown that if there are at least $n+o(n)$ edges of each colour in a proper $n$-edge-colouring of a bipartite multigraph then there is a rainbow matching using every colour.


title: Stabilizing Embedology: Geometry-Preserving Delay-Coordinate Maps

abstract: Delay-coordinate mapping is an effective and widely used technique for reconstructing and analyzing the dynamics of a nonlinear system based on time-series outputs. The efficacy of delay-coordinate mapping has long been supported by Takens' embedding theorem, which guarantees that delay-coordinate maps use the time-series output to provide a reconstruction of the hidden state space that is a one-to-one embedding of the system's attractor. While this topological guarantee ensures that distinct points in the reconstruction correspond to distinct points in the original state space, it does not characterize the quality of this embedding or illuminate how the specific parameters affect the reconstruction. In this paper, we extend Takens' result by establishing conditions under which delay-coordinate mapping is guaranteed to provide a stable embedding of a system's attractor. Beyond only preserving the attractor topology, a stable embedding preserves the attractor geometry by ensuring that distances between points in the state space are approximately preserved. In particular, we find that delay-coordinate mapping stably embeds an attractor of a dynamical system if the stable rank of the system is large enough to be proportional to the dimension of the attractor. The stable rank reflects the relation between the sampling interval and the number of delays in delay-coordinate mapping. Our theoretical findings give guidance to choosing system parameters, echoing the trade-off between irrelevancy and redundancy that has been heuristically investigated in the literature. Our initial result is stated for attractors that are smooth submanifolds of Euclidean space, with extensions provided for the case of strange attractors.


title: The Minimum Halo Mass for Star Formation at z = 6 - 8

abstract: Recent analysis of strongly-lensed sources in the Hubble Frontier Fields indicates that the rest-frame UV luminosity function of galaxies at $z=$6--8 rises as a power law down to $M_\mathrm{UV}=-15$, and possibly as faint as -12.5. We use predictions from a cosmological radiation hydrodynamic simulation to map these luminosities onto physical space, constraining the minimum dark matter halo mass and stellar mass that the Frontier Fields probe. While previously-published theoretical studies have suggested or assumed that early star formation was suppressed in halos less massive than $10^9$--$10^{11} M_\odot$, we find that recent observations demand vigorous star formation in halos at least as massive as (3.1, 5.6, 10.5)$\times10^9 M_\odot$ at $z=(6,7,8)$. Likewise, we find that Frontier Fields observations probe down to stellar masses of (8.1, 18, 32)$\times10^6 M_\odot$; that is, they are observing the likely progenitors of analogues to Local Group dwarfs such as Pegasus and M32. Our simulations yield somewhat different constraints than two complementary models that have been invoked in similar analyses, emphasizing the need for further observational constraints on the galaxy-halo connection.


title: A review of matrix scaling and Sinkhorn's normal form for matrices and positive maps

abstract: Given a nonnegative matrix $A$, can you find diagonal matrices $D_1,~D_2$ such that $D_1AD_2$ is doubly stochastic? The answer to this question is known as Sinkhorn's theorem. It has been proved with a wide variety of methods, each presenting a variety of possible generalisations. Recently, generalisations such as to positive maps between matrix algebras have become more and more interesting for applications. This text gives a review of over 70 years of matrix scaling. The focus lies on the mathematical landscape surrounding the problem and its solution as well as the generalisation to positive maps and contains hardly any nontrivial unpublished results.


title: Extension of the osp(m|n)~ so(m-n) Correspondence to the Infinite-Dimensional Chiral Spinors and Self Dual Tensors

abstract: The spinor representations of the orthosymplectic Lie superalgebras osp(m|n) are considered and constructed. These are infinite-dimensional irreducible representations, of which the superdimension coincides with the dimension of the spinor representation of so(m-n). Next, we consider the self dual tensor representations of osp(m|n) and their generalizations: these are also infinite-dimensional and correspond to the highest irreducible component of the $p^{th}$ power of the spinor representation. We determine the character of these representations, and deduce a superdimension formula. From this, it follows that also for these representations the osp(m|n)~ so(m-n) correspondence holds.


title: Cosmology with Independently Varying Neutrino Temperature and Number

abstract: We consider Big Bang nucleosynthesis and the cosmic microwave background when both the neutrino temperature and neutrino number are allowed to vary from their standard values. The neutrino temperature is assumed to differ from its standard model value by a fixed factor from Big Bang nucleosynthesis up to the present. In this scenario, the effective number of relativistic degrees of freedom, $N_{\rm eff}^{\rm CMB}$, derived from observations of the cosmic microwave background is not equal to the true number of neutrinos, $N_\nu$. We determine the element abundances predicted by Big Bang nucleosynthesis as a function of the neutrino number and temperature, converting the latter to the equivalent value of $N_{\rm eff}^{\rm CMB}$. We find that a value of $N_{\rm eff}^{\rm CMB} \approx 3$ can be made consistent with $N_\nu = 4$ with a decrease in the neutrino temperature of $\sim 5\%$, while $N_\nu = 5$ is excluded for any value of $N_{\rm eff}^{\rm CMB}$. No observationally-allowed values for $N_{\rm eff}^{\rm CMB}$ and $N_\nu$ can solve the lithium problem.


title: Magnetic field dependence of the internal quality factor and noise performance of lumped-element kinetic inductance detectors

abstract: We present a technique for increasing the internal quality factor of kinetic inductance detectors (KIDs) by nulling ambient magnetic fields with a properly applied magnetic field. The KIDs used in this study are made from thin-film aluminum, they are mounted inside a light-tight package made from bulk aluminum, and they are operated near $150 \, \mathrm{mK}$. Since the thin-film aluminum has a slightly elevated critical temperature ($T_\mathrm{c} = 1.4 \, \mathrm{K}$), it therefore transitions before the package ($T_\mathrm{c} = 1.2 \, \mathrm{K}$), which also serves as a magnetic shield. On cooldown, ambient magnetic fields as small as approximately $30 \, \mathrm{\mu T}$ can produce vortices in the thin-film aluminum as it transitions because the bulk aluminum package has not yet transitioned and therefore is not yet shielding. These vortices become trapped inside the aluminum package below $1.2 \, \mathrm{K}$ and ultimately produce low internal quality factors in the thin-film superconducting resonators. We show that by controlling the strength of the magnetic field present when the thin film transitions, we can control the internal quality factor of the resonators. We also compare the noise performance with and without vortices present, and find no evidence for excess noise beyond the increase in amplifier noise, which is expected with increasing loss.


title: Degraded Broadcast Channel with Secrecy Outside a Bounded Range

abstract: The $K$-receiver degraded broadcast channel with secrecy outside a bounded range is studied, in which a transmitter sends $K$ messages to $K$ receivers, and the channel quality gradually degrades from receiver $K$ to receiver 1. Each receiver $k$ is required to decode message $W_1,\ldots,W_k$, for $1\leq k\leq K$, and to be kept ignorant of $W_{k+2},\ldots,W_K$, for $k=1,\ldots, K-2$. Thus, each message $W_k$ is kept secure from receivers with at least two-level worse channel quality, i.e., receivers 1, $\ldots$, $k-2$. The secrecy capacity region is fully characterized. The achievable scheme designates one superposition layer to each message with binning employed for each layer. Joint embedded coding and binning are employed to protect all upper-layer messages from lower-layer receivers. Furthermore, the scheme allows adjacent layers to share rates so that part of the rate of each message can be shared with its immediate upper-layer message to enlarge the rate region. More importantly, an induction approach is developed to perform Fourier-Motzkin elimination of $2K$ variables from the order of $K^2$ bounds to obtain a close-form achievable rate region. An outer bound is developed that matches the achievable rate region, whose proof involves recursive construction of the rate bounds and exploits the intuition gained from the achievable scheme.


title: Recognizing Detailed Human Context In-the-Wild from Smartphones and Smartwatches

abstract: The ability to automatically recognize a person's behavioral context can contribute to health monitoring, aging care and many other domains. Validating context recognition in-the-wild is crucial to promote practical applications that work in real-life settings. We collected over 300k minutes of sensor data with context labels from 60 subjects. Unlike previous studies, our subjects used their own personal phone, in any way that was convenient to them, and engaged in their routine in their natural environments. Unscripted behavior and unconstrained phone usage resulted in situations that are harder to recognize. We demonstrate how fusion of multi-modal sensors is important for resolving such cases. We present a baseline system, and encourage researchers to use our public dataset to compare methods and improve context recognition in-the-wild.


title: Outlaw distributions and locally decodable codes

abstract: Locally decodable codes (LDCs) are error correcting codes that allow for decoding of a single message bit using a small number of queries to a corrupted encoding. Despite decades of study, the optimal trade-off between query complexity and codeword length is far from understood. In this work, we give a new characterization of LDCs using distributions over Boolean functions whose expectation is hard to approximate (in~$L_\infty$~norm) with a small number of samples. We coin the term `outlaw distributions' for such distributions since they `defy' the Law of Large Numbers. We show that the existence of outlaw distributions over sufficiently `smooth' functions implies the existence of constant query LDCs and vice versa. We give several candidates for outlaw distributions over smooth functions coming from finite field incidence geometry, additive combinatorics and from hypergraph (non)expanders. We also prove a useful lemma showing that (smooth) LDCs which are only required to work on average over a random message and a random message index can be turned into true LDCs at the cost of only constant factors in the parameters.


title: Review on carrier multiplication in graphene

abstract: The remarkable gapless and linear band structure of graphene opens up new carrier relaxation channels bridging the valence and the conduction band. These Auger scattering processes change the number of charge carriers and can give rise to a significant multiplication of optically excited carriers in graphene. This is an ultrafast many-particle phenomenon that is of great interest both for fundamental many-particle physics as well as technological applications. Here, we review the research on carrier multiplication in graphene and Landau-quantized graphene including theoretical modelling and experimental demonstration.


title: Sparse Blind Deconvolution and Demixing Through $\ell_{1,2}$-Minimization

abstract: This paper concerns solving the sparse deconvolution and demixing problem using $\ell_{1,2}$-minimization. We show that under a certain structured random model, robust and stable recovery is possible. The results extend results of Ling and Strohmer [Self Calibration and Biconvex Compressive Sensing, Inverse Problems, 2015], and in particular theoretically explain certain experimental findings from that paper. Our results do not only apply to the deconvolution and demixing problem, but to recovery of column-sparse matrices in general.


title: Breakdown of the equivalence between active gravitational mass and energy for a quantum body

abstract: We determine active gravitational mass operator of the simplest composite quantum body - a hydrogen atom - within the semiclassical approach to the Einstein equation for a gravitational field. We show that the expectation value of the mass is equivalent to energy for stationary quantum states. On the other hand, it occurs that, for quantum superpositions of stationary states with constant expectation values of energy, the expectation values of the gravitational mass exhibit time-dependent oscillations. This breaks the equivalence between active gravitational mass and energy and can be observed as a macroscopic effect for a macroscopic ensemble of coherent quantum states of the atoms. The corresponding experiment could be the first direct observation of quantum effects in General Relativity.


title: Game among Interdependent Networks: The Impact of Rationality on System Robustness

abstract: Many real-world systems are composed of interdependent networks that rely on one another. Such networks are typically designed and operated by different entities, who aim at maximizing their own payoffs. There exists a game among these entities when designing their own networks. In this paper, we study the game investigating how the rational behaviors of entities impact the system robustness. We first introduce a mathematical model to quantify the interacting payoffs among varying entities. Then we study the Nash equilibrium of the game and compare it with the optimal social welfare. We reveal that the cooperation among different entities can be reached to maximize the social welfare in continuous game only when the average degree of each network is constant. Therefore, the huge gap between Nash equilibrium and optimal social welfare generally exists. The rationality of entities makes the system inherently deficient and even renders it extremely vulnerable in some cases. We analyze our model for two concrete systems with continuous strategy space and discrete strategy space, respectively. Furthermore, we uncover some factors (such as weakening coupled strength of interdependent networks, designing suitable topology dependency of the system) that help reduce the gap and the system vulnerability.


title: Grassmann phase-space methods for fermions: uncovering classical probability structure

abstract: The phase-space description of bosonic quantum systems has numerous applications in such fields as quantum optics, trapped ultracold atoms, and transport phenomena. Extension of this description to the case of fermionic systems leads to formal Grassmann phase-space quasiprobability distributions and master equations. The latter are usually considered as not possessing probabillistic interpretation and as not directly computationally accessible. Here, we describe how to construct $c$-number interpretations of Grassmann phase space representations and their master equations. As a specific example, the Grassmann $B$ representation is considered. We disscuss how to introduce $c$-number probability distributions on Grassmann algebra and how to integrate them. A measure of size and proximity is defined for Grassmann numbers, and the Grassmann derivatives are introduced which are based on infinitesimal variations of function arguments. An example of $c$-number interpretation of formal Grassmann equations is presented.

categories: quant-ph cond-mat.quant-gas math-ph math.MP


title: Prethermalization from a low-density Holstein-Primakoff expansion

abstract: We consider the non-equilibrium dynamics arising after a quench of the transverse magnetic field of a quantum Ising chain, together with the sudden switch-on of a long-range interaction term. The dynamics after the quantum quench is mapped onto a fully-connected model of hard-core bosons, after a suitable combination of a Holstein-Primakoff transformation and of a low-density expansion in the quasi-particles injected by the quench. This mapping holds for a broad class of initial states and for quenches which do not cross the critical point of the transverse field Ising model. We then study the algebraic relaxation in time of a number of observables towards a metastable, pre-thermal state, which becomes the asymptotic steady state in the thermodynamic limit.


title: Impact of the tidal p-g instability on the gravitational wave signal from coalescing binary neutron stars

abstract: Recent studies suggest that coalescing neutron stars are subject to a fluid instability involving the nonlinear coupling of the tide to $p$-modes and $g$-modes. Its influence on the inspiral dynamics and thus the gravitational wave signal is, however, uncertain because we do not know precisely how the instability saturates. Here we construct a simple, physically motivated model of the saturation that allows us to explore the instability's impact as a function of the model parameters. We find that for plausible assumptions about the saturation, current gravitational wave detectors might miss $> 70\%$ of events if only point particle waveforms are used. Parameters such as the chirp mass, component masses, and luminosity distance might also be significantly biased. On the other hand, we find that relatively simple modifications to the point particle waveform can alleviate these problems and enhance the science that emerges from the detection of binary neutron stars.


title: Stationary averaging for multi-scale continuous time Markov chains using parallel replica dynamics

abstract: We propose two algorithms for simulating continuous time Markov chains in the presence of metastability. We show that the algorithms correctly estimate, under the ergodicity assumption, stationary averages of the process. Both algorithms, based on the idea of the parallel replica method, use parallel computing in order to explore metastable sets more efficiently. The algorithms require no assumptions on the Markov chains beyond ergodicity and the presence of identifiable metastability. In particular, there is no assumption on reversibility. For simpler illustration of the algorithms, we assume that a synchronous architecture is used throughout of the paper. We present error analyses, as well as numerical simulations on multi-scale stochastic reaction network models in order to demonstrate consistency of the method and its efficiency.


title: Sparse Bounds for Oscillatory and Random Singular Integrals

abstract: Let $ T_{P } f (x) = \int e ^{i P (y)} K (y) f (x-y) \; dy $, where $ K (y)$ is a smooth Calder\'on-Zygmund kernel on $ \mathbb R ^{n}$, and $ P$ be a polynomial. We show that there is a sparse bound for the bilinear form $ \langle T_P f, g \rangle$. This in turn easily implies $ A_p $ inequalities. The method of proof is applied in a random discrete setting, yielding the first weighted inequalities for operators defined on sparse sets of integers.


title: Phenomenological Implications of the $p_T$ Spectra of $\phi$ and $\Omega$ produced at LHC and RHIC

abstract: The data on the $p_T$ spectra of $\phi$ and $\Omega$ at LHC can be presented in a format that shows exponential behavior up to $p_T\approx 6$ GeV/c with the same slope for both particles and for nearly all centralities. They are empirical properties that are shared at lower energies with the inverse slope showing a power-law dependence on $\sqrt{s_{NN}}$. The shared properties of the spectra are shown to emerge naturally from the recombination model. No flow is needed. We find experimental hints for the possibility that $\phi$ and $\Omega$ are mostly produced in the ridge that are generated by minijets. Appropriate experimental test is suggested.


title: The JCMT Gould Belt Survey: Dense Core Clusters in Orion A

abstract: The Orion A molecular cloud is one of the most well-studied nearby star-forming regions, and includes regions of both highly clustered and more dispersed star formation across its full extent. Here, we analyze dense, star-forming cores identified in the 850 {\mu}m and 450 {\mu}m SCUBA-2 maps from the JCMT Gould Belt Legacy Survey. We identify dense cores in a uniform manner across the Orion A cloud and analyze their clustering properties. Using two independent lines of analysis, we find evidence that clusters of dense cores tend to be mass segregated, suggesting that stellar clusters may have some amount of primordial mass segregation already imprinted in them at an early stage. We also demonstrate that the dense core clusters have a tendency to be elongated, perhaps indicating a formation mechanism linked to the filamentary structure within molecular clouds.


title: Anomalous Hall hysteresis in Tm3Fe5O12/Pt with strain-induced perpendicular magnetic anisotropy

abstract: We demonstrate robust interface strain-induced perpendicular magnetic anisotropy in atomically flat ferrimagnetic insulator Tm3Fe5O12 films grown with pulsed laser deposition on substituted-Gd3Ga5O12 substrate which maximizes the tensile strain at the interface. In bilayers consisting of Pt and TIG, we observe large squared Hall hysteresis loops over a wide range of thicknesses of Pt at room temperature. When a thin Cu layer is inserted between Pt and TIG, the Hall hysteresis magnitude decays but stays finite as the thickness of Cu increases up to 5 nm. However, if the Cu layer is placed atop Pt instead, the Hall hysteresis magnitude is consistently larger than when the Cu layer with the same thickness is inserted in between for all Cu thicknesses. These results suggest that both the proximity-induced ferromagnetism and spin current contribute to the anomalous Hall effect.


title: Privacy-preserving Average Consensus: Privacy Analysis and Optimal Algorithm Design

abstract: Privacy-preserving average consensus aims to guarantee the privacy of initial states and asymptotic consensus on the exact average of the initial value. In existing work, it is achieved by adding and subtracting variance decaying and zero-sum random noises to the consensus process. However, there is lack of theoretical analysis to quantify the degree of the privacy protection. In this paper, we introduce the maximum disclosure probability that the other nodes can infer one node's initial state within a given small interval to quantify the privacy. We develop a novel privacy definition, named $(\epsilon, \delta)$-data-privacy, to depict the relationship between maximum disclosure probability and estimation accuracy. Then, we prove that the general privacy-preserving average consensus (GPAC) provides $(\epsilon, \delta)$-data-privacy, and provide the closed-form expression of the relationship between $\epsilon$ and $\delta$. Meanwhile, it is shown that the added noise with uniform distribution is optimal in terms of achieving the highest $(\epsilon, \delta)$-data-privacy. We also prove that when all information used in the consensus process is available, the privacy will be compromised. Finally, an optimal privacy-preserving average consensus (OPAC) algorithm is proposed to achieve the highest $(\epsilon, \delta)$-data-privacy and avoid the privacy compromission. Simulations are conducted to verify the results.


title: Generalized Kalman Smoothing: Modeling and Algorithms

abstract: State-space smoothing has found many applications in science and engineering. Under linear and Gaussian assumptions, smoothed estimates can be obtained using efficient recursions, for example Rauch-Tung-Striebel and Mayne-Fraser algorithms. Such schemes are equivalent to linear algebraic techniques that minimize a convex quadratic objective function with structure induced by the dynamic model. These classical formulations fall short in many important circumstances. For instance, smoothers obtained using quadratic penalties can fail when outliers are present in the data, and cannot track impulsive inputs and abrupt state changes. Motivated by these shortcomings, generalized Kalman smoothing formulations have been proposed in the last few years, replacing quadratic models with more suitable, often nonsmooth, convex functions. In contrast to classical models, these general estimators require use of iterated algorithms, and these have received increased attention from control, signal processing, machine learning, and optimization communities. In this survey we show that the optimization viewpoint provides the control and signal processing community great freedom in the development of novel modeling and inference frameworks for dynamical systems. We discuss general statistical models for dynamic systems, making full use of nonsmooth convex penalties and constraints, and providing links to important models in signal processing and machine learning. We also survey optimization techniques for these formulations, paying close attention to dynamic problem structure. Modeling concepts and algorithms are illustrated with numerical examples.


title: Automorphic cohomology, motivic cohomology, and the adjoint $L$-function

abstract: We propose a relationship between the cohomology of arithmetic groups, and the motivic cohomology of certain (Langlands-)attached motives. The motivic cohomology group in question is that related, by Beilinson's conjecture, to the adjoint $L$-function at $s=1$. We present evidence for the conjecture using the theory of periods of automorphic forms, and using analytic torsion.


title: Robust Estimation of Multiple Inlier Structures

abstract: The robust estimator presented in this paper processes each structure independently. The scales of the structures are estimated adaptively and no threshold is involved in spite of different objective functions. The user has to specify only the number of elemental subsets for random sampling. After classifying all the input data, the segmented structures are sorted by their strengths and the strongest inlier structures come out at the top. Like any robust estimators, this algorithm also has limitations which are described in detail. Several synthetic and real examples are presented to illustrate every aspect of the algorithm.


title: On top dimensional Lyubeznik numbers in mixed characteristic

abstract: We prove that the top mixed characteristic Lyubeznik number of a ring $S$ that is a quotient of a complete unramified regular local ring of mixed characteristc with algebraically closed residue field is $1$ provided that depth $S \geq 2$ and dim $S \geq 3$ using a second vanishing theorem in mixed characteristic proved in arXiv:1609.05846


title: Optimal Circuit-Level Decoding for Surface Codes

abstract: Surface codes exploit topological protection to increase error resilience in quantum computing devices and can in principle be implemented in existing hardware. They are one of the most promising candidates for active error correction, not least due to a polynomial-time decoding algorithm which admits one of the highest predicted error thresholds. We consider the dependence of this threshold on underlying assumptions including different noise models, and analyze the performance of a minimum weight perfect matching (MWPM) decoding compared to a mathematically optimal maximum likelihood (ML) decoding. Our ML algorithm tracks the success probabilities for all possible corrections over time and accounts for individual gate failure probabilities and error propagation due to the syndrome measurement circuit. We present the very first evidence for the true error threshold of an optimal circuit level decoder, allowing us to draw conclusions about what kind of improvements are possible over standard MWPM.


title: A Consumer BCI for Automated Music Evaluation Within a Popular On-Demand Music Streaming Service - Taking Listener's Brainwaves to Extremes

abstract: We investigated the possibility of using a machine-learning scheme in conjunction with commercial wearable EEG-devices for translating listener's subjective experience of music into scores that can be used for the automated annotation of music in popular on-demand streaming services. Based on the established -neuroscientifically sound- concepts of brainwave frequency bands, activation asymmetry index and cross-frequency-coupling (CFC), we introduce a Brain Computer Interface (BCI) system that automatically assigns a rating score to the listened song. Our research operated in two distinct stages: i) a generic feature engineering stage, in which features from signal-analytics were ranked and selected based on their ability to associate music induced perturbations in brainwaves with listener's appraisal of music. ii) a personalization stage, during which the efficiency of ex- treme learning machines (ELMs) is exploited so as to translate the derived pat- terns into a listener's score. Encouraging experimental results, from a pragmatic use of the system, are presented.


title: A Theory of Interactive Debugging of Knowledge Bases in Monotonic Logics

abstract: A broad variety of knowledge-based applications such as recommender, expert, planning or configuration systems usually operate on the basis of knowledge represented by means of some logical language. Such a logical knowledge base (KB) enables intelligent behavior of such systems by allowing them to automatically reason, answer queries of interest or solve complex real-world problems. Nowadays, where information acquisition comes at low costs and often happens automatically, the applied KBs are continuously growing in terms of size, information content and complexity. These developments foster the emergence of errors in these KBs and thus pose a significant challenge on all people and tools involved in KB evolution, maintenance and application. If some minimal quality criteria such as logical consistency are not met by some KB, it becomes useless for knowledge-based applications. To guarantee the compliance of KBs with given requirements, (non-interactive) KB debuggers have been proposed. These however often cannot localize all potential faults, suggest too large or incorrect modifications of the faulty KB or suffer from poor scalability due to the inherent complexity of the KB debugging problem. As a remedy to these issues, based on a well-founded theoretical basis this work proposes complete, sound and optimal methods for the interactive debugging of KBs that suggest the one (minimally invasive) error correction of the faulty KB that yields a repaired KB with exactly the intended semantics. Users, e.g. domain experts, are involved in the debugging process by answering automatically generated queries whether some given statements must or must not hold in the domain that should be modeled by the problematic KB at hand.


title: Tidally Induced Pulsations in Kepler Eclipsing Binary KIC 3230227

abstract: KIC 3230227 is a short period ($P\approx 7.0$ days) eclipsing binary with a very eccentric orbit ($e=0.6$). From combined analysis of radial velocities and {\it Kepler} light curves, this system is found to be composed of two A-type stars, with masses of $M_1=1.84\pm 0.18M_{\odot}$, $M_2=1.73\pm 0.17M_{\odot}$ and radii of $R_1=2.01\pm 0.09R_{\odot}$, $R_2=1.68\pm 0.08 R_{\odot}$ for the primary and secondary, respectively. In addition to an eclipse, the binary light curve shows a brightening and dimming near periastron, making this a somewhat rare eclipsing heartbeat star system. After removing the binary light curve model, more than ten pulsational frequencies are present in the Fourier spectrum of the residuals, and most of them are integer multiples of the orbital frequency. These pulsations are tidally driven, and both the amplitudes and phases are in agreement with predictions from linear tidal theory for $l=2, m=-2$ prograde modes.


title: Geometry-Based Next Frame Prediction from Monocular Video

abstract: We consider the problem of next frame prediction from video input. A recurrent convolutional neural network is trained to predict depth from monocular video input, which, along with the current video image and the camera trajectory, can then be used to compute the next frame. Unlike prior next-frame prediction approaches, we take advantage of the scene geometry and use the predicted depth for generating the next frame prediction. Our approach can produce rich next frame predictions which include depth information attached to each pixel. Another novel aspect of our approach is that it predicts depth from a sequence of images (e.g. in a video), rather than from a single still image. We evaluate the proposed approach on the KITTI dataset, a standard dataset for benchmarking tasks relevant to autonomous driving. The proposed method produces results which are visually and numerically superior to existing methods that directly predict the next frame. We show that the accuracy of depth prediction improves as more prior frames are considered.


title: Linear-time string indexing and analysis in small space

abstract: The field of succinct data structures has flourished over the last 16 years. Starting from the compressed suffix array (CSA) by Grossi and Vitter (STOC 2000) and the FM-index by Ferragina and Manzini (FOCS 2000), a number of generalizations and applications of string indexes based on the Burrows-Wheeler transform (BWT) have been developed, all taking an amount of space that is close to the input size in bits. In many large-scale applications, the construction of the index and its usage need to be considered as one unit of computation. Efficient string indexing and analysis in small space lies also at the core of a number of primitives in the data-intensive field of high-throughput DNA sequencing. We report the following advances in string indexing and analysis. We show that the BWT of a string $T\in \{1,\ldots,\sigma\}^n$ can be built in deterministic $O(n)$ time using just $O(n\log{\sigma})$ bits of space, where $\sigma \leq n$. Within the same time and space budget, we can build an index based on the BWT that allows one to enumerate all the internal nodes of the suffix tree of $T$. Many fundamental string analysis problems can be mapped to such enumeration, and can thus be solved in deterministic $O(n)$ time and in $O(n\log{\sigma})$ bits of space from the input string. We also show how to build many of the existing indexes based on the BWT, such as the CSA, the compressed suffix tree (CST), and the bidirectional BWT index, in randomized $O(n)$ time and in $O(n\log{\sigma})$ bits of space. The previously fastest construction algorithms for BWT, CSA and CST, which used $O(n\log{\sigma})$ bits of space, took $O(n\log{\log{\sigma}})$ time for the first two structures, and $O(n\log^{\epsilon}n)$ time for the third, where $\epsilon$ is any positive constant. Contrary to the state of the art, our bidirectional BWT index supports every operation in constant time per element in its output.


title: Global Caching for the Alternation-free $\mu$-Calculus

abstract: We present a sound, complete, and optimal single-pass tableau algorithm for the alternation-free $\mu$-calculus. The algorithm supports global caching with intermediate propagation and runs in time $2^{\mathcal{O}(n)}$. In game-theoretic terms, our algorithm integrates the steps for constructing and solving the B\"uchi game arising from the input tableau into a single procedure; this is done on-the-fly, i.e. may terminate before the game has been fully constructed. This suggests a slogan to the effect that global caching = game solving on-the-fly. A prototypical implementation shows promising initial results.


title: Recognizing Implicit Discourse Relations via Repeated Reading: Neural Networks with Multi-Level Attention

abstract: Recognizing implicit discourse relations is a challenging but important task in the field of Natural Language Processing. For such a complex text processing task, different from previous studies, we argue that it is necessary to repeatedly read the arguments and dynamically exploit the efficient features useful for recognizing discourse relations. To mimic the repeated reading strategy, we propose the neural networks with multi-level attention (NNMA), combining the attention mechanism and external memories to gradually fix the attention on some specific words helpful to judging the discourse relations. Experiments on the PDTB dataset show that our proposed method achieves the state-of-art results. The visualization of the attention weights also illustrates the progress that our model observes the arguments on each level and progressively locates the important words.


title: Consensus-based Privacy-preserving Data Aggregation

abstract: Privacy-preserving data aggregation in ad hoc networks is a challenging problem, considering the distributed communication and control requirement, dynamic network topology, unreliable communication links, etc. Different from the widely used cryptographic approaches, in this paper, we address this challenging problem by exploiting the distributed consensus technique. We first propose a secure consensus-based data aggregation (SCDA) algorithm that guarantees an accurate sum aggregation while preserving the privacy of sensitive data. Then, we prove that the proposed algorithm converges accurately and is $(\epsilon, \sigma)$-data-privacy, and the mathematical relationship between $\epsilon$ and $\sigma$ is provided. Extensive simulations have shown that the proposed algorithm has high accuracy and low complexity, and they are robust against network dynamics.


title: Finding Model-Checkable Needles in Large Source Code Haystacks: Modular Bug-Finding via Static Analysis and Dynamic Invariant Discovery

abstract: In this paper, we present a novel marriage of static and dynamic analysis. Given a large code base with many functions and a mature test suite, we propose using static analysis to find functions 1) with assertions or other evident correctness properties (e.g., array bounds requirements or pointer access) and 2) with simple enough control flow and data use to be amenable to predicate-abstraction based or bounded model checking without human intervention. Because most such functions in realistic software systems in fact rely on many input preconditions not specified by the language's type system (or annotated in any way), we propose using dynamically discovered invariants based on a program's test suite to characterize likely preconditions, in order to reduce the problem of false positives. While providing little in the way of verification, this approach may provide an additional quick and highly scalable bug-finding method for programs that are usually considered "too large to model check." We present a simple example showing that the technique can be useful for a more typically "model-checkable" code base, even in the presence of a poorly designed test suite and bad invariants.


title: Omnidirectional Gravitational Wave Detector with a Laser-Interferometric Gravitational Compass

abstract: Based on the Szekeres-Pirani gravitational compass we suggest the addition of a fourth, non-coplanar mass/mirror to the presently existing laser based gravitational wave observatories, enabling them to operate omnidirectionally, to filter out ambiguous interpretations and to point out the direction of the gravitational wave source.


title: Yarkovsky V-shape identification of asteroid families

abstract: There are only a few known main belt (MB) asteroid families with ages greater than 2 Gyr (Bro\v{z} et al. 2013, Spoto et al. 2015). Estimates based on the family producing collision rate suggest that the lack of >2 Gyr-old families may be due to a selection bias in current techniques used to identify families. Family fragments disperse in their orbital elements, semi-major axis, $a$, eccentricity, $e$, and inclination, $i$, due to secular resonances, close encounters with massive asteroids and the non-gravitational Yarkovsky force. This causes the family fragments to be indistinguishable from the background of the main belt making them more difficult to identify with the hierarchical clustering method (HCM) with increasing family age. The discovery of the Eulalia and new Polana families in the inner belt relied on new techniques because Yarkovsky spreading made them too disperse to be identified using the classical HCM. The technique used to discover the new Polana and Eulalia families is modified to identify asteroid families by searching for correlations between $a$ and asteroid diameter, $D$, or absolute magnitude, $H$. A group of asteroids is identified as a collisional family if its boundary in the $a$ vs. 1/$D$ or $a$ vs. $H$ planes has a characteristic V-shape which is due to the size dependent Yarkovsky spreading. The V-shape boundary is identified with two separate techniques. The first identifies a border by measuring the drop between the number of objects inside and outside of the border. The second identifies the V-shape border by measuring the number density of objects in $a$ vs. 1/$D$,$H$ space. The V-shape techniques are demonstrated on known families Erigone, Vesta, Koronis, and families difficult to identify by HCM such as Flora, Baptistina, new Polana, Eulalia and Karin. Application of the technique in a search for >2 Gyr-old families throughout the MB is discussed.


title: Multiclass Classification Calibration Functions

abstract: In this paper we refine the process of computing calibration functions for a number of multiclass classification surrogate losses. Calibration functions are a powerful tool for easily converting bounds for the surrogate risk (which can be computed through well-known methods) into bounds for the true risk, the probability of making a mistake. They are particularly suitable in non-parametric settings, where the approximation error can be controlled, and provide tighter bounds than the common technique of upper-bounding the 0-1 loss by the surrogate loss. The abstract nature of the more sophisticated existing calibration function results requires calibration functions to be explicitly derived on a case-by-case basis, requiring repeated efforts whenever bounds for a new surrogate loss are required. We devise a streamlined analysis that simplifies the process of deriving calibration functions for a large number of surrogate losses that have been proposed in the literature. The effort of deriving calibration functions is then surmised in verifying, for a chosen surrogate loss, a small number of conditions that we introduce. As case studies, we recover existing calibration functions for the well-known loss of Lee et al. (2004), and also provide novel calibration functions for well-known losses, including the one-versus-all loss and the logistic regression loss, plus a number of other losses that have been shown to be classification-calibrated in the past, but for which no calibration function had been derived.


title: Cat codes with optimal decoherence suppression for a lossy bosonic channel

abstract: We investigate cat codes that can correct multiple excitation losses and identify two types of logical errors: bit-flip errors due to excessive excitation loss and dephasing errors due to quantum back-action from the environment. We show that selected choices of logical subspace and coherent amplitude can efficiently reduce dephasing errors. The trade-off between the two major errors enables optimized performance of cat codes in terms of minimized decoherence. With high coupling efficiency, we show that one-way quantum repeaters with cat codes feature drastically boosted secure communication rate per mode compared with conventional encoding schemes, and thus showcase the promising potential of quantum information processing with continuous variable quantum codes.


title: $\gamma$ Doradus Pulsations in the Eclipsing Binary Star KIC 6048106

abstract: We present the ${\it Kepler}$ photometry of KIC 6048106 exhibiting O'Connell effect and multiperiodic pulsations. Including a starspot on either of the components, light-curve synthesis indicates that this system is a semi-detached Algol with a mass ratio of 0.211, an orbital inclination of 73.9 deg, and a large temperature difference of 2,534 K. To examine in detail both spot variations and pulsations, we separately analyzed the {\it Kepler} time-series data at the interval of an orbital period by an iterative way. The results reveal that the variable asymmetries of the light maxima can be interpreted as the changes of a magnetic cool spot on the secondary component with time. Multiple frequency analyses were performed in the outside-eclipse light residuals after removal of the binarity effects from the observed {\it Kepler} data. We detected 30 frequencies with signal to noise amplitude ratios larger than 4.0, of which six ($f_2$--$f_6$ and $f_{10}$) can be identified as high-order (17 $\le n \le$ 25) low-degree ($\ell$ = 2) gravity-mode pulsations that were stable during the observing run of 200 d. In contrast, the other frequencies may be harmonic and combination terms. For the six frequencies, the pulsation periods and pulsation constants are in the ranges of 0.352$-$0.506 d and 0.232$-$0.333 d, respectively. These values and the position on the HR diagram demonstrate that the primary star is a $\gamma$ Dor variable. The evolutionary status and the pulsation nature of KIC 6048106 are discussed.


title: Enhancing stellar spectroscopy with extreme adaptive optics and photonics

abstract: Extreme adaptive optics systems are now in operation across the globe. These systems, capable of high order wavefront correction, deliver Strehl ratios of 90% in the near-infrared. Originally intended for the direct imaging of exoplanets, these systems are often equipped with advanced coronagraphs that suppress the on-axis-star, interferometers to calibrate wavefront errors, and low order wavefront sensors to stabilize any tip/tilt residuals to a degree never seen before. Such systems are well positioned to facilitate the detailed spectroscopic characterization of faint substellar companions at small angular separations from the host star. Additionally, the increased light concentration of the point-spread function and the unprecedented stability create opportunities in other fields of astronomy as well, including spectroscopy. With such Strehl ratios, efficient injection into single-mode fibers or photonic lanterns becomes possible. With diffraction-limited components feeding the instrument, calibrating a spectrograph's line profile becomes considerably easier, as modal noise or imperfect scrambling of the fiber output are no longer an issue. It also opens up the possibility of exploiting photonic technologies for their advanced functionalities, inherent replicability, and small, lightweight footprint to design and build future instrumentation. In this work, we outline how extreme adaptive optics systems will enable advanced photonic and diffraction-limited technologies to be exploited in spectrograph design and the impact it will have on spectroscopy. We illustrate that the precision of an instrument based on these technologies, would be limited by the spectral content and stellar noise on cool stars and capable of achieving a radial velocity precision of several m/s; the level required for detecting an exo-Earth in the HZ of a nearby M-dwarf.


title: Quantum Locality in Game Strategy

abstract: Game theory is a well established branch of mathematics whose formalism has a vast range of applications from the social sciences, biology, to economics. Motivated by quantum information science, there has been a leap in the formulation of novel game strategies that lead to new (quantum Nash) equilibrium points whereby players in some classical games are always outperformed if sharing and processing joint information ruled by the laws of quantum physics is allowed. We show that, for a bipartite non zero-sum game, input local quantum correlations, and separable states in particular, suffice to achieve an advantage over any strategy that uses classical resources, thus dispensing with quantum nonlocality, entanglement, or even discord between the players' input states. This highlights the remarkable key role played by pure quantum coherence at powering some protocols. Finally, we propose an experiment that uses separable states and basic photon interferometry to demonstrate the locally-correlated quantum advantage.


title: Learning HMMs with Nonparametric Emissions via Spectral Decompositions of Continuous Matrices

abstract: Recently, there has been a surge of interest in using spectral methods for estimating latent variable models. However, it is usually assumed that the distribution of the observations conditioned on the latent variables is either discrete or belongs to a parametric family. In this paper, we study the estimation of an $m$-state hidden Markov model (HMM) with only smoothness assumptions, such as H\"olderian conditions, on the emission densities. By leveraging some recent advances in continuous linear algebra and numerical analysis, we develop a computationally efficient spectral algorithm for learning nonparametric HMMs. Our technique is based on computing an SVD on nonparametric estimates of density functions by viewing them as \emph{continuous matrices}. We derive sample complexity bounds via concentration results for nonparametric density estimation and novel perturbation theory results for continuous matrices. We implement our method using Chebyshev polynomial approximations. Our method is competitive with other baselines on synthetic and real problems and is also very computationally efficient.


title: A projective variety with discrete, non-finitely generated automorphism group

abstract: We construct a projective variety with discrete, non-finitely generated automorphism group. As an application, we show that there exists a complex projective variety with infinitely many non-isomorphic real forms.


title: Ultrafast Magnetization Reversal by Picosecond Electrical Pulses

abstract: The field of spintronics involves the study of both spin and charge transport in solid state devices with a view toward increasing their functionality and efficiency. Alternatively, the field of ultrafast magnetism focuses on the use of femtosecond laser pulses to excite electrons in magnetic materials, which allows the magnetic order to be dramatically changed on unprecedented sub-picosecond time-scales. Here, we unite these two distinct research activities by using picosecond electrical pulses to rapidly excite electrons in a magnetic metal. We are able to deterministically and repetitively reverse the magnetization of a GdFeCo film with sub-10 picosecond electrical pulses. The magnetization reverses in ~10ps, which is more than an order of magnitude faster than any other electrically controlled magnetic switching. We attribute the deterministic switching of the magnetization to ultrafast excitation of the electrons, a fundamentally different mechanism from other current driven switching mechanisms such as spin-transfer-torque (STT) or spin-orbit-torque (SOT). The energy density required for switching is measured and the process is found to be efficient, projecting to only 4 fJ needed to switch a (20 nm)^3 cell, which is comparable to other state-of-the-art STT-MRAM memory devices. This discovery will launch a new field of research into picosecond spintronic phenomena and devices.


title: Measurement of $CP$ asymmetry in the $D^{0} \to K^0_S K^0_S$ decay at Belle

abstract: We report a measurement of the time-integrated $CP$ asymmetry in the neutral charm meson decay $D^0 \to K^0_S K^0_S$ using 921~fb$^{-1}$ data collected at the $\Upsilon(4S)$ and $\Upsilon(5S)$ resonances with the Belle detector at the KEKB asymmetric-energy $e^+e^-$ collider. The observed asymmetry is $$ A_{CP}(D^0 \to K^0_S K^0_S) = (-0.02 \pm 1.53 \pm 0.17) \%, $$ where the first uncertainty is statistical and the second systematic. This latter uncertainty is dominated by the error of the normalisation channel. The result is consistent with Standard Model expectations and improves the uncertainty with respect to previous measurement of this quantity by more than a factor of three.


title: Existence and nonexistence of solutions for the heat equation with a superlinear source term

abstract: Classification theory on the existence and non-existence of local in time solutions for initial value problems of nonlinear heat equations are investigated. Without assuming a concrete growth rate on a nonlinear term, we reveal the threshold integrability of initial data which classify existence and nonexistence of solutions via a quasi-scaling and its invariant integral. Typical nonlinear terms, for instance polynomial type, exponential type and its sum, product and composition, can be treated as applications.


title: Analyzing Dependent Placements of Small Cells in a Two-Layer Heterogeneous Network with a Rate Coverage Constraint

abstract: We consider the downlink of a two-layer heterogeneous network, comprising macro cells (MCs) and small cells (SCs). The existing literature generally assumes independent placements of the access points (APs) in different layers; in contrast, we analyze a dependent placement where SC APs are placed at locations with poor service from the MC layer. Our goal is to obtain an estimate of the number of SCs required to maintain a target outage rate. Such an analysis is trivial if the MCs are located according to a Poisson point process (PPP), which provides a lower bound on performance. Here, we consider MCs placed on a hexagonal grid, which complements the PPP model by providing an upper bound on performance. We first provide accurate bounds for the average interference within a MC when SCs are not used. Then, by obtaining the outage areas, we estimate the number of SCs required within an MC to overcome outage. If resource allocation amongst SCs is not used, we show that the problem of outage is not solved completely, and the residual outage area depends on whether co-channel or orthogonal SCs are used. Simulations show that a much smaller residual outage area is obtained with orthogonal SCs.


title: On the ideal generated by all squarefree monomials of a given degree

abstract: An explicit construction is given of a minimal free resolution of the ideal generated by all squarefree monomials of a given degree. The construction relies upon and exhibits the natural action of the symmetric group on the syzygy modules. The resolution is obtained over an arbitrary coefficient ring; in particular, it is characteristic free. Two applications are given: an equivariant resolution of De Concini-Procesi rings indexed by hook partitions, and a resolution of FI-modules.


title: Universal self-field critical current for thin-film superconductors

abstract: For any practical superconductor the magnitude of the critical current density, $J_\textrm{c}$, is crucially important. It sets the upper limit for current in the conductor. Usually $J_\textrm{c}$ falls rapidly with increasing external magnetic field but even in zero external field the current flowing in the conductor generates a self-field which limits $J_\textrm{c}$. Here we show for thin films of thickness less than the London penetration depth, $\lambda$, this limiting $J_\textrm{c}$ adopts a universal value for all superconductors - metals, oxides, cuprates, pnictides, borocarbides and heavy Fermions. For type I superconductors, it is $H_{\textrm{c}}/\lambda$ where $H_\textrm{c}$ is the thermodynamic critical field. But surprisingly for type II superconductors we find the self-field $J_\textrm{c}$ is $H_{\textrm{c}1}/\lambda$ where $H_{\textrm{c}1}$ is the lower critical field. $J_\textrm{c}$ is thus fundamentally determined and this provides a simple means to extract absolute values of $\lambda(T)$ and, from its temperature dependence, the symmetry and magnitude of the superconducting gap.


title: Accretion of Supersonic Winds on Boson Stars

abstract: We present the evolution of a supersonic wind interacting with a Boson Star (BS) and compare the resulting wind density profile with that of the shock cone formed when the wind is accreted by a non-rotating Black Hole (BH) of the same mass. The physical differences between these accretors are that a BS, unlike a BH has no horizon, it does not have a mechanical surface either and thus the wind is expected to trespass the BS. Despite these conditions, on the BS space-time the gas achieves a stationary flux with the gas accumulating in a high density elongated structure comparable to the shock cone formed behind a BH. The highest density resides in the center of the BS whereas in the case of the BH it is found on the downstream part of the BH near the event horizon. The maximum density of the gas is smaller in the BS than in the BH case. Our results indicate that the highest density of the wind is more similar on the BS to that on the BH when the BS has high self-interaction, when it is more compact and when the wind velocity is higher. We expect this and similar analyses help to know whether BSs can still be Black Hole mimickers or can be ruled out.


title: Stochastic dynamics and the predictability of big hits in online videos

abstract: The competition for the attention of users is a central element of the Internet. Crucial issues are the origin and predictability of big hits, the few items that capture a big portion of the total attention. We address these issues analyzing 10 million time series of videos' views from YouTube. We find that the average gain of views is linearly proportional to the number of views a video already has, in agreement with usual rich-get-richer mechanisms and Gibrat's law, but this fails to explain the prevalence of big hits. The reason is that the fluctuations around the average views are themselves heavy tailed. Based on these empirical observations, we propose a stochastic differential equation with L\'evy noise as a model of the dynamics of videos. We show how this model is substantially better in estimating the probability of an ordinary item becoming a big hit, which is considerably underestimated in the traditional proportional-growth models.


title: Efficient silicon metasurfaces for visible light

abstract: Dielectric metasurfaces require high refractive index contrast materials for optimum performance. This requirement imposes a severe restraint; devices have either been demonstrated at wavelengths of 700nm and above using high-index semiconductors such as silicon, or they use lower index dielectric materials such as TiO$_{2}$ or Si$_{3}$N$_{4}$ and operate in the visible wavelength regime. Here, we show that the high refractive index of silicon can be exploited at wavelengths as short as 532 nm by demonstrating a silicon metasurface with a transmission efficiency of 47% at this wavelength. The metasurface consists of a graded array of silicon posts arranged in a square lattice on a quartz substrate. We show full 2{\pi} phase control and we experimentally demonstrate polarization-independent beam deflection at 532nm wavelength. The crystalline silicon is placed on a quartz substrate by a bespoke layer transfer technique and we note that an efficiency >70% may be achieved for a further optimized structure in the same material. Our results open a new way for realizing efficient metasurfaces based on silicon in the visible wavelength regime.


title: Data-driven discovery of partial differential equations

abstract: We propose a sparse regression method capable of discovering the governing partial differential equation(s) of a given system by time series measurements in the spatial domain. The regression framework relies on sparsity promoting techniques to select the nonlinear and partial derivative terms terms of the governing equations that most accurately represent the data, bypassing a combinatorially large search through all possible candidate models. The method balances model complexity and regression accuracy by selecting a parsimonious model via Pareto analysis. Time series measurements can be made in an Eulerian framework where the sensors are fixed spatially, or in a Lagrangian framework where the sensors move with the dynamics. The method is computationally efficient, robust, and demonstrated to work on a variety of canonical problems of mathematical physics including Navier-Stokes, the quantum harmonic oscillator, and the diffusion equation. Moreover, the method is capable of disambiguating between potentially non-unique dynamical terms by using multiple time series taken with different initial data. Thus for a traveling wave, the method can distinguish between a linear wave equation or the Korteweg-deVries equation, for instance. The method provides a promising new technique for discovering governing equations and physical laws in parametrized spatio-temporal systems where first-principles derivations are intractable.


title: Exact Sampling of the Infinite Horizon Maximum of a Random Walk Over a Non-linear Boundary

abstract: We present the first algorithm that samples $\max_{n\geq0}\{S_{n}-n^{\alpha}\},$ where $S_n$ is a mean zero random walk, and $n^{\alpha}$ with $\alpha\in(1/2,1)$ defines a nonliner boundary. We show that our algorithm has finite expected running time. We also apply the algorithm to construct the first exact simulation method for the steady-state departure process of a $GI/GI/\infty$ queue where the service time distribution has infinite mean.


title: Exact rule-based stochastic simulations for systems with unlimited number of molecular species

abstract: We introduce expandable partial propensity direct method (EPDM) - a new exact stochastic simulation algorithm suitable for systems involving many interacting molecular species. The algorithm is especially efficient for sparsely populated systems, where the number of species that may potentially be generated is much greater than the number of species actually present in the system at any given time. The number of operations per reaction scales linearly with the number of species, but only those which have one or more molecules. To achieve this kind of performance we are employing a data structure which allows to add and remove species and their interactions on the fly. When a new specie is added, its interactions with every other specie are generated dynamically by a set of user-defined rules. By removing the records involving the species with zero molecules, we keep the number of species as low as possible. This enables simulations of systems for which listing all species is not practical. The algorithm is based on partial propensities direct method (PDM) by Ramaswamy et al. for sampling trajectories of the chemical master equation.


title: KU-ISPL Language Recognition System for NIST 2015 i-Vector Machine Learning Challenge

abstract: In language recognition, the task of rejecting/differentiating closely spaced versus acoustically far spaced languages remains a major challenge. For confusable closely spaced languages, the system needs longer input test duration material to obtain sufficient information to distinguish between languages. Alternatively, if languages are distinct and not acoustically/linguistically similar to others, duration is not a sufficient remedy. The solution proposed here is to explore duration distribution analysis for near/far languages based on the Language Recognition i-Vector Machine Learning Challenge 2015 (LRiMLC15) database. Using this knowledge, we propose a likelihood ratio based fusion approach that leveraged both score and duration information. The experimental results show that the use of duration and score fusion improves language recognition performance by 5% relative in LRiMLC15 cost.


title: A Logic of Knowing Why

abstract: When we say "I know why he was late", we know not only the fact that he was late, but also an explanation of this fact. We propose a logical framework of "knowing why" inspired by the existing formal studies on why-questions, scientific explanation, and justification logic. We introduce the Ky_i operator into the language of epistemic logic to express "agent i knows why phi" and propose a Kripke-style semantics of such expressions in terms of knowing an explanation of phi. We obtain two sound and complete axiomatizations w.r.t. two different model classes depending on different assumptions about introspection.


title: New views of the distant stellar halo

abstract: Currently only a small number of Milky Way (MW) stars are known to exist beyond 100 kpc from the Galactic center. Though the distribution of these stars in the outer halo is believed to be sparse, they can provide evidence of more recent accretion events than in the inner halo and help map out the MW's dark matter halo to its virial radius. We have re-examined the outermost regions of 11 existing stellar halo models with two synthetic surveys: one mimicking present-day searches for distant M giants and another mimicking RR Lyrae (RRLe) projections for LSST. Our models suggest that color and proper motion cuts currently used to select M giant candidates for follow-up successfully remove nearly all halo dwarf self-contamination and are useful for focusing observations on distant M giants, of which there are thousands to tens of thousands beyond 100 kpc in our models. We likewise expect that LSST will identify comparable numbers of RRLe at these distances. We demonstrate that several observable properties of both tracers, such as proximity of neighboring stars, proper motions, and distances (for RRLe) could help us separate different accreted structures from one another. We also discuss prospects for using ratios of M giants to RRLe as a proxy for accretion time, which in the future could provide new constraints on the recent accretion history of our Galaxy.


title: Majorana modes in solid state systems and its dynamics

abstract: We review the properties of Majorana fermions in particle physics and point out that Majorana modes in solid state systems are significantly different. The key reason is the concept of anti-particle in solid state systems is different from its counterpart in particle physics. We define Majorana modes as the eigenstates of Majorana operators and find that they can exist both at edges and in the bulk. According to our definition, only one single Majorana mode can exist in a system no matter at edges or in the bulk. Kitaev's spinless $p$-wave superconductor is used to illustrate our results and the dynamical behavior of the Majorana modes.


title: Control Barrier Function Based Quadratic Programs for Safety Critical Systems

abstract: Safety critical systems involve the tight coupling between potentially conflicting control objectives and safety constraints. As a means of creating a formal framework for controlling systems of this form, and with a view toward automotive applications, this paper develops a methodology that allows safety conditions -- expressed as control barrier functions -- to be unified with performance objectives -- expressed as control Lyapunov functions -- in the context of real-time optimization-based controllers. Safety conditions are specified in terms of forward invariance of a set, and are verified via two novel generalizations of barrier functions; in each case, the existence of a barrier function satisfying Lyapunov-like conditions implies forward invariance of the set, and the relationship between these two classes of barrier functions is characterized. In addition, each of these formulations yields a notion of control barrier function (CBF), providing inequality constraints in the control input that, when satisfied, again imply forward invariance of the set. Through these constructions, CBFs can naturally be unified with control Lyapunov functions (CLFs) in the context of a quadratic program (QP); this allows for the achievement of control objectives (represented by CLFs) subject to conditions on the admissible states of the system (represented by CBFs). The mediation of safety and performance through a QP is demonstrated on adaptive cruise control and lane keeping, two automotive control problems that present both safety and performance considerations coupled with actuator bounds.


title: Why are pulsar planets rare?

abstract: Pulsar timing observations have revealed planets around only a few pulsars. We suggest that the rarity of these planets is due mainly to two effects. First, we show that the most likely formation mechanism requires the destruction of a companion star. Only pulsars with a suitable companion (with an extreme mass ratio) are able to form planets. Second, while a dead zone (a region of low turbulence) in the disk is generally thought to be essential for planet formation, it is most probably rare in disks around pulsars because of the irradiation from the pulsar. The irradiation strongly heats the inner parts of the disk pushing the inner boundary of the dead zone out. We suggest that the rarity of pulsar planets can be explained by the low probability for these two requirements - a very low-mass companion and a dead zone - to be satisfied.