title: OpenMP, OpenMP/MPI, and CUDA/MPI C programs for solving the time-dependent dipolar Gross-Pitaevskii equation

abstract: We present new versions of the previously published C and CUDA programs for solving the dipolar Gross-Pitaevskii equation in one, two, and three spatial dimensions, which calculate stationary and non-stationary solutions by propagation in imaginary or real time. Presented programs are improved and parallelized versions of previous programs, divided into three packages according to the type of parallelization. First package contains improved and threaded version of sequential C programs using OpenMP. Second package additionally parallelizes three-dimensional variants of the OpenMP programs using MPI, allowing them to be run on distributed-memory systems. Finally, previous three-dimensional CUDA-parallelized programs are further parallelized using MPI, similarly as the OpenMP programs. We also present speedup test results obtained using new versions of programs in comparison with the previous sequential C and parallel CUDA programs. The improvements to the sequential version yield a speedup of 1.1 to 1.9, depending on the program. OpenMP parallelization yields further speedup of 2 to 12 on a 16-core workstation, while OpenMP/MPI version demonstrates a speedup of 11.5 to 16.5 on a computer cluster with 32 nodes used. CUDA/MPI version shows a speedup of 9 to 10 on a computer cluster with 32 nodes.

categories: cond-mat.quant-gas cs.MS nlin.PS physics.comp-ph


title: Cosmology with a light ghost

abstract: We study the creation and evolution of cosmological perturbations in renormalizable quadratic gravity with a Weyl term. We adopt a prescription that implies the stability of the vacuum at the price of introducing a massive spin-two ghost state, leading to the loss of unitarity. The theory may still be predictive regardless the interpretation of non-unitary processes provided that their rate is negligible compared to the Universe expansion rate. This implies that the ghost is effectively stable. In such a setup, there are two scalar degrees of freedom excited during inflation. The first one is the usual curvature perturbation whose power spectrum appears to coincide with that of single-field inflation. The second one is a scalar component of the ghost encoded in the shift vector of the metric in the uniform inflaton gauge. The amplitudes of primordial tensor and vector perturbations are strongly suppressed. After inflation the ghost field starts to oscillate and its energy density shortly becomes dominant in the Universe. For all ghost masses allowed by laboratory constraints ghosts should have "overclosed" the Universe at temperatures higher than that of primordial nucleosynthesis. Thus, the model with the light Weyl ghost is ruled out.


title: Self-similar real trees defined as fixed-points and their geometric properties

abstract: We consider fixed-point equations for probability measures charging measured compact metric spaces that naturally yield continuum random trees. On the one hand, we study the existence/uniqueness of the fixed-points and the convergence of the corresponding iterative schemes. On the other hand, we study the geometric properties of the random measured real trees that are fixed-points, in particular their fractal properties. We obtain bounds on the Minkowski and Hausdorff dimension, that are proved tight in a number of applications, including the very classical continuum random tree, but also for the dual trees of random recursive triangulations of the disk introduced by Curien and Le Gall [Ann Probab, vol. 39, 2011]. The method happens to be especially efficient to treat cases for which the mass measure on the real tree induced by natural encodings only provides weak estimates on the Hausdorff dimension.


title: Vector magnetic field measurements along a cooled stereo-imaged coronal loop

abstract: The variation of the vector magnetic field along structures in the solar corona remains unmeasured. Using a unique combination of spectropolarimetry and stereoscopy, we infer and compare the vector magnetic field structure and three-dimensional morphology of an individuated coronal loop structure undergoing a thermal instability. We analyze spectropolarimetric data of the He I 10830 {\AA} triplet ($1s2s{\ }^{3}S_{1} - 1s2p{\ }^{3}P_{2,1,0}$) obtained at the Dunn Solar Telescope with the Facility Infrared Spectropolarimeter on 19 September 2011. Cool coronal loops are identified by their prominent drainage signatures in the He I data (redshifts up to 185 km sec$^{-1}$). Extinction of EUV background radiation along these loops is observed by both the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory and the Extreme Ultraviolet Imager onboard spacecraft A of the Solar Terrestrial Relations Observatory, and is used to stereoscopically triangulate the loop geometry up to heights of 70 Mm ($0.1$ $R_{sun}$) above the solar surface. The He I polarized spectra along this loop exhibit signatures indicative of atomic-level polarization as well as magnetic signatures through the Hanle and Zeeman effects. Spectropolarimetric inversions indicate that the magnetic field is generally oriented along the coronal loop axis, and provide the height dependence of the magnetic field intensity. The technique we demonstrate is a powerful one that may help better understand the thermodynamics of coronal fine structure magnetism.


title: Matrix Thermalization

abstract: Matrix quantum mechanics offers an attractive environment for discussing gravitational holography, in which both sides of the holographic duality are well-defined. Similarly to higher-dimensional implementations of holography, collapsing shell solutions in the gravitational bulk correspond in this setting to thermalization processes in the dual quantum mechanical theory. We construct an explicit, fully nonlinear supergravity solution describing a generic collapsing dilaton shell, specify the holographic renormalization prescriptions necessary for computing the relevant boundary observables, and apply them to evaluating thermalizing two-point correlation functions in the dual matrix theory.


title: Limits of Femtosecond Fiber Amplification by Parabolic Pre-Shaping

abstract: We explore parabolic pre-shaping as a means of generating and amplifying ultrashort pulses. We develop a theoretical framework for modeling the technique and use its conclusions to design a femtosecond fiber amplifier. Starting from 9 ps pulses, we obtain 4.3 $\mu$J, nearly transform-limited pulses 275 fs in duration, simultaneously achieving over 40 dB gain and 33-fold compression. Finally, we show that this amplification scheme is limited by Raman scattering, and outline a method by which the pulse duration and energy may be further improved and tailored for a given application.


title: Bounding averages rigorously using semidefinite programming: mean moments of the Lorenz system

abstract: We describe methods for proving bounds on infinite-time averages in differential dynamical systems. The methods rely on the construction of nonnegative polynomials with certain properties, similarly to the way nonlinear stability can be proved using Lyapunov functions. Nonnegativity is enforced by requiring the polynomials to be sums of squares, a condition which is then formulated as a semidefinite program (SDP) that can be solved computationally. Although such computations are subject to numerical error, we demonstrate two ways to obtain rigorous results: using interval arithmetic to control the error of an approximate SDP solution, and finding exact analytical solutions to relatively small SDPs. Previous formulations are extended to allow for bounds depending analytically on parametric variables. These methods are illustrated using the Lorenz equations, a system with three state variables $(x,y,z)$ and three parameters $(\beta,\sigma,r)$. Bounds are reported for infinite-time averages of all eighteen moments $x^ly^mz^n$ up to quartic degree that are symmetric under $(x,y)\mapsto(-x,-y)$. These bounds apply to all solutions regardless of stability, including chaotic trajectories, periodic orbits, and equilibrium points. The analytical approach yields two novel bounds that are sharp: the mean of $z^3$ can be no larger than its value of $(r-1)^3$ at the nonzero equilibria, and the mean of $xy^3$ must be nonnegative. The interval arithmetic approach is applied at the standard chaotic parameters to bound eleven average moments that all appear to be maximized on the shortest periodic orbit. Our best upper bound on each such average exceeds its value on the maximizing orbit by less than 1%. Many bounds reported here are much tighter than would be possible without computer assistance.


title: Absorption of Direct Factors With Respect to the Minimal Faithful Permutation Degree of a Finite Group

abstract: The minimal faithful permutation degree $\mu(G)$ of a finite group $G$ is the least nonnegative integer $n$ such that $G$ embeds in the symmetric group $\Sym(n)$. We prove that if $H$ is a group then $\mu(G)=\mu(G\times H)$ for some group $G$ then $H$ embeds in $A\times Q^k$ for some abelian group of odd order, some generalised quaternion $2$-group and some nonnegative integer $k$. As a consequence, $\mu(G^{n+1})=\mu(G^n)$ for some nonnegative integer $n$ if and only if $G$ is trivial.


title: Observation of semi-localized dispersive states in the strongly correlated electron-doped ferromagnet Eu$_{1-x}$Gd$_{x}$O

abstract: Chemical substitution plays a key role in controlling the electronic and magnetic properties of complex materials. For instance, in EuO, carrier doping can induce a spin-polarized metallic state, colossal magnetoresistance, and significantly enhance the Curie temperature. Here, we employ a combination of molecular-beam epitaxy, angle-resolved photoemission spectroscopy, and an effective model calculation to investigate and understand how semi-localized states evolve in lightly electron doped Eu$_{1-x}$Gd$_{x}$O above the ferromagnetic Curie temperature. Our studies reveal a characteristic length scale for the spatial extent of the donor wavefunctions which remains constant as a function of doping, consistent with recent tunneling studies of doped EuO. Our work sheds light on the nature of the semiconductor-to-metal transition in Eu$_{1-x}$Gd$_{x}$O and should be generally applicable for doped complex oxides.


title: The software package SpectralSequences

abstract: We describe the computer algebra software package SpectralSequences for the computer algebra system Macaulay2. This package implements many data types, objects and algorithms which pertain to, among other things, filtered complexes, spectral sequence pages and maps therein. We illustrate some of the syntax and capabilities of SpectralSequences by way of several examples.


title: Random Holographic "Large Worlds" with Emergent Dimensions

abstract: I propose a random network model governed by a Gaussian weight corresponding to Ising link antiferromagnetism as a model for emergent quantum space-time. In this model, discrete space is fundamental, not a regularization, its spectral dimension $d_s$ is not a model input but is, rather, completely determined by the antiferromagnetic coupling constant. Perturbative terms suppressing triangles and favouring squares lead to locally Euclidean ground states that are Ricci flat "large worlds" with power-law extension. I then consider the quenched graphs of lowest energy for $d_s=2$ and $d_s=3$ and I show how quenching leads to the spontaneous emergence of embedding spaces of Hausdorff dimension $d_H=4$ and $d_H=5$, respectively. One of the additional, spontaneous dimensions can be interpreted as time, causality being an emergent property that arises in the large $N$ limit (with $N$ the number of vertices). For $d_s=2$, the quenched graphs constitute a discrete version of a 5D-space-filling surface with a number of fundamental degrees of freedom scaling like $N^{2/5}$, a graph version of the holographic principle. These holographic degrees of freedom can be identified with the squares of the quenched graphs, which, being these triangle-free, are the fundamental area (or loop) quanta.


title: Quadratic sequences of powers and Mohanty's Conjecture

abstract: We prove under the Bombieri-Lang conjecture for surfaces that there is an absolute bound on the length of sequences of integer squares with constant second differences, for sequences which are not formed by the squares of integers in arithmetic progression. This answers a question proposed in 2010 by J. Browkin and J. Brzezinski, and independently by E. Gonzalez-Jimenez and X. Xarles. We also show that under the Bombieri-Lang conjecture for surfaces, for every $k\geq 3$ there is an absolute bound on the length of sequences formed by $k$-th powers with constant second differences. This gives a conditional result on one of Mohanty's conjectures on arithmetic progressions in Mordell's elliptic curves $y^2=x^3+b$. Moreover, we obtain an unconditional result regarding infinite families of such arithmetic progressions. We also study the case of hyperelliptic curves of the form $y^2=x^k+b$. These results are proved by unconditionally finding all curves of genus $0$ or $1$ on certain surfaces of general type. Moreover, we prove the unconditional analogues of these arithmetic results for function fields by finding all the curves of low genus on these surfaces.


title: Hydrodynamic Overview at Hot Quarks 2016

abstract: This presents an overview of relativistic hydrodynamic modeling in heavy-ion collisions prepared for Hot Quarks 2016, at South Padre Island, TX, USA. The influence of the initial state and viscosity on various experimental observables are discussed. Specific problems that arise in the hydrodynamical modeling at the Beam Energy Scan are briefly discussed.


title: Tracing the Hercules stream with Gaia and LAMOST: new evidence for a fast bar in the Milky Way

abstract: The length and pattern speed of the Milky Way bar are still controversial. Photometric and spectroscopic surveys of the inner Galaxy, as well as gas kinematics, favour a long and slowly rotating bar, with corotation around a Galactocentric radius of 6 kpc. On the other hand, the existence of the Hercules stream in local velocity space favours a short and fast bar with corotation around 4 kpc. This follows from the fact that the Hercules stream looks like a typical signature of the outer Lindblad resonance of the bar. As we showed recently, reconciling this local stream with a slow bar would need to find a yet unknown alternative explanation, based for instance on the effect of spiral arms. Here, by combining the TGAS catalogue of the Gaia DR1 with LAMOST radial velocities, we show that the position of Hercules in velocity space as a function of radius in the outer Galaxy indeed varies exactly as predicted by fast bar models with a pattern speed no less than 1.8 times the circular frequency at the Sun's position.


title: Secondary Upsilon invariants of knots

abstract: The knot invariant Upsilon, defined by Ozsvath, Stipsicz, and Szabo, induces a homomorphism from the smooth knot concordance group to the group of piecewise linear functions on the interval [0,2]. Here we define a set of related secondary invariants, each of which assigns to a knot a piecewise linear function on [0,2]. These secondary invariants provide bounds on the genus and concordance genus of knots. Examples of knots for which Upsilon vanishes but which are detected by these secondary invariants are presented.


title: A variant of Bombieri-Vinogradov theorem with explicit constants

abstract: In this paper we improve the result of Akbary and Hambrook by a factor of log x by obtaining a better version of Vaughan's inequality and using the explicit variant of an inequality connected to the M\"obius function, derived by Helfgott in his work on ternary Goldbach conjecture.


title: Spin and hyperelliptic structures of log twisted differentials

abstract: Using stable log maps, we introduce log twisted differentials extending the notion of abelian differentials to the Deligne-Mumford boundary of stable curves. The moduli stack of log twisted differentials provides a compactification of the strata of abelian differentials. The open strata can have up to three connected components, due to spin and hyperelliptic structures. We prove that the spin parity can be distinguished on the boundary of the log compactification. Moreover, combining the techniques of log geometry and admissible covers, we introduce log twisted hyperelliptic differentials, and prove that their moduli stack provides a toroidal compactification of the hyperelliptic loci in the open strata.


title: A $L^{2}$ to $L^{\infty}$ approach for the Landau Equation

abstract: Consider the Landau equation with Coulomb potential in a periodic box. We develop a new $L^{2}\rightarrow L^{\infty }$ framework to construct global unique solutions near Maxwellian with small $L^{\infty }\ $norm. The first step is to establish global $L^{2}$ estimates with strong velocity weight and time decay, under the assumption of $L^{\infty }$ bound, which is further controlled by such $L^{2}$ estimates via De Giorgi's method \cite{golse2016harnack} and \cite{mouhot2015holder}. The second step is to employ estimates in $S_{p}$ spaces to control velocity derivatives to ensure uniqueness, which is based on Holder estimates via De Giorgi's method \cite{golse2016harnack}, \cite{golse2015holder}, and \cite{mouhot2015holder}.


title: Link Prediction in evolving networks based on the popularity of nodes

abstract: Link prediction aims to uncover the underlying relationship behind networks, which could be utilized to predict the missing edges or identify the spurious edges, and attracts much attention from various fields. The key issue of link prediction is to estimate the likelihood of two nodes in networks. Most current approaches of link prediction base on static structural analysis and ignore the temporal aspects of evolving networks. Unlike previous work, in this paper, we propose a popularity based structural perturbation method (PBSPM) that characterizes the similarity of an edge not only from existing connections of networks, but also from the popularity of its two endpoints, since popular nodes have much more probability to form links between themselves. By taking popularity of nodes into account, PBSPM could suppress nodes that have high importance, but gradually become inactive. Therefore the proposed method is inclined to predict potential edges between active nodes, rather than edges between inactive nodes. Experimental results on four real networks show that the proposed method outperforms the state-of-the-art methods both in accuracy and robustness in evolving networks.


title: Markov State Models of Gene Regulatory Networks

abstract: Gene regulatory networks with dynamics characterized by multiple stable states underlie cell fate-decisions. Quantitative models that can link molecular-level knowledge of gene regulation to a global understanding of network dynamics have the potential to guide cell-reprogramming strategies. Networks are often modeled by the stochastic Chemical Master Equation, but methods for systematic identification of key properties of the global dynamics are currently lacking. We present a method for analyzing global dynamics of gene networks using the Markov State Model (MSM) framework, which utilizes a separation-of-timescales based clustering method to obtain a coarse-grained state transition graph that approximates global gene network dynamics. The method identifies the number, phenotypes, and lifetimes of long-lived network states. Application of transition path theory to the constructed MSM decomposes global dynamics into a set of dominant transition paths and associated relative probabilities for stochastic state-switching. In this proof-of-concept study, we found that the MSM provides a general framework for analyzing and visualizing stochastic multistability and state-transitions in gene networks. Our results suggest that the MSM framework, adopted from the field of atomistic Molecular Dynamics, can be a useful tool for quantitative Systems Biology at the network scale.


title: Kazhdan-Lusztig polynomials of thagomizer matroids

abstract: We introduce thagomizer matroids and compute the Kazhdan-Lusztig polynomial of a rank $n+1$ thagomizer matroid by showing that the coefficient of $t^k$ is equal to the number of Dyck paths of semilength $n$ with $k$ long ascents. We also give a conjecture for the $S_n$-equivariant Kazhdan-Lusztig polynomial of a thagomizer matroid.


title: How Well Do Local Algorithms Solve Semidefinite Programs?

abstract: Several probabilistic models from high-dimensional statistics and machine learning reveal an intriguing --and yet poorly understood-- dichotomy. Either simple local algorithms succeed in estimating the object of interest, or even sophisticated semi-definite programming (SDP) relaxations fail. In order to explore this phenomenon, we study a classical SDP relaxation of the minimum graph bisection problem, when applied to Erd\H{o}s-Renyi random graphs with bounded average degree $d>1$, and obtain several types of results. First, we use a dual witness construction (using the so-called non-backtracking matrix of the graph) to upper bound the SDP value. Second, we prove that a simple local algorithm approximately solves the SDP to within a factor $2d^2/(2d^2+d-1)$ of the upper bound. In particular, the local algorithm is at most $8/9$ suboptimal, and $1+O(1/d)$ suboptimal for large degree. We then analyze a more sophisticated local algorithm, which aggregates information according to the harmonic measure on the limiting Galton-Watson (GW) tree. The resulting lower bound is expressed in terms of the conductance of the GW tree and matches surprisingly well the empirically determined SDP values on large-scale Erd\H{o}s-Renyi graphs. We finally consider the planted partition model. In this case, purely local algorithms are known to fail, but they do succeed if a small amount of side information is available. Our results imply quantitative bounds on the threshold for partial recovery using SDP in this model.


title: Spline surfaces with T-junctions

abstract: This paper develops a new way to create smooth piecewise polynomial free-form spline surfaces from quad- meshes that include T-junctions, where surface strips start or terminate. All mesh nodes can be interpreted as control points of geometrically-smooth, piecewise polynomials that we call GT-splines. GT-splines are B-spline-like and cover T-junctions by two or four patches of degree bi-4. They complement multi-sided surface constructions in generating free-form surfaces with adaptive layout. Since GT-splines do not require a global coordination of knot intervals, GT-constructions are easy to deploy and can provide smooth surfaces with T-junctions where T-splines can not have a smooth parameterization. GT-constructions display a uniform highlight line distribution on input meshes where alternatives, such as Catmull-Clark subdivision, exhibit oscillations.


title: Effective electro-optic modulation in low-loss graphene-plasmonic slot waveguides

abstract: Surface plasmon polaritons enable light concentration within subwavelength regions, opening thereby new avenues for miniaturizing the device and strengthening light-matter interactions. Here we realize effective electro-optic modulation in low-loss plasmonic waveguides with the aid of graphene, and the devices are fully integrated in the silicon-on-insulator platform. By advantageously exploiting low-loss plasmonic slot-waveguide modes, which weakly leak into a substrate while feature strong fields within the two-layer-graphene covered slots in metal, we successfully achieve a tunability of 0.13 dB/um for our fabricated graphene-plasmonic waveguide devices with extremely low insertion loss, which outperforms previously reported graphene-plasmonic devices. Our results highlight the potential of graphene plasmonic leaky-mode hybrid waveguides to realized active ultra-compact devices for optoelectronic applications.


title: Classification of Homogeneous Fourier Matrices

abstract: Modular data are commonly studied in mathematics and physics. A modular datum defines a finite-dimensional representation of the modular group $SL_2(\mathbb{Z})$. In this paper, we show that there is a one-to-one correspondence between Fourier matrices associated to modular data and self-dual $C$-algebras that satisfy a certain condition. Also, we prove that a homogenous $C$-algebra arising from a Fourier matrix has all the degrees equal to $1$.


title: Reconstructing Partitions from their Multisets of $k$-Minors

abstract: For non-negative integers $n$ and $k$ with $n \ge k$, a {\em $k$-minor} of a partition $\lambda = [\lambda_1, \lambda_2, \dots]$ of $n$ is a partition $\mu = [\mu_1, \mu_2, \dots]$ of $n-k$ such that $\mu_i \le \lambda_i$ for all $i$. The multiset $\widehat{M}_k(\lambda)$ of $k$-minors of $\lambda$ is defined as the multiset of $k$-minors $\mu$ with multiplicity of $\mu$ equal to the number of standard Young tableaux of skew shape $\lambda / \mu$. We show that there exists a function $G(n)$ such that the partitions of $n$ can be reconstructed from their multisets of $k$-minors if and only if $k \le G(n)$. Furthermore, we prove that $\lim_{n \rightarrow \infty} G(n)/n = 1$ with $n-G(n) = O(n/\log n)$. As a direct consequence of this result, the irreducible representations of the symmetric group $S_n$ can be reconstructed from their restrictions to $S_{n-k}$ if and only if $k \le G(n)$ for the same function $G(n)$. For a minor $\mu$ of the partition $\lambda$, we study the excitation factor $E_\mu (\lambda)$, which appears as a crucial part in Naruse's Skew-Shape Hook Length Formula. We observe that certain excitation factors of $\lambda$ can be expressed as a $\mathbb{Q}[k]$-linear combination of the elementary symmetric polynomials of the hook lengths in the first row of $\lambda$ where $k = \lambda_1$ is the number of cells in the first row of $\lambda$.


title: A simple finite element method for the Stokes equations

abstract: The goal of this paper is to introduce a simple finite element method to solve the Stokes and the Navier-Stokes equations. This method is in primal velocity-pressure formulation and is so simple such that both velocity and pressure are approximated by piecewise constant functions. Implementation issues as well as error analysis are investigated. A basis for a divergence free subspace of the velocity field is constructed so that the original saddle point problem can be reduced to a symmetric and positive definite system with much fewer unknowns. The numerical experiments indicate that the method is accurate and robust.


title: An alternative view on dissipation in turbulent flows

abstract: An original experimental setup has been elaborated in order to get a better view of turbulent flows in a von Karman geometry. The availability of a very fast camera allowed to follow in time the evolution of the flows. A surprising finding is that the development of smaller whorls ceases earlier than expected and the aspect of the flows remains the same above Reynolds number of a few thousand. This fact provides an explanation of the constancy of the reduced dissipation in the same range without the need of singularity. Its cause could be in relation with the same type of behavior observed in a rotating frame.


title: The Standard-Model Extension and Gravitational Tests

abstract: The Standard-Model Extension (SME) provides a comprehensive effective field-theory framework for the study of CPT and Lorentz symmetry. This work reviews the structure and philosophy of the SME and provides some intuitive examples of symmetry violation. The results of recent gravitational tests performed within the SME are summarized including analysis of results from the Laser Interferometer Gravitational-Wave Observatory (LIGO), sensitivities achieved in short-range gravity experiments, constraints from cosmic-ray data, and results achieved by studying planetary ephemerids. Some proposals and ongoing efforts will also be considered including gravimeter tests, tests of the Weak Equivalence Principle, and antimatter experiments. Our review of the above topics is augmented by several original extensions of the relevant work. We present new examples of symmetry violation in the SME and use the cosmic-ray analysis to place first-ever constraints on 81 additional operators.


title: Psycho-Demographic Analysis of the Facebook Rainbow Campaign

abstract: Over the past decade, online social media has had a tremendous impact on the way people engage in social activism. For instance, about 26M Facebook users expressed their support in upholding the cause of marriage equality by overlaying their profile pictures with rainbow-colored filters. Similarly, hundreds of thousands of users changed their profile pictures to a black dot condemning incidents of sexual violence in India. This act of demonstrating support for social causes by changing online profile pictures is being referred to as pictivism. In this paper, we analyze the psycho-demographic profiles, social networking behavior, and personal interests of users who participated in the Facebook Rainbow campaign. Our study is based on a sample of about 800K detailed profiles of Facebook users combining questionnaire-based psychological scores with Facebook profile data. Our analysis provides detailed insights into psycho-demographic profiles of the campaign participants. We found that personality traits such as openness and neuroticism are both positively associated with the likelihood of supporting the campaign, while conscientiousness exhibited a negative correlation. We also observed that females, religious disbelievers, democrats and adults in the age group of 20 to 30 years are more likely to be a part of the campaign. Our research further confirms the findings of several previous studies which suggest that a user is more likely to participate in an online campaign if a large fraction of his/her friends are already doing so. We also developed machine learning models for predicting campaign participation. Users' personal interests, approximated by Facebook user like activity, turned out to be the best indicator of campaign participation. Our results demonstrated that a predictive model which leverages the aforementioned features accurately identifies campaign participants (AUC=0.76).


title: A Machine Learns to Predict the Stability of Tightly Packed Planetary Systems

abstract: The requirement that planetary systems be dynamically stable is often used to vet new discoveries or set limits on unconstrained masses or orbital elements. This is typically carried out via computationally expensive N-body simulations. We show that characterizing the complicated and multi-dimensional stability boundary of tightly packed systems is amenable to machine learning methods. We find that training an XGBoost machine learning algorithm on physically motivated features yields an accurate classifier of stability in packed systems. On the stability timescale investigated ($10^7$ orbits), it is 3 orders of magnitude faster than direct N-body simulations. Optimized machine learning classifiers for dynamical stability may thus prove useful across the discipline, e.g., to characterize the exoplanet sample discovered by the upcoming Transiting Exoplanet Survey Satellite (TESS). This proof of concept motivates investing computational resources to train algorithms capable of predicting stability over longer timescales and over broader regions of phase space.


title: Universality of the nodal length of bivariate random trigonometric polynomials

abstract: We consider random trigonometric polynomials of the form \[ f_n(x,y)=\sum_{1\le k,l \le n} a_{k,l} \cos(kx) \cos(ly), \] where the entries $(a_{k,l})_{k,l\ge 1}$ are i.i.d. random variables that are centered with unit variance. We investigate the length $\ell_K(f_n)$ of the nodal set $Z_K(f_n)$ of the zeros of $f_n$ that belong to a compact set $K \subset \mathbb R^2$. We first establish a local universality result, namely we prove that, as $n$ goes to infinity, the sequence of random variables $n\, \ell_{K/n}(f_n)$ converges in distribution to a universal limit which does not depend on the particular law of the entries. We then show that at a macroscopic scale, the expectation of $\ell_{[0,\pi]^2}(f_n)/n$ also converges to an universal limit. Our approach provides two main byproducts: (i) a general result regarding the continuity of the volume of the nodal sets with respect to $C^1$-convergence which refines previous findings of Rusakov et al., Iksanov et al. and Aza\"is et al., and (ii) a new strategy for proving small ball estimates in random trigonometric models, providing in turn uniform local controls of the nodal volumes.


title: End-to-end attention-based distant speech recognition with Highway LSTM

abstract: End-to-end attention-based models have been shown to be competitive alternatives to conventional DNN-HMM models in the Speech Recognition Systems. In this paper, we extend existing end-to-end attention-based models that can be applied for Distant Speech Recognition (DSR) task. Specifically, we propose an end-to-end attention-based speech recognizer with multichannel input that performs sequence prediction directly at the character level. To gain a better performance, we also incorporate Highway long short-term memory (HLSTM) which outperforms previous models on AMI distant speech recognition task.


title: The Observable Signatures of GRB Cocoons

abstract: As a long GRB jet propagates within the surrounding stellar atmosphere it creates a cocoon composed of an outer Newtonian shocked stellar material and an inner (possibly relativistic) shocked jet material. The jet deposits $10^{51}-10^{52}$ erg into this cocoon. This energy is comparable to the GRB's energy and to the energy of the accompanying supernova, yet its signature has been largely neglected so far. A fraction of the cocoon energy is released during its expansion following the breakout from the star and later as it interacts with the surrounding matter. We explore here the possible signatures of the cocoon emission and outline a framework to calculate them from the conditions of the cocoon at the time of the jet breakout. We show that the cocoon signature depends strongly on the level of mixing between the shocked jet and shocked stellar material that fills it, which is currently unknown. We find that if there is no mixing at all then the $\gamma$-ray emission from the cocoon is so bright that it should have been already detected, and the lack of such detections indicates that mixing at some level must take place. We calculate also the expected signal for partial and full mixing. While the typical signals are weaker than GRBs' afterglows, the latter are highly beamed while the former have wide angles. We predict that future optical, UV and X-ray transient searches, like LSST, ZTF, ULTRASAT, ISS-Lobster and others will most likely detect such signals, providing a wealth of information on the progenitors and jets of GRBs. While we focus on long GRBs, we note that analogous (but weaker) cocoons may arise in short GRBs as well. Their signatures might be the most promising electromagnetic counterparts for gravitational waves merger's signals.


title: All-dielectric nanophotonics: fundamentals, fabrication, and applications

abstract: In this Article, we review a novel, rapidly developing field of modern light science named all-dielectric nanophotonics. This branch of nanophotonics is based on the properties of high-index dielectric nanoparticles which allow for controlling both magnetic and electric responses of a nanostructured matter. Here, we discuss optical properties of high-index dielectric nanoparticles, methods of their fabrication, and recent advances in practical applications, including the quantum source emission engineering, Fano resonances in all-dielectric nanoclusters, surface enhanced spectroscopy and sensing, coupled-resonator optical waveguides, metamaterials and metasurfaces, and nonlinear nanophotonics.


title: Topological quantum phase transition and superconductivity induced by pressure in the bismuth tellurohalide BiTeI

abstract: A pressure-induced topological quantum phase transition has been theoretically predicted for the semiconductor BiTeI with giant Rashba spin splitting. In this work, the evolution of the electrical transport properties in BiTeI and BiTeBr is investigated under high pressure. The pressure-dependent resistivity in a wide temperature range passes through a minimum at around 3 GPa, indicating the predicted transition in BiTeI. Superconductivity is observed in both BiTeI and BiTeBr while the resistivity at higher temperatures still exhibits semiconducting behavior. Theoretical calculations suggest that the superconductivity may develop from the multi-valley semiconductor phase. The superconducting transition temperature Tc increases with applied pressure and reaches a maximum value of 5.2 K at 23.5 GPa for BiTeI (4.8 K at 31.7 GPa for BiTeBr), followed by a slow decrease. Our results demonstrate that BiTeX (X = I, Br) compounds with non-trivial topology of electronic states display new ground states upon compression.


title: The structure of almost Abelian Lie algebras

abstract: An almost Abelian Lie algebra is a non-Abelian Lie algebra with a codimension 1 Abelian ideal. Most 3-dimensional real Lie algebras are almost Abelian, and they appear in every branch of physics that deals with anisotropic media - cosmology, crystallography etc. In differential geometry and theoretical physics, almost Abelian Lie groups have given rise to some of the simplest solvmanifolds on which various geometric structures such as symplectic, K\"ahler, spin etc., are currently studied in explicit terms. However, a systematic study of almost Abelian Lie groups and algebras from mathematics perspective has not been carried out yet, and the present paper is the first step in addressing this wide and diverse class of groups and algebras. The present paper studies the structure and important algebraic properties of almost Abelian Lie algebras of arbitrary dimension over any field of scalars. A classification of almost Abelian Lie algebras is given. All Lie subalgebras and ideals, automorphisms and derivations, Lie orthogonal operators and quadratic Casimir elements are described exactly.


title: Foot force models of crowd dynamics on a wobbly bridge

abstract: Modern pedestrian and suspension bridges are designed using industry-standard packages, yet disastrous resonant vibrations are observed, necessitating multi-million dollar repairs. Recent examples include pedestrian induced vibrations during the openings of the Solf\'erino Bridge in Paris in 1999 and the increased bouncing of the Squibb Park Bridge in Brooklyn in 2014. The most prominent example of an unstable lively bridge is the London Millennium Bridge which started wobbling as a result of pedestrian-bridge interactions. Pedestrian phase-locking due to footstep phase adjustment, is suspected to be the main cause of its large lateral vibrations; however, its role in the initiation of wobbling was debated. In this paper, we develop foot force models of pedestrians' response to bridge motion and detailed, yet analytically tractable models of crowd phase-locking. We use bio-mechanically inspired models of crowd lateral movement to investigate to what degree pedestrian synchrony must be present for a bridge to wobble significantly and what is a critical crowd size. Our results can be used as a safety guideline for designing pedestrian bridges or limiting the maximum occupancy of an existing bridge. The pedestrian models can be used as "crash test dummies" when numerically probing a specific bridge design. This is particularly important because the US code for designing pedestrian bridges does not contain explicit guidelines that account for the collective pedestrian behavior.


title: The random matrix hard edge: rare events and a transition

abstract: We study probabilities of various rare events for the limiting point process that appears at the random matrix hard edge. We also show a transition from hard edge to bulk behavior. Asymptotic events studied include a central limit theorem and large deviation result for the number of points in a growing interval $[0,\lambda]$ as $\lambda \to \infty$. We study these results for the square root of the hard edge process. In this setting many of these behaviors mimic those of the Sine$_\beta$ process.


title: Monolithic CMOS-compatible zero-index metamaterials

abstract: Zero-index materials exhibit exotic optical properties that can be utilized for integrated-optics applications. However, practical implementation requires compatibility with complementary metallic-oxide-semiconductor (CMOS) technologies. We demonstrate a CMOS-compatible zero-index metamaterial consisting of a square array of air holes in a 220-nm-thick silicon-on-insulator (SOI) wafer. This design is achieved through a Dirac-cone dispersion. The metamaterial is entirely composed of silicon and offers compatibility through low-aspect-ratio structures that can be simply fabricated in a standard device layer. This platform enables mass adoption and exploration of zero-index-based photonic devices at low cost and high fidelity.


title: The Final SDSS High-Redshift Quasar Sample of 52 Quasars at z>5.7

abstract: We present the discovery of nine quasars at $z\sim6$ identified in the Sloan Digital Sky Survey (SDSS) imaging data. This completes our survey of $z\sim6$ quasars in the SDSS footprint. Our final sample consists of 52 quasars at $5.7<z\le6.4$, including 29 quasars with $z_{\rm AB}\le20$ mag selected from 11,240 deg$^2$ of the SDSS single-epoch imaging survey (the main survey), 10 quasars with $20\le z_{\rm AB}\le20.5$ selected from 4223 deg$^2$ of the SDSS overlap regions (regions with two or more imaging scans), and 13 quasars down to $z_{\rm AB}\approx22$ mag from the 277 deg$^2$ in Stripe 82. They span a wide luminosity range of $-29.0\le M_{1450}\le-24.5$. This well-defined sample is used to derive the quasar luminosity function (QLF) at $z\sim6$. After combining our SDSS sample with two faint ($M_{1450}\ge-23$ mag) quasars from the literature, we obtain the parameters for a double power-law fit to the QLF. The bright-end slope $\beta$ of the QLF is well constrained to be $\beta=-2.8\pm0.2$. Due to the small number of low-luminosity quasars, the faint-end slope $\alpha$ and the characteristic magnitude $M_{1450}^{\ast}$ are less well constrained, with $\alpha=-1.90_{-0.44}^{+0.58}$ and $M^{\ast}=-25.2_{-3.8}^{+1.2}$ mag. The spatial density of luminous quasars, parametrized as $\rho(M_{1450}<-26,z)=\rho(z=6)\,10^{k(z-6)}$, drops rapidly from $z\sim5$ to 6, with $k=-0.72\pm0.11$. Based on our fitted QLF and assuming an IGM clumping factor of $C=3$, we find that the observed quasar population cannot provide enough photons to ionize the $z\sim6$ IGM at $\sim90$\% confidence. Quasars may still provide a significant fraction of the required photons, although much larger samples of faint quasars are needed for more stringent constraints on the quasar contribution to reionization.


title: An accurate metric for the spacetime around neutron stars

abstract: The problem of having an accurate description of the spacetime around neutron stars is of great astrophysical interest. For astrophysical applications, one needs to have a metric that captures all the properties of the spacetime around a neutron star. Furthermore, an accurate appropriately parameterised metric, i.e., a metric that is given in terms of parameters that are directly related to the physical structure of the neutron star, could be used to solve the inverse problem, which is to infer the properties of the structure of a neutron star from astrophysical observations. In this work we present such an approximate stationary and axisymmetric metric for the exterior of neutron stars, which is constructed using the Ernst formalism and is parameterised by the relativistic multipole moments of the central object. This metric is given in terms of an expansion on the Weyl-Papapetrou coordinates with the multipole moments as free parameters and is shown to be extremely accurate in capturing the physical properties of a neutron star spacetime as they are calculated numerically in general relativity. Because the metric is given in terms of an expansion, the expressions are much simpler and easier to implement, in contrast to previous approaches. For the parameterisation of the metric in general relativity, the recently discovered universal 3-hair relations are used to produce a 3-parameter metric. Finally, a straightforward extension of this metric is given for scalar-tensor theories with a massless scalar field, which also admit a formulation in terms of an Ernst potential.


title: (2,2) Superconformal Bootstrap in Two Dimensions

abstract: We find a simple relation between two-dimensional BPS N=2 superconformal blocks and bosonic Virasoro conformal blocks, which allows us to analyze the crossing equations for BPS 4-point functions in unitary (2,2) superconformal theories numerically with semidefinite programming. We constrain gaps in the non-BPS spectrum through the operator product expansion of BPS operators, in ways that depend on the moduli of exactly marginal deformations through chiral ring coefficients. In some cases, our bounds on the spectral gaps are observed to be saturated by free theories, by N=2 Liouville theory, and by certain Landau-Ginzburg models.


title: The Inviscid Limit and Boundary Layers for Navier-Stokes Flows

abstract: The validity of the vanishing viscosity limit, that is, whether solutions of the Navier-Stokes equations modeling viscous incompressible flows converge to solutions of the Euler equations modeling inviscid incompressible flows as viscosity approaches zero, is one of the most fundamental issues in mathematical fluid mechanics. The problem is classified into two categories: the case when the physical boundary is absent, and the case when the physical boundary is present and the effect of the boundary layer becomes significant. The aim of this article is to review recent progress on the mathematical analysis of this problem in each category.


title: The SOFIA Massive (SOMA) Star Formation Survey. I. Overview and First Results

abstract: We present an overview and first results of the Stratospheric Observatory For Infrared Astronomy Massive (SOMA) Star Formation Survey, which is using the FORCAST instrument to image massive protostars from $\sim10$--$40\:\rm{\mu}\rm{m}$. These wavelengths trace thermal emission from warm dust, which in Core Accretion models mainly emerges from the inner regions of protostellar outflow cavities. Dust in dense core envelopes also imprints characteristic extinction patterns at these wavelengths, causing intensity peaks to shift along the outflow axis and profiles to become more symmetric at longer wavelengths. We present observational results for the first eight protostars in the survey, i.e., multiwavelength images, including some ancillary ground-based MIR observations and archival {\it{Spitzer}} and {\it{Herschel}} data. These images generally show extended MIR/FIR emission along directions consistent with those of known outflows and with shorter wavelength peak flux positions displaced from the protostar along the blueshifted, near-facing sides, thus confirming qualitative predictions of Core Accretion models. We then compile spectral energy distributions and use these to derive protostellar properties by fitting theoretical radiative transfer models. Zhang and Tan models, based on the Turbulent Core Model of McKee and Tan, imply the sources have protostellar masses $m_*\sim10$--50$\:M_\odot$ accreting at $\sim10^{-4}$--$10^{-3}\:M_\odot\:{\rm{yr}}^{-1}$ inside cores of initial masses $M_c\sim30$--500$\:M_\odot$ embedded in clumps with mass surface densities $\Sigma_{\rm{cl}}\sim0.1$--3$\:{\rm{g\:cm}^{-2}}$. Fitting Robitaille et al. models typically leads to slightly higher protostellar masses, but with disk accretion rates $\sim100\times$ smaller. We discuss reasons for these differences and overall implications of these first survey results for massive star formation theories.


title: Stable computations with flat radial basis functions using vector-valued rational approximations

abstract: One commonly finds in applications of smooth radial basis functions (RBFs) that scaling the kernels so they are `flat' leads to smaller discretization errors. However, the direct numerical approach for computing with flat RBFs (RBF-Direct) is severely ill-conditioned. We present an algorithm for bypassing this ill-conditioning that is based on a new method for rational approximation (RA) of vector-valued analytic functions with the property that all components of the vector share the same singularities. This new algorithm (RBF-RA) is more accurate, robust, and easier to implement than the Contour-Pad\'e method, which is similarly based on vector-valued rational approximation. In contrast to the stable RBF-QR and RBF-GA algorithms, which are based on finding a better conditioned base in the same RBF-space, the new algorithm can be used with any type of smooth radial kernel, and it is also applicable to a wider range of tasks (including calculating Hermite type implicit RBF-FD stencils). We present a series of numerical experiments demonstrating the effectiveness of this new method for computing RBF interpolants in the flat regime. We also demonstrate the flexibility of the method by using it to compute implicit RBF-FD formulas in the flat regime and then using these for solving Poisson's equation in a 3D spherical shell.


title: Compact Linearization for Binary Quadratic Problems subject to Assignment Constraints

abstract: We prove new necessary and sufficient conditions to carry out a compact linearization approach for a general class of binary quadratic problems subject to assignment constraints as it has been proposed by Liberti in 2007. The new conditions resolve inconsistencies that can occur when the original method is used. We also present a mixed-integer linear program to compute a minimally-sized linearization. When all the assignment constraints have non-overlapping variable support, this program is shown to have a totally unimodular constraint matrix. Finally, we give a polynomial-time combinatorial algorithm that is exact in this case and can still be used as a heuristic otherwise.


title: Probabilistic Safety Programs

abstract: Achieving safe control under uncertainty is a key problem that needs to be tackled for enabling real-world autonomous robots and cyber-physical systems. This paper introduces Probabilistic Safety Programs (PSP) that embed both the uncertainty in the environment as well as invariants that determine safety parameters. The goal of these PSPs is to evaluate future actions or trajectories and determine how likely it is that the system will stay safe under uncertainty. We propose to perform these evaluations by first compiling the PSP to a graphical model then using a fast variational inference algorithm. We highlight the efficacy of the framework on the task of safe control of quadrotors and autonomous vehicles in dynamic environments.


title: Optimizing Open-Ended Crowdsourcing: The Next Frontier in Crowdsourced Data Management

abstract: Crowdsourcing is the primary means to generate training data at scale, and when combined with sophisticated machine learning algorithms, crowdsourcing is an enabler for a variety of emergent automated applications impacting all spheres of our lives. This paper surveys the emerging field of formally reasoning about and optimizing open-ended crowdsourcing, a popular and crucially important, but severely understudied class of crowdsourcing---the next frontier in crowdsourced data management. The underlying challenges include distilling the right answer when none of the workers agree with each other, teasing apart the various perspectives adopted by workers when answering tasks, and effectively selecting between the many open-ended operators appropriate for a problem. We describe the approaches that we've found to be effective for open-ended crowdsourcing, drawing from our experiences in this space.


title: Redundant interdependencies boost the robustness of multilayer networks

abstract: In the analysis of the robustness of multiplex networks, it is commonly assumed that a node is functioning only if its interdependent nodes are simultaneously functioning. According to this model, a multiplex network becomes more and more fragile as the number of layers increases. In this respect, the addition of a new layer of interdependent nodes to a preexisting multiplex network will never improve its robustness. Whereas such a model seems appropriate to understand the effect of interdependencies in the simplest scenario of a network composed of only two layers, it may seem not suitable to characterize the robustness of real systems formed by multiple network layers. It seems in fact unrealistic that a real system, evolved, through the development of multiple layers of interactions, towards a fragile structure. In this paper, we introduce a model of percolation where the condition that makes a node functional is that the node is functioning in at least two of the layers of the network. The model reduces to the commonly adopted percolation model for multiplex networks when the number of layers equals two. For larger number of layers, however, the model describes a scenario where the addition of new layers boosts the robustness of the system by creating redundant interdependencies among layers. We prove this fact thanks to the development of a message-passing theory able to characterize the model in both synthetic and real-world multiplex graphs.


title: 3D Higher spin gravity and the fractional quantum Hall effect

abstract: This article is based on the talk "Fractional Spin Gravity" presented in the 31st International Colloquium on Group Theoretical Methods in Physics, Rio de Janeiro, 19-25th June 2016. There we emphasised an implication of the works [1,2] by N. Boulanger, P. Sundell and the author on fractional spin extensions of 2+1D higher spin gravity. This is that higher spin gravity may govern interactions of pseudo-particles excitations in the (fractional) quantum Hall effect. More generally, fractional spin currents in 2+1D source higher spin gravity curvatures.


title: Cancellation in the Additive Twists of Fourier Coefficients for $\mathrm{GL}_2$ and $\mathrm{GL}_3$ over Number Fields

abstract: In this article, we study the sum of additively twisted Fourier coefficients of an irreducible cuspidal automorphic representation of $\mathrm{GL}_2$ or $\mathrm{GL}_3$ over an arbitrary number field. When the representation is unramified at all non-archimedean places, we prove the Wilton type bound for $\mathrm{GL}_2$ and the Miller type bound for $\mathrm{GL}_3$ which are uniform in terms of the additive character.


title: Fully Hodge-Newton decomposable Shimura varieties

abstract: The motivation for this paper is the study of arithmetic properties of Shimura varieties, in particular the Newton stratification of the special fiber of a suitable integral model at a prime with parahoric level structure. This is closely related to the structure of Rapoport-Zink spaces and of affine Deligne-Lusztig varieties. We prove a Hodge-Newton decomposition for affine Deligne-Lusztig varieties and for the special fibres of Rapoport-Zink spaces, relating these spaces to analogous ones defined in terms of Levi subgroups, under a certain condition (Hodge-Newton decomposability) which can be phrased in combinatorial terms. Second, we study the Shimura varieties in which every non-basic $\sigma$-isogeny class is Hodge-Newton decomposable. We show that (assuming the axioms of \cite{HR}) this condition is equivalent to nice conditions on either the basic locus, or on all the non-basic Newton strata of the Shimura varieties. We also give a complete classification of Shimura varieties satisfying these conditions. While previous results along these lines often have restrictions to hyperspecial (or at least maximal parahoric) level structure, and/or quasi-split underlying group, we handle the cases of arbitrary parahoric level structure, and of possibly non-quasi-split underlying groups. This results in a large number of new cases of Shimura varieties where a simple description of the basic locus can be expected. As a striking consequence of the results, we obtain that this property is independent of the parahoric subgroup chosen as level structure. We expect that our conditions are closely related to the question whether the weakly admissible and admissible loci coincide.


title: Towards an automated tool to evaluate the impact of the nuclear modification of the gluon density on quarkonium, D and B meson production in proton-nucleus collisions

abstract: We propose a simple and model-independent procedure to account for the impact of the nuclear modification of the gluon density as encoded in nuclear collinear PDF sets on two-to-two partonic hard processes in proton-nucleus collisions. This applies to a good approximation to quarkonium, D and B meson production, generically referred to H. Our procedure consists in parametrising the square of the parton scattering amplitude, A_{gg -> H X} and constraining it from the proton-proton data. Doing so, we have been able to compute the corresponding nuclear modification factors for J/psi, Upsilon and D^0 as a function of y and P_T at sqrt(s_NN)=5 and 8 TeV in the kinematics of the various LHC experiments in a model independent way. It is of course justified since the most important ingredient in such evaluations is the probability of each kinematical configuration. Our computations for D mesons can also be extended to B meson production. To further illustrate the potentiality of the tool, we provide --for the first time-- predictions for the nuclear modification factor for eta_c production in pPb collisions at the LHC.


title: Detection of intensity bursts using Hawkes processes: an application to high frequency financial data

abstract: Given a stationary point process, an intensity burst is defined as a short time period during which the number of counts is larger than the typical count rate. It might signal a local non-stationarity or the presence of an external perturbation to the system. In this paper we propose a novel procedure for the detection of intensity bursts within the Hawkes process framework. By using a model selection scheme we show that our procedure can be used to detect intensity bursts when both their occurrence time and their total number is unknown. Moreover, the initial time of the burst can be determined with a precision given by the typical inter-event time. We apply our methodology to the mid-price change in FX markets showing that these bursts are frequent and that only a relatively small fraction is associated to news arrival. We show lead-lag relations in intensity burst occurrence across different FX rates and we discuss their relation with price jumps.


title: A Short Guide to Anyons and Modular Functors

abstract: To the working physicist, anyon theory is meant to describe certain quasi-particle excitations occurring in two dimensional topologically ordered systems. A typical calculation using this theory will involve operations such as $\otimes$ to combine anyons, $F^{abc}_{d}$ to re-associate such combinations, and $R^{ab}_c$ to commute or braid these anyons. Although there is a powerful string-diagram notation that greatly assists these manipulations, we still appear to be operating on particles arranged on a one-dimensional line, algebraically ordered from left to right. The obvious question is, where is the other dimension? The topological framework for considering these anyons as truly living in a two dimensional space is known as a modular functor, or topological quantum field theory. In this work we show how the apparently one-dimensional algebraic anyon theory is secretly the theory of anyons living in a fully two-dimensional system. The mathematical literature covering this secret is vast, and we try to distill this down into something more manageable.


title: Laser remote magnetometry using mesospheric sodium

abstract: We have demonstrated a remote magnetometer based on sodium atoms in the Earth's mesosphere, at a 106-kilometer distance from our instrument. A 1.33-watt laser illuminated the atoms, and the magnetic field was inferred from back-scattered light collected by a telescope with a 1.55-meter-diameter aperture. The measurement sensitivity was 162 nT/$\sqrt{Hz}$. The value of magnetic field inferred from our measurement is consistent with an estimate based on the Earth's known field shape to within a fraction of a percent. Projected improvements in optics could lead to sensitivity of 20 nT/$\sqrt{Hz}$, and the use of advanced lasers or a large telescope could approach 1-nT/$\sqrt{Hz}$ sensitivity. All experimental and theoretical sensitivity values are based on a 60$^\circ$ angle between the laser beam axis and the magnetic field vector; at the optimal 90$^\circ$ angle sensitivity would be improved by about a factor of two.

categories: astro-ph.IM astro-ph.EP physics.atom-ph physics.ins-det physics.space-ph


title: Squeezing giant spin states via geometric phase control in cavity-assisted Raman transitions

abstract: Squeezing ensemble of spins provides a way to surpass the standard quantum limit (SQL) in quantum metrology and test the fundamental physics as well, and therefore attracts broad interest. Here we propose an experimentally accessible protocol to squeeze a giant ensemble of spins via the geometric phase control. Using the cavity-assisted Raman transitions in a double $\Lambda$-type system, we realize an effective Dicke model. Under the condition of vanishing effective spin transition frequency, we find a particular evolution time where the cavity decouples from the spins and the spin ensemble is squeezed considerably. Our scheme has the potential to improve the sensitivity in quantum metrology with spins by about two orders.


title: Conjectures involving a generalization of the sums of powers of integers

abstract: Some class of sums which naturally include the sums of powers of integers is considered. A number of conjectures concerning a representation of these sums is made.


title: Quantum estimation in neutrino oscillations

abstract: Neutrino oscillations are at the forefront of advances in Physics beyond the Standard Model. Increasing accuracy in measurements of the neutrino mixing matrix is an important challenge in current experiments. It depends on parameters that do not directly correspond to observables of the neutrino system. This type of estimation problem is handled by Quantum Estimation Theory (QET) via the Fisher Information (FI) and the Quantum Fisher Information (QFI). In this work, we analyze two-flavor neutrino oscillations within the framework of QET. We compute the QFI for the mixing angle $\theta$ and show that mass measurements are the ones that achieve optimal precision. We also study the FI associated with flavor measurements and show that they are optimized at specific neutrino times-of-flight. Therefore, although the usual population measurement does not realize the precision limit set by the QFI, it can in principle be implemented with the best possible sensitivity to $\theta $. We study how these quantifiers relate to the single-particle, mode entanglement. We demonstrate that this form of entanglement does not enhance neither of them. In particular, this shows that in single-particle settings, entanglement is not directly connected with the optimal precision in metrological tasks.


title: Singe-photon multi-ports router based on the coupled cavity optomechanical system

abstract: A scheme of single-photon multi-port router is put forward by coupling two optomechanical cavities with waveguides. It is shown that the coupled two optomechanical cavities can exhibit photon blockade effect, which is generated from interference of three mode interaction. A single-photon travel along the system is calculated. The results show that the single photon can be controlled in the multi-port system because of the radiation pressure, which should be useful for constructing quantum network.


title: Explicit symplectic algorithms based on generating functions for relativistic charged particle dynamics in time-dependent electromagnetic field

abstract: Relativistic dynamics of a charged particle in time-dependent electromagnetic fields has theoretical significance and a wide range of applications. It is often multi-scale and requires accurate long-term numerical simulations using symplectic integrators. For modern large-scale particle simulations in complex, time-dependent electromagnetic field, explicit symplectic algorithms are much more preferable. In this paper, we treat the relativistic dynamics of a particle as a Hamiltonian system on the cotangent space of the space-time, and construct for the first time explicit symplectic algorithms for relativistic charged particles of order 2 and 3 using the sum-split technique and generating functions.


title: Transport exponents of states with large support

abstract: We extend results of Damanik and Tcheremchantsev on estimating transport exponents to initial states supported on more than one site. These general results for upper and lower bounds are then applied to several classes of models, including Sturmian, quiasi-periodic and substitution-generated potentials, and random polymer model.


title: AutoGP: Exploring the Capabilities and Limitations of Gaussian Process Models

abstract: We investigate the capabilities and limitations of Gaussian process models by jointly exploring three complementary directions: (i) scalable and statistically efficient inference; (ii) flexible kernels; and (iii) objective functions for hyperparameter learning alternative to the marginal likelihood. Our approach outperforms all previously reported GP methods on the standard MNIST dataset; performs comparatively to previous kernel-based methods using the RECTANGLES-IMAGE dataset; and breaks the 1% error-rate barrier in GP models using the MNIST8M dataset, showing along the way the scalability of our method at unprecedented scale for GP models (8 million observations) in classification problems. Overall, our approach represents a significant breakthrough in kernel methods and GP models, bridging the gap between deep learning approaches and kernel machines.


title: On the Stability of Satellite Planes I: Effects of Mass, Velocity, Halo Shape and Alignment

abstract: The recently discovered vast thin plane of dwarf satellites orbiting the Andromeda Galaxy (M31) adds to the mystery of the small scale distribution of the Local Group's galaxy population. Such well defined planar structures are apparently rare occurrences in cold dark matter cosmological simulations, and we lack a coherent explanation of their formation and existence. In this paper, we explore the long-term survivability of thin planes of dwarfs in galactic halos, focusing, in particular, on systems mimicking the observed Andromeda distribution. The key results show that, in general, planes of dwarf galaxies are fragile, sensitive to the shape of the dark matter halo and other perturbing effects. In fact, long lived planes of satellites only exist in polar orbits in spherical dark matter halos, presenting a challenge to the observed Andromeda plane which is significantly tilted with respect to the optical disk. Our conclusion is that, in standard cosmological models, planes of satellites are generally short lived, and hence we must be located at a relatively special time in the evolution of the Andromeda Plane, lucky enough to see its coherent pattern.


title: Sequential Learning without Feedback

abstract: In many security and healthcare systems a sequence of features/sensors/tests are used for detection and diagnosis. Each test outputs a prediction of the latent state, and carries with it inherent costs. Our objective is to {\it learn} strategies for selecting tests to optimize accuracy \& costs. Unfortunately it is often impossible to acquire in-situ ground truth annotations and we are left with the problem of unsupervised sensor selection (USS). We pose USS as a version of stochastic partial monitoring problem with an {\it unusual} reward structure (even noisy annotations are unavailable). Unsurprisingly no learner can achieve sublinear regret without further assumptions. To this end we propose the notion of weak-dominance. This is a condition on the joint probability distribution of test outputs and latent state and says that whenever a test is accurate on an example, a later test in the sequence is likely to be accurate as well. We empirically verify that weak dominance holds on real datasets and prove that it is a maximal condition for achieving sublinear regret. We reduce USS to a special case of multi-armed bandit problem with side information and develop polynomial time algorithms that achieve sublinear regret.


title: Diffusive stability of spatially periodic patterns with a conservation law

abstract: Applying the Lyapunov-Schmidt reduction approach introduced by Mielke and Schneider in their analysis of the fourth-order scalar Swift-Hohenberg equation, we carry out a rigorous small-amplitude stability analysis of Turing patterns for the model introduced by Matthews and Cox for pattern formation with a conservation law. Our results confirm that stability is accurately predicted in the small-amplitude limit by the formal modified Ginzburg-Lanadau system (mGL) consisting of a coupled Ginzburg-Landau equation and mean mode equation derived by Matthews and Cox, rigorously validating the standard weakly unstable approximation.


title: Star Formation Law in the Milky Way

abstract: The Schmidt law (SF law) in the Milky Way was investigated using 3D distribution maps of HII regions, HI and molecular (\Htwo) gases with spatial resolutions of $\sim 1$ kpc in the Galactic plane and a few tens of pc in the vertical direction. HII regions were shown to be distributed in a star-forming (SF) disk with nearly constant vertical full thickness 92 pc in spatial coincidence with the molecular gas disk. The vertically averaged volume star formation rate (SFR) $\rho_{\rm SFR}$ in the SF disk is related to the surface SFR $\Sigma_{\rm SFR}$ by $\rho_{\rm SFR} /[{\rm M_\odot y^{-1} kpc^{-3}}] =9.26\times \Sigma_{\rm SFR}/[{\rm M_\odot y^{-1} kpc^{-2}}]$. The SF law fitted by a single power law of gas density in the form of $\Sigma_{\rm SFR} \propto \rho_{\rm SFR} \propto \rho_{\rm gas}^\alpha$ and $\propto \Sigma_{\rm gas}^\beta$ showed indices of $\alpha=0.78 \pm 0.05$ for $\rho_{\rm H_2}$ and $2.15 \pm 0.08$ for $\rho_{\rm total}$, and $\beta=1.14\pm 0.23$ for $\Sigma_{\rm total}$, where $\rho$ and $\Sigma$ denote volume and surface densities, respectively. The star formation rate is shown to be directly related to the molecular gas, but indirectly to HI and total gas densities. The dependence of the SF law on the gaseous phase is explained by the phase transition theory between HI and H$_2$ gases.


title: An application of continuous logic to fixed point theory

abstract: In aiming to apply to a broader class of examples the Avigad-Iovino "ultraproducts and metastability" approach to obtaining uniformity for convergence of sequences, we construct a framework using continuous logic that in particular is able to handle discontinuous functions in its domain of discourse. This setup weakens the usual continuity requirements for functions, but compensates for the loss of control by introducing a notion of "linear structure" that captures in a quite general way the situation of having geodesics between every pair of points, and has as a special case the vector space structure of Banach spaces. We use this to apply the Avigad-Iovino method to specific convergence results from functional analysis involving iterations of discontinuous functions, and so obtain uniform metastable convergence in those results.


title: Localization of weakly interacting bosons in two dimensions: disorder vs lattice geometry effects

abstract: We investigate the effects of disorder and lattice geometry against localisation phenomena in a weakly interacting ultracold bosonic gas confined in a 2D optical lattice. The behaviour of the quantum fluid is studied at the mean-field level performing computational experiments, as a function of disorder strength for lattices of sizes similar to current experiments. Quantification of localisation, away from the Bose glass phase, was obtained directly from the stationary density profiles through a robust statistical analysis of the condensate component, as a function of the disorder amplitude. Our results show a smooth transition, or crossover, to localisation induced by disorder in square and triangular lattices. In contrast, associated to its larger tunneling amplitude, honeycomb lattices show absence of localisation for the same range of disorder strengths and same lattice amplitude, while also exhibiting partial localisation for large disorder amplitudes. We also conclude that the coordination number z have a partial influence on how fast this smooth transition occurs as the system size increases. Signatures of disorder are also found in the ground state energy spectrum, where a continuous distribution emerges instead of a distribution of sharp peaks proper to the system in the absence of disorder.


title: Blue Phosphorene Oxide: Strain-tunable Quantum Phase Transitions and Novel 2D Emergent Fermions

abstract: Tunable quantum phase transitions and novel emergent fermions in solid state materials are fascinating subjects of research. Here, we propose a new stable two-dimensional (2D) material, the blue phosphorene oxide (BPO), which exhibits both. Based on first-principles calculations, we show that its equilibrium state is a narrow-bandgap semiconductor with three bands at low energy. Remarkably, a moderate strain can drive a semiconductor-to-semimetal quantum phase transition in BPO. At the critical transition point, the three bands cross at a single point at Fermi level, around which the quasiparticles are a novel type of 2D pseudospin-1 fermions. Going beyond the transition, the system becomes a symmetry-protected semimetal, for which the conduction and valence bands touch quadratically at a single Fermi point that is protected by symmetry, and the low-energy quasiparticles become another novel type of 2D double Weyl fermions. We construct effective models characterizing the phase transition and these novel emergent fermions, and we point out several exotic effects, including super Klein tunneling, supercollimation, and universal optical absorbance. Our result reveals BPO as an intriguing platform for the exploration of fundamental properties of quantum phase transitions and novel emergent fermions, and also suggests its great potential in nanoscale device applications.


title: Going off the Grid: Iterative Model Selection for Biclustered Matrix Completion

abstract: We consider the problem of performing matrix completion with side information on row-by-row and column-by-column similarities. We build upon recent proposals for matrix estimation with smoothness constraints with respect to row and column graphs. We present a novel iterative procedure for directly minimizing an information criterion in order to select an appropriate amount row and column smoothing, namely perform model selection. We also discuss how to exploit the special structure of the problem to scale up the estimation and model selection procedure via the Hutchinson estimator. We present simulation results and an application to predicting associations in imaging-genomics studies.


title: Uniquely solvable and energy stable decoupled schemes for Cahn-Hilliard-Stokes-Darcy system for two-phase flows in karstic geometry

abstract: We propose and analyze two novel decoupled numerical schemes for solving the Cahn-Hilliard-Stokes-Darcy (CHSD) model for two-phase flows in karstic geometry. In the first numerical scheme, we explore a fractional step method (operator splitting) to decouple the phase-field (Cahn-Hilliard equation) from the velocity field (Stokes-Darcy fluid equations). To further decouple the Stokes-Darcy system, we introduce a first order pressure stabilization term in the Darcy solver in the second numerical scheme so that the Stokes system is decoupled from the Darcy system and hence the CHSD system can be solved in a fully decoupled manner. We show that both decoupled numerical schemes are uniquely solvable, energy stable, and mass conservative. Ample numerical results are presented to demonstrate the accuracy and efficiency of our schemes.


title: VRPBench: A Vehicle Routing Benchmark Tool

abstract: The number of optimization techniques in the combinatorial domain is large and diversified. Nevertheless, there is still a lack of real benchmarks to validate optimization algorithms. In this work we introduce VRPBench, a tool to create instances and visualize solutions to the Vehicle Routing Problem (VRP) in a planar graph embedded in the Euclidean 2D space. We use VRPBench to model a real-world mail delivery case of the city of Artur Nogueira. Such scenarios were characterized as a multi-objective optimization of the VRP. We extracted a weighted graph from a digital map of the city to create a challenging benchmark for the VRP. Each instance models one generic day of mail delivery with hundreds to thousands of delivery points, thus allowing both the comparison and validation of optimization algorithms for routing problems.


title: Runaway gas accretion and gap opening versus type~I migration

abstract: Growing planets interact with their natal protoplanetary disc, which exerts a torque onto them allowing them to migrate in the disc. Small mass planets do not affect the gas profile and migrate in the fast type~I migration. Although type~I migration can be directed outwards for planets smaller than $20-30 M_\oplus$, planets above this mass should be lost into the central star long before the disc disperses. Massive planets open a gap and subsequently migrate in the slower, type~II migration, which could save them from migrating all the way to the star. Hence, growing giant planets can be saved if and only if they can reach the gap opening mass, because this extends their migration time-scale, allowing them to eventually survive at large orbits until the disc itself disperses. However, most of the previous studies only measured the torques on planets with fixed masses and orbits to determine the migration rate. Additionally, the transition between type~I and type~II migration itself is not well studied. Here we use isothermal 2D disc simulations with FARGO-2D1D to study the migration behaviour of gas accreting protoplanets in discs. We find that migrating giant planets always open gaps in the disc. We further show analytically and numerically that planets growing in the runaway gas accretion regime will reach gap opening masses before migrating all the way to the central star if the disc is not extremely viscous and/or thick. An accretion rate limited to the radial gas flow in the disc, in contrast, is not fast enough. When gas accretion by the planet is taken into account, the gap opening process is accelerated because the planet accretes material originating from its horseshoe region. This allows an accreting planet to transition to type~II migration before being lost even if gas fails to be provided for a rapid enough growth and the gap opening mass is not reached.


title: An Alternative Approach to Mean Field Game with Major and Minor Players, and Applications to Herders Impacts

abstract: The goal of the paper is to introduce a formulation of the mean field game with major and minor players as a fixed point on a space of controls. This approach emphasizes naturally the role played by McKean-Vlasov dynamics in some of the players optimization problems. We apply this approach to linear quadratic models for which we recover the existing solutions for open loop equilibria, and we show that we can also provide solutions for closed loop versions of the game. Finally, we implement numerically our theoretical results on a simple model of flocking.


title: Homomorphism reductions on Polish groups

abstract: In an earlier paper, we introduced the following pre-order on the subgroups of a given Polish group: if $G$ is a Polish group and $H,L \subseteq G$ are subgroups, we say $H$ is {\em homomorphism reducible} to $L$ iff there is a continuous group homomorphism $\varphi : G \rightarrow G$ such that $H = \varphi^{-1} (L)$. We previously showed that there is a $K_\sigma$ subgroup, $L$, of the countable power of any locally compact Polish group, $G$, such that every $K_\sigma$ subgroup of $G^\omega$ is homomorphism reducible to $L$. In the present work, we show that this fails in the countable power of the group of increasing homeomorphisms of the unit interval.


title: Strong edge-colorings of sparse graphs with large maximum degree

abstract: A {\em strong $k$-edge-coloring} of a graph $G$ is a mapping from $E(G)$ to $\{1,2,\ldots,k\}$ such that every two adjacent edges or two edges adjacent to the same edge receive distinct colors. The {\em strong chromatic index} $\chi_s'(G)$ of a graph $G$ is the smallest integer $k$ such that $G$ admits a strong $k$-edge-coloring. We give bounds on $\chi_s'(G)$ in terms of the maximum degree $\Delta(G)$ of a graph $G$. when $G$ is sparse, namely, when $G$ is $2$-degenerate or when the maximum average degree ${\rm Mad}(G)$ is small. We prove that the strong chromatic index of each $2$-degenerate graph $G$ is at most $5\Delta(G) +1$. Furthermore, we show that for a graph $G$, if ${\rm Mad}(G)< 8/3$ and $\Delta(G)\geq 9$, then $\chi_s'(G)\leq 3\Delta(G) -3$ (the bound $3\Delta(G) -3$ is sharp) and if ${\rm Mad}(G)<3$ and $\Delta(G)\geq 7$, then $\chi_s'(G)\leq 3\Delta(G)$ (the restriction ${\rm Mad}(G)<3$ is sharp).


title: Contribution of the first galaxies to the cosmic far-infrared/sub-millimeter background - I. Mean background level

abstract: We study the contribution of the first galaxies to the far-infrared/sub-millimeter (FIR/sub-mm) extragalactic background light (EBL) by implementing an analytical model for dust emission. We explore different dust models, assuming different grain size distributions and chemical compositions. According to our findings, observed re-radiated emission from dust in dwarf-size galaxies at $z \sim 10$ would peak at a wavelength of $\sim 500 \mu {\rm m}$ with observed fluxes of $\sim 10^{-3} - 10^{-2}$ nJy, which is below the capabilities of current observatories. In order to be detectable, model sources at these high redshifts should exhibit luminosities of $\gtrsim 10^{12} L_{\odot}$, comparable to that of local ultra-luminous systems. The FIR/sub-mm EBL generated by primeval galaxies peaks at $\sim 500 \mu {\rm m}$, with an intensity ranging from $\sim 10^{-4}$ to $10^{-3} {\rm nW \ m^{-2} \ sr^{-1}}$, depending on dust properties. These values are $\sim 3 - 4$ orders of magnitude below the absolute measured cosmic background level, suggesting that the first galaxies would not contribute significantly to the observed FIR/sub-mm EBL. Our model EBL exhibits a strong correlation with the dust-to-metal ratio, where we assume a fiducial value of $D = 0.005$, increasing almost proportionally to it. Thus, measurements of the FIR/sub-mm EBL could provide constraints on the amount of dust in the early Universe. Even if the absolute signal from primeval dust emission may be undetectable, it might still be possible to obtain information about it by exploring angular fluctuations at $\sim 500 \mu {\rm m}$, close to the peak of dust emission from the first galaxies.


title: Finite State Mean Field Games with Major and Minor Players

abstract: The goal of the paper is to develop the theory of finite state mean field games with major and minor players when the state space of the game is finite. We introduce the finite player games and derive a mean field game formulation in the limit when the number of minor players tends to infinity. In this limit, we prove that the value functions of the optimization problems are viscosity solutions of PIDEs of the HJB type, and we construct the best responses for both types of players. From there, we prove existence of Nash equilibria under reasonable assumptions. Finally we prove that a form of propagation of chaos holds in the present context and use this result to prove existence of approximate Nash equilibria for the finite player games from the solutions of the mean field games. this vindicate our formulation of the mean field game problem.


title: Split Nash Equilibria for Related Games and Applications to Economics

abstract: In this paper we introduce the concept of split Nash equilibrium problems associated with two related noncooperative strategic games. Then we apply the Fan-KKM theorem to prove the existence of solutions to split Nash equilibrium problems of related noncooperative strategic games, in which the strategy sets of the players are nonempty closed and convex subsets in Banach spaces. As application of this existence to economics, an example is provided to study the existence of split Nash equilibrium of utilities of two related economies. As applications, we study the existence of split Nash equilibrium in the dual extended Bertrant duopoly model of price competition


title: Two-scale analysis of solar magnetic helicity

abstract: We develop a two-scale formalism to determine global magnetic helicity spectra in systems where the local magnetic helicity has opposite signs on both sides of the equator, giving rise to cancelation with conventional methods. We verify this approach using first a synthetic one-dimensional magnetic field and then two-dimensional slices from a three-dimensional alpha effect-type dynamo-generated magnetic field with forced turbulence of opposite helicity above and below the midplane of the domain. We then apply this formalism to global solar synoptic vector magnetograms. To improve the statistics, data from three consecutive Carrington rotations (2161--2163) are combined into a single map. We find that the spectral magnetic helicity representative of the northern hemisphere is negative at all wavenumbers and peaks at ~ 0.06 Mm^{-1} (scales around 100 Mm). There is no evidence of bihelical magnetic fields that are found in three-dimensional turbulence simulations of helicity-driven alpha effect-type dynamos.


title: Gate-independent energy gap in non-covalently intercalated bilayer graphene on SiC(0001)

abstract: Our first-principles calculations show that an energy gap around 0.12-0.25 eV can be engineered in epitaxial graphene on SiC(0001) through the non-covalent intercalation of transition- or alkali-metals, yet originated from the distinct mechanisms. The former is attributed to the combined effects of metal induced perpendicular electric field and interaction while the latter is solely attributed to the built-in electric field. A great advantage of this scheme is that the gap size is almost independent of the gate voltage up to 1 V/nm, thus reserving the electric means to tune the Fermi level of graphene when configured as field-effect transistors. Given the recent progress in experimental techniques on the intercalated graphene, our findings provide a practical way to incorporate the graphene with the current semiconductor industry.


title: Almost-Sure Reachability in Stochastic Multi-Mode System

abstract: A constant-rate multi-mode system is a hybrid system that can switch freely among a finite set of modes, and whose dynamics is specified by a finite number of real-valued variables with mode-dependent constant rates. We introduce and study a stochastic extension of a constant-rate multi-mode system where the dynamics is specified by mode-dependent compactly supported probability distributions over a set of constant rate vectors. Given a tolerance $\varepsilon > 0$, the almost-sure reachability problem for stochastic multi-mode systems is to decide the existence of a control strategy that steers the system almost-surely from an arbitrary start state to an $\varepsilon$-neighborhood of an arbitrary target state while staying inside a pre-specified safety set. We prove a necessary and sufficient condition to decide almost-sure reachability and, using this condition, we show that almost-sure reachability can be decided in polynomial time. Our algorithm can be used as a path-following algorithm in combination with any off-the-shelf path-planning algorithm to make a robot or an autonomous vehicle with noisy low-level controllers follow a given path with arbitrary precision.


title: Relative quantum coherence, incompatibility, and quantum correlations of states

abstract: Quantum coherence, incompatibility, and quantum correlations are fundamental features of quantum physics. A unified view of those features is crucial for revealing quantitatively their intrinsic connections. We define the relative quantum coherence of two states as the coherence of one state in the reference basis spanned by the eigenvectors of another one and establish its quantitative connections with the extent of mutual incompatibility of two states. We also show that the proposed relative quantum coherence, which can take any form of measures such as $l_1$ norm and relative entropy, can be interpreted as or connected to various quantum correlations such as quantum discord, symmetric discord, entanglement of formation, and quantum deficits. Our results reveal conceptual implications and basic connections of quantum coherence, mutual incompatibility, and quantum correlations.


title: The rigidity of hypersurface in Euclidean space

abstract: In the present paper, we revisit the rigidity of hypersurfaces in Euclidean space. We highlight Darboux equation and give new proof of rigidity of hypersurfaces by energy method and maximal principle.


title: Weak Convergence (IA). Sequences of Random Vectors

abstract: (English) This monograph aims at presenting the core weak convergence theory for sequences of random vectors with values in $\mathbb{R}^k$. In some places, a more general formulation in metric spaces is provided. It lays out the necessary foundation that paves the way to applications in particular subfields of the theory. In particular, the needs of Asymptotic Statistics are addressed. A whole chapter is devoted to weak convergence in $\mathbb{R}$ where specific tools, for example for handling weak convergence of sequences using independent and indentically distributed random variables such that the Renyi's representations by means of standard uniform or exponential random variables, are stated. The function empirical process is presented as a powerful tool for solving a considerable number of asymptotic problems in Statistics. The text is written in a self-contained approach whith the proofs of all used results at the exception of the general Skorohod-Wichura Theorem.


title: Refined Shapley-Folkman Lemma and Its Application in Duality Gap Estimation

abstract: Based on concepts like kth convex hull and finer characterization of nonconvexity of a function, we propose a refinement of the Shapley-Folkman lemma and derive a new estimate for the duality gap of nonconvex optimization problems with separable objective functions. We apply our result to a network flow problem and the dynamic spectrum management problem in communication systems as examples to demonstrate that the new bound can be qualitatively tighter than the existing ones. The idea is also applicable to cases with general nonconvex constraints.


title: Numeric Solution of Advection-Diffusion Equations by a Discrete Time Random Walk Scheme

abstract: Explicit numerical finite difference schemes for partial differential equations are well known to be easy to implement but they are particularly problematic for solving equations whose solutions admit shocks, blowups and discontinuities. Here we present an explicit numerical scheme for solving non-linear advection-diffusion equations admitting shock solutions that is both easy to implement and stable. The numerical scheme is obtained by considering the continuum limit of a discrete time and space stochastic process for non-linear advection-diffusion. The stochastic process is well posed and this guarantees the stability of the scheme. Several examples are provided to highlight the importance of the formulation of the stochastic process in obtaining a stable and accurate numerical scheme.


title: Analysis of the semileptonic decay \Lambda_c->ne^+\nu_e

abstract: The semileptonic weak decay process of the $\Lambda_c$ baryon to the neutron $\Lambda_c\rightarrow ne^+\nu_e$ is examined. The transition form factors are investigated with light-cone QCD sum rules. The differential decay width is obtained in the dynamical region by fitting the sum rules-allowed results with the dipole formula. The total decay width and the branching ratio are estimated to be $\Gamma(\Lambda_c\rightarrow ne^+\nu_e)=(8.57\pm0.41)\times10^{-15}\,\mbox{GeV}$ and $\mbox{Br}(\Lambda_c\rightarrow ne^+\nu_e)=0.26\pm0.01\%$, respectively.


title: A Joint Indoor WLAN Localization and Outlier Detection Scheme Using LASSO and Elastic-Net Optimization Techniques

abstract: In this paper, we introduce two indoor Wireless Local Area Network (WLAN) positioning methods using augmented sparse recovery algorithms. These schemes render a sparse user's position vector, and in parallel, minimize the distance between the online measurement and radio map. The overall localization scheme for both methods consists of three steps: 1) coarse localization, obtained from comparing the online measurements with clustered radio map. A novel graph-based method is proposed to cluster the offline fingerprints. In the online phase, a Region Of Interest (ROI) is selected within which we search for the user's location; 2) Access Point (AP) selection; and 3) fine localization through the novel sparse recovery algorithms. Since the online measurements are subject to inordinate measurement readings, called outliers, the sparse recovery methods are modified in order to jointly estimate the outliers and user's position vector. The outlier detection procedure identifies the APs whose readings are either not available or erroneous. The proposed localization methods have been tested with Received Signal Strength (RSS) measurements in a typical office environment and the results show that they can localize the user with significantly high accuracy and resolution which is superior to the results from competing WLAN fingerprinting localization methods.


title: Collective cell migration without proliferation: density determines cell velocity and wave velocity

abstract: Collective cell migration contributes to embryogenesis, wound healing and tumor metastasis. Cell monolayer migration experiments help understanding what determines the movement of cells far from the leading edge. Inhibiting cell proliferation limits cell density increase and prevents jamming; we observe long-duration migration and quantify space-time characteristics of the velocity profile over large length- and time-scales. Velocity waves propagate backwards and their frequency depends only on cell density at the moving front. Both cell average velocity and wave velocity increase linearly with the cell effective radius regardless of the distance to the front. Inhibiting lamellipodia decreases cell velocity while waves either disappear or have a lower frequency. Our model combines conservation laws, monolayer mechanical properties and a phenomenological coupling between strain and polarity: advancing cells pull on their followers which then become polarized. With reasonable values of parameters, this model agrees with several of our experimental observations. Together, our experiments and model disantangle the respective contributions of active velocity and of proliferation in monolayer migration, explain how cells maintain their polarity far from the moving front, and highlight the importance of strain-polarity coupling and density in long-range information propagation.


title: Structured Group Sparsity: A Novel Indoor WLAN Localization, Outlier Detection, and Radio Map Interpolation Scheme

abstract: This paper introduces novel schemes for indoor localization, outlier detection, and radio map interpolation using Wireless Local Area Networks (WLANs). The localization method consists of a novel multicomponent optimization technique that minimizes the squared $\ell_{2}$-norm of the residuals between the radio map and the online Received Signal Strength (RSS) measurements, the $\ell_{1}$-norm of the user's location vector, and weighted $\ell_{2}$-norms of layered groups of Reference Points (RPs). RPs are grouped using a new criterion based on the similarity between the so-called Access Point (AP) coverage vectors. In addition, since AP readings are prone to containing inordinate readings, called outliers, an augmented optimization problem is proposed to detect the outliers and localize the user with cleaned online measurements. Moreover, a novel scheme to record fingerprints from a smaller number of RPs and estimate the radio map at RPs without recorded fingerprints is developed using sparse recovery techniques. All localization schemes are tested on RSS fingerprints collected from a real environment. The overall scheme has comparable complexity with competing approaches, while it performs with high accuracy under a small number of APs and finer granularity of RPs.


title: On subgroups of first homology

abstract: We prove several new theorems regarding first homology. Some dichotomies for first homology of Peano continua are presented, as well as a notion of strong abelianization for arbitrary path connected metric spaces. We also show that the fundamental group of the Hawaiian earring has Borel subgroups of almost all multiplicative and additive types.