title: A no-go theorem for theories that decohere to quantum mechanics

abstract: To date, there has been no experimental evidence that invalidates quantum theory. Yet it may only be an effective description of the world, in the same way that classical physics is an effective description of the quantum world. We ask whether there exists an operationally-defined theory superseding quantum theory, but which reduces to it via a decoherence-like mechanism. We prove that no such post-quantum theory exists if it is demanded that it satisfy two natural physical principles: causality and purification. Causality formalises the statement that information propagates from present to future, and purification that each state of incomplete information arises in an essentially unique way due to lack of information about an environment. Hence, our result can either be viewed as evidence that the fundamental theory of Nature is quantum, or as showing in a rigorous manner that any post-quantum theory must abandon causality, purification, or both.


title: X-ray Binaries

abstract: This chapter discusses the implications of X-ray binaries on our knowledge of Type Ibc and Type II supernovae. X-ray binaries contain accreting neutron stars and stellar--mass black holes which are the end points of massive star evolution. Studying these remnants thus provides clues to understanding the evolutionary processes that lead to their formation. We focus here on the distributions of dynamical masses, space velocities and chemical anomalies of their companion stars. These three observational features provide unique information on the physics of core collapse and supernovae explosions within interacting binary systems. There is suggestive evidence for a gap between ~2-5 Msun in the observed mass distribution. This might be related to the physics of the supernova explosions although selections effects and possible systematics may be important. The difference between neutron star mass measurements in low-mass X-ray binaries (LMXBs) and pulsar masses in high-mass X-ray binaries (HMXBs) reflect their different accretion histories, with the latter presenting values close to birth masses. On the other hand, black holes in LMXBs appear to be limited to <~12 Msun because of strong mass-loss during the wind Wolf-Rayet phase. Detailed studies of a limited sample of black-hole X-ray binaries suggest that the more massive black holes have a lower space velocity, which could be explained if they formed through direct collapse. Conversely, the formation of low-mass black holes through a supernova explosion implies that large escape velocities are possible through ensuing natal and/or Blaauw kicks. Finally, chemical abundance studies of the companion stars in seven X-ray binaries indicate they are metal-rich (all except GRO J1655-40) and possess large peculiar abundances of alpha-elements (Abridged)


title: Stability interchanges in a curved Sitnikov problem

abstract: We consider a curved Sitnikov problem, in which an infinitesimal particle moves on a circle under the gravitational influence of two equal masses in Keplerian motion within a plane perpendicular to that circle. There are two equilibrium points, whose stability we are studying. We show that one of the equilibrium points undergoes stability interchanges as the semi-major axis of the Keplerian ellipses approaches the diameter of that circle. To derive this result, we first formulate and prove a general theorem on stability interchanges, and then we apply it to our model. The motivation for our model resides with the $n$-body problem in spaces of constant curvature.


title: WOMBAT: A Scalable and High Performance Astrophysical MHD Code

abstract: We present a new code for astrophysical magneto-hydrodynamics specifically designed and optimized for high performance and scaling on modern and future supercomputers. We describe a novel hybrid OpenMP/MPI programming model that emerged from a collaboration between Cray, Inc. and the University of Minnesota. This design utilizes MPI-RMA optimized for thread scaling, which allows the code to run extremely efficiently at very high thread counts ideal for the latest generation of the multi-core and many-core architectures. Such performance characteristics are needed in the era of "exascale" computing. We describe and demonstrate our high-performance design in detail with the intent that it may be used as a model for other, future astrophysical codes intended for applications demanding exceptional performance.


title: Iterative methods for solving factorized linear systems

abstract: Stochastic iterative algorithms such as the Kaczmarz and Gauss-Seidel methods have gained recent attention because of their speed, simplicity, and the ability to approximately solve large-scale linear systems of equations without needing to access the entire matrix. In this work, we consider the setting where we wish to solve a linear system in a large matrix X that is stored in a factorized form, X = UV; this setting either arises naturally in many applications or may be imposed when working with large low-rank datasets for reasons of space required for storage. We propose a variant of the randomized Kaczmarz method for such systems that takes advantage of the factored form, and avoids computing X. We prove an exponential convergence rate and supplement our theoretical guarantees with experimental evidence demonstrating that the factored variant yields significant acceleration in convergence.


title: Unifying microscopic and continuum treatments of van der Waals and Casimir interactions

abstract: We present an approach for computing long-range van der Waals (vdW) interactions between complex molecular systems and arbitrarily shaped macroscopic bodies, melding atomistic treatments of electronic fluctuations based on density functional theory in the former, with continuum descriptions of strongly shape-dependent electromagnetic fields in the latter, thus capturing many-body and multiple scattering effects to all orders. Such a theory is especially important when considering vdW interactions at mesoscopic scales, i.e. between molecules and structured surfaces with features on the scale of molecular sizes, in which case the finite sizes, complex shapes, and resulting nonlocal electronic excitations of molecules are strongly influenced by electromagnetic retardation and wave effects that depend crucially on the shapes of surrounding macroscopic bodies. We show that these effects together can modify vdW interactions by orders of magnitude compared to previous treatments based on Casimir--Polder or non-retarded approximations, which are valid only at macroscopically large or atomic-scale separations, respectively.


title: Finite volume calculation of $K$-theory invariants

abstract: Odd index pairings of $K_1$-group elements with Fredholm modules are of relevance in index theory, differential geometry and applications such as to topological insulators. For the concrete setting of operators on a Hilbert space over a lattice, it is shown how to calculate the resulting index as the signature of a suitably constructed finite-dimensional matrix, more precisely the finite volume restriction of what we call the spectral localizer. In presence of real symmetries, secondary $\mathbb{Z}_2$-invariants can be obtained as the sign of the Pfaffian of the spectral localizer. These results reconcile two complementary approaches to invariants of topological insulators.


title: Control Allocation for Wide Area Coordinated Damping

abstract: In this work, a modal-based sparse control allocation (CA) is proposed for coordinated and fault-tolerant wide-area damping controllers (WADCs). In our proposed method, the supervisory CA only communicates with necessary actuators to achieve the required damping performance and in case of actuator failures (e.g., due to loss of communication or scheduling), capabilities of the remaining actuators are fully used before the nominal performance is degraded. This method offers the advantages of modular design where WADC is initially designed to achieve satisfactory damping without the detailed knowledge of actuators. In the next step, CA is designed to manage actuator failures and limitations without the need to redesign the nominal WADC. The proposed approach is applied to a modified $286$-bus Western Electricity Coordinating Council (WECC) system to verify the feasibility on a complex power system. Simulation results indicate the effectiveness of the proposed method in coordinating multiple actuators and building resiliency.


title: A Search for Fast Radio Bursts with the GBNCC Pulsar Survey

abstract: We report on a search for Fast Radio Bursts (FRBs) with the Green Bank Northern Celestial Cap (GBNCC) Pulsar Survey at 350 MHz. Pointings amounting to a total on-sky time of 61 days were searched to a DM of 3000 pc cm$^{-3}$ while the rest (23 days; 29% of the total time) were searched to a DM of 500 pc cm$^{-3}$. No FRBs were detected in the pointings observed through May 2016. We estimate a 95% confidence upper limit on the FRB rate of $3.6\times 10^3$ FRBs sky$^{-1}$ day$^{-1}$ above a peak flux density of 0.63 Jy at 350 MHz for an intrinsic pulse width of 5 ms. We place constraints on the spectral index $\alpha$ by running simulations for different astrophysical scenarios and cumulative flux density distributions. The non-detection with GBNCC is consistent with the 1.4-GHz rate reported for the Parkes surveys for $\alpha > +0.35 $ in the absence of scattering and free-free absorption and $\alpha > -0.3$ in the presence of scattering, for a Euclidean flux distribution. The constraints imply that FRBs exhibit either a flat spectrum or a spectral turnover at frequencies above 400 MHz. These constraints also allow estimation of the number of bursts that can be detected with current and upcoming surveys. We predict that CHIME may detect anywhere from several to $\sim$50 FRBs a day (depending on model assumptions), making it well suited for interesting constraints on spectral index, the log $N$-log $S$ slope and pulse profile evolution across its bandwidth (400-800 MHz).


title: Multiparameter quantum Pfaffians

abstract: The multiparameter quantum Pfaffian of the $(p, \lambda)$-quantum group is introduced and studied together with the quantum determinant, and an identity relating the two invariants is given. Generalization to the multiparameter hyper-Pfaffian and relationship with the quantum minors are also considered.


title: Emergence of a stellar cusp by a dark matter cusp in a low-mass compact ultra-faint dwarf galaxy

abstract: Recent observations have been discovering new ultra-faint dwarf galaxies as small as $\sim20~{\rm pc}$ in half-light radius and $\sim3~{\rm km~s^{-1}}$ in line-of-sight velocity dispersion. In these galaxies, dynamical friction on a star against dark matter can be significant and alter their stellar density distribution. The effect can strongly depend on a central density profile of dark matter, i.e. cusp or core. In this study, I perform computations using a classical and a modern analytic formulae and $N$-body simulations to study how dynamical friction changes a stellar density profile and how different it is between cuspy and cored dark matter haloes. This study shows that, if a dark matter halo has a cusp, dynamical friction can cause shrivelling instability which results in emergence of a stellar cusp in the central region $\simeq2~{\rm pc}$. On the other hand, if it has a constant-density core, dynamical friction is significantly weaker and does not generate a stellar cusp even if the galaxy has the same line-of-sight velocity dispersion. In such a compact and low-mass galaxy, since the shrivelling instability by dynamical friction is inevitable if it has a dark matter cusp, absence of a stellar cusp implies that the galaxy has a dark-matter core. I expect that this could be used to diagnose a dark matter density profile in these compact ultra-faint dwarf galaxies.


title: Sums of squares and products of Bessel functions

abstract: Let $r_k(n)$ denote the number of representations of the positive integer $n$ as the sum of $k$ squares. We rigorously prove for the first time a Voronoi summation formula for $r_k(n), k\geq2,$ proved incorrectly by A. I. Popov and later rediscovered by A. P. Guinand, but without proof and without conditions on the functions associated in the transformation. Using this summation formula we establish a new transformation between a series consisting of $r_k(n)$ and a product of two Bessel functions, and a series involving $r_k(n)$ and the Gaussian hypergeometric function. This transformation can be considered as a massive generalization of well-known results of G. H. Hardy, and of A. L. Dixon and W. L. Ferrar, as well as of a classical result of A. I. Popov that was completely forgotten. An analytic continuation of this transformation yields further useful results that generalize those obtained earlier by Dixon and Ferrar.


title: Lower bounds on the quantum Fisher information based on the variance and various types of entropies

abstract: We examine important properties of the difference between the variance and the quantum Fisher information over four, i.e., $(\Delta A)^2-F_{\rm Q}[\varrho,A]/4.$ We find that it is equal to a generalized variance defined in Petz [J. Phys. A 35, 929 (2002)] and Gibilisco, Hiai, and Petz [IEEE Trans. Inf. Theory 55, 439 (2009)]. We present an upper bound on this quantity that is proportional to the linear entropy. As expected, our relation shows that for states that are close to being pure, the quantum Fisher information over four is close to the variance. We also obtain the variance and the quantum Fisher information averaged over all Hermitian operators, and examine its relation to the von Neumann entropy. Apart from the usual quantum Fisher information, we also consider the Kubo-Mori-Bogoliubov quantum Fisher information.


title: Quartic time-dependent oscillatons

abstract: In this paper, we will study some properties of oscillaton, spherically symmetric object made of a real time-dependent scalar field, Using a self- interaction quartic scalar potential instead of a quadratic or exponential ones discussed in previous works. Since the oscillatons can be regarded as models for astrophysical objects which play the role of dark matter, there- fore investigation of their properties has more importance place in present time of physics; research. Therefore we investigate the properties of these objects by Solving the system of differential equations obtained from the Einstein Klein Gordon (EKG) equations and will show their importance as new candidates for the role of dark matter in the galactic scales.


title: Sorting by Reversals and the Theory of 4-Regular Graphs

abstract: We show that the theory of sorting by reversals fits into the well-established theory of circuit partitions of 4-regular multigraphs (which also involves the combinatorial structures of circle graphs and delta-matroids). In this way, we expose strong connections between the two theories that have not been fully appreciated before. We also discuss a generalization of sorting by reversals involving the double-cut-and-join (DCJ) operation. Finally, we also show that the theory of sorting by reversals is closely related to that of gene assembly in ciliates.


title: Tunnelling in Dante's Inferno

abstract: We study quantum tunnelling in Dante's Inferno model of large field inflation. Such a tunnelling process, which will terminate inflation, becomes problematic if the tunnelling rate is rapid compared to the Hubble time scale at the time of inflation. Consequently, we constrain the parameter space of Dante's Inferno model by demanding a suppressed tunnelling rate during inflation. The constraints are derived and explicit numerical bounds are provided for representative examples. Our considerations are at the level of an effective field theory; hence, the presented constraints have to hold regardless of any UV completion.


title: Relating the finite-volume spectrum and the two-and-three-particle $S$ matrix for relativistic systems of identical scalar particles

abstract: Working in relativistic quantum field theory, we derive the quantization condition satisfied by coupled two- and three-particle systems of identical scalar particles confined to a cubic spatial volume with periodicity $L$. This gives the relation between the finite-volume spectrum and the infinite-volume $\textbf 2 \to \textbf 2$, $\textbf 2 \to \textbf 3$ and $\textbf 3 \to \textbf 3$ scattering amplitudes for such theories. The result holds for relativistic systems composed of scalar particles with nonzero mass $m$, whose center of mass energy lies below the four-particle threshold, and for which the two-particle $K$ matrix has no singularities below the three-particle threshold. The quantization condition is exact up to corrections of the order $\mathcal{O}(e^{-mL})$ and holds for any choice of total momenta satisfying the boundary conditions.


title: The Gaussian Stiffness of Graphene deduced from a Continuum Model based on Molecular Dynamics Potentials

abstract: We consider a discrete model of a graphene sheet with atomic interactions governed by a harmonic approximation of the 2nd-generation Brenner potential that depends on bond lengths, bond angles, and two types of dihedral angles. A continuum limit is then deduced that fully describes the bending behavior. In particular, we deduce for the first time an analytical expression of the Gaussian stiffness, a scarcely investigated parameter ruling the rippling of graphene, for which contradictory values have been proposed in the literature. We disclose the atomic-scale sources of both bending and Gaussian stiffnesses and provide for them quantitative evaluations.


title: Isometries of Clifford algebras II

abstract: Let F be a field of characteristic different from 2, and let $F^{n}$ denote the vector space of n-tuples of elements in F. Let ${e_{1}, ... , e_{n}}$ denote the canonical basis of $F^{n}$. Let r and s be nonnegative integers such that r + s = n, and let Q denote the nondegenerate bilinear form on $F^{n}$ such that $Q(e_{i}, e_{j}) = 0$ if i,j are distinct, $Q(e_{i},e_{i}) = 1$ if $1 \leq i \leq r$ and $Q(e_{r+j},e_{r+j}) = -1$ if $1 \leq j \leq s$. Let $C\ell(r,s)$ denote the Clifford algebra determined by Q and $F^{n}$. There is a canonical extension of Q to a nondegenerate, symmetric, bilinear form $\bar{Q}$ on $C\ell(r,s)$. An element g of $C\ell(r,s)$ will be called an isometry of $C\ell(r,s)$ if left and right translations by g preserve $\bar{Q}$. Let $G_{r,s}$ denote the group of all isometries of $C\ell(r,s)$. We construct a Lie algebra $LG_{r,s}$ over F that equals the Lie algebra of $G_{r,s}$ in the case that F = R or C. The Lie algebra $LG_{r,s}$ admits an involutive automorphism whose +1 and -1 eigenspaces determine a Cartan decomposition $LG_{r,s} = K_{r,s} \oplus P_{r,s}$. We compute the bracket relations for a natural system of generators of $LG_{r,s}$. Finally, we determine $LG_{r,s}$ in the case that F = R.


title: Symplectic resolutions for Higgs moduli spaces

abstract: In this paper, we study the algebraic symplectic geometry of the singular moduli spaces of Higgs bundles of degree $0$ and rank $n$ on a compact Riemann surface $X$ of genus $g$. In particular, we prove that such moduli spaces are symplectic singularities, in the sense of Beauville [Bea00], and admit a projective symplectic resolution if and only if $g=1$ or $(g, n)=(2,2)$. These results are an application of a recent paper by Bellamy and Schedler [BS16] via the so-called Isosingularity Theorem.


title: Numerical Approximations for a three components Cahn-Hilliard phase-field Model based on the Invariant Energy Quadratization method

abstract: How to develop efficient numerical schemes while preserving the energy stability at the discrete level is a challenging issue for the three component Cahn-Hilliard phase-field model. In this paper, we develop first and second order temporal approximation schemes based on the "Invariant Energy Quadratization" approach, where all nonlinear terms are treated semi-explicitly. Consequently, the resulting numerical schemes lead to a well-posed linear system with the symmetric positive definite operator to be solved at each time step. We rigorously prove that the proposed schemes are unconditionally energy stable. Various 2D and 3D numerical simulations are presented to demonstrate the stability and the accuracy of the schemes.


title: Decidability, Complexity, and Expressiveness of First-Order Logic Over the Subword Ordering

abstract: We consider first-order logic over the subword ordering on finite words, where each word is available as a constant. Our first result is that the $\Sigma_1$ theory is undecidable (already over two letters). We investigate the decidability border by considering fragments where all but a certain number of variables are alternation bounded, meaning that the variable must always be quantified over languages with a bounded number of letter alternations. We prove that when at most two variables are not alternation bounded, the $\Sigma_1$ fragment is decidable, and that it becomes undecidable when three variables are not alternation bounded. Regarding higher quantifier alternation depths, we prove that the $\Sigma_2$ fragment is undecidable already for one variable without alternation bound and that when all variables are alternation bounded, the entire first-order theory is decidable.


title: Mixed Higgs-Radion States at the LHC -- a Detailed Study

abstract: Light radions constitute one of the few surviving possibilities for observable new particle states at the sub-TeV level which arise in models with extra spacetime dimensions. It is already known that the 125 GeV state discovered at CERN is unlikely to be a pure radion state, since its decays resemble those of the Standard Model Higgs boson too closely. However, due to experimental errors in the measured decay widths, the possibility still remains that it could be a mixture of the radion with one (or more) Higgs states. We use the existing LHC data at 8 and 13 TeV to make a thorough investigation of this possibility. Not surprisingly, it turns out that this model is already constrained quite effectively by direct LHC searches for an additional scalar heavier than 125 GeV. We then make a detailed study of the so-called 'conformal point', where this heavy state practically decouples from (most of) the Standard Model fields. Some projections for the future are also included.


title: The maximum number of cliques in graphs without long cycles

abstract: The Erd\H{o}s--Gallai Theorem states that for $k\geq 3$ every graph on $n$ vertices with more than $\frac{1}{2}(k-1)(n-1)$ edges contains a cycle of length at least $k$. Kopylov proved a strengthening of this result for 2-connected graphs with extremal examples $H_{n,k,t}$ and $H_{n,k,2}$. In this note, we generalize the result of Kopylov to bound the number of $s$-cliques in a graph with circumference less than $k$. Furthermore, we show that the same extremal examples that maximize the number of edges also maximize the number of cliques of any fixed size. Finally, we obtain the extremal number of $s$-cliques in a graph with no path on $k$-vertices.


title: Implementation of Tetris as a Model Counter

abstract: Solving SharpSAT problems is an important area of work. In this paper, we discuss implementing Tetris, an algorithm originally designed for handling natural joins, as an exact model counter for the SharpSAT problem. Tetris uses a simple geometric framework, yet manages to achieve the fractional hypertree-width bound. Its design allows it to handle complex problems involving extremely large numbers of clauses on which other state-of-the-art model counters do not perform well, yet still performs strongly on standard SAT benchmarks. We have achieved the following objectives. First, we have found a natural set of model counting benchmarks on which Tetris outperforms other model counters. Second, we have constructed a data structure capable of efficiently handling and caching all of the data Tetris needs to work on over the course of the algorithm. Third, we have modified Tetris in order to move from a theoretical, asymptotic-time-focused environment to one that performs well in practice. In particular, we have managed to produce results keeping us within a single order of magnitude as compared to other solvers on most benchmarks, and outperform those solvers by multiple orders of magnitude on others.


title: Exploiting Convolutional Neural Network for Risk Prediction with Medical Feature Embedding

abstract: The widespread availability of electronic health records (EHRs) promises to usher in the era of personalized medicine. However, the problem of extracting useful clinical representations from longitudinal EHR data remains challenging. In this paper, we explore deep neural network models with learned medical feature embedding to deal with the problems of high dimensionality and temporality. Specifically, we use a multi-layer convolutional neural network (CNN) to parameterize the model and is thus able to capture complex non-linear longitudinal evolution of EHRs. Our model can effectively capture local/short temporal dependency in EHRs, which is beneficial for risk prediction. To account for high dimensionality, we use the embedding medical features in the CNN model which hold the natural medical concepts. Our initial experiments produce promising results and demonstrate the effectiveness of both the medical feature embedding and the proposed convolutional neural network in risk prediction on cohorts of congestive heart failure and diabetes patients compared with several strong baselines.


title: Projected Primal-Dual Gradient Flow of Augmented Lagrangian with Application to Distributed Maximization of the Algebraic Connectivity of a Network

abstract: In this paper, a projected primal-dual gradient flow of augmented Lagrangian is presented to solve convex optimization problems that are not necessarily strictly convex. The optimization variables are restricted by a convex set with computable projection operation on its tangent cone as well as equality constraints. As a supplement of the analysis in \cite{niederlander2016distributed}, we show that the projected dynamical system converges to one of the saddle points and hence finding an optimal solution. Moreover, the problem of distributedly maximizing the algebraic connectivity of an undirected network by optimizing the port gains of each nodes (base stations) is considered. The original semi-definite programming (SDP) problem is relaxed into a nonlinear programming (NP) problem that will be solved by the aforementioned projected dynamical system. Numerical examples show the convergence of the aforementioned algorithm to one of the optimal solutions. The effect of the relaxation is illustrated empirically with numerical examples. A methodology is presented so that the number of iterations needed to reach the equilibrium is suppressed. Complexity per iteration of the algorithm is illustrated with numerical examples.


title: The Smirnov class for spaces with the complete Pick property

abstract: We show that every function in a reproducing kernel Hilbert space with a normalized complete Pick kernel is the quotient of a multiplier and a cyclic multiplier. This extends a theorem of Alpay, Bolotnikov and Kaptano\u{g}lu. We explore various consequences of this result regarding zero sets, spaces on compact sets and Gleason parts. In particular, using a construction of Salas, we exhibit a rotationally invariant complete Pick space of analytic functions on the unit disc for which the corona theorem fails.


title: Group Testing using left-and-right-regular sparse-graph codes

abstract: We consider the problem of non-adaptive group testing of $N$ items out of which $K$ or less items are known to be defective. We propose a testing scheme based on left-and-right-regular sparse-graph codes and a simple iterative decoder. We show that for any arbitrarily small $\epsilon>0$ our scheme requires only $m=c_\epsilon K\log \frac{c_1N}{K}$ tests to recover $(1-\epsilon)$ fraction of the defective items with high probability (w.h.p) i.e., with probability approaching $1$ asymptotically in $N$ and $K$, where the value of constants $c_\epsilon$ and $\ell$ are a function of the desired error floor $\epsilon$ and constant $c_1=\frac{\ell}{c_\epsilon}$ (observed to be approximately equal to 1 for various values of $\epsilon$). More importantly the iterative decoding algorithm has a sub-linear computational complexity of $\mathcal{O}(K\log \frac{N}{K})$ which is known to be optimal. Also for $m=c_2 K\log K\log \frac{N}{K}$ tests our scheme recovers the \textit{whole} set of defective items w.h.p. These results are valid for both noiseless and noisy versions of the problem as long as the number of defective items scale sub-linearly with the total number of items, i.e., $K=o(N)$. The simulation results validate the theoretical results by showing a substantial improvement in the number of tests required when compared to the testing scheme based on left-regular sparse-graphs.


title: Third Law of Thermodynamics as a Single Inequality

abstract: The third law of thermodynamics in the form of the unattainability principle states that exact ground-state cooling requires infinite resources. Here we investigate the amount of non-equilibrium resources needed for approximate cooling. We consider as resource any system out of equilibrium, allowing for resources beyond the i.i.d. assumption and including the input of work as a particular case. We establish in full generality a sufficient and a necessary condition for cooling and show that for a vast class of non-equilibrium resources these two conditions coincide, providing a single necessary and sufficient criterion. Such conditions are expressed in terms of a single function playing a similar role for the third law to the one of the free energy for the second law. From a technical point of view we provide new results about concavity/convexity of certain Renyi-divergences, which might be of independent interest.


title: Epidemiological modeling of the 2005 French riots: a spreading wave and the role of contagion

abstract: As a large-scale instance of dramatic collective behaviour, the 2005 French riots started in a poor suburb of Paris, then spread in all of France, lasting about three weeks. Remarkably, although there were no displacements of rioters, the riot activity did travel. Access to daily national police data has allowed us to explore the dynamics of riot propagation. Here we show that an epidemic-like model, with just a few parameters and a single sociological variable characterizing neighbourhood deprivation, accounts quantitatively for the full spatio-temporal dynamics of the riots. This is the first time that such data-driven modelling involving contagion both within and between cities (through geographic proximity or media) at the scale of a country, and on a daily basis, is performed. Moreover, we give a precise mathematical characterization to the expression "wave of riots", and provide a visualization of the propagation around Paris, exhibiting the wave in a way not described before. The remarkable agreement between model and data demonstrates that geographic proximity played a major role in the propagation, even though information was readily available everywhere through media. Finally, we argue that our approach gives a general framework for the modelling of the dynamics of spontaneous collective uprisings.


title: Numerical Approximations for the Cahn-Hilliard phase field model of the binary fluid-surfactant system

abstract: In this paper, we consider the numerical approximations for the commonly used binary fluid-surfactant phase field model that consists two nonlinearly coupled Cahn-Hilliard equations. The main challenge in solving the system numerically is how to develop easy-to-implement time stepping schemes while preserving the unconditional energy stability. We solve this issue by developing two linear and decoupled, first order and a second order time-stepping schemes using the so-called "Invariant Energy Quadratization" approach for the double well potentials and a subtle explicit-implicit technique for the nonlinear coupling potential. Moreover, the resulting linear system is well-posed and the linear operator is symmetric positive definite. We rigorously prove the first order scheme is unconditionally energy stable. Various numerical simulations are presented to demonstrate the stability and the accuracy thereafter.


title: Learning Word-Like Units from Joint Audio-Visual Analysis

abstract: Given a collection of images and spoken audio captions, we present a method for discovering word-like acoustic units in the continuous speech signal and grounding them to semantically relevant image regions. For example, our model is able to detect spoken instances of the word 'lighthouse' within an utterance and associate them with image regions containing lighthouses. We do not use any form of conventional automatic speech recognition, nor do we use any text transcriptions or conventional linguistic annotations. Our model effectively implements a form of spoken language acquisition, in which the computer learns not only to recognize word categories by sound, but also to enrich the words it learns with semantics by grounding them in images.


title: Switched control for quantized feedback systems: invariance and limit cycle analysis

abstract: We study feedback control for discrete-time linear time-invariant systems in the presence of quantization both in the control action and in the measurement of the controlled variable. While in some application the quantization effects can be neglected, when high-precision control is needed, they have to be explicitly accounted for in control design. In this paper we propose a switched control solution for minimizing the effect of quantization of both the control and controlled variables in the case of a simple integrator with unitary delay, a model that is quite common in the computing systems domain, for example in thread scheduling, clock synchronization, and resource allocation. We show that the switched solution outperforms the one without switching, designed by neglecting quantization, and analyze necessary and sufficient conditions for the controlled system to exhibit periodic solutions in the presence of an additive constant disturbance affecting the control input. Simulation results provide evidence of the effectiveness of the approach.


title: A Model-based Projection Technique for Segmenting Customers

abstract: We consider the problem of segmenting a large population of customers into non-overlapping groups with similar preferences, using diverse preference observations such as purchases, ratings, clicks, etc. over subsets of items. We focus on the setting where the universe of items is large (ranging from thousands to millions) and unstructured (lacking well-defined attributes) and each customer provides observations for only a few items. These data characteristics limit the applicability of existing techniques in marketing and machine learning. To overcome these limitations, we propose a model-based projection technique, which transforms the diverse set of observations into a more comparable scale and deals with missing data by projecting the transformed data onto a low-dimensional space. We then cluster the projected data to obtain the customer segments. Theoretically, we derive precise necessary and sufficient conditions that guarantee asymptotic recovery of the true customer segments. Empirically, we demonstrate the speed and performance of our method in two real-world case studies: (a) 84% improvement in the accuracy of new movie recommendations on the MovieLens data set and (b) 6% improvement in the performance of similar item recommendations algorithm on an offline dataset at eBay. We show that our method outperforms standard latent-class and demographic-based techniques.


title: Monitoring and Intervention: Concepts and Formal Models

abstract: Our machines, products, utilities, and environments have long been monitored by embedded software systems. Our professional, commercial, social and personal lives are also subject to monitoring as they are mediated by software systems. Data on nearly everything now exists, waiting to be collected and analysed for all sorts of reasons. Given the rising tide of data we pose the questions: What is monitoring? Do diverse and disparate monitoring systems have anything in common? We attempt answer these questions by proposing an abstract conceptual framework for studying monitoring. We argue that it captures a structure common to many different monitoring practices, and that from it detailed formal models can be derived, customised to applications. The framework formalises the idea that monitoring is a process that observes the behaviour of people and objects in a context. The entities and their behaviours are represented by abstract data types and the observable attributes by logics. Since monitoring usually has a specific purpose, we extend the framework with protocols for detecting attributes or events that require interventions and, possibly, a change in behaviour. Our theory is illustrated by a case study from criminal justice, that of electronic tagging.


title: Relay-Assisted Mixed FSO/RF Systems over M\'alaga-$\mathcal{M}$ and $\kappa$-$\mu$ Shadowed Fading Channels

abstract: This letter presents a unified analytical framework for relay-assisted mixed FSO/RF transmission. In addition to accounting for different FSO detection techniques, the mathematical model offers a twofold unification of mixed FSO/RF systems by considering mixed M\'alaga-$\mathcal{M}$/$\kappa$-$\mu$ shadowed fading, which includes as special cases nearly all linear turbulence/fading models adopted in the open literature.


title: Efficient computation of multidimensional theta functions

abstract: An important step in the efficient computation of multi-dimensional theta functions is the construction of appropriate symplectic transformations for a given Riemann matrix assuring a rapid convergence of the theta series. An algorithm is presented to approximately map the Riemann matrix to the Siegel fundamental domain. The shortest vector of the lattice generated by the Riemann matrix is identified exactly, and the algorithm ensures that its length is larger than $\sqrt{3}/2$. The approach is based on a previous algorithm by Deconinck et al. using the LLL algorithm for lattice reductions. Here, the LLL algorithm is replaced by exact Minkowski reductions for small genus and an exact identification of the shortest lattice vector for larger values of the genus.


title: A Linear, Decoupled and Energy stable scheme for smectic-A Liquid Crystal Flows

abstract: In this paper, we consider numerical approximations for the model of smectic-A liquid crystal flows. The model equation, that is derived from the variational approach of the de Gennes free energy, is a highly nonlinear system that couples the incompressible Navier-Stokes equations, and two nonlinear coupled second-order elliptic equations. Based on some subtle explicit--implicit treatments for nonlinear terms, we develop a unconditionally energy stable, linear and decoupled time marching numerical scheme. We also rigorously prove that the proposed scheme obeys the energy dissipation law at the discrete level. Various numerical simulations are presented to demonstrate the accuracy and the stability thereafter.


title: Joint Power Allocation and Beamforming for Energy-Efficient Two-Way Multi-Relay Communications

abstract: This paper considers the joint design of user power allocation and relay beamforming in relaying communications, in which multiple pairs of single-antenna users exchange information with each other via multiple-antenna relays in two time slots. All users transmit their signals to the relays in the first time slot while the relays broadcast the beamformed signals to all users in the second time slot. The aim is to maximize the system's energy efficiency (EE) subject to quality-of-service (QoS) constraints in terms of exchange throughput requirements. The QoS constraints are nonconvex with many nonlinear cross-terms, so finding a feasible point is already computationally challenging. The sum throughput appears in the numerator while the total consumption power appears in the denominator of the EE objective function. The former is a nonconcave function and the latter is a nonconvex function, making fractional programming useless for EE optimization. Nevertheless, efficient iterations of low complexity to obtain its optimized solutions are developed. The performances of the multiple-user and multiple-relay networks under various scenarios are evaluated to show the merit of the paper development.


title: Free fermions on a piecewise linear four-manifold. II: Pachner moves

abstract: This is the second in a series of papers where we construct an invariant of a four-dimensional piecewise linear manifold $M$ with a given middle cohomology class $h\in H^2(M,\mathbb C)$. This invariant is the square root of the torsion of unusual chain complex introduced in Part I (arXiv:1605.06498) of our work, multiplied by a correcting factor. Here we find this factor by studying the behavior of our construction under all four-dimensional Pachner moves, and show that it can be represented in a multiplicative form: a product of same-type multipliers over all 2-faces, multiplied by a product of same-type multipliers over all pentachora.


title: What Are People Tweeting about Zika? An Exploratory Study Concerning Symptoms, Treatment, Transmission, and Prevention

abstract: The purpose of this study was to do a dataset distribution analysis, a classification performance analysis, and a topical analysis concerning what people are tweeting about four disease characteristics: symptoms, transmission, prevention, and treatment. A combination of natural language processing and machine learning techniques were used to determine what people are tweeting about Zika. Specifically, a two-stage classifier system was built to find relevant tweets on Zika, and then categorize these into the four disease categories. Tweets in each disease category were then examined using latent dirichlet allocation (LDA) to determine the five main tweet topics for each disease characteristic. Results 1,234,605 tweets were collected. Tweets by males and females were similar (28% and 23% respectively). The classifier performed well on the training and test data for relevancy (F=0.87 and 0.99 respectively) and disease characteristics (F=0.79 and 0.90 respectively). Five topics for each category were found and discussed with a focus on the symptoms category. Through this process, we demonstrate how misinformation can be discovered so that public health officials can respond to the tweets with misinformation.


title: On Lipschitz continuous optimal stopping boundaries

abstract: We obtain a probabilistic proof of the local Lipschitz continuity for the optimal stopping boundary of a class of problems with state space $[0,T]\times\mathbb{R}^d$, $d\ge 1$. To the best of our knowledge this is the only existing proof that relies exclusively upon stochastic calculus, all the other proofs making use of PDE techniques and integral equations. Thanks to our approach we obtain our result for a class of diffusions whose associated second order differential operator is not necessarily uniformly elliptic. The latter condition is normally assumed in the related PDE literature.


title: Path abstraction

abstract: Given the set of paths through a digraph, the result of uniformly deleting some vertices and identifying others along each path is coherent in such a way as to yield the set of paths through another digraph, called a \emph{path abstraction} of the original digraph. The construction of path abstractions is detailed and relevant basic results are established; generalizations are also discussed. Connections with random digraphs are also illustrated.


title: Ballistic, diffusive, and arrested transport in disordered momentum-space lattices

abstract: Ultracold atoms in optical lattices offer a unique platform for investigating disorder-driven phenomena. While static disordered site potentials have been explored in a number of optical lattice experiments, a more general control over site-energy and off-diagonal tunneling disorder has been lacking. The use of atomic quantum states as "synthetic dimensions" has introduced the spectroscopic, site-resolved control necessary to engineer new, more tailored realizations of disorder. Here, by controlling laser-driven dynamics of atomic population in a momentum-space lattice, we extend the range of synthetic-dimension-based quantum simulation and present the first explorations of dynamical disorder and tunneling disorder in an atomic system. By applying static tunneling phase disorder to a one-dimensional lattice, we observe ballistic quantum spreading as in the case of uniform tunneling. When the applied disorder fluctuates on timescales comparable to intersite tunneling, we instead observe diffusive atomic transport, signaling a crossover from quantum to classical expansion dynamics. We compare these observations to the case of static site-energy disorder, where we directly observe quantum localization in the momentum-space lattice.


title: Non-Stoquastic Interactions in Quantum Annealing via the Aharonov-Anandan Phase

abstract: We argue that a complete description of quantum annealing (QA) implemented with continuous variables must take into account the non-adiabatic Aharonov-Anandan geometric phase that arises when the system Hamiltonian changes during the anneal. We show that this geometric effect leads to the appearance of non-stoquastic terms in the effective quantum Ising Hamiltonians that are typically used to describe QA with flux-qubits. We explicitly demonstrate the effect of these geometric interactions when QA is performed with a system of one and two coupled flux-qubits. The realization of non-stoquastic Hamiltonians has important implications from a computational complexity perspective, since it is believed that in many cases QA with stoquastic Hamiltonians can be efficiently simulated via classical algorithms such as Quantum Monte Carlo. It is well-known that the direct implementation of non-stoquastic interactions with flux-qubits is particularly challenging. Our results suggest an alternative path for the implementation of non-stoquastic interactions via geometric phases that can be exploited for computational purposes.


title: Two-Party Function Computation on the Reconciled Data

abstract: In this paper, we initiate a study of a new problem termed function computation on the reconciled data, which generalizes a set reconciliation problem in the literature. Assume a distributed data storage system with two users $A$ and $B$. The users possess a collection of binary vectors $S_{A}$ and $S_{B}$, respectively. They are interested in computing a function $\phi$ of the reconciled data $S_{A} \cup S_{B}$. It is shown that any deterministic protocol, which computes a sum and a product of reconciled sets of binary vectors represented as nonnegative integers, has to communicate at least $2^n + n - 1$ and $2^n + n - 2$ bits in the worst-case scenario, respectively, where $n$ is the length of the binary vectors. Connections to other problems in computer science, such as set disjointness and finding the intersection, are established, yielding a variety of additional upper and lower bounds on the communication complexity. A protocol for computation of a sum function, which is based on use of a family of hash functions, is presented, and its characteristics are analyzed.


title: Phylogenetic Factor Analysis

abstract: Phylogenetic comparative methods explore the relationships between quantitative traits adjusting for shared evolutionary history. This adjustment often occurs through a Brownian diffusion process along the branches of the phylogeny that generates model residuals or the traits themselves. For high-dimensional traits, inferring all pair-wise correlations within the multivariate diffusion is limiting. To circumvent this problem, we propose phylogenetic factor analysis (PFA) that assumes a small unknown number of independent evolutionary factors arise along the phylogeny and these factors generate clusters of dependent traits. Set in a Bayesian framework, PFA provides measures of uncertainty on the factor number and groupings, combines both continuous and discrete traits, integrates over missing measurements and incorporates phylogenetic uncertainty with the help of molecular sequences. We develop Gibbs samplers based on dynamic programming to estimate the PFA posterior distribution, over three-fold faster than for multivariate diffusion and a further order-of-magnitude more efficiently in the presence of latent traits. We further propose a novel marginal likelihood estimator for previously impractical models with discrete data and find that PFA also provides a better fit than multivariate diffusion in evolutionary questions in columbine flower development, placental reproduction transitions and triggerfish fin morphometry.


title: A Type B Analogue to Ribbon Tableaux

abstract: We introduce a shifted analogue of the ribbon tableaux defined by James and Kerber. For any positive integer $k$, we give a bijection between the $k$-ribbon fillings of a shifted shape and regular fillings of a $\lfloor k/2\rfloor$-tuple of shapes called its $k$-quotient. We also define the corresponding generating functions, and prove that they are symmetric, Schur positive and Schur $Q$-positive.


title: On rescheduling due to machine disruption while to minimize the total weighted completion time

abstract: We investigate a single machine rescheduling problem that arises from an unexpected machine unavailability, after the given set of jobs has already been scheduled to minimize the total weighted completion time. Such a disruption is represented as an unavailable time interval and is revealed to the production planner before any job is processed; the production planner wishes to reschedule the jobs to minimize the alteration to the originally planned schedule, which is measured as the maximum time deviation between the original and the new schedules for all the jobs. The objective function in this rescheduling problem is to minimize the sum of the total weighted completion time and the weighted maximum time deviation, under the constraint that the maximum time deviation is bounded above by a given value. That is, the maximum time deviation is taken both as a constraint and as part of the objective function. We present a pseudo-polynomial time exact algorithm and a fully polynomial time approximation scheme, the latter of which is the best possible given that the general problem is NP-hard.


title: Geometric Stability of Pseudo-plane Ideal Flows

abstract: Geometric analysis of steady pseudo-plane ideal flow reveals a fundamental relation between vertical coherence and streamline topology. It shows vertical alignment only exists in straightline jet and circular vortex. A geometric stability theory is then developed as a structural stability in contrast to traditional dynamical stability, suggesting all vertical-aligned pseudo-plane flows are geometrically stable and all non-aligned non-straightline flows are unstable. It provides a topological explanation for the tendency of axisymmetrization and vertical alignment observed in vortex experiments.


title: Scalable Architecture for Anomaly Detection and Visualization in Power Generating Assets

abstract: Power-generating assets (e.g., jet engines, gas turbines) are often instrumented with tens to hundreds of sensors for monitoring physical and performance degradation. Anomaly detection algorithms highlight deviations from predetermined benchmarks with the goal of detecting incipient faults. We are developing an integrated system to address three key challenges within analyzing sensor data from power-generating assets: (1) difficulty in ingesting and analyzing data from large numbers of machines; (2) prevalence of false alarms generated by anomaly detection algorithms resulting in unnecessary downtime and maintenance; and (3) lack of an integrated visualization that helps users understand and explore the flagged anomalies and relevant sensor context in the energy domain. We present preliminary results and our key findings in addressing these challenges. Our system's scalable event ingestion framework, based on OpenTSDB, ingests nearly 400,000 sensor data samples per seconds using a 30 machine cluster. To reduce false alarm rates, we leverage the False Discovery Rate (FDR) algorithm which significantly reduces the number of false alarms. Our visualization tool presents the anomalies and associated content flagged by the FDR algorithm to inform users and practitioners in their decision making process. We believe our integrated platform will help reduce maintenance costs significantly while increasing asset lifespan. We are working to extend our system on multiple fronts, such as scaling to more data and more compute nodes (70 in total).


title: Locality and Availability of Array Codes Constructed from Subspaces

abstract: Ever-increasing amounts of data are created and processed in internet-scale companies such as Google, Facebook, and Amazon. The efficient storage of such copious amounts of data has thus become a fundamental and acute problem in modern computing. No single machine can possibly satisfy such immense storage demands. Therefore, distributed storage systems (DSS), which rely on tens of thousands of storage nodes, are the only viable solution. Such systems are broadly used in all modern internet-scale systems. However, the design of a DSS poses a number of crucial challenges, markedly different from single-user storage systems. Such systems must be able to reconstruct the data efficiently, to overcome failure of servers, to correct errors, etc. Lots of research was done in the last few years to answer these challenges and the research is increasing in parallel to the increasing amount of stored data. The main goal of this paper is to consider codes which have two of the most important features of distributed storage systems, namely, locality and availability. Our codes are array codes which are based on subspaces of a linear space over a finite field. We present several constructions of such codes which are $q$-analog to some of the known block codes. Some of these codes possess independent intellectual merit. We examine the locality and availability of the constructed codes. In particular we distinguish between two types of locality and availability, node vs.~symbol, locality and availability. To our knowledge this is the first time that such a distinction is given in the literature.


title: Global aspects of polarization optics and coset space geometry

abstract: We use group theoretic ideas and coset space methods to deal with problems in polarization optics of a global nature. These include the possibility of a globally smooth phase convention for electric fields for all points on the Poincar\'{e} sphere, and a similar possibility of real or complex bases of transverse electric vectors for all possible propagation directions. It is shown that these methods help in understanding some known results in an effective manner, and in answering new questions as well. We find that apart from the groups $SU(2)$ and $SO(3)$ which occur naturally in these problems, the group $SU(3)$ also plays an important role.


title: Angular distribution of single photon superradiance in a dilute and cold atomic ensemble

abstract: On the basis of a quantum microscopic approach we study the dynamics of the afterglow of a dilute Gaussian atomic ensemble excited by pulsed radiation. Taking into account the vector nature of the electromagnetic field we analyze in detail the angular and polarization distribution of single-photon superradiance of such an ensemble. The dependence of the angular distribution of superradiance on the length of the pulse and its carrier frequency as well as on the size and the shape of the atomic clouds is studied. We show that there is substantial dependence of the superradiant emission on the polarization and the direction of fluorescence. We observe essential peculiarities of superradiance in the region of the forward diffraction zone and in the area of the coherent backscattering cone. We demonstrate that there are directions for which the rate of fluorescence is several times more than the decay rate of the timed-Dicke state. We show also that single-photon superradiance can be excited by incoherent excitation when atomic polarization in the ensemble is absent. Besides a quantum microscopic approach, we analyze single-photon superradiance on the basis of the theory of incoherent multiple scattering in optically thick media (random walk theory). In the case of very short resonant and long nonresonant pulses we derive simple analytical expressions for the decay rate of single-photon superradiance for incoherent fluorescence in an arbitrary direction.


title: Nuclear fusion enhances cancer cell killing efficacy in a protontherapy model

abstract: Protontherapy is hadrontherapy fastest-growing modality and a pillar in the battle against cancer. Hadrontherapy superiority lies in its inverted depth-dose profile, hence tumour-confined irradiation. Protons, however, lack distinct radiobiological advantages over photons or electrons. Higher LET (Linear Energy Transfer) $^{12}$C ions can overcome cancer radioresistance: DNA lesion complexity increases with LET, resulting in efficient cell killing, i.e. higher Relative Biological Effectiveness (RBE). However, economic and radiobiological issues hamper $^{12}$C-ion clinical amenability. Thus, enhancing proton RBE is desirable. To this end, we exploited the p + $^{11}$B $\rightarrow$3$\alpha$ reaction to generate high-LET alpha particles with a clinical proton beam. To maximize the reaction rate, we used sodium borocaptate (BSH) with natural boron content. Boron-Neutron Capture Therapy (BNCT) uses $^{10}$B-enriched BSH for neutron irradiation-triggered alpha-particles. We recorded significantly increased cellular lethality and chromosome aberration complexity. A strategy combining protontherapy ballistic precision with the higher RBE promised by BNCT and $^{12}$C-ion therapy is thus demonstrated.


title: Extracting single crystal diffraction pattern from powder diffraction by intensity correlation functions

abstract: We applied the analysis of x-ray intensity angular correlation function to dilute ensembles of identical spinel crystals. Firstly, we show that the angular correlation from measured diffraction patterns with many crystals per shot converges to the correlation for single crystal. Secondly, we determined the diffraction pattern from single crystal from diffraction patterns of many crystals. Finally, we discussed the full 3D reciprocal space retrieval and its potential application in solving structure with dilute powder diffraction data, where the intensity fluctuation could be observed on each ring.


title: Bayesian Hierarchical Models with Conjugate Full-Conditional Distributions for Dependent Data from the Natural Exponential Family

abstract: We introduce a Bayesian approach for analyzing (possibly) high-dimensional dependent data that are distributed according to a member from the natural exponential family of distributions. This problem requires extensive methodological advancements, as jointly modeling high-dimensional dependent data leads to the so-called "big n problem." The computational complexity of the "big n problem" is further exacerbated when allowing for non-Gaussian data models, as is the case here. Thus, we develop new computationally efficient distribution theory for this setting. In particular, we introduce the "conjugate multivariate distribution," which is motivated by the univariate distribution introduced in Diaconis and Ylvisaker (1979). Furthermore, we provide substantial theoretical and methodological development including: results regarding conditional distributions, an asymptotic relationship with the multivariate normal distribution, conjugate prior distributions, and full-conditional distributions for a Gibbs sampler. To demonstrate the wide-applicability of the proposed methodology, we provide two simulation studies and three applications based on an epidemiology dataset, a federal statistics dataset, and an environmental dataset, respectively.


title: The Fifth Moment of modular L-functions

abstract: We prove a sharp bound on the fifth moment of modular L-functions of fixed small weight, and large prime level.


title: Some estimates for $\theta$-type Calder\'on-Zygmund operators and linear commutators on certain weighted amalgam spaces

abstract: In this paper, we first introduce some new kinds of weighted amalgam spaces. Then we discuss the strong type and weak type estimates for a class of Calder\'on--Zygmund type operators $T_\theta$ in these new weighted spaces. Furthermore, the strong type estimate and endpoint estimate of linear commutators $[b,T_{\theta}]$ formed by $b$ and $T_{\theta}$ are established. Also we study related problems about two-weight, weak type inequalities for $T_{\theta}$ and $[b,T_{\theta}]$ in the weighted amalgam spaces and give some results.


title: Thin films with precisely engineered nanostructures

abstract: Synthesis of rationally designed nanostructured materials with optimized mechanical properties, e.g., high strength with considerable ductility, requires rigorous control of diverse microstructural parameters including the mean size, size dispersion and spatial distribution of grains. However, currently available synthesis techniques can seldom satisfy these requirements. Here, we report a new methodology to synthesize thin films with unprecedented microstructural control via systematic, in-situ seeding of nanocrystals into amorphous precursor films. When the amorphous films are subsequently crystallized by thermal annealing, the nanocrystals serve as preferential grain nucleation sites and control their microstructure. We demonstrate the capability of this approach by precisely tailoring the size, geometry and spatial distribution of nanostructured grains in structural (TiAl) as well as functional (TiNi) thin films. The approach, which is applicable to a broad class of metallic alloys and ceramics, enables explicit microstructural control of thin film materials for a wide spectrum of applications.


title: Nitrogen Fractionation in Protoplanetary Disks from the H13CN/HC15N Ratio

abstract: Nitrogen fractionation is commonly used to assess the thermal history of Solar System volatiles. With ALMA it is for the first time possible to directly measure 14N/15N ratios in common molecules during the assembly of planetary systems. We present ALMA observations of the H13CN and HC15N J=3-2 lines at 0".5 angular resolution, toward a sample of six protoplanetary disks, selected to span a range of stellar and disk structure properties. Adopting a typical 12C/13C ratio of 70, we find comet-like 14N/15N ratios of 80-160 in 5/6 of the disks (3 T Tauri and 2 Herbig Ae disks) and lack constraints for one of the T Tauri disks (IM Lup). There are no systematic differences between T Tauri and Herbig Ae disks, or between full and transition disks within the sample. In addition, no correlation is observed between disk-averaged D/H and 14N/15N ratios in the sample. One of the disks, V4046 Sgr, presents unusually bright HCN isotopologue emission, enabling us to model the radial profiles of H13CN and HC15N. We find tentative evidence of an increasing 14N/15N ratio with radius, indicating that selective photodissociation in the inner disk is important in setting the 14N/15N ratio during planet formation.


title: Mini-BFSS in Silico

abstract: We study a mass-deformed $\mathcal{N}=4$ version of the BFSS matrix model with three matrices and gauge group $SU(2)$. This model has zero Witten index. Despite this, we give numerical evidence for the existence of four supersymmetric ground states, two bosonic and two fermionic, in the limit where the mass deformation is tuned to zero.


title: ACIA, not ACID: Conditions, Properties and Challenges

abstract: Although ACID is the previous golden rule for transaction support, durability is now not a basic requirement for data storage. Rather, high availability is becoming the first-class property required by online applications. We show that high availability of data is almost surely a stronger property than durability. We thus propose ACIA (Atomicity, Consistency, Isolation, Availability) as the new standard for transaction support. Essentially, the shift from ACID to ACIA is due to the change of assumed conditions for data management. Four major condition changes exist. With ACIA transactions, more diverse requirements can be flexibly supported for applications through the specification of consistency levels, isolation levels and fault tolerance levels. Clarifying the ACIA properties enables the exploitation of techniques used for ACID transactions, as well as bringing about new challenges for research.


title: Clouds in the atmospheres of extrasolar planets. V. The impact of CO2 ice clouds on the outer boundary of the habitable zone

abstract: Clouds have a strong impact on the climate of planetary atmospheres. The potential scattering greenhouse effect of CO2 ice clouds in the atmospheres of terrestrial extrasolar planets is of particular interest because it might influence the position and thus the extension of the outer boundary of the classic habitable zone around main sequence stars. Here, the impact of CO2 ice clouds on the surface temperatures of terrestrial planets with CO2 dominated atmospheres, orbiting different types of stars is studied. Additionally, their corresponding effect on the position of the outer habitable zone boundary is evaluated. For this study, a radiative-convective atmospheric model is used the calculate the surface temperatures influenced by CO2 ice particles. The clouds are included using a parametrised cloud model. The atmospheric model includes a general discrete ordinate radiative transfer that can describe the anisotropic scattering by the cloud particles accurately. A net scattering greenhouse effect caused by CO2 clouds is only obtained in a rather limited parameter range which also strongly depends on the stellar effective temperature. For cool M-stars, CO2 clouds only provide about 6 K of additional greenhouse heating in the best case scenario. On the other hand, the surface temperature for a planet around an F-type star can be increased by 30 K if carbon dioxide clouds are present. Accordingly, the extension of the habitable zone due to clouds is quite small for late-type stars. Higher stellar effective temperatures, on the other hand, can lead to outer HZ boundaries about 0.5 au farther out than the corresponding clear-sky values.


title: On the existence of connecting orbits for critical values of the energy

abstract: We consider an open connected set $\Omega$ and a smooth potential $U$ which is positive in $\Omega$ and vanishes on $\partial\Omega$. We study the existence of orbits of the mechanical system \[ \ddot{u}=U_x(u), \] that connect different components of $\partial\Omega$ and lie on the zero level of the energy. We allow that $\partial\Omega$ contains a finite number of critical points of $U$. The case of symmetric potential is also considered.


title: Q-analogues of the Fibo-Stirling numbers

abstract: Let $F_n$ denote the $n^{th}$ Fibonacci number relative to the initial conditions $F_0=0$ and $F_1=1$. Bach, Paudyal, and Remmel introduced Fibonacci analogues of the Stirling numbers called Fibo-Stirling numbers of the first and second kind. These numbers serve as the connection coefficients between the Fibo-falling factorial basis $\{(x)_{\downarrow_{F,n}}:n \geq 0\}$ and the Fibo-rising factorial basis $\{(x)_{\uparrow_{F,n}}:n \geq 0\}$ which are defined by $(x)_{\downarrow_{F,0}} = (x)_{\uparrow_{F,0}} = 1$ and for $k \geq 1$, $(x)_{\downarrow_{F,k}} = x(x-F_1) \cdots (x-F_{k-1})$ and $(x)_{\uparrow_{F,k}} = x(x+F_1) \cdots (x+F_{k-1})$. We gave a general rook theory model which allowed us to give combinatorial interpretations of the Fibo-Stirling numbers of the first and second kind. There are two natural $q$-analogues of the falling and rising Fibo-factorial basis. That is, let $[x]_q = \frac{q^x-1}{q-1}$. Then we let $[x]_{\downarrow_{q,F,0}} = \overline{[x]}_{\downarrow_{q,F,0}} = [x]_{\uparrow_{q,F,0}} = \overline{[x]}_{\uparrow_{q,F,0}}=1$ and, for $k > 0$, we let $[x]_{\downarrow_{q,F,k}} = [x]_q [x-F_1]_q \cdots [x-F_{k-1}]_q$, $\overline{[x]}_{\downarrow_{q,F,k}}= [x]_q ([x]_q-[F_1]_q) \cdots ([x]_q-[F_{k-1}]_q)$, $[x]_{\uparrow_{q,F,k}}= [x]_q [x+F_1]_q \cdots [x+F_{k-1}]_q$, and $\overline{[x]}_{\uparrow_{q,F,k}}= [x]_q ([x]_q+[F_1]_q) \cdots ([x]_q+[F_{k-1}]_q)$. In this paper, we show we can modify the rook theory model of Bach, Paudyal, and Remmel to give combinatorial interpretations for the two different types $q$-analogues of the Fibo-Stirling numbers which arise as the connection coefficients between the two different $q$-analogues of the Fibonacci falling and rising factorial bases. \end{abstract}


title: Non-colocated Time-Reversal MUSIC: High-SNR Distribution of Null Spectrum

abstract: We derive the asymptotic distribution of the null spectrum of the well-known Multiple Signal Classification (MUSIC) in its computational Time-Reversal (TR) form. The result pertains to a single-frequency non-colocated multistatic scenario and several TR-MUSIC variants are here investigated. The analysis builds upon the 1st-order perturbation of the singular value decomposition and allows a simple characterization of null-spectrum moments (up to the 2nd order). This enables a comparison in terms of spectrums stability. Finally, a numerical analysis is provided to confirm the theoretical findings.


title: Hardware Translation Coherence for Virtualized Systems

abstract: To improve system performance, modern operating systems (OSes) often undertake activities that require modification of virtual-to-physical page translation mappings. For example, the OS may migrate data between physical frames to defragment memory and enable superpages. The OS may migrate pages of data between heterogeneous memory devices. We refer to all such activities as page remappings. Unfortunately, page remappings are expensive. We show that translation coherence is a major culprit and that systems employing virtualization are especially badly affected by their overheads. In response, we propose hardware translation invalidation and coherence or HATRIC, a readily implementable hardware mechanism to piggyback translation coherence atop existing cache coherence protocols. We perform detailed studies using KVM-based virtualization, showing that HATRIC achieves up to 30% performance and 10% energy benefits, for per-CPU area overheads of 2%. We also quantify HATRIC's benefits on systems running Xen and find up to 33% performance improvements.


title: On The Compound MIMO Wiretap Channel with Mean Feedback

abstract: Compound MIMO wiretap channel with double sided uncertainty is considered under channel mean information model. In mean information model, channel variations are centered around its mean value which is fed back to the transmitter. We show that the worst case main channel is anti-parallel to the channel mean information resulting in an overall unit rank channel. Further, the worst eavesdropper channel is shown to be isotropic around its mean information. Accordingly, we provide the capacity achieving beamforming direction. We show that the saddle point property holds under mean information model, and thus, compound secrecy capacity equals to the worst case capacity over the class of uncertainty. Moreover, capacity achieving beamforming direction is found to require matrix inversion, thus, we derive the null steering (NS) beamforming as an alternative suboptimal solution that does not require matrix inversion. NS beamformer is in the direction orthogonal to the eavesdropper mean channel that maintains the maximum possible gain in mean main channel direction. Extensive computer simulation reveals that NS performs very close to the optimal solution. It also verifies that, NS beamforming outperforms both maximum ratio transmission (MRT) and zero forcing (ZF) beamforming approaches over the entire SNR range. Finally, An equivalence relation with MIMO wiretap channel in Rician fading environment is established.


title: SO*(2N) coherent states for loop quantum gravity

abstract: A SU(2) intertwiner with N legs can be interpreted as the quantum state of a convex polyhedron with N faces (when working in 3d). We show that the intertwiner Hilbert space carries a representation of the non-compact group SO*(2N). This group can be viewed as the subgroup of the symplectic group Sp(4N,R) which preserves the SU(2) invariance. We construct the associated Perelomov coherent states and discuss the notion of semi-classical limit, which is more subtle that we could expect. Our work completes the work by Freidel and Livine which focused on the U(N) subgroup of SO*(2N).


title: Reaching for the quantum limits in the simultaneous estimation of phase and phase diffusion

abstract: Phase diffusion invariably accompanies all phase estimation strategies -- quantum or classical. A precise estimation of the former can often provide valuable understanding of the physics of the phase generating phenomena itself. We theoretically examine the performance of fixed-particle number probe states in the simultaneous estimation of phase and collective phase diffusion. We derive analytical quantum limits associated with the simultaneous local estimation of phase and phase diffusion within the quantum Cramer-Rao bound framework in the regimes of large and small phase diffusive noise. The former is for a general fixed-particle number state and the latter for Holland Burnett states, for which we show quantum-enhanced estimation of phase as well as phase diffusion. We next investigate the simultaneous attainability of these quantum limits using projective measurements acting on a single copy of the state in terms of a trade-off relation. In particular, we are interested how this trade-off varies as a function of the dimension of the state. We derive an analytical bound for this trade-off in the large phase diffusion regime for a particular form of the measurement, and show that the maximum of 2, set by the quantum Cramer-Rao bound, is attainable. Further, we show numerical evidence that as diffusion approaches zero, the optimal trade-off relation approaches 1 for Holland-Burnett states. These numerical results are valid in the small particle number regime and suggest that the trade-off for estimating one parameter with quantum-limited precision leads to a complete lack of precision for the other parameter as the diffusion strength approaches zero. Finally, we provide numerical results showing behaviour of the trade-off for a general value of phase diffusion when using Holland-Burnett probe states.


title: Floor Scale Modulo Lifting for QC-LDPC codes

abstract: In the given paper we present a novel approach for constructing a QC-LDPC code of smaller length by lifting a given QC-LDPC code. The proposed method can be considered as a generalization of floor lifting. Also we prove several probabilistic statements concerning a theoretical improvement of the method with respect to the number of small cycles. Making some offline calculation of scale parameter it is possible to construct a sequence of QC-LDPC codes with different circulant sizes generated from a single exponent matrix using only floor and scale operations. The only parameter we store in memory is a constant needed for scaling.


title: Joint Uplink-Downlink Cell Associations for Interference Networks with Local Connectivity

abstract: We study information theoretic models of interference networks that consist of K Base Station (BS) - Mobile Terminal (MT) pairs. Each BS is connected to the MT carrying the same index as well as L following MTs, where the connectivity parameter L >= 1. We fix the value of L and study large networks as K goes to infinity. We assume that each MT can be associated with Nc BSs, and these associations are determined by a cloud-based controller that has a global view of the network. An MT has to be associated with a BS, in order for the BS to transmit its message in the downlink, or decode its message in the uplink. In previous work, the cell associations that maximize the average uplink-downlink per user degrees of freedom (puDoF) were identified for the case when L=1. Further, when only the downlink is considered, the problem was settled for all values of L when we are restricted to use only zero-forcing interference cancellation schemes. In this work, we first propose puDoF inner bounds for arbitrary values of L when only the uplink is considered, and characterize the uplink puDoF value when only zero-forcing schemes are allowed. We then introduce new achievable average uplink-downlink puDoF values. We show that the new scheme is optimal for the range when Nc <= L/2 when we restrict our attention to zero forcing schemes. Additionally we conjecture that the having unity puDoF during uplink is optimal when Nc >= L.


title: Searching for Large-gap Quantum Spin Hall Insulators: Boron-Nitride/(Pb, Sn)/{\alpha}-Al2O3 Sandwich Structures

abstract: Topological insulators hold great potential for efficient information processing and storage. Using density functional theory calculations, we predict that a honeycomb lead monolayer can be stabilized on the Al2O3 (0001) substrate and becomes topologically non-trivial with a sizeable band gap (~0.27 eV). Furthermore, we propose to use hexagonal boron-nitride (h-BN) monolayer as a protection for the topological states of Pb/Al2O3 and Sn/Al2O3. Our findings suggest new possibilities for designing and protecting two-dimensional TIs for practical applications.


title: The Role of Transmitter Cooperation in Linear Interference Networks with Block Erasures

abstract: In this work, we explore the potential and optimal use of transmitter cooperation in wireless interference networks with deep fading conditions. We consider a linear interference network with K transmitter-receiver pairs, where each transmitter can be connected to two neighboring receivers. Long-term fluctuations (shadow fading) in the wireless channel can lead to any link being erased with probability p. Each receiver is interested in one unique message that can be available at two transmitters. The considered rate criterion is the average per user degrees of freedom (puDoF) as K goes to infinity. Prior to this work, the optimal assignment of messages to transmitters were identified in the two limits as p goes to 0 and as p goes to 1. We identify new schemes that achieve average puDoF values that are higher than the state of the art for a significant part of the range 0 < p < 1. The key idea to our results is to understand that the role of cooperation shifts from increasing the probability of delivering a message to its intended destination at high values of p, to interference cancellation at low values of p. Our schemes are based on an algorithm that achieves the optimal DoF value in any network realization, when restricted to a given message assignment as well as the use of zero-forcing schemes.


title: Khovanov complexes of rational tangles

abstract: We show that the Khovanov complex of a rational tangle has a very simple representative whose backbone of non-zero morphisms forms a zig-zag. Furthermore, this minimal complex can be computed quickly by an inductive algorithm. (For example, we calculate $Kh(8_2)$ by hand.) We find that the bigradings of the subobjects in these minimal complexes can be described by matrix actions, which after a change of basis is the reduced Burau representation of $B_3$.


title: Possibility of Realizing Quantum Spin Hall Effect at Room Temperature in Stanene/Al2O3(0001)

abstract: Two-dimensional quantum spin Hall (QSH) insulators with reasonably wide band gaps are imperative for the development of various innovative technologies. Through systematic density functional calculations and tight-binding simulations, we found that stanene on {\alpha}-alumina surface may possess a sizeable topologically nontrivial band gap (~0.25 eV) at the {\Gamma} point. Furthermore, stanene is atomically bonded to but electronically decoupled from the substrate, providing high structural stability and isolated QSH states to a large extent. The underlying physical mechanism is rather general, and this finding may lead to the opening of a new vista for the exploration of QSH insulators for room temperature device applications.


title: Formation and assembly history of stellar components in galaxies as a function of stellar and halo mass

abstract: Galaxy mass assembly is an end product of structure formation in the $\Lambda$CDM cosmology. As an extension of Lee \& Yi (2013), we investigate the assembly history of stellar components in galaxies as a function of halo environments and stellar mass using semi-analytic approaches. In our fiducial model, halo mass intrinsically determines the formation and assembly of the stellar mass. Overall, the ex situ fraction slowly increases in central galaxies with increasing halo mass but sharply increases for $\log M_{*}/M_{\odot}\gtrsim11$. A similar trend is also found in satellite galaxies, which implies that mergers are essential to build stellar masses above $\log M_{*}/M_{\odot}\sim11$. We also examine the time evolution of the contribution of mass growth channels. Mergers become the primary channel in the mass growth of central galaxies when their host halo mass begins to exceed $\log M_{200}/M_{\odot}\sim13$. However, satellite galaxies seldom reach the merger-dominant phase despite their reduced star formation activities due to environmental effects.


title: Visibility of 4-covers of elliptic curves

abstract: Let $C$ be a $4$-cover of an elliptic curve $E$, written as a quadric intersection in $\mathbb{P}^3$. Let $E'$ be another elliptic curve with $4$-torsion isomorphic to that of $E$. We show how to write down the $4$-cover $C'$ of $E'$ with the property that $C$ and $C'$ are represented by the same cohomology class on the $4$-torsion. In fact we give equations for $C'$ as a curve of degree $8$ in $\mathbb{P}^5$. We also study the K3-surfaces fibred by the curves $C'$ as we vary $E'$. In particular we show how to write down models for these surfaces as complete intersections of quadrics in $\mathbb{P}^5$ with exactly $16$ singular points. This allows us to give examples of elliptic curves over $\mathbb{Q}$ that have elements of order $4$ in their Tate-Shafarevich group that are not visible in a principally polarized abelian surface.


title: Transport reversal for model reduction of hyperbolic partial differential equations

abstract: Snapshot matrices built from solutions to hyperbolic partial differential equations exhibit slow decay in singular values, whereas fast decay is crucial for the success of projection- based model reduction methods. To overcome this problem, we build on previous work in symmetry reduction [Rowley and Marsden, Physica D (2000), pp. 1-19] and propose an iterative algorithm that decomposes the snapshot matrix into multiple shifting profiles, each with a corresponding speed. Its applicability to typical hyperbolic problems is demonstrated through numerical examples, and other natural extensions that modify the shift operator are considered. Finally, we give a geometric interpretation of the algorithm.


title: Actively Calibrated Line Mountable Capacitive Voltage Transducer For Power Systems Applications

abstract: A class of Actively Calibrated Line Mounted Capacitive Voltage Transducers (LMCVT) are introduced as a viable line mountable instrumentation option for deploying large numbers of voltage transducers onto the medium and high voltage systems. Active Calibration is shown to reduce the error of line mounted voltage measurements by an order of magnitude from previously published techniques. The instrument physics and sensing method is presented and the performance is evaluated in a laboratory setting. Finally, a roadmap to a fully deployable prototype is shown.


title: Design of Improved Quasi-Cyclic Protograph-Based Raptor-Like LDPC Codes for Short Block-Lengths

abstract: Protograph-based Raptor-like low-density parity-check codes (PBRL codes) are a recently proposed family of easily encodable and decodable rate-compatible LDPC (RC-LDPC) codes. These codes have an excellent iterative decoding threshold and performance across all design rates. PBRL codes designed thus far, for both long and short block-lengths, have been based on optimizing the iterative decoding threshold of the protograph of the RC code family at various design rates. In this work, we propose a design method to obtain better quasi-cyclic (QC) RC-LDPC codes with PBRL structure for short block-lengths (of a few hundred bits). We achieve this by maximizing an upper bound on the minimum distance of any QC-LDPC code that can be obtained from the protograph of a PBRL ensemble. The obtained codes outperform the original PBRL codes at short block-lengths by significantly improving the error floor behavior at all design rates. Furthermore, we identify a reduction in complexity of the design procedure, facilitated by the general structure of a PBRL ensemble.


title: Imaging electron flow and quantum dot formation in MoS2 nanostructures

abstract: Among newly discovered two-dimensional (2D) materials, semiconducting ultrathin sheets of MoS2 show potential for nanoelectronics. However, the carrier mobility in MoS2 is limited by scattering from surface impurities and the substrate. To probe the sources of scattering, we use a cooled scanning probe microscope (SPM) to image the flow of electrons in a MoS2 Hall bar sample at 4.2 K. Capacitive coupling to the SPM tip changes the electron density below and scatters electrons flowing nearby; an image of flow can be obtained by measuring the change in resistance between two contacts as the tip is raster scanned across the sample. We present images of current flow through a large contact that decay exponentially away from the sample edge. In addition, the images show the characteristic "bullseye" pattern of Coulomb blockade conductance rings around a quantum dot as the density is depleted with a back gate. We estimate the size and position of these quantum dots using a capacitive model.


title: Constructing tame supercuspidal representations

abstract: A new approach to Jiu-Kang Yu's construction of tame supercuspidal representations of $p$-adic reductive groups is presented. Connections with the theory of cuspidal Deligne-Lusztig representations of finite groups of Lie type are also discussed.


title: A homotopy method for computing the largest eigenvalue of an irreducible nonnegative tensor

abstract: In this paper we propose a homotopy method to compute the largest eigenvalue and a corresponding eigenvector of a nonnegative tensor. We prove that it converges to the desired eigenpair when the tensor is irreducible. We also implement the method using an prediction-correction approach for path following. Some numerical results are provided to illustrate the efficiency of the method.


title: Stratified Splitting for Efficient Monte Carlo Integration

abstract: The efficient evaluation of high-dimensional integrals is of importance in both theoretical and practical fields of science, such as data science, statistical physics, and machine learning. However, exact computation methods suffer from the curse of dimensionality. However, due to the curse of dimensionality, deterministic numerical methods are inefficient in high-dimensional settings. Consequentially, for many practical problems, one must resort to Monte Carlo estimation. In this paper, we introduce a novel Sequential Monte Carlo technique called Stratified Splitting. The method provides unbiased estimates and can handle various integrand types including indicator functions, which are used in rare-event probability estimation problems. Moreover, we demonstrate that a variant of the algorithm can achieve polynomial complexity. The results of our numerical experiments suggest that the Stratified Splitting method is capable of delivering accurate results for a variety of integration problems.


title: A homotopy method for solving multilinear systems with M-tensors

abstract: Multilinear systems of equations arise in various applications, such as numerical partial differential equations, data mining, and tensor complementarity problems. In this paper, we propose a homotopy method for finding the unique positive solution to a multilinear system with a nonsingular M-tensor and a positive right side vector. We analyze the method and prove its convergence to the desired solution. We report some numerical results based on an implementation of the proposed method using a prediction-correction approach for path following.


title: Radial length, radial John disks and $K$-quasiconformal harmonic mappings

abstract: In this article, we continue our investigations of the boundary behavior of harmonic mappings. We first discuss the classical problem on the growth of radial length and obtain a sharp growth theorem of the radial length of $K$-quasiconformal harmonic mappings. Then we present an alternate characterization of radial John disks. In addition, we investigate the linear measure distortion and the Lipschitz continuity on $K$-quasiconformal harmonic mappings of the unit disk onto a radial John disk. Finally, using Pommerenke interior domains, we characterize certain differential properties of $K$-quasiconformal harmonic mappings.


title: The join construction

abstract: In homotopy type theory we can define the join of maps as a binary operation on maps with a common co-domain. This operation is commutative, associative, and the unique map from the empty type into the common codomain is a neutral element. Moreover, we show that the idempotents of the join of maps are precisely the embeddings, and we prove the `join connectivity theorem', which states that the connectivity of the join of maps equals the join of the connectivities of the individual maps. We define the image of a map $f:A\to X$ in $U$ via the join construction, as the colimit of the finite join powers of $f$. The join powers therefore provide approximations of the image inclusion, and the join connectivity theorem implies that the approximating maps into the image increase in connectivity. A modified version of the join construction can be used to show that for any map $f:A\to X$ in which $X$ is only assumed to be locally small, the image is a small type. We use the modified join construction to give an alternative construction of set-quotients, the Rezk completion of a precategory, and we define the $n$-truncation for any $n:\mathbb{N}$. Thus we see that each of these are definable operations on a univalent universe for Martin-L\"of type theory with a natural numbers object, that is moreover closed under homotopy coequalizers.


title: When heterodyning beats homodyning: an assessment with quadrature moments

abstract: We examine the moment-reconstruction performance of both the homodyne and heterodyne (double-homodyne) measurement schemes for arbitrary quantum states and introduce moment estimators that optimize the respective schemes for any given data. In the large-data limit, these estimators are as efficient as the maximum-likelihood estimators. We then illustrate the superiority of the heterodyne measurement for the reconstruction of the first and second moments by analyzing Gaussian states and many other significant non-classical states.


title: Sample Complexity of the Boolean Multireference Alignment Problem

abstract: The Boolean multireference alignment problem consists in recovering a Boolean signal from multiple shifted and noisy observations. In this paper we obtain an expression for the error exponent of the maximum A posteriori decoder. This expression is used to characterize the number of measurements needed for signal recovery in the low SNR regime, in terms of higher order autocorrelations of the signal. The characterization is explicit for various signal dimensions, such as prime and even dimensions.


title: An ALMA and MagAO Study of the Substellar Companion GQ Lup B

abstract: Multi-wavelength observations provide a complementary view of the formation of young directly-imaged planet-mass companions. We report the ALMA 1.3 mm and Magellan adaptive optics (MagAO) H-alpha, i', z', and Ys observations of the GQ Lup system, a classical T Tauri star with a 10-40 Mjup substellar companion at ~110 AU projected separation. We estimate the accretion rates for both components from the observed H-alpha fluxes. In our 0.05 arcsec resolution ALMA map, we resolve GQ Lup A's disk in dust continuum, but no signal is found from the companion. The disk is compact, with a radius of ~22 AU, a dust mass of ~6 Earth masses, an inclination angle of ~56 deg, and a very flat surface density profile indicative of a radial variation in dust grain sizes. No gaps or inner cavity are found in the disk, so there is unlikely a massive inner companion to scatter GQ Lup B outward. Thus, GQ Lup B might have formed in situ via disk fragmentation or prestellar core collapse. We also show that GQ Lup A's disk is misaligned with its spin axis, and possibly with GQ Lup B's orbit. Our analysis on the tidal truncation radius of GQ Lup A's disk suggests that GQ Lup B's orbit might have a low eccentricity.


title: Lepton identification at particle flow oriented detector for the future $e^{+}e^{-}$ Higgs factories

abstract: The lepton identification is essential for the physics programs at high-energy frontier, especially for the precise measurement of the Higgs boson. For this purpose, a Toolkit for Multivariate Data Analysis (TMVA) based lepton identification (LICH) has been developed for detectors using high granularity calorimeters. Using the conceptual detector geometry for the Circular Electron-Positron Collider (CEPC) and single charged particle samples with energy larger than 2 GeV, LICH identifies electrons/muons with efficiencies higher than 99.5% and controls the mis-identification rate of hadron to muons/electrons to better than 1%/0.5%. Reducing the calorimeter granularity by 1-2 orders of magnitude, the lepton identification performance is stable for particles with E > 2 GeV. Applied to fully simulated eeH/$\mu\mu$H events, the lepton identification performance is consistent with the single particle case: the efficiency of identifying all the high energy leptons in an event, is 95.5-98.5%.


title: FPGA Architecture for Deep Learning and its application to Planetary Robotics

abstract: Autonomous control systems onboard planetary rovers and spacecraft benefit from having cognitive capabilities like learning so that they can adapt to unexpected situations in-situ. Q-learning is a form of reinforcement learning and it has been efficient in solving certain class of learning problems. However, embedded systems onboard planetary rovers and spacecraft rarely implement learning algorithms due to the constraints faced in the field, like processing power, chip size, convergence rate and costs due to the need for radiation hardening. These challenges present a compelling need for a portable, low-power, area efficient hardware accelerator to make learning algorithms practical onboard space hardware. This paper presents a FPGA implementation of Q-learning with Artificial Neural Networks (ANN). This method matches the massive parallelism inherent in neural network software with the fine-grain parallelism of an FPGA hardware thereby dramatically reducing processing time. Mars Science Laboratory currently uses Xilinx-Space-grade Virtex FPGA devices for image processing, pyrotechnic operation control and obstacle avoidance. We simulate and program our architecture on a Xilinx Virtex 7 FPGA. The architectural implementation for a single neuron Q-learning and a more complex Multilayer Perception (MLP) Q-learning accelerator has been demonstrated. The results show up to a 43-fold speed up by Virtex 7 FPGAs compared to a conventional Intel i5 2.3 GHz CPU. Finally, we simulate the proposed architecture using the Symphony simulator and compiler from Xilinx, and evaluate the performance and power consumption.


title: Spherical Planetary Robot for Rugged Terrain Traversal

abstract: Wheeled planetary rovers such as the Mars Exploration Rovers (MERs) and Mars Science Laboratory (MSL) have provided unprecedented, detailed images of the Mars surface. However, these rovers are large and are of high-cost as they need to carry sophisticated instruments and science laboratories. We propose the development of low-cost planetary rovers that are the size and shape of cantaloupes and that can be deployed from a larger rover. The rover named SphereX is 2 kg in mass, is spherical, holonomic and contains a hopping mechanism to jump over rugged terrain. A small low-cost rover complements a larger rover, particularly to traverse rugged terrain or roll down a canyon, cliff or crater to obtain images and science data. While it may be a one-way journey for these small robots, they could be used tactically to obtain high-reward science data. The robot is equipped with a pair of stereo cameras to perform visual navigation and has room for a science payload. In this paper, we analyze the design and development of a laboratory prototype. The results show a promising pathway towards development of a field system.


title: Radiometric Actuators for Spacecraft Attitude Control

abstract: CubeSats and small satellites are emerging as low-cost tools to perform astronomy, exoplanet searches and earth observation. These satellites can be dedicated to pointing at targets for weeks or months at a time. This is typically not possible on larger missions where usage is shared. Current satellites use reaction wheels and where possible magneto-torquers to control attitude. However, these actuators can induce jitter due to various sources. In this work, we introduce a new class of actuators that exploit radiometric forces induced by gasses on surface with a thermal gradient. Our work shows that a CubeSat or small spacecraft mounted with radiometric actuators can achieve precise pointing of few arc-seconds or less and avoid the jitter problem. The actuator is entirely solid-state, containing no moving mechanical components. This ensures high-reliability and long-life in space. A preliminary design for these actuators is proposed, followed by feasibility analysis of the actuator performance.