title: The strong rotation of M5 (NGC 5904) as seen from the MIKiS Survey of Galactic Globular Clusters

abstract: In the context of the ESO-VLT Multi-Instrument Kinematic Survey (MIKiS) of Galactic globular clusters, we present the line-of-sight rotation curve and velocity dispersion profile of M5 (NGC 5904), as determined from the radial velocity of more than 800 individual stars observed out to 700" (~ 5 half-mass radii) from the center. We find one of the cleanest and most coherent rotation patterns ever observed for globular clusters, with a very stable rotation axis (having constant position angle of 145^o at all surveyed radii) and a well-defined rotation curve. The density distribution turns out to be flattened in the direction perpendicular to the rotation axis, with a maximum ellipticity of 0.15. The rotation velocity peak (~3 km/s in projection) is observed at ~0.6 half-mass radii, and its ratio with respect to the central velocity dispersion (~0.3-0.4 at 4 projected half-mass radii) indicates that ordered motions play a significant dynamical role. This result strengthens the growing empirical evidence of the kinematic complexity of Galactic globular clusters and motivates the need of fundamental investigations of the role of angular momentum in collisional stellar dynamics.


title: Tidal breakup of quadruple stars in the Galactic Centre

abstract: The most likely origin of hypervelocity stars (HVSs) is the tidal disruption of a binary star by the supermassive black hole (MBH) in the Galactic Centre (GC). However, HE0437-5439, a $9$ M$_\odot$ B-type main-sequence star moving with a heliocentric radial velocity of about $720$ km s$^{-1}$ at a distance of $\sim 60$ kpc, and the recent discovered hypervelocity binary candidate (HVB), traveling at $\sim 570$ km s$^{-1}$, challenge this standard scenario. Recently, Fragione & Gualandris (2018) have demonstrated that the tidal breakup of a triple star leads to an insufficient rate. Observations show that quadruple stars made up of two binaries orbiting their common center of mass (the so-called 2+2 quadruples) are $\approx 4\%$ of the stars in the solar neighborhood. Although rarer than triples, 2+2 quadruple stars may have a role in ejecting HVBs as due to their larger energy reservoir. We present a numerical study of 2+2 quadruple disruptions by the MBH in the GC and find that the production of HVBs has a probability $\lesssim 2-4\%$, which translates into an ejection rate of $\lesssim 1$ Gyr$^{-1}$, comparable to the triple disruption scenario. Given the low ejection rate, we suggest that alternative mechanisms are responsible for the origin of HVBs, as the ejection from the interaction of a young star cluster with the MBH in the GC and the origin in the Large Magellanic Cloud.


title: Complex saddle trajectories for multidimensional quantum wave packet/coherent state propagation: application to a many-body system

abstract: A practical search technique for finding the complex saddle points used in wave packet or coherent state propagation is developed which works for a large class of Hamiltonian dynamical systems with many degrees of freedom. The method can be applied to problems in atomic, molecular, and optical physics, and other domains. A Bose-Hubbard model is used to illustrate the application to a many-body system where discrete symmetries play an important and fascinating role. For multidimensional wave packet propagation, locating the necessary saddles involves the seemingly insurmountable difficulty of solving a boundary value problem in a high-dimensional complex space, followed by determining whether each particular saddle found actually contributes. In principle, this must be done for each propagation time considered. The method derived here identifies a real search space of minimal dimension, which leads to a complete set of contributing saddles up to intermediate times much longer than the Ehrenfest time scale for the system. The analysis also gives a powerful tool for rapidly identifying the various dynamical regimes of the system.


title: Probabilistic Verification for Obviously Strategyproof Mechanisms

abstract: Obviously strategyproof (OSP) mechanisms maintain the incentive compatibility of agents that are not fully rational. They have been object of a number of studies since their recent definition. A research agenda, initiated in [Ferraioli&Ventre, AAAI 2017], is to find a small (possibly, the smallest) set of conditions allowing to implement an OSP mechanism. To this aim, we define a model of probabilistic verification wherein agents are caught misbehaving with a certain probability, and show how OSP mechanisms can implement every social choice function at the cost of either imposing very large fines for lies or verifying a linear number of agents.


title: 2-Parts of real class sizes

abstract: We investigate the structure of finite groups whose non-central real class sizes have the same $2$-part. In particular, we prove that such groups are solvable and have $2$-length one. As a consequence, we show that a finite group is solvable if it has two real class sizes. This confirms a conjecture due to G. Navarro, L. Sanus and P. Tiep.


title: The Markov-quantile process attached to a family of Marginals

abstract: Let $\mu$ = ($\mu$t)t$\in$R be any 1-parameter family of probability measures on R. Its quantile process (Gt)t$\in$R : ]0, 1[ $\rightarrow$ RR, given by Gt($\alpha$) = inf{x $\in$ R : $\mu$t(]--$\infty$, x]) > $\alpha$}, is not Markov in general. We modify it to build the Markov process we call "Markov-quantile".We first describe the discrete analogue: if ($\mu$n)n$\in$Z is a family of probability measures on R, a Markov process Y = (Yn)n$\in$Z such that Law(Yn) = $\mu$n is given by the data of its couplings from n to n + 1, i.e. Law((Yn, Yn+1)), and the process Y is the inhomogeneous Markov chain having those couplings as transitions. Therefore, there is a canonical Markov process with marginals $\mu$n and as similar as possible to the quantile process: the chain whose transitions are the quantile couplings. We show that an analogous process exists for a continuous parameter t: there is a unique Markov process X with the measures $\mu$t as marginals, and being a limit for the finite dimensional topology of quantile processes where the past is made independent of the future at finitely many times (many non-Markovian limits exist in general). The striking fact is that the construction requires no regularity for the family $\mu$. We rely on order arguments, which seems to be completely new for the purpose.We also prove new results the Markov-quantile process yields in two contemporary frameworks:-- In case $\mu$ is increasing for the stochastic order, X has increasing trajectories. This is an analogue of a result of Kellerer dealing with the convex order, peacocks and martingales. Modifiying Kellerer's proof, we also prove simultaneously his result and ours in this case.-- If $\mu$ is absolutely continuous in Wasserstein space P2(R) then X is solution of a Benamou--Brenier transport problem with marginals $\mu$t. Itprovides a Markov probabilistic representation of the continuity equation, unique in a certain sense.


title: Multipoint Cauchy problem for Schr\"odinger type equations with general elliptic part

abstract: In this paper, the existence, uniqueness and regularity properties, Strichartz type estimates for solution of multipoint Cauchy problem for linear and nonlinear Schr\"odinger equations with general elliptic leading part is obtained.


title: Rate-Splitting Multiple Access for Coordinated Multi-Point Joint Transmission

abstract: As a promising downlink multiple access scheme, Rate-Splitting Multiple Access (RSMA) has been shown to achieve superior spectral and energy efficiencies compared with Space-Division Multiple Access (SDMA) and Non-Orthogonal Multiple Access (NOMA) in downlink single-cell systems. By relying on linearly precoded rate-splitting at the transmitter and successive interference cancellation at the receivers, RSMA has the capability of partially decoding the interference and partially treating the interference as noise, and therefore copes with a wide range of user deployments and network loads. In this work, we further study RSMA in downlink Coordinated Multi-Point (CoMP) Joint Transmission (JT) networks by investigating the optimal beamformer design to maximize the Weighted Sum-Rate (WSR) of all users subject to individual Quality of Service (QoS) rate constraints and per base station power constraints. Numerical results show that, in CoMP JT, RSMA achieves significant WSR improvement over SDMA and NOMA in a wide range of inter-user and inter-cell channel strength disparities. Specifically, SDMA (resp. NOMA) is more suited to deployments with little (resp. large) inter-user channel strength disparity and large (resp. little) inter-cell channel disparity, while RSMA is suited to any deployment. We conclude that RSMA provides rate, robustness and QoS enhancements over SDMA and NOMA in CoMP JT networks.


title: Fourth order nematic elasticity and modulated nematic phases: a poor man's approach

abstract: We propose an extension of Frank-Oseen's elastic energy for bulk nematic liquid crystals which is based on the hypothesis that the fundamental deformations allowed in nematic liquid crystals are splay, twist and bend. The extended elastic energy is a fourth order form in the fundamental deformations. The existence of bulk spontaneous modulated or deformed nematic liquid crystal ground states is investigated. The analysis is limited to bulk nematic liquid crystals in the absence of limiting surfaces and/or external fields. The non deformed ground state is stable only when Frank-Oseen's elastic constants are positive. In case where at least one of them is negative, the ground state becomes deformed. The analysis of the stability of the deformed states in the space of the elastic parameters allows to characterize different types of deformed nematic phases. Some of them are new nematic phases, for instance a twist -- splay nematic phase is predicted. Inequalities between second order elastic constants which govern the stability of the twist--bend and splay--bend state are obtained.


title: Biology is a constructive physics

abstract: Yuri Manin's approach to Zipf's law (Kolmogorov complexity as energy) is applied to investigation of biological evolution. Model of constructive statistical mechanics where complexity is a contribution to energy is proposed to model genomics. Scaling laws in genomics are discussed in relation to Zipf's law. This gives a model of Eugene Koonin's Third Evolutionary Synthesis (physical model which should describe scaling in genomics).

categories: q-bio.PE cond-mat.stat-mech cs.IT math.IT


title: Energy-dependent topological anti-de Sitter black holes in Gauss-Bonnet Born-Infeld gravity

abstract: Employing higher curvature corrections to Einstein--Maxwell gravity has garnered a great deal of attention motivated by the high energy regime in quantum nature of black hole physics. In addition, one may employ gravity's rainbow to encode quantum gravity effects into the black hole solutions. In this paper, we regard an energy dependent static spacetime with various topologies and study its black hole solutions in the context of Gauss--Bonnet Born--Infeld (GB--BI) gravity. We study thermodynamic properties and examine the first law of thermodynamics. Using suitable local transformation, we endow the Ricci--flat black hole solutions with a global rotation and study the effects of rotation on thermodynamic quantities. We also investigate thermal stability in canonical ensemble through calculating the heat capacity. We obtain the effects of various parameters on the horizon radius of stable black holes. Finally, we discuss second order phase transition in the extended phase space thermodynamics and investigate the critical behavior.


title: Using Machine Learning to Improve Cylindrical Algebraic Decomposition

abstract: Cylindrical Algebraic Decomposition (CAD) is a key tool in computational algebraic geometry, best known as a procedure to enable Quantifier Elimination over real-closed fields. However, it has a worst case complexity doubly exponential in the size of the input, which is often encountered in practice. It has been observed that for many problems a change in algorithm settings or problem formulation can cause huge differences in runtime costs, changing problem instances from intractable to easy. A number of heuristics have been developed to help with such choices, but the complicated nature of the geometric relationships involved means these are imperfect and can sometimes make poor choices. We investigate the use of machine learning (specifically support vector machines) to make such choices instead. Machine learning is the process of fitting a computer model to a complex function based on properties learned from measured data. In this paper we apply it in two case studies: the first to select between heuristics for choosing a CAD variable ordering; the second to identify when a CAD problem instance would benefit from Groebner Basis preconditioning. These appear to be the first such applications of machine learning to Symbolic Computation. We demonstrate in both cases that the machine learned choice outperforms human developed heuristics.


title: An Integrative Introduction to Human Augmentation Science

abstract: Human Augmentation (HA) spans several technical fields and methodological approaches, including Experimental Psychology, Human-Computer Interaction, Psychophysiology, and Artificial Intelligence. Augmentation involves various strategies for optimizing and controlling cognitive states, which requires an understanding of biological plasticity, dynamic cognitive processes, and models of adaptive systems. As an instructive lesson, we will explore a few HA-related concepts and outstanding issues. Next, we focus on inducing and controlling HA using experimental methods by introducing three techniques for HA implementation: learning augmentation, augmentation using physical media, and extended phenotype modeling. To conclude, we will review integrative approaches to augmentation, which transcend specific functions.


title: Axion mass and quantum information

abstract: I present a novel way of linking the shift in the axion mass due to non-perturbative gravitational effects and the quantum information interpretation of the holographic principle. The main result of this article is a novel way of expressing the shift in the axion mass by means of properties of the tensor network used in the construction of the quantum information interpretation of the holographic principle. Indeed it appears that the axion mass shift can be seen as a geometric property of the holographic bulk entanglement. Moreover, the explicit breaking of the continuous Peccei-Quinn symmetry by gravitational instantons can be regarded as a result of including tensors encoding gauge fields on the holographic bulk tensor network.


title: On the principle of linearized stability in interpolation spaces for quasilinear evolution equations

abstract: We give a proof for the asymptotic exponential stability of equilibria of quasilinear parabolic evolution equations in admissible interpolation spaces.


title: Non-convex regularization of bilinear and quadratic inverse problems by tensorial lifting

abstract: Considering the question: how non-linear may a non-linear operator be in order to extend the linear regularization theory, we introduce the class of dilinear mappings, which covers linear, bilinear, and quadratic operators between Banach spaces. The corresponding dilinear inverse problems cover blind deconvolution, deautoconvolution, parallel imaging in MRI, and the phase retrieval problem. Based on the universal property of the tensor product, the central idea is here to lift the non-linear mappings to linear representatives on a suitable topological tensor space. At the same time, we extend the class of usually convex regularization functionals to the class of diconvex functionals, which are likewise defined by a tensorial lifting. Generalizing the concepts of subgradients and Bregman distances from convex analysis to the new framework, we analyse the novel class of dilinear inverse problems and establish convergence rates under similar conditions than in the linear setting. Considering the deautoconvolution problem as specific application, we derive satisfiable source conditions and validate the theoretical convergence rates numerically.


title: Bootstrapping holographic warped CFTs or: how I learned to stop worrying and tolerate negative norms

abstract: We use modular invariance to derive constraints on the spectrum of warped conformal field theories (WCFTs) --- nonrelativistic quantum field theories described by a chiral Virasoro and $U(1)$ Kac-Moody algebra. We focus on holographic WCFTs and interpret our results in the simplest holographic set up: three dimensional gravity with Compere-Song-Strominger boundary conditions. Holographic WCFTs feature a negative $U(1)$ level that is responsible for negative norm descendant states. Despite the violation of unitarity we show that the modular bootstrap is still viable provided the (Virasoro-Kac-Moody) primaries carry positive norm. In particular, we show that holographic WCFTs must feature either primary states with negative norm or states with imaginary $U(1)$ charge, the latter of which have a natural holographic interpretation. For large central charge and arbitrary level, we show that the first excited primary state in any WCFT satisfies the Hellerman bound. Moreover, when the level is positive we point out that known bounds for CFTs with internal $U(1)$ symmetries readily apply to unitary WCFTs.


title: A Strong Stability Preserving Analysis for Multistage Two-Derivative Time-Stepping Schemes Based on Taylor Series Conditions

abstract: High order strong stability preserving (SSP) time discretizations are often needed to ensure the nonlinear (and sometimes non-inner-product) strong stability properties of spatial discretizations specially designed for the solution of hyperbolic PDEs. Multiderivative time-stepping methods have recently been increasingly used for evolving hyperbolic PDEs, and the strong stability properties of these methods are of interest. In our prior work we explored time discretizations that preserve the strong stability properties of spatial discretizations coupled with forward Euler and a second derivative formulation. However, many spatial discretizations do not satisfy strong stability properties when coupled with this second derivative formulation, but rather with a more natural Taylor series formulation. In this work we demonstrate sufficient conditions for an explicit two-derivative multistage method to preserve the strong stability properties of spatial discretizations in a forward Euler and Taylor series formulation. We call these strong stability preserving Taylor series (SSP-TS) methods. We also prove that the maximal order of SSP-TS methods is p = 6, and define an optimization procedure that allows us to find such SSP methods. Several types of these methods are presented and their efficiency compared. Finally, these methods are tested on several PDEs to demonstrate the benefit of SSP-TS methods, the need for the SSP property, and the sharpness of the SSP time-step in many cases.


title: Detecting and modeling worst-case dependence structures between random inputs of computational reliability models

abstract: Uncertain information on input parameters of reliability models is usually modeled by considering these parameters as random, and described by marginal distributions and a dependence structure of these variables. In numerous real-world applications, while information is mainly provided by marginal distributions, typically from samples , little is really known on the dependence structure itself. Faced with this problem of incomplete or missing information, risk studies are often conducted by considering independence of input variables, at the risk of including irrelevant situations. This approach is especially used when reliability functions are considered as black-box computational models. Such analyses remain weakened in absence of in-depth model exploration, at the possible price of a strong risk misestimation. Considering the frequent case where the reliability output is a quantile, this article provides a methodology to improve risk assessment, by exploring a set of pessimistic dependencies using a copula-based strategy. In dimension greater than two, a greedy algorithm is provided to build input regular vine copulas reaching a minimum quantile to which a reliability admissible limit value can be compared, by selecting pairwise components of sensitive influence on the result. The strategy is tested over toy models and a real industrial case-study. The results highlight that current approaches can provide non-conservative results, and that a nontrivial dependence structure can be exhibited to define a worst-case scenario.


title: OH$^+$ emission from cometary knots in planetary nebulae

abstract: We model the molecular emission from cometary knots in planetary nebulae (PNe) using a combination of photoionization and photodissociation region (PDR) codes, for a range of central star properties and gas densities. Without the inclusion of ionizing extreme ultraviolet (EUV) radiation, our models require central star temperatures $T_*$ to be near the upper limit of the range investigated in order to match observed H$_2$ and OH$^+$ surface brightnesses consistent with observations - with the addition of EUV flux, our models reproduce observed OH$^+$ surface brightnesses for $T_* \ge 100 \, {\rm kK}$. For $T_* < 80 \, {\rm kK}$, the predicted OH$^+$ surface brightness is much lower, consistent with the non-detection of this molecule in PNe with such central star temperatures. Our predicted level of H$_2$ emission is somewhat weaker than commonly observed in PNe, which may be resolved by the inclusion of shock heating or fluorescence due to UV photons. Some of our models also predict ArH$^+$ and HeH$^+$ rotational line emission above detection thresholds, despite neither molecule having been detected in PNe, although the inclusion of photodissociation by EUV photons, which is neglected by our models, would be expected to reduce their detectability.


title: The mre11 A470 alleles influence the hereditability and the segregation of telosomes in Saccharomyces cerevisiae (original with corrected raw data)

abstract: Telomeres, the nucleoprotein complexes at the termini of linear chromosomes, are essential for the processes of end replication, end protection, and chromatin segregation. The Mre11 complex is involved in multiple cellular roles in DNA repair and structure in the regulation and function of telomere size homeostasis. In this study, we characterize yeast telomere chromatin structure, phenotypic heritability, and chromatin segregation in both wild-type [MRE11] and A470 motif alleles. MRE11 strains confer a telomere size of 300 base pairs of G+T irregular simple sequence repeats. This DNA and a portion of subtelomeric DNA is embedded in a telosome: a MNase-resistant non-nucleosomal particle. Chromatin immunoprecipitation shows a three to four-fold lower occupancy of Mre11A470T proteins than wild-type proteins in telosomes. Telosomes containing the Mre11A470T protein confer a greater resistance to MNase digestion than wild-type telosomes. The integration of a wild-type MRE11 allele into an ectopic locus in the genome of an mre11A470T mutant and the introduction of an mre11A470T allele at an ectopic site in a wild-type strain lead to unexpectedly differing results. In each case, the replicated sister chromatids inherit telosomes containing only the protein encoded by the genomic mre11 locus, even in the presence of protein encoded by the opposing ectopic allele. We hypothesize that the telosome segregates by a conservative mechanism. These data support a mechanism for the linkage between sister chromatid replication and maintenance of either identical mutant or identical wild-type telosomes after replication of sister chromatids. These data suggest the presence of an active mechanism for chromatin segregation in yeast.


title: Abstract Mining

abstract: We have developed an application that will take a "MEDLINE" output from the PubMed database and allows the user to cluster all non-trivial words of the abstracts of the PubMed output. The number of clusters to use can be selected by the user. A specific cluster may be selected, and the PMIDs and dates for all publications in the selected cluster are displayed underneath. See figure 2, where cluster 12 is selected. The application also has an "Abstracts" tab, where the abstracts for the selected cluster can be perused. Here, it is also possible to download a HTML file containing the PMID, date, title, and abstract for each publication in the selected cluster. A third tab is called "Titles", where all the titles for the selected cluster are displayed. Via a "Use Cluster" button, the selected Cluster can itself be clustered. A "Back" button allows the user to return to any previous state. Finally, it is also possible to exclude documents whose abstracts contain certain words (see figure 3). The application will allow researchers to enter general search terms in the PubMed search engine, then use the application to search for publications of special interest within those search terms.


title: Generative Steganography by Sampling

abstract: In this paper, a novel data-driven information hiding scheme called generative steganography by sampling (GSS) is proposed. Unlike in traditional modification-based steganography, in our method the stego image is directly sampled by a powerful generator: no explicit cover is used. Both parties share a secret key used for message embedding and extraction. The Jensen-Shannon divergence is introduced as a new criterion for evaluating the security of generative steganography. Based on these principles, we propose a simple practical generative steganography method that uses semantic image inpainting. The message is written in advance to an uncorrupted region that needs to be retained in the corrupted image. Then, the corrupted image with the secret message is fed into a Generator trained by a generative adversarial network (GAN) for semantic completion. Message loss and prior loss terms are proposed for penalizing message extraction error and unrealistic stego image. In our design, we first train a generator whose training target is the generation of new data samples from the same distribution as that of existing training data. Next, for the trained generator, backpropagation to the message and prior loss are introduced to optimize the coding of the input noise data for the generator. The presented experiments demonstrate the potential of the proposed framework based on both qualitative and quantitative evaluations of the generated stego images.


title: A family of monomial ideals with the persistence property

abstract: In this paper we introduce a family of monomial ideals with the persistence property. Given positive integers $n$ and $t$, we consider the monomial ideal $I=Ind_t(P_n)$ generated by all monomials $\textbf{x} ^F$, where $F$ is an independent set of vertices of the path graph $P_n$ of size $t$, which is indeed the facet ideal of the $t$-th skeleton of the independence complex of $P_n$. We describe the set of associated primes of all powers of $I$ explicitly. It turns out that any such ideal $I$ has the persistence property. Moreover the index of stability of $I$ and the stable set of associated prime ideals of $I$ are determined.


title: Plasmon-enhanced spin-orbit interaction of light in graphene

abstract: We develop a novel theoretical framework describing polariton-enhanced spin-orbit interaction of light on the surface of two-dimensional media. Starting from the integral formulation of electromagnetic scattering, we exploit the reduced dimensionality of the system to introduce a quantum-like formalism particularly suitable to fully take advantage of rotational invariance. Our description is closely related to that of a fictitious spin one quantum particle living in the atomically thin medium, whose orbital, spin and total angular momenta play a key role in the scattering process. Conservation of total angular momentum upon scattering enables to physically unveil the interaction between radiation and the two-dimensional material along with the detailed exchange processes among orbital and spin components. In addition, we specialize our model to doped extended graphene, finding such spin-orbit interaction to be dramatically enhanced by the excitation of surface plasmon polaritons propagating radially along the graphene sheet. We provide several examples of the enormous possibilities offered by plasmon-enhanced spin-orbit interaction of light including vortex generation, mixing, and engineering of tunable deep subwavelength arrays of optical traps in the near field. Our results hold great potential for the development of nano-scaled quantum active elements and logic gates for the manipulation of hyper-entangled photon states as well as for the design of artificial media imprinted by engineered photonic lattices tweezing cold atoms into the desired patterns.


title: Mather theory and symplectic rigidity

abstract: Using methods from symplectic topology, we prove existence of invariant variational measures associated to the flow $\phi_H$ of a Hamiltonian $H\in C^{\infty}(M)$ on a symplectic manifold $(M,\omega)$. These measures coincide with Mather measures (from Aubry-Mather theory) in the Tonelli case. We compare properties of the supports of these measures to classical Mather measures and we construct an example showing that their support can be extremely unstable when $H$ fails to be convex, even for nearly integrable $H$. Parts of these results extend work by Viterbo and Vichery. Using ideas due to Entov-Polterovich we also detect interesting invariant measures for $\phi_H$ by studying a generalization of the symplectic shape of sublevel sets of $H$. This approach differs from the first one in that it works also for $(M,\omega)$ in which every compact subset can be displaced. We present applications to Hamiltonian systems on $\mathbb{R}^{2n}$ and twisted cotangent bundles.


title: Learning Non-Stationary Space-Time Models for Environmental Monitoring

abstract: One of the primary aspects of sustainable development involves accurate understanding and modeling of environmental phenomena. Many of these phenomena exhibit variations in both space and time and it is imperative to develop a deeper understanding of techniques that can model space-time dynamics accurately. In this paper we propose NOSTILL-GP - NOn-stationary Space TIme variable Latent Length scale GP, a generic non-stationary, spatio-temporal Gaussian Process (GP) model. We present several strategies, for efficient training of our model, necessary for real-world applicability. Extensive empirical validation is performed using three real-world environmental monitoring datasets, with diverse dynamics across space and time. Results from the experiments clearly demonstrate general applicability and effectiveness of our approach for applications in environmental monitoring.


title: Critical phenomena of charged Einstein-Gauss-Bonnet black holes with charged scalar hair

abstract: Einstein-Gauss-Bonnet-gravity (EGB) coupled minimally to a $U(1)$ gauged, massive scalar field possesses -- for appropriate choices of the $U(1)$ charge -- black hole solutions that carry charged scalar hair if the frequency of the harmonic time-dependence of the scalar field is equal to the upper bound on the superradiant frequency. The existence of these solutions has first been discussed in \cite{Grandi:2017zgz}. In this paper, we demonstrate that the critical value of the scalar charge results from the requirement of non-extremality of the charged black hole solutions and the fact that the scalar field should not escape to infinity. Moreover, we investigate the hairy black holes in more detail and demonstrate that the branch of these solutions joins the branch of the corresponding charged EGB black hole for vanishing scalar field, but is {\it not} connected to the branch of boson stars in the limit of vanishing horizon radius. This indicates that it is unlikely that these black holes appear from the collapse of the corresponding boson stars. Finally, we prove a No-hair theorem for charged scalar fields with harmonic time-dependence for static, spherically symmetric, asymptotically flat electro-vacuum black holes in $d$ space-time dimensions and hence demonstrate that the GB term is crucial for the existence of the hairy black holes discussed in this paper.


title: HD66051: the first eclipsing binary hosting an early-type magnetic star

abstract: Early-type magnetic stars are rarely found in close binary systems. No such objects were known in eclipsing binaries prior to this study. Here we investigated the eclipsing, spectroscopic double-lined binary HD66051, which exhibits out-of-eclipse photometric variations suggestive of surface brightness inhomogeneities typical of early-type magnetic stars. Using a new set of high-resolution spectropolarimetric observations, we discovered a weak magnetic field on the primary and found intrinsic, element-dependent variability in its spectral lines. The magnetic field structure of the primary is dominated by a nearly axisymmetric dipolar component with a polar field strength $B_{\rm d}\approx600$ G and an inclination with respect to the rotation axis of $\beta_{\rm d}=13^{\rm o}$. A weaker quadrupolar component is also likely to be present. We combined the radial velocity measurements derived from our spectra with archival optical photometry to determine fundamental masses (3.16 and 1.75 $M_\odot$) and radii (2.78 and 1.39 $R_\odot$) with a 1-3% precision. We also obtained a refined estimate of the effective temperatures (13000 and 9000 K) and studied chemical abundances for both components with the help of disentangled spectra. We demonstrate that the primary component of HD66051 is a typical late-B magnetic chemically peculiar star with a non-uniform surface chemical abundance distribution. It is not an HgMn-type star as suggested by recent studies. The secondary is a metallic-line star showing neither a strong, global magnetic field nor intrinsic spectral variability. Fundamental parameters provided by our work for this interesting system open unique possibilities for probing interior structure, studying atomic diffusion, and constraining binary star evolution.


abstract: In this paper, we discuss the integer decomposition property for Cayley sums and Minkowski sums of lattice polytopes. In fact, we characterize when Cayley sums have the integer decomposition property in terms of Minkowski sums. Moreover, by using this characterization, we consider when Cayley sums and Minkowski sums of $2$-convex-normal lattice polytopes have the integer decomposition property. Finally, we also discuss the level property for Minkowski sums and Cayley sums.


title: Effect of uniaxial stress on the magnetic phases of CeAuSb$_2$

abstract: We present results of measurements of resistivity of \CAS{} under the combination of $c$-axis magnetic field and in-plane uniaxial stress. In unstressed \CAS{} there are two magnetic phases. The low-field A phase is a single-component spin-density wave (SDW), with $\mathbf{q} = (\eta, \pm \eta, 1/2)$, and the high-field B phase consists of microscopically coexisting $(\eta, \eta, 1/2)$ and $(\eta, -\eta, 1/2)$ spin-density waves. Pressure along a $\langle 100 \rangle$ lattice direction is a transverse field to both of these phases, and so initially has little effect, however eventually induces new low- and high-field phases in which the principal axes of the SDW components appear to have rotated to the $\langle 100 \rangle$ directions. Under this strong $\langle 100 \rangle$ compression, the field evolution of the resistivity is much smoother than at zero strain: In zero strain, there is a strong first-order transition, while under strong $\langle 100 \rangle$ it becomes much broader. We hypothesize that this is a consequence of the uniaxial stress lifting the degeneracy between the (100) and (010) directions.


title: A theory of linear typings as flows on 3-valent graphs

abstract: Building on recently established enumerative connections between lambda calculus and the theory of embedded graphs (or "maps"), this paper develops an analogy between typing (of lambda terms) and coloring (of maps). Our starting point is the classical notion of an abelian group-valued "flow" on an abstract graph (Tutte, 1954). Typing a linear lambda term may be naturally seen as constructing a flow (on an embedded 3-valent graph with boundary) valued in a more general algebraic structure consisting of a preordered set equipped with an "implication" operation and unit satisfying composition, identity, and unit laws. Interesting questions and results from the theory of flows (such as the existence of nowhere-zero flows) may then be re-examined from the standpoint of lambda calculus and logic. For example, we give a characterization of when the local flow relations (across vertices) may be categorically lifted to a global flow relation (across the boundary), proving that this holds just in case the underlying map has the orientation of a lambda term. We also develop a basic theory of rewriting of flows that suggests topological meanings for classical completeness results in combinatory logic, and introduce a polarized notion of flow, which draws connections to the theory of proof-nets in linear logic and to bidirectional typing.


title: A matrix-free approach to parallel and memory-efficient deformable image registration

abstract: We present a novel computational approach to fast and memory-efficient deformable image registration. In the variational registration model, the computation of the objective function derivatives is the computationally most expensive operation, both in terms of runtime and memory requirements. In order to target this bottleneck, we analyze the matrix structure of gradient and Hessian computations for the case of the normalized gradient fields distance measure and curvature regularization. Based on this analysis, we derive equivalent matrix-free closed-form expressions for derivative computations, eliminating the need for storing intermediate results and the costs of sparse matrix arithmetic. This has further benefits: (1) matrix computations can be fully parallelized, (2) memory complexity for derivative computation is reduced from linear to constant, and (3) overall computation times are substantially reduced. In comparison with an optimized matrix-based reference implementation, the CPU implementation achieves speedup factors between 3.1 and 9.7, and we are able to handle substantially higher resolutions. Using a GPU implementation, we achieve an additional speedup factor of up to 9.2. Furthermore, we evaluated the approach on real-world medical datasets. On ten publicly available lung CT images from the DIR-Lab 4DCT dataset, we achieve the best mean landmark error of 0.93 mm compared to other submissions on the DIR-Lab website, with an average runtime of only 9.23 s. Complete non-rigid registration of full-size 3D thorax-abdomen CT volumes from oncological follow-up is achieved in 12.6 s. The experimental results show that the proposed matrix-free algorithm enables the use of variational registration models also in applications which were previously impractical due to memory or runtime restrictions.


title: Spindle Nodal Chain in Three-Dimensional alpha' Boron

abstract: Topological metal/semimetals (TMs) have emerged as a new frontier in the field of quantum materials. A few two-dimensional (2D) boron sheets have been suggested as Dirac materials, however, to date TMs made of three-dimensional (3D) boron structures have not been found. Herein, by means of systematic first principles computations, we discovered that a rather stable 3D boron allotrope, namely 3D-alpha' boron, is a nodal-chain semimetal. In the momentum space, six nodal lines and rings contact each other and form a novel spindle nodal chain. This 3D-alpha' boron can be formed by stacking 2D wiggle alpha' boron sheets, which are also nodal-ring semimetals. In addition, our chemical bond analysis revealed that the topological properties of the 3D and 2D boron structures are related to the pi bonds between boron atoms, however, the bonding characteristics are different from those in the 2D and 3D carbon structures.


title: Quantum chaos and entanglement in ergodic and non-ergodic systems

abstract: We study entanglement entropy (EE) as a signature of quantum chaos in ergodic and non-ergodic systems. In particular we look at the quantum kicked top and kicked rotor as multi-qubit systems, and investigate the single qubit EE which characterizes bipartite entanglement of this qubit with the rest of the system. We study the correspondence of the Kolmogorov-Sinai entropy of the classical kicked systems with the EE of their quantum counterparts. We find that EE is a signature of global chaos in ergodic systems, and local chaos in non-ergodic systems. In particular, we show that EE can be maximised even when systems are highly non-ergodic, when the corresponding classical system is locally chaotic. In contrast, we find evidence that the quantum analogue of Kolmogorov-Arnol'd-Noser (KAM) tori are tori of low entanglement entropy. We conjecture that entanglement should play an important role in any quantum KAM theory.


title: Persistent Monitoring of Stochastic Spatio-temporal Phenomena with a Small Team of Robots

abstract: This paper presents a solution for persistent monitoring of real-world stochastic phenomena, where the underlying covariance structure changes sharply across time, using a small number of mobile robot sensors. We propose an adaptive solution for the problem where stochastic real-world dynamics are modeled as a Gaussian Process (GP). The belief on the underlying covariance structure is learned from recently observed dynamics as a Gaussian Mixture (GM) in the low-dimensional hyper-parameters space of the GP and adapted across time using Sequential Monte Carlo methods. Each robot samples a belief point from the GM and locally optimizes a set of informative regions by greedy maximization of the submodular entropy function. The key contributions of this paper are threefold: adapting the belief on the covariance using Markov Chain Monte Carlo (MCMC) sampling such that particles survive even under sharp covariance changes across time; exploiting the belief to transform the problem of entropy maximization into a decentralized one; and developing an approximation algorithm to maximize entropy on a set of informative regions in the continuous space. We illustrate the application of the proposed solution through extensive simulations using an artificial dataset and multiple real datasets from fixed sensor deployments, and compare it to three competing state-of-the-art approaches.


title: $\Delta L \ge 4$ processes

abstract: We discuss the experimental prospects for observing processes which violate lepton number ($\Delta L$) in four units (or more). First, we reconsider neutrinoless quadruple beta decay, deriving a model independent and very conservative lower limit on its half-life of the order of $10^{41}$ ys for $^{150}$Nd. This renders quadruple beta decay unobservable for any feasible experiment. We then turn to a more general discussion of different possible low-energy processes with values of $\Delta L \ge 4$. A simple operator analysis leads to rather pessimistic conclusions about the observability at low-energy experiments in all cases we study. However, the situation looks much brighter for accelerator experiments. For two example models with $\Delta L=4$ and another one with $\Delta L=5$, we show how the LHC or a hypothetical future pp-collider, such as the FCC, could probe multi-lepton number violating operators at the TeV scale.


title: On $\mathrm{C}^*$-algebras associated to product systems

abstract: Let $P$ be a unital subsemigroup of a group $G$. We propose an approach to $\mathrm{C}^*$-algebras associated to product systems over $P$. We call the $\mathrm{C}^*$-algebra of a given product system $\mathcal{E}$ its covariance algebra and denote it by $A\times_{\mathcal{E}}P$, where $A$ is the coefficient $\mathrm{C}^*$-algebra. We prove that our construction does not depend on the embedding $P\hookrightarrow G$ and that a representation of $A\times_{\mathcal{E}}P$ is faithful on the fixed-point algebra for the canonical coaction of $G$ if and only if it is faithful on $A$. We compare this with other constructions in the setting of irreversible dynamical systems, such as Cuntz--Nica--Pimsner algebras, Fowler's Cuntz--Pimsner algebra, semigroup $\mathrm{C}^*$-algebras of Xin Li and Exel's crossed products by interaction groups.


title: Numerical comparison of mass-conservative schemes for the Gross-Pitaevskii equation

abstract: In this paper we present a numerical comparison of various mass-conservative discretizations for the time-dependent Gross-Pitaevskii equation. We have three main objectives. First, we want to clarify how purely mass-conservative methods perform compared to methods that are additionally energy-conservative or symplectic. Second, we shall compare the accuracy of energy-conservative and symplectic methods among each other. Third, we will investigate if a linearized energy-conserving method suffers from a loss of accuracy compared to an approach which requires to solve a full nonlinear problem in each time-step. In order to obtain a representative comparison, our numerical experiments cover different physically relevant test cases, such as traveling solitons, stationary multi-solitons, Bose-Einstein condensates in an optical lattice and vortex pattern in a rapidly rotating superfluid. We shall also consider a computationally severe test case involving a pseudo Mott insulator. Our space discretization is based on finite elements throughout the paper. We will also give special attention to long time behavior and possible coupling conditions between time-step sizes and mesh sizes. The main observation of this paper is that mass conservation alone will not lead to a competitive method in complex settings. Furthermore, energy-conserving and symplectic methods are both reliable and accurate, yet, the energy-conservative schemes achieve a visibly higher accuracy in our test cases. Finally, the scheme that performs best throughout our experiments is an energy-conserving relaxation scheme with linear time-stepping proposed by C. Besse (SINUM,42(3):934--952,2004).


title: Size and shape of skyrmions for variable Dzyaloshinskii-Moriya interaction and uniaxial anisotropy

abstract: We have performed micromagnetic simulations to study the formation of skyrmions in ferromagnetic elements with different shapes having perpendicular anisotropy. The strength of Dzyaloshinskii-Moriya interaction (D) and uniaxial anisotropy (K) are varied to elucidate the regime in which skyrmion formation can take place. It is found that for a certain combination of D and K skyrmion formation does not happen. Further we also observed that for large D and small K values, finite size effect dominates which in turn hinders formation of typical N'eel (spherical) skyrmions. However the resulting magnetic phase is skyrmionic in nature and has different shape. We also have found that the shape of the magnetic nano element has a significant role in determining the final magnetic state in addition to the competing D and K values.


title: Maximal discrete sparsity in parabolic optimal control with measures

abstract: We consider variational discretization of a parabolic optimal control problem governed by space-time measure controls. For the state discretization we use a Petrov-Galerkin method employing piecewise constant states and piecewise linear and continuous test functions in time. For the space discretization we use piecewise linear and continuous functions. As a result the controls are composed of Dirac measures in space-time, centered at points on the discrete space-time grid. We prove that the optimal discrete states and controls converge strongly in $L^q$ and weakly-$*$ in $\mathcal{M}$, respectively, to their smooth counterparts, where $q \in (1,\min\{2,1+2/d\}]$ is the spatial dimension. Furthermore, we compare our approach to a approach by Casas, E. and Kunisch, K., where the corresponding control problem is discretized employing a discontinuous Galerkin method for the state discretization and where the discrete controls are piecewise constant in time and Dirac measures in space. Numerical experiments highlight the features of our discrete approach.


title: Horizon temperature on the real line

abstract: We illustrate the analogue of the Unruh effect for a quantum system on the real line. Our derivation relies solely on basic elements of representation theory of the group of affine transformations without a notion of time or metric. Our result shows that a thermal distribution naturally emerges in connecting quantum states belonging to representations associated to distinct notions of translational symmetry.


title: Quantum gravity effect on the Hawking radiation of charged rotating BTZ black hole

abstract: In this study, the quantum gravity effect on the tunnelling radiation of charged massive spin-0 scalar particle from 2+1 dimensional charged rotating Banados-Teitelboim-Zanelli (BTZ) black hole is looked into by using the Hamilton-Jacobi approach. For this, we calculate the modified Hawking temperature of the black hole by using the modified Klein-Gordon equation based on the Generalized Uncertainty Principle (GUP), and we noticed that the modified Hawking temperature of the black hole depends not only on the black hole properties, but also on the angular momentum, energy, charge and mass of the tunnelling scalar particle. Using the modified Hawking temperature, we discussed the stability of the black hole in the context of the modified heat capacity, and observed that it might undergo both first and second-type phase transitions in the presence of the quantum gravity effect, but just a first-type transition in the absence of the quantum gravity effect. Furthermore, we investigated the modified Hawking temperature of the black hole by using the tunnelling processes of the charged massive Dirac and vector boson particles. We observed that scalar, Dirac and vector particles are tunnelled from the black hole completely differently from each other in the presence of the quantum gravity effect.


title: Analysis of a time-stepping scheme for time fractional diffusion problems with nonsmooth data

abstract: This paper establishes the convergence of a time-steeping scheme for time fractional diffusion problems with nonsmooth data. We first analyze the regularity of the model problem with nonsmooth data, and then prove that the time-steeping scheme possesses optimal convergence rates in $ L^2(0,T;L^2(\Omega)) $-norm and $ L^2(0,T;H_0^1(\Omega)) $-norm with respect to the regularity of the solution. Finally, numerical results are provided to verify the theoretical results.


title: Speckle interferometry at SOAR in 2016 and 2017

abstract: The results of speckle interferometric observations at the 4.1 m SOAR telescope in 2016 and 2017 are given, totaling 2483 measurements of 1570 resolved pairs and 609 non-resolutions. We describe briefly recent changes in the instrument and observing method and quantify the accuracy of the pixel scale and position angle calibration. Comments are given on 44 pairs resolved here for the first time. Orbital motion of the newly resolved subsystem BU~83 Aa,Ab roughly agrees with its 36 year astrometric orbit proposed by J.~Dommanget. Most Tycho binaries examined here turned out to be spurious.


title: Random Asynchronous Iterations in Distributed Coordination Algorithms

abstract: Distributed coordination algorithms (DCA) carry out information processing processes among a group of networked agents without centralized information fusion. Though it is well known that DCA characterized by an SIA (stochastic, indecomposable, aperiodic) matrix generate consensus asymptotically via synchronous iterations, the dynamics of DCA with asynchronous iterations have not been studied extensively, especially when viewed as stochastic processes. This paper aims to show that for any given irreducible stochastic matrix, even non-SIA, the corresponding DCA lead to consensus successfully via random asynchronous iterations under a wide range of conditions on the transition probability. Particularly, the transition probability is neither required to be independent and identically distributed, nor characterized by a Markov chain.


title: A New Quantum Operator for Distance

abstract: We introduce a new semi-relativistic quantum operator for the length of the worldline a particle traces out as it moves. In this article the operator is constructed in a heuristic way and some of its elementary properties are explored. The operator ends up depending in a very complicated way on the potential of the system it is to act on so as a proof of concept we use it to analyze the expected distance traveled by a free Gaussian wavepacket with some initial momentum. It is shown in this case that the distance such a particle travels becomes light-like as its mass vanishes and agrees with the classical result for macroscopic masses. This preliminary result has minor implications for the Weak Equivalence Principle (WEP) in quantum mechanics. In particular it shows that the logical relationship between two formulations of the WEP in classical mechanics extends to quantum mechanics. That our result is qualitatively consistent with the work of others emboldens us to start the task of evaluating the new operator in non-zero potentials. However, we readily acknowledge that the looseness in the definition of our operator means that all of our so-called results are highly speculative. Plans for future work with the new operator are discussed in the last section.


title: Convergence and Concentration of Empirical Measures under Wasserstein Distance in Unbounded Functional Spaces

abstract: We provide upper bounds of the expected Wasserstein distance between a probability measure and its empirical version, generalizing recent results for finite dimensional Euclidean spaces and bounded functional spaces. Such a generalization can cover Euclidean spaces with large dimensionality, with the optimal dependence on the dimensionality. Our method also covers the important case of Gaussian processes in separable Hilbert spaces, with rate-optimal upper bounds for functional data distributions whose coordinates decay geometrically or polynomially. Moreover, our bounds of the expected value can be combined with mean-concentration results to yield improved exponential tail probability bounds for the Wasserstein error of empirical measures under Bernstein-type or log Sobolev-type conditions.


title: Net baryon diffusion in fluid dynamic simulations of relativistic heavy-ion collisions

abstract: A hybrid (hydrodynamics + hadronic transport) theoretical framework is assembled to model the bulk dynamics of relativistic heavy-ion collisions at energies accessible in the Beam Energy Scan (BES) program at the Relativistic Heavy-Ion Collider (RHIC) and the NA61/SHINE experiment at CERN. The system's energy-momentum tensor and net baryon current are evolved according to relativistic hydrodynamics with finite shear viscosity and non-zero net baryon diffusion. Our hydrodynamic description is matched to a hadronic transport model in the dilute region. With this fully integrated theoretical framework, we present a pilot study of the hadronic chemistry, particle spectra, and anisotropic flow. Phenomenological effects of a non-zero net-baryon current and its diffusion on hadronic observables are presented for the first time. The importance of the hadronic transport phase is also investigated.


title: Optimal storage of a single photon by a single intra-cavity atom

abstract: We theoretically analyse the efficiency of a quantum memory for single photons. The photons propagate along a transmission line and impinge on one of the mirrors of a high-finesse cavity. The quantum memory is constituted by a single atom within the optical resonator. Photon storage is realised by the controlled transfer of the photonic excitation into a metastable state of the atom and occurs via a Raman transition with a suitably tailored laser pulse, which drives the atom. Our study is supported by numerical simulations, in which we include the modes of the transmission line and we use the experimental parameters of existing experimental setups. It reproduces the results derived using input-output theory in the corresponding regime and can be extended to compute dynamics where the input-output formalism cannot be straightforwardly applied. Our analysis determines the maximal storage efficiency, namely, the maximal probability to store the photon in a stable atomic excitation, in the presence of spontaneous decay and cavity parasitic losses. It further delivers the form of the laser pulse that achieves the maximal efficiency by partially compensating parasitic losses. We numerically assess the conditions under which storage based on adiabatic dynamics is preferable to non-adiabatic pulses. Moreover, we systematically determine the shortest photon pulse that can be efficiently stored as a function of the system parameters.


title: The time-scales of major mergers from simulations of isolated binary galaxy collisions

abstract: A six-dimensional parameter space based on high-resolution numerical simulations of isolated binary galaxy collisions has been constructed to investigate the dynamical friction timescales, $\tau_{\rm mer}$, for major mergers. Our experiments follow the gravitational encounters between $\sim 600$ pairs of similarly massive late- and early-type galaxies with orbital parameters compliant to the predictions of the LambdaCDM cosmology. We analyze the performance of different schemes for tracking the secular evolution of mergers, finding that the product of the intergalactic distance and velocity is best suited to identify the time of coalescence. In contrast, a widely used merger time estimator such as the exhaustion of the orbital spin is shown to systematically underpredict $\tau_ {\rm mer}$, resulting in relative errors that can reach 60% for nearly radial encounters. Regarding the internal spins of the progenitors, we find that they can lead to total variations in the merger times above 30% in highly circular encounters, whereas only that of the principal halo is capable of modulating the strength of the interaction prevailing throughout a merger. The comparison of our simulated merger times with predictions from different variants of a well-known fitting formula has revealed an only partially satisfactory agreement, which has led us to recalculate the values of the coefficients of these expressions to obtain relations that fit perfectly major mergers. The observed biases between data and predictions, that do not only apply to the present work, are inconsistent with expectations from differences in the degree of idealization of the collisions, their metric, spin-related biases, or the simulation set-up. This hints to a certain lack of accuracy of the dynamical friction modelling, arising perhaps from a still not quite complete identification of the parameters governing orbital decay.


title: Quantum Walk Search on Kronecker Graphs

abstract: Kronecker graphs, obtained by repeatedly performing the Kronecker product of the adjacency matrix of an "initiator" graph with itself, have risen in popularity in network science due to their ability to generate complex networks with real-world properties. In this paper, we explore spatial search by continuous-time quantum walk on Kronecker graphs. Specifically, we give analytical proofs for quantum search on first-, second-, and third-order Kronecker graphs with the complete graph as the initiator, showing that search takes Grover's $O(\sqrt{N})$ time. Numerical simulations indicate that higher-order Kronecker graphs with the complete initiator also support optimal quantum search.


title: A variant of the Mordell-Lang conjecture

abstract: The Mordell-Lang conjecture (proven by Faltings, Vojta and McQuillan) states that the intersection of a subvariety $V$ of a semiabelian variety $G$ defined over an algebraically closed field $\mathbb{k}$ of characteristic $0$ with a finite rank subgroup $\Gamma \le G(\mathbb{k})$ is a finite union of cosets of subgroups of $\Gamma$. We explore a variant of this conjecture when $G$ is a product of an abelian variety $A$ defined over $\mathbb{k}$ with the additive group $\mathbb{G}_a$.


title: Contribution of phase-mixing of Alfv\'en waves to coronal heating in multi-harmonic loop oscillations

abstract: Kink oscillations of a coronal loop are observed and studied in detail because they provide a unique probe into the structure of coronal loops through MHD seismology and a potential test of coronal heating through the phase-mixing of Alfv\'en waves. In particular, recent observations show that standing oscillations of loops often involve also higher harmonics, beside the fundamental mode. The damping of these kink oscillations is explained by mode coupling with Alfv\'en waves. We investigate the consequences for wave-based coronal heating of higher harmonics and what coronal heating observational signatures we may use to infer the presence of higher harmonic kink oscillations. We perform a set of non-ideal MHD simulations where the damping of the kink oscillation of a flux tube via mode coupling is modelled. Our MHD simulation parameters are based on the seismological inversion of an observation for which the first three harmonics are detected. We study the phase-mixing of Alfv\'en waves that leads to the deposition of heat in the system, and we apply the seismological inversion techniques to the MHD simulation output. We find that the heating due to phase-mixing of the Alfv\'en waves triggered by the damping of the kink oscillation is relatively small, however we can illustrate i) how the heating location drifts due to the subsequent damping of lower order harmonics. We also address the role of the higher order harmonics and the width of the boundary shell in the energy deposition. We conclude that the coronal heating due to phase-mixing seems not to provide enough energy to maintain the thermal structure of the solar corona even when multi-harmonics oscillations are included, and these oscillations play an inhibiting role in the development of smaller scale structures.


title: Intermediate Data Caching Optimization for Multi-Stage and Parallel Big Data Frameworks

abstract: In the era of big data and cloud computing, large amounts of data are generated from user applications and need to be processed in the datacenter. Data-parallel computing frameworks, such as Apache Spark, are widely used to perform such data processing at scale. Specifically, Spark leverages distributed memory to cache the intermediate results, represented as Resilient Distributed Datasets (RDDs). This gives Spark an advantage over other parallel frameworks for implementations of iterative machine learning and data mining algorithms, by avoiding repeated computation or hard disk accesses to retrieve RDDs. By default, caching decisions are left at the programmer's discretion, and the LRU policy is used for evicting RDDs when the cache is full. However, when the objective is to minimize total work, LRU is woefully inadequate, leading to arbitrarily suboptimal caching decisions. In this paper, we design an algorithm for multi-stage big data processing platforms to adaptively determine and cache the most valuable intermediate datasets that can be reused in the future. Our solution automates the decision of which RDDs to cache: this amounts to identifying nodes in a direct acyclic graph (DAG) representing computations whose outputs should persist in the memory. Our experiment results show that our proposed cache optimization solution can improve the performance of machine learning applications on Spark decreasing the total work to recompute RDDs by 12%.


title: Periodic shearing motions in the Jovian magnetosphere causing a localized peak in the main auroral emission close to noon

abstract: Recently, a transient localized brightness enhancement has been observed in Jupiter's main auroral emission close to noon by Palmaerts et al. (2014). We use results from three-dimensional global MHD simulations to understand what is causing this localized peak in the main emission. In the simulations, the peak occurs every rotation period and is due to shearing motions in the magnetodisk. These shearing motions are caused by heavy flux-tubes being accelerated to large azimuthal speeds at dawn. The centrifugal force acting on these flux-tubes is then so high that they rapidly move away from the planet. When they reach noon, their azimuthal velocity decreases, thus reducing the centrifugal force, and allowing the flux-tubes to move back closer to Jupiter. The shearing motions associated with this periodic phenomenon locally increase the field aligned currents in the simulations, thus causing a transient brightness enhancement in the main auroral emission, similar to the one observed by Palmaerts et al. (2014).


title: Certified Graph View Maintenance with Regular Datalog

abstract: We employ the Coq proof assistant to develop a mechanically-certified framework for evaluating graph queries and incrementally maintaining materialized graph instances, also called views. The language we use for defining queries and views is Regular Datalog (RD) -- a notable fragment of non-recursive Datalog that can express complex navigational queries, with transitive closure as native operator. We first design and encode the theory of RD and then mechanize a RD-specific evaluation algorithm capable of fine-grained, incremental graph view computation, which we prove sound with respect to the declarative RD semantics. By using the Coq extraction mechanism, we test an Ocaml version of the verified engine on a set of preliminary benchmarks. Our development is particularly focused on leveraging existing verification and notational techniques to: a) define mechanized properties that can be easily understood by logicians and database researchers and b) attain formal verification with limited effort. Our work is the first step towards a unified, machine-verified, formal framework for dynamic graph query languages and their evaluation engines. This paper is under consideration for acceptance in TPLP.


title: Beyond the random phase approximation with a local exchange vertex

abstract: With the aim of constructing an electronic structure approach that systematically goes beyond the GW and random phase approximation (RPA) we introduce a vertex correction based on the exact-exchange (EXX) potential of time-dependent density functional theory. The EXX vertex function is constrained to be local but is expected to capture similar physics as the Hartree-Fock vertex. With the EXX vertex we then consistently unify different beyond-RPA approaches such as the various re-summations of RPA with exchange (RPAx) and the second order screened exchange (SOSEX) approximation. The theoretical analysis is supported by numerical studies on the hydrogen dimer and the electron gas, and we discuss the importance of including the vertex correction in both the screened interaction and the self energy. Finally, we give details on our implementation within the plane-wave pseudo potential framework and demonstrate the excellent performance of the different RPAx methods in describing the energetics of hydrogen and van der Waals bonds.


title: UPC condition with parameter for subanalytic sets

abstract: In 1986 Paw{\l}ucki and Ple\'sniak introduced the notion of {\sl uniformly polynomially cuspidal} (UPC) sets and proved that every relatively compact and fat subanalytic subset of ${\Rz}^n$ satisfies the UPC condition. Herein we investigate the UPC property of the sections of a relatively compact open subanalytic set $E\subset{\Rz}^k\times{\Rz}^n$ and we show that two of the three parameters in the UPC condition can be chosen independently of the section, while the third one depends generally on the point defining the section.


title: Symmetry-protected topological phases in lattice gauge theories: topological QED$_2$

abstract: The interplay of symmetry, topology, and many-body effects in the classification of possible phases of matter poses a formidable challenge that is attracting great attention in condensed-matter physics. Such many-body effects are typically induced by inter-particle interactions involving an action at a distance, such as the Coulomb interaction between electrons in a symmetry-protected topological (SPT) phase. In this work, we show that similar phenomena also appear in high-energy physics, where inter-particle interactions are mediated by gauge bosons, and constrained by a local gauge symmetry. We introduce a variant of the so-called Schwinger model, which describes quantum electrodynamics in (1+1) dimensions (QED$_2$), and show that it can host SPT phases with a topological contribution to the vacuum {\theta} angle, which leads to a new type of topological QED$_2.$ We use bosonization and density-matrix renormalization group techniques to study its rich phase diagram in great detail, and present a scheme for its realization in experiments of ultra-cold atoms in optical lattices.

categories: cond-mat.quant-gas cond-mat.str-el hep-lat hep-th quant-ph


title: Eigenvalue variation under moving mixed Dirichlet-Neumann boundary conditions and applications

abstract: We deal with the sharp asymptotic behaviour of eigenvalues of elliptic operators with varying mixed Dirichlet-Neumann boundary conditions. In case of simple eigenvalues, we compute explicitly the constant appearing in front of the expansion's leading term. This allows inferring some remarkable consequences for Aharonov-Bohm eigenvalues when the singular part of the operator has two coalescing poles.


title: Coordinate space calculation of two- and three-loop sunrise-type diagrams, elliptic functions and truncated Bessel integral identities

abstract: We integrate three-loop sunrise-type vacuum diagrams in $D_0=4$ dimensions with four different masses using configuration space techniques. The finite parts of our results are in numerical agreement with corresponding three-loop calculations in momentum space. Using some of the closed form results of the momentum space calculation we arrive at new integral identities involving truncated integrals of products of Bessel functions. For the non-degenerate finite two-loop sunrise-type vacuum diagram in $D_0=2$ dimensions we make use of the known closed form $p$-space result to express the moment of a product of three Bessel functions in terms of a sum of Claussen polylogarithms. Using results for the nondegenerate two-loop sunrise diagram from the literature in $D_0=2$ dimensions we obtain a Bessel function integral identity in terms of elliptic functions.


title: Improving timing sensitivity in the microhertz frequency regime: limits from PSR J1713$+$0747 on gravitational waves produced by super-massive black-hole binaries

abstract: We search for continuous gravitational waves (CGWs) produced by individual super-massive black-hole binaries (SMBHBs) in circular orbits using high-cadence timing observations of PSR J1713$+$0747. We observe this millisecond pulsar using the telescopes in the European Pulsar Timing Array (EPTA) with an average cadence of approximately 1.6 days over the period between April 2011 and July 2015, including an approximately daily average between February 2013 and April 2014. The high-cadence observations are used to improve the pulsar timing sensitivity across the GW frequency range of $0.008-5$ $\mu$Hz. We use two algorithms in the analysis, including a spectral fitting method and a Bayesian approach. For an independent comparison, we also use a previously published Bayesian algorithm. We find that the Bayesian approaches provide optimal results and the timing observations of the pulsar place a 95 per cent upper limit on the sky-averaged strain amplitude of CGWs to be $\lesssim3.5\times10^{-13}$ at a reference frequency of 1 $\mu$Hz. We also find a 95 per cent upper limit on the sky-averaged strain amplitude of low-frequency CGWs to be $\lesssim1.4\times10^{-14}$ at a reference frequency of 20~nHz.


title: Activity-dependent self-regulation of viscous length scales in biological systems

abstract: Cellular cortex, which is a highly viscous thin cytoplasmic layer just below the cell membrane, controls the cell's mechanical properties, which can be characterized by a hydrodynamic length scale $\ell$. Cells actively regulate $\ell$ via the activity of force generating molecules, such as myosin II. Here we develop a general theory for such systems through coarse-grained hydrodynamic approach including activity in the static description of the system providing an experimentally accessible parameter and elucidate the detailed mechanism of how a living system can actively self-regulate its hydrodynamic length scale, controlling the rigidity of the system. Remarkably, we find that $\ell$, as a function of activity, behaves universally and roughly inversely proportional to the activity of the system. Our theory rationalizes a number of experimental findings on diverse systems and comparison of our theory with existing experimental data show good agreement.


title: Geometry and temperature chaos in mixed spherical spin glasses at low temperature - the perturbative regime

abstract: We study the Gibbs measure of mixed spherical $p$-spin glass models at low temperature, in (part of) the 1-RSB regime, including, in particular, models close to pure in an appropriate sense. We show that the Gibbs measure concentrates on spherical bands around deep critical points of the (extended) Hamiltonian restricted to the sphere of radius $\sqrt N q_\star$, where $q_\star^2$ is the rightmost point in the support of the overlap distribution. We also show that the relevant critical points are pairwise orthogonal for two different low temperatures. This allows us to explain why temperature chaos occurs for those models, in contrast to the pure spherical models.


title: Decoupled Parallel Backpropagation with Convergence Guarantee

abstract: Backpropagation algorithm is indispensable for the training of feedforward neural networks. It requires propagating error gradients sequentially from the output layer all the way back to the input layer. The backward locking in backpropagation algorithm constrains us from updating network layers in parallel and fully leveraging the computing resources. Recently, several algorithms have been proposed for breaking the backward locking. However, their performances degrade seriously when networks are deep. In this paper, we propose decoupled parallel backpropagation algorithm for deep learning optimization with convergence guarantee. Firstly, we decouple the backpropagation algorithm using delayed gradients, and show that the backward locking is removed when we split the networks into multiple modules. Then, we utilize decoupled parallel backpropagation in two stochastic methods and prove that our method guarantees convergence to critical points for the non-convex problem. Finally, we perform experiments for training deep convolutional neural networks on benchmark datasets. The experimental results not only confirm our theoretical analysis, but also demonstrate that the proposed method can achieve significant speedup without loss of accuracy.


title: The Estimation Lie Algebra Associated with Quantum Filters

abstract: We introduce the Lie algebra of super-operators associated with a quantum filter, specifically emerging from the Stratonovich calculus. In classical filtering, the analogue algebra leads to a geometric theory of nonlinear filtering which leads to well-known results by Brockett and by Mitter characterizing potential models where the curse-of-dimensionality may be avoided, and finite dimensional filters obtained. We discuss the quantum analogue to these results. In particular, we see that, in the case where all outputs are subjected to homodyne measurement, the Lie algebra of super-operators is isomorphic to a Lie algebra of system operators from which one may approach the question of the existence of finite-dimensional filters.


title: Free energy landscapes in spherical spin glasses

abstract: We introduce and analyze free energy landscapes defined by associating to any point inside the sphere a free energy calculated on a thin spherical band around it, using many orthogonal replicas. This allows us to reinterpret, rigorously prove and extend for general spherical models the main ideas of the Thouless-Anderson-Palmer (TAP) approach originally introduced in the 70s for the Sherrington-Kirkpatrick model. In particular, we establish a TAP representation for the free energy, valid for any overlap value which can be sampled as many times as we wish in an appropriate sense. We call such overlaps multi-samplable. The correction to the Hamiltonian in the TAP representation arises in our analysis as the free energy of a certain model on an overlap dependent band. For the largest multi-samplable overlap it coincides with the Onsager reaction term from physics. For smaller multi-samplable overlaps the formula we obtain is new. We also derive the corresponding TAP equation for critical points. We prove all the above without appealing to the celebrated Parisi formula or the ultrametricity property. We prove that any overlap value in the support of the Parisi measure is multi-samplable. For generic models, we further show that the set of multi-samplable overlaps coincides with a certain set that arises in the characterization for the Parisi measure by Talagrand. The ultrametric tree of pure states can be embedded in the interior of the sphere in a natural way. For this embedding, we show that the points on the tree uniformly maximize the free energies we define. From this we conclude that the Hamiltonian at each point on the tree is approximately maximal over the sphere of same radius, and that points on the tree approximately solve the TAP equations for critical points.


title: Alleviating partisan gerrymandering: can math and computers help to eliminate wasted votes?

abstract: Partisan gerrymandering is a major cause for voter disenfranchisement in United States. However, convincing US courts to adopt specific measures to quantify gerrymandering has been of limited success to date. Recently, Stephanopoulos and McGhee introduced a new and precise measure of partisan gerrymandering via the so-called "efficiency gap" that computes the absolutes difference of wasted votes between two political parties in a two-party system. Quite importantly from a legal point of view, this measure was found legally convincing enough in a US appeals court in a case that claims that the legislative map of the state of Wisconsin was gerrymandered; the case is now pending in US Supreme Court. In this article, we show the following: (a) We provide interesting mathematical and computational complexity properties of the problem of minimizing the efficiency gap measure. To the best of our knowledge, these are the first non-trivial theoretical and algorithmic analyses of this measure of gerrymandering. (b) We provide a simple and fast algorithm that can "un-gerrymander" the district maps for the states of Texas, Virginia, Wisconsin and Pennsylvania by bring their efficiency gaps to acceptable levels from the current unacceptable levels. Our work thus shows that, notwithstanding the general worst case approximation hardness of the efficiency gap measure as shown by us,finding district maps with acceptable levels of efficiency gaps is a computationally tractable problem from a practical point of view. Based on these empirical results, we also provide some interesting insights into three practical issues related the efficiency gap measure. We believe that, should the US Supreme Court uphold the decision of lower courts, our research work and software will provide a crucial supporting hand to remove partisan gerrymandering.


title: Planetary Migration in Protoplanetary Disks

abstract: The known exoplanet population displays a great diversity of orbital architectures, and explaining the origin of this is a major challenge for planet formation theories. The gravitational interaction between young planets and their protoplanetary disks provides one way in which planetary orbits can be shaped during the formation epoch. Disk-planet interactions are strongly influenced by the structure and physical processes that drive the evolution of the protoplanetary disk. In this review we focus on how disk-planet interactions drive the migration of planets when different assumptions are made about the physics of angular momentum transport, and how it drives accretion flows in protoplanetary disk models. In particular, we consider migration in discs where: (i) accretion flows arise because turbulence diffusively transports angular momentum; (ii) laminar accretion flows are confined to thin, ionised layers near disk surfaces and are driven by the launching of magneto-centrifugal winds, with the midplane being completely inert; (iii) laminar accretion flows pervade the full column density of the disc, and are driven by a combination of large scale horizontal and vertical magnetic fields.


title: Supernova remnants in clumpy media: particle propagation and gamma-ray emission

abstract: Observations from the radio to the gamma-ray wavelengths indicate that supernova remnant (SNR) shocks are sites of effective particle acceleration. It has been proposed that the presence of dense clumps in the environment where supernovae explode might have a strong impact in shaping the hadronic gamma-ray spectrum. Here we present a detailed numerical study about the penetration of relativistic protons into clumps which are engulfed by a SNR shock, taking into account the magneto-hydrodynamical properties of the background plasma. We show that the spectrum of protons inside clumps is much harder than that in the diffuse inter-clump medium and we discuss the implications for the formation of the spectrum of hadronic gamma rays, which does not reflect anymore the acceleration spectrum of protons, resulting substantially modified inside the clumps due to propagation effects. For the Galactic SNR RX J1713.7-3946, we show that a hadronic scenario including dense clumps inside the remnant shell is able to reproduce the broadband gamma-ray spectrum from GeV to TeV energies. Moreover, we argue that small clumps crossed by the shock could provide a natural explanation to the non-thermal X-ray variability observed in some hot spots of RX J1713.7-3946. Finally we discuss the detectability of gamma-ray emission from clumps with the upcoming Cherenkov Telescope Array and the possible detection of the clumps themselves through molecular lines.


title: On the interior geometry of a typical black hole microstate

abstract: We argue that the region behind the horizon of a one-sided black hole can be probed by an analogue of the double-trace deformation protocol of Gao-Jafferis-Wall. This is achieved via a deformation of the CFT Hamiltonian by a term of the form ${\cal O} \widetilde{\cal O}$, where $\widetilde{\cal O}$ denote the state-dependent "mirror operators". We argue that this deformation creates negative energy shockwaves in the bulk, which allow particles inside the horizon to escape and to get directly detected in the CFT. This provides evidence for the smoothness of the horizon of black holes dual to typical states. We argue that the mirror operators allow us to perform an analogue of the Hayden-Preskill decoding protocol. Our claims rely on a technical conjecture about the chaotic behavior of out-of-time-order correlators on typical pure states at scrambling time.


title: Lack of null-controllability for the fractional heat equation and related equations

abstract: We consider the equation $(\partial_t + \rho(\sqrt{-\Delta}))f(t,x) = \mathbf 1_\omega u(t,x)$, $x\in \mathbb R$ or $\mathbb T$. We prove it is not null-controllable if $\rho$ is analytic on a conic neighborhood of $\mathbb R_+$ and $\rho(\xi) = o(|\xi|)$. The proof relies essentially on geometric optics, i.e.\ estimates for the evolution of semiclassical coherent states. The method also applies to other equations. The most interesting example might be the Kolmogorov-type equation $(\partial_t -\partial_v^2 + v^2\partial_x)f(t,x,v) = \mathbf 1_\omega u(t,x,v)$ for $(x,v)\in \Omega_x\times \Omega_v$ with $\Omega_x = \mathbb R$ or $\mathbb T$ and $\Omega_v = \mathbb R$ or $(-1,1)$. We prove it is not null-controllable in any time if $\omega$ is a vertical band $\omega_x\times \Omega_v$. The idea is to remark that, for some families of solutions, the Kolmogorov equation behaves like the rotated fractional heat equation $(\partial_t + \sqrt i(-\Delta)^{1/4})g(t,x) = \mathbf 1_\omega u(t,x)$, $x\in \mathbb T$.


title: Dusk Over Dawn O$_2$ Asymmetry in Europa's Near-Surface Atmosphere

abstract: The evolution of Europa's water-product exosphere over its 85-hour day, based on current models, has not been shown to exhibit any diurnal asymmetries. Here we simulate Europa's exosphere using a 3-D Monte Carlo routine including, for the first time, the role of Europa's rotation on the evolution of exospheric molecules throughout the orbit. We focus on O$_2$, sputtered by a trailing hemisphere source with a temperature-dependence under isotropic plasma conditions. We find that the O$_2$ component, while global, is not homogenous in Europa local time. Rather, the O$_2$ accumulates at the dusk hemisphere. When rotation is explicitly excluded, no diurnal asymmetries exist. We find that the assumed thermal-dependence on the O$_2$ source is critical for a diurnal asymmetry: the diurnal surface temperature profile is imprinted on to the near-surface O$_2$ atmosphere, due to the small hop times of the non-adsorbing O$_2$ effectively rotating with Europa. Altogether, our simulations conclude that the dusk-over-dawn asymmetry is driven by Europa's day-night O$_2$ cycle synchronized with Europa's orbital period based on our model assumptions on O$_2$ production and loss. This conclusion is in agreement with the recent understanding that a non-adsorbing, rotating O$_2$ source peaking at noon will naturally accumulate from dawn-to-dusk, should the O$_2$ lifetime be sufficiently long compared to the orbital period. Lastly we compare hemispherically-averaged dusk-over-dawn ratios to the recently observed oxygen emission data by the Hubble Space Telescope. We find that while the simulations are globally consistent with brighter oxygen emission at dusk than at dawn, the orbital evolution of the asymmetries in our simulations can be improved by ameliorating the O$_2$ source & loss rates, and possibly adsorption onto the regolith.


title: An analytical closed-form solution for free vibration of stepped circular/annular Mindlin functionally graded plate

abstract: An exact solution based on a unique procedure is presented for free vibration of stepped circular and annular functionally graded (FG) plates via first-order shear deformation plate theory of Mindlin. A power-law distribution of the volume fraction of the components is considered for the Young's Modulus and Poisson's ratio of the studied FG plate. Free vibration of the plate is solved by introducing some new potential functions and the use of separation of variables method. Finally, several comparisons of the developed model were presented with the FEA analysis, to demonstrate the accuracy of the proposed exact procedure. The effect of the geometrical parameters such as step thickness ratios and step locations on the natural frequencies of FG plates is also investigated.


title: End-to-end correlations in the Kitaev chain

abstract: The interdependence between long range correlations and topological signatures in fermionic arrays is examined. End-to-end correlations, in particular classical correlations, maintain a characteristic pattern in the presence of delocalized excitations and this behavior can be used as an operational criterion to identify Majorana fermions in one-dimensional systems. The study discusses how to obtain the chain eigenstates in tensor-state representation together with the proposed assessment of correlations. Outstandingly, the final result can be written as a simple analytical expression that underlines the link with the system's topological phases.


title: Coherence and decoherence in quantum absorption refrigerators

abstract: Absorption refrigerators transfer thermal energy from a cold reservoir to a hot reservoir using input energy from a third, so-called work reservoir. We examine the operation of quantum absorption refrigerators when coherences between eigenstates survive in the steady state limit. In our model, the working medium comprises a discrete, four-level system. We manifest that eigenbasis quantum coherences within this system generally suppress the cooling current in the refrigerator, while minimally affecting the coefficient of performance (cooling efficiency). We rationalize the behavior of the four-level refrigerator by studying two, three-level model systems for energy transport and refrigeration. Our calculations further illuminate the shortcomings of secular quantum master equations, and the necessity of employing dynamical equations of motion that retain couplings between population and coherences.


title: Approximating Nash Equilibria for Black-Box Games: A Bayesian Optimization Approach

abstract: Game theory has emerged as a powerful framework for modeling a large range of multi-agent scenarios. Many algorithmic solutions require discrete, finite games with payoffs that have a closed-form specification. In contrast, many real-world applications require modeling with continuous action spaces and black-box utility functions where payoff information is available only in the form of empirical (often expensive and/or noisy) observations of strategy profiles. To the best of our knowledge, few tools exist for solving the class of expensive, black-box continuous games. In this paper, we develop a method to find equilibria for such games in a sequential decision-making framework using Bayesian Optimization. The proposed approach is validated on a collection of synthetic game problems with varying degree of noise and action space dimensions. The results indicate that it is capable of improving the maximum regret in noisy and high dimensions to a greater extent than hierarchical or discretized methods.


title: An improvement of the convergence proof of the ADAM-Optimizer

abstract: A common way to train neural networks is the Backpropagation. This algorithm includes a gradient descent method, which needs an adaptive step size. In the area of neural networks, the ADAM-Optimizer is one of the most popular adaptive step size methods. It was invented in \cite{Kingma.2015} by Kingma and Ba. The $5865$ citations in only three years shows additionally the importance of the given paper. We discovered that the given convergence proof of the optimizer contains some mistakes, so that the proof will be wrong. In this paper we give an improvement to the convergence proof of the ADAM-Optimizer.


title: Green functions of conormal derivative problems for stationary Stokes system

abstract: We study Green functions for stationary Stokes systems satisfying the conormal derivative boundary condition. We establish existence, uniqueness, and various estimates for the Green function under the assumption that weak solutions of the Stokes system are continuous in the interior of the domain. Also, we establish the global pointwise bound for the Green function under the additional assumption that weak solutions of the conormal derivative problem for the Stokes system are locally bounded up to the boundary. We provide some examples satisfying such continuity and boundedness properties.


title: The origin and fate of O$_2$ in Europa's ice: an atmospheric perspective

abstract: The early prediction and subsequent detection of an O$_2$ atmosphere on Europa, coupled with the discovery that Europa has an ocean under its ice mantle, has made this moon a prime astrobiologic target, soon to be visited by the JUICE and Europa Clipper spacecraft. In spite of the considerable number of observational, modeling, and laboratory efforts, understanding the physics leading to the observed morphology of Europa's near surface O$_2$ atmosphere has been problematic. This is the case as the observed emissions depend on the local incident plasma ion flux, the local temperature and composition of the regolith, as well as on the near surface electron temperature and density. Here we rely heavily on earlier reviews briefly summarizing the observational, laboratory and simulation efforts. Although it is agreed that radiolysis of the surface ice by the incident Jovian plasma is the ultimate source of observed O$_2$, a recent, simple model of thermal desorption from a regolith permeated with O$_2$ has changed the usual paradigm. This suggests that the observed orbital dependence of the local source of the near-surface O$_2$ atmosphere is due to thermal release of O$_2$ likely trapped on the ice grains at dangling bonds with a smaller contribution due to direct sputtering. This could also impact our understanding of the suggestion that the radiolytic products in Europa's regolith might be a source of oxidants for its underground ocean.


title: Queuing Theoretic Models for Multicast and Coded-Caching in Downlink Wireless Systems

abstract: We consider a server connected to $L$ users over a shared finite capacity link. Each user is equipped with a cache. File requests at the users are generated as independent Poisson processes according to a popularity profile from a library of $M$ files. The server has access to all the files in the library. Users can store parts of the files or full files from the library in their local caches. The server should send missing parts of the files requested by the users. The server attempts to fulfill the pending requests with minimal transmissions exploiting multicasting and coding opportunities among the pending requests. We study the performance of this system in terms of queuing delays for the naive multicasting and several coded multicasting schemes proposed in the literature. We also provide approximate expressions for the mean queuing delay for these models and establish their effectiveness with simulations.


title: Quantum Algorithms for Connectivity and Related Problems

abstract: An important family of span programs, st-connectivity span programs, have been used to design quantum algorithms in various contexts, including a number of graph problems and formula evaluation problems. The complexity of the resulting algorithms depends on the largest positive witness size of any 1-input, and the largest negative witness size of any 0-input. Belovs and Reichardt first showed that the positive witness size is exactly characterized by the effective resistance of the input graph, but only rough upper bounds were known previously on the negative witness size. We show that the negative witness size in an st-connectivity span program is exactly characterized by the capacitance of the input graph. This gives a tight analysis for algorithms based on st-connectivity span programs on any set of inputs. We use this analysis to give a new quantum algorithm for estimating the capacitance of a graph. We also describe a new quantum algorithm for deciding if a graph is connected, which improves the previous best quantum algorithm for this problem if we're promised that either the graph has at least kappa > 1 components, or the graph is connected and has small average resistance, which is upper bounded by the diameter. We also give an alternative algorithm for deciding if a graph is connected that can be better than our first algorithm when the maximum degree is small. Finally, using ideas from our second connectivity algorithm, we give an algorithm for estimating the algebraic connectivity of a graph, the second largest eigenvalue of the Laplacian.


title: On the orbital variability of Ganymede's atmosphere

abstract: Ganymede's atmosphere is produced by radiative interactions with its surface, sourced by the Sun and the Jovian plasma. The sputtered and thermally desorbed molecules are tracked in our Exospheric Global Model (EGM), a 3-D parallelized collisional model. This program was developed to reconstruct the formation of the upper atmosphere/exosphere of planetary bodies interacting with solar photon flux and magnetospheric and/or the solar wind plasmas. Here, we describe the spatial distribution of the H$_2$O and O$_2$ components of Ganymede's atmosphere, and their temporal variability along Ganymede's rotation around Jupiter. In particular, we show that Ganymede's O$_2$ atmosphere is characterized by time scales of the order of Ganymede's rotational period with Jupiter's gravity being a significant driver of the spatial distribution of the heaviest exospheric components. Both the sourcing and the Jovian gravity are needed to explain some of the characteristics of the observed aurora emissions. As an example, the O$_2$ exosphere should peak at the equator with systematic maximum at the dusk equator terminator. The sputtering rate of the H$_2$O exosphere should be maximum on the leading hemisphere because of the shape of the open/close field lines boundary and displays some significant variability with longitude.