title: Sparre-Andersen identity and the last passage time

abstract: It is shown that the celebrated result of Sparre Andersen for random walks and L\'evy processes has intriguing consequences when the last time of the process in $(-\infty,0]$, say $\sigma$, is added to the picture. In the case of no positive jumps this leads to six random times, all of which have the same distribution - the uniform distribution on $[0,\sigma]$. Surprisingly, this result does not appear in the literature, even though it is based on some classical observations concerning exchangeable increments.


title: A test for monomial containment

abstract: We present an algorithm to decide whether a given ideal in the polynomial ring contains a monomial without using Gr\"obner bases, factorization or sub-resultant computations.


title: Superconductivity at 52 K in hydrogen-substituted LaFeAsO1-xHx under high pressure

abstract: The 1111-type iron-based superconductor LnFeAsO1-xFx (Ln stands for lanthanide) is the first material with a Tc above 50 K, other than cuprate superconductors. Electron doping into LaFeAsO by H, rather than F, revealed a double-dome-shaped Tc-x diagram, with a first dome (SC1, 0.05<x<0.20) and a second dome (SC2, 0.2<x<0.5). Here, we report the Tc for the whole hydrogen-doping range in LaFeAsO1-xHx under pressures of up to 19 GPa. Tc rises to 52 K at 6 GPa for the Tc-valley composition between the two Tc domes. This is the first instance of the Tc exceeding 50 K in La-1111-type iron-based superconductors. On the other hand, the Tc of SmFeAsO1-xHx decreased continually, keeping its single-dome structure up to 15 GPa. The present findings strongly suggest that the main reason for realization of the Tc >50 K observed in RE-1111 compounds (RE: Pr, Sm, and Gd) at ambient pressure is the merging of SC1 and SC2.


title: Chordal generators and the hydrodynamic normalization for the unit ball

abstract: Let $c\geq0$ and denote by $\mathcal{K}(\mathbb{H},c)$ the set of all infinitesimal generators $G:\mathbb{H}\to\mathbb{C}$ on the upper half-plane $\mathbb{H}$ such that $\limsup_{y\to\infty}y\cdot |G(iy)|\leq c.$ This class is related to univalent functions $f:\mathbb{H}\to\mathbb{H}$ with hydrodynamic normalization and appears in the so called chordal Loewner equation. In this paper, we generalize the class $\mathcal{K}(\mathbb{H},c)$ and the hydrodynamic normalization to the Euclidean unit ball in $\mathbb{C}^n$. The generalization is based on the observation that $G\in\mathcal{K}(\mathbb{H},c)$ can be characterized by an inequality for the hyperbolic length of $G(z).$


title: Small scales structures and neutrino masses

abstract: We review the impact of massive neutrinos on cosmological observables at the linear order. By means of N-body simulations we investigate the signatures left by neutrinos on the fully non-linear regime. We present the effects induced by massive neutrinos on the matter power spectrum, the halo mass function and on the halo-matter bias in massive neutrino cosmologies. We also investigate the clustering of cosmic neutrinos within galaxy clusters.


title: A multiplicity result for double singularly perturbed elliptic systems

abstract: We show that the number of low energy solutions of a double singularly perturbed Schroedinger Maxwell system type on a smooth 3 dimensional manifold (M,g) depends on the topological properties of the manifold. The result is obtained via Lusternik Schnirelmann category theory.


title: Effect of Volatility Clustering on Indifference Pricing of Options by Convex Risk Measures

abstract: In this article, we look at the effect of volatility clustering on the risk indifference price of options described by Sircar and Sturm in their paper (Sircar, R., & Sturm, S. (2012). From smile asymptotics to market risk measures. Mathematical Finance. Advance online publication. doi:10.1111/mafi.12015). The indifference price in their article is obtained by using dynamic convex risk measures given by backward stochastic differential equations. Volatility clustering is modelled by a fast mean-reverting volatility in a stochastic volatility model for stock price. Asymptotics of the indifference price of options and their corresponding implied volatility are obtained in this article, as the mean-reversion time approaches zero. Correction terms to the asymptotic option price and implied volatility are also obtained.


title: Existence and stability results on a class of Non Linear Schroedinger Equations in bounded domains with Dirichlet boundary conditions

abstract: Existence of solution and stability results on a class of Non Linear Schroedinger type equations with a bounded nonlinearity are obtained, for a bounded domain and with Dirichlet boundary conditions. The kind of stability under discussion shows that the corresponding solution exhibits features of a solitary wave type.


title: Magneto-convective models of red dwarfs: constraints imposed by the lithium abundance

abstract: Magnetic fields impede the onset of convection, thereby altering the thermal structure of a convective envelope in a low mass star: this has an effect on the amount of lithium depletion in a magnetized star. In order to quantify this effect, we have applied a magneto-convective model to two low mass stars for which lithium abundances and precise structural parameters are known: YY Gem and CU Cnc. For both stars, we have obtained models which satisfy empirical constraints on the following parameters: R, L, surface magnetic field strength, and Li abundance. In the case of YY Gem, we have obtained a model which satisfies the empirical constraints with an internal magnetic field of several megagauss: such a field strength is within the range of a dynamo where the field energy is in equipartition with rotational energy deep inside the convection zone. However, in the case of CU Cnc, the Li requires an internal magnetic field which is probably too strong for a dynamo origin: we suggest possible alternatives which might account for the reported Li abundance in CU Cnc.


title: Relativistic Lagrangians for the Lorentz-Dirac equation

abstract: We present two types of relativistic Lagrangians for the Lorentz-Dirac equation written in terms of an arbitrary world-line parameter. One of the Lagrangians contains an exponential damping function of the proper time and explicitly depends on the world-line parameter. Another Lagrangian includes additional cross-terms consisting of auxiliary dynamical variables and does not depend explicitly on the world-line parameter. We demonstrate that both the Lagrangians actually yield the Lorentz-Dirac equation with a source-like term.


title: Solving the Klein-Gordon equation using Fourier spectral methods: A benchmark test for computer performance

abstract: The cubic Klein-Gordon equation is a simple but non-trivial partial differential equation whose numerical solution has the main building blocks required for the solution of many other partial differential equations. In this study, the library 2DECOMP&FFT is used in a Fourier spectral scheme to solve the Klein-Gordon equation and strong scaling of the code is examined on thirteen different machines for a problem size of 512^3. The results are useful in assessing likely performance of other parallel fast Fourier transform based programs for solving partial differential equations. The problem is chosen to be large enough to solve on a workstation, yet also of interest to solve quickly on a supercomputer, in particular for parametric studies. Unlike other high performance computing benchmarks, for this problem size, the time to solution will not be improved by simply building a bigger supercomputer.


title: A Heuristic EDF Uplink Scheduler for Real Time Application in WiMAX Communication

abstract: WiMAX, Worldwide Interoperability for Microwave Access, is a developing wireless communication scheme that can provide broadband access to large-scale coverage. WiMAX belongs to the family of standards of IEEE-802.16. To satisfy user demands and support a new set of real time services and applications, a realistic and dynamic resource allocation algorithm is mandatory. One of the most efficient algorithm is EDF (earliest deadline first). But the problem is that when the difference between deadlines is large enough, then lower priority queues have to starve. So in this paper, we present a heuristic earliest deadline first (H-EDF) approach of the uplink scheduler of the WiMAX real time system. This H-EDF presents a way for efficient allocation of the bandwidth for uplink, so that bandwidth utilization is proper and appropriate fairness is provided to the system. We use Opnet simulator for implementing the WiMAX network, which uses this H-EDF scheduling algorithm. We will analysis the performance of the H-EDF algorithm in consideration with throughput as well as involvement of delay.


title: Noise Robustness of the Incompatibility of Quantum Measurements

abstract: The existence of incompatible measurements is a fundamental phenomenon having no explanation in classical physics. Intuitively, one considers given measurements to be incompatible within a framework of a physical theory, if their simultaneous implementation on a single physical device is prohibited by the theory itself. In the mathematical language of quantum theory, measurements are described by POVMs (positive operator valued measures), and given POVMs are by definition incompatible if they cannot be obtained via coarse-graining from a single common POVM; this notion generalizes noncommutativity of projective measurements. In quantum theory, incompatibility can be regarded as a resource necessary for manifesting phenomena such as Clauser-Horne-Shimony-Holt (CHSH) Bell inequality violations or Einstein-Podolsky-Rosen (EPR) steering which do not have classical explanation. We define operational ways of quantifying this resource via the amount of added classical noise needed to render the measurements compatible, i.e., useless as a resource. In analogy to entanglement measures, we generalize this idea by introducing the concept of incompatibility measure, which is monotone in local operations. In this paper, we restrict our consideration to binary measurements, which are already sufficient to explicitly demonstrate nontrivial features of the theory. In particular, we construct a family of incompatibility monotones operationally quantifying violations of certain scaled versions of the CHSH Bell inequality, prove that they can be computed via a semidefinite program, and show how the noise-based quantities arise as special cases. We also determine maximal violations of the new inequalities, demonstrating how Tsirelson's bound appears as a special case. The resource aspect is further motivated by simple quantum protocols where our incompatibility monotones appear as relevant figures of merit.


title: X-Ray Spectral Parameters for a Sample of 95 Active Galactic Nuclei

abstract: We present a broadband X-ray analysis of a new homogeneous sample of 95 active galactic nuclei (AGN) from the 22-month Swift/BAT all-sky survey. For this sample we treated jointly the X-ray spectra observed by XMM-Newton and INTEGRAL missions for the total spectral range of 0.5-250 keV. Photon index \Gamma, relative reflection R, equivalent width of Fe $K_{\alpha}$ line (EW Fe $K_{\alpha}$), hydrogen column density $N_{H}$, exponential cut-off energy $E_{c}$ and intrinsic luminosity $L_{corr}$ are determined for all objects of the sample. We investigated correlations \Gamma - R, EW Fe $K_{\alpha}$ - $L_{corr}$, \Gamma - $E_{c}$, EW Fe $K_{\alpha}$ - $N_{H}$. Dependence \Gamma - R for Seyfert 1 and 2 type of galaxies has been investigated separately. We found that the relative reflection parameter at low power-law indexes for Seyfert 2 galaxies is systematically higher than for Seyfert 1 ones. This can be related to an increasing contribution of the reflected radiation from the gas-dust torus. Our data show that there exists some anticorrelation between EW Fe $K_{\alpha}$ and $L_{corr}$, but it is not strong. We have not found statistically significant deviations from the AGN Unified Model.


title: Fragmentation functions of polarized heavy quarkonium

abstract: Study of the polarized heavy quarkonium production in recently proposed QCD factorization formalism requires knowledge of a large number of input fragmentation functions (FFs) from a single parton or a heavy quark-antiquark pair to a polarized heavy quarkonium. In this paper, we calculate these FFs at the input scale in terms of nonrelativistic QCD (NRQCD) factorization. We derive all relevant polarized NRQCD long-distance matrix elements based symmetries and propose a self-consistent scheme to define them in arbitrary $d$ dimensions. We then calculate polarized input FFs contributed from all $S$-wave and $P$-wave NRQCD intermediate states. With our calculation of the polarized input FFs, and the partonic hard part available in literatures, the QCD factorization formalism is ready to be applied to polarized heavy quarkonium production.


title: Precision-Aware application execution for Energy-optimization in HPC node system

abstract: Power consumption is a critical consideration in high performance computing systems and it is becoming the limiting factor to build and operate Petascale and Exascale systems. When studying the power consumption of existing systems running HPC workloads, we find that power, energy and performance are closely related which leads to the possibility to optimize energy consumption without sacrificing (much or at all) the performance. In this paper, we propose a HPC system running with a GNU/Linux OS and a Real Time Resource Manager (RTRM) that is aware and monitors the healthy of the platform. On the system, an application for disaster management runs. The application can run with different QoS depending on the situation. We defined two main situations. Normal execution, when there is no risk of a disaster, even though we still have to run the system to look ahead in the near future if the situation changes suddenly. In the second scenario, the possibilities for a disaster are very high. Then the allocation of more resources for improving the precision and the human decision has to be taken into account. The paper shows that at design time, it is possible to describe different optimal points that are going to be used at runtime by the RTOS with the application. This environment helps to the system that must run 24/7 in saving energy with the trade-off of losing precision. The paper shows a model execution which can improve the precision of results by 65% in average by increasing the number of iterations from 1e3 to 1e4. This also produces one order of magnitude longer execution time which leads to the need to use a multi-node solution. The optimal trade-off between precision vs. execution time is computed by the RTOS with the time overhead less than 10% against a native execution.


title: From complexity to algebra and back: digraph classes, collapsibility and the PGP

abstract: Inspired by computational complexity results for the quantified constraint satisfaction problem, we study the clones of idempotent polymorphisms of certain digraph classes. Our first results are two algebraic dichotomy, even "gap", theorems. Building on and extending [Martin CP'11], we prove that partially reflexive paths bequeath a set of idempotent polymorphisms whose associated clone algebra has: either the polynomially generated powers property (PGP); or the exponentially generated powers property (EGP). Similarly, we build on [DaMM ICALP'14] to prove that semicomplete digraphs have the same property. These gap theorems are further motivated by new evidence that PGP could be the algebraic explanation that a QCSP is in NP even for unbounded alternation. Along the way we also effect a study of a concrete form of PGP known as collapsibility, tying together the algebraic and structural threads from [Chen Sicomp'08], and show that collapsibility is equivalent to its $\Pi_2$-restriction. We also give a decision procedure for $k$-collapsibility from a singleton source of a finite structure (a form of collapsibility which covers all known examples of PGP for finite structures). Finally, we present a new QCSP trichotomy result, for partially reflexive paths with constants. Without constants it is known these QCSPs are either in NL or Pspace-complete [Martin CP'11], but we prove that with constants they attain the three complexities NL, NP-complete and Pspace-complete.


title: Random gauge field effects on the conductivity of graphene sheets with disordered ripples

abstract: We study the effect of disordered ripples on the conductivity of monolayer graphene flakes. We calculate the relaxation times and the Boltzmann conductivities associated with two mechanisms. First, we study the conductivity correction due to an external in-plane magnetic field $B_\parallel$. Due to the irregular local curvature found at graphene sheets deposited over a substrate, $B_\parallel$ can be mapped into an effective random magnetic field perpendicular to the graphene surface. Second, we study the electron momentum relaxation due to intrinsic pseudo magnetic fields originated from deformations and strain. We find that the competition between these mechanisms gives rise to a strong anisotropy in the conductivity tensor. This result provides a new strategy to quantitatively infer the strength of pseudo-magnetic fields in rippled graphene flakes.


title: Transductive Multi-view Zero-Shot Learning

abstract: Most existing zero-shot learning approaches exploit transfer learning via an intermediate-level semantic representation shared between an annotated auxiliary dataset and a target dataset with different classes and no annotation. A projection from a low-level feature space to the semantic representation space is learned from the auxiliary dataset and is applied without adaptation to the target dataset. In this paper we identify two inherent limitations with these approaches. First, due to having disjoint and potentially unrelated classes, the projection functions learned from the auxiliary dataset/domain are biased when applied directly to the target dataset/domain. We call this problem the projection domain shift problem and propose a novel framework, transductive multi-view embedding, to solve it. The second limitation is the prototype sparsity problem which refers to the fact that for each target class, only a single prototype is available for zero-shot learning given a semantic representation. To overcome this problem, a novel heterogeneous multi-view hypergraph label propagation method is formulated for zero-shot learning in the transductive embedding space. It effectively exploits the complementary information offered by different semantic representations and takes advantage of the manifold structures of multiple representation spaces in a coherent manner. We demonstrate through extensive experiments that the proposed approach (1) rectifies the projection shift between the auxiliary and target domains, (2) exploits the complementarity of multiple semantic representations, (3) significantly outperforms existing methods for both zero-shot and N-shot recognition on three image and video benchmark datasets, and (4) enables novel cross-view annotation tasks.


title: The combined equilibrium of business land-use and its congestion pricing principles (working paper)

abstract: This working paper is divided into two parts. Firstly, we develop a new combined equilibrium model of business land-use, which puts travelers' traffic equilibrium and business companies' competitive location equilibrium into a unified framework. A variational inequality is presented for the combined equilibrium and the properties of equilibrium solution are investigated. Secondly, the congestion pricing principles associated with the combined equilibrium are studied. From the mathematical point of view, we prove that there exists an optimal road pricing scheme that can minimize the social cost of travelers. This road pricing scheme generalizes the traditional link-based optimal road pricing scheme. Furthermore, when simultaneously imposing charges on travelers and companies is allowed, we prove that there exists an optimal congestion pricing scheme that can derive a combined equilibrium toward an overall system optimum. The economic meaning of every pricing scheme proposed in this paper is discussed in detail. At last, a simple numerical example is used to demonstrate that the optimal congestion pricing scheme may indeed reduce the social cost.


title: Planetary systems based on a quantum-like model

abstract: Planetary systems have their origin in the gravitational collapse of a cloud of gas and dust. Through a process of accretion, is formed a massive star and a disk of planetesimals orbiting the star. Using a formalism analogous to quantum mechanics (quantum-like model), the star-planetesimal system is described and the flow quantizing the gravitational field theoretical model parameters are obtained. Goodness of fit (chi-square) of the observed data with model quantum-like, to the solar system, satellites, exoplanets and protoplanetary disk around HL Tauri is determined. Shows that the radius, eccentricity, energy, angular momentum and orbital inclination of planetary objects formed take discrete values depending only on the mass star.


title: Autonomous Brownian motor driven by nonadiabatic variation of internal parameters

abstract: We discuss an autonomous motor based on a Brownian particle driven from thermal equilibrium by periodic in time variation of the internal potential through which the particle interacts with molecules of the surrounding thermal bath. We demonstrate for such a motor the absence of a linear response regime: The average driving force and drift velocity are shown to be quadratic in both the frequency and amplitude of the variation. The adiabatic approximation (of an infinitely slow variation) and the leading correction to it (linear in the variation's frequency) both lead to zero drift and are insufficient to describe the motor's operation.


title: Analyzing the Impact of Inter Cooperation Region Interference in Coordinated Multi-Point Uplink Networks

abstract: We analyze the uplink of coordinated multi-point (CoMP) networks in which cooperation can be amongst 2 or 3 base stations (BSs). We consider a 2D network of BSs on a regular hexagonal lattice wherein the cooperation tessellates the 2D plane into cooperation regions (CRs); specifically, we analyze the impact of the interference between the CRs in the network. Our model accounts realistic propagation conditions, particularly including shadowing. We obtain accurate, closed-form, approximations for the user capacity coverage probability (CCP) and the ergodic capacity at each point within the CR. To provide a network-level analysis, we focus on the locations within each CR with the minimum CCP, the worstcase point(s). The worst-case CCP and or ergodic capacity can be used in parametric studies for network design. Here, the analysis is applied to obtain the relationship between cell size and CCP and, thereby, the required density of BSs to achieve a chosen target capacity coverage. The analysis also allows for a comparison between different orders of BS cooperation, quantifying the reduced required BS density from higher orders of cooperation. Comprehensive simulations are used to illustrate the accuracy of our analysis, including the approximations used for analytic tractability.


title: Residual-driven online Generalized Multiscale Finite Element Methods

abstract: The construction of local reduced-order models via multiscale basis functions has been an area of active research. In this paper, we propose online multiscale basis functions which are constructed using the offline space and the current residual. Online multiscale basis functions are constructed adaptively in some selected regions based on our error indicators. We derive an error estimator which shows that one needs to have an offline space with certain properties to guarantee that additional online multiscale basis function will decrease the error. This error decrease is independent of physical parameters, such as the contrast and multiple scales in the problem. The offline spaces are constructed using Generalized Multiscale Finite Element Methods (GMsFEM). We show that if one chooses a sufficient number of offline basis functions, one can guarantee that additional online multiscale basis functions will reduce the error independent of contrast. We note that the construction of online basis functions is motivated by the fact that the offline space construction does not take into account distant effects. Using the residual information, we can incorporate the distant information provided the offline approximation satisfies certain properties. In the paper, theoretical and numerical results are presented. Our numerical results show that if the offline space is sufficiently large (in terms of the dimension) such that the coarse space contains all multiscale spectral basis functions that correspond to small eigenvalues, then the error reduction by adding online multiscale basis function is independent of the contrast. We discuss various ways computing online multiscale basis functions which include a use of small dimensional offline spaces.


title: Boltzmann vs Gibbs: a finite-size match

abstract: The long standing contrast between Boltzmann's and Gibbs' approach to statistical thermodynamics has been recently rekindled by Dunkel and Hilbert [1], who criticize the notion of negative absolute temperature (NAT), as a misleading consequence of Boltzmann's definition of entropy. A different definition, due to Gibbs, has been proposed, which forbids NAT and makes the energy equipartition rigorous in arbitrary sized systems. The two approaches, however, are shown to converge to the same results in the thermodynamical limit. A vigorous debate followed ref.[1], with arguments against [2,3] and in favor [4,5,6,7] of Gibbs' entropy. In an attempt to leave the speculative level and give the discussion some deal of concreteness, we analyze the practical consequences of Gibbs' definition in two finite-size systems: a non interacting gas of N atoms with two-level internal spectrum, and an Ising model of N interacting spins. It is shown that for certain measurable quantities, the difference resulting from Boltzmann's and Gibbs' approach vanishes as the inverse square rrot of N, much less rapidly than the 1/N slope expected. As shown by numerical estimates, this makes the experimental solution of the controversy a feasible task.


title: Henize 2-10: the ongoing formation of a nuclear star cluster around a massive black hole

abstract: The central region of the galaxy Henize 2-10 has a central black hole (BH) with a mass of about $2\times 10^6$ M$_\odot$. While this black hole does not appear to coincide with any central stellar over density, it is surrounded by 11 young massive clusters with masses above $10^5$ M$_\odot$. The availability of high quality data on the structure of the galaxy and the age and mass of the clusters provides excellent initial conditions for studying the dynamical evolution of Henize 2-10's nucleus. Here we present a set of $N$-body simulations of the central clusters and black hole to understand whether and how they will merge to form a nuclear star cluster. Nuclear star clusters (NSCs) are present in a majority of galaxies with stellar mass similar to Henize 2-10. Despite the results depend on the choice of initial conditions, we find that a NSC with mass $M_{NSC}\simeq 4-6\times 10^6$ M$_\odot$ and effective radius $r_{NSC}\simeq 2.6-4.1$ pc will form within $0.2$ Gyr. This work is the first showing, in a realistic realization of the host galaxy and its star cluster system, that the formation of a bright nucleus is a process that can happen after the formation of a central massive BH leading to a composite NSC+BH central system. The merging process of the clusters does not affect significantly the kinematics of the BH, whose motion, after the globular cluster merger, is limited to a $\sim 1$ pc oscillation at less than $2$ kms$^{-1}$ speed.


title: Wormholes in Horava gravity with cosmological constant

abstract: By combining analytical and numerical methods we find that the solutions of the complete Horava theory with negative cosmological constant that satisfy the conditions of staticity, spherical symmetry and vanishing of the shift function are two kinds of geometry: (i) a wormhole-like solution with two sides joined by a throat and (ii) a single side with a naked singularity at the origin. We study the second-order effective action. We consider the case when the coupling constant of the (partial ln N)^2 term, which is the unique deviation from general relativity in the effective action, is small. At one side the wormhole acquires a kind of deformed AdS asymptotia and at the other side there is an asymptotic essential singularity. The deformation of AdS essentially means that the lapse function N diverges asymptotically a bit faster than AdS. This can also be interpreted as an anisotropic Lifshitz scaling that the solutions acquire asymptotically.


title: On the contribution of the double Drell-Yan process to WW and ZZ production at the LHC

abstract: In this paper we investigate consequences of an assumption that the discrepancy of the predicted and observed W+W- production cross sections at the LHC is caused by the missing contribution of the double Drell-Yan process (DDYP). Using our simple model of DDYP of Ref. [1] we show that inclusion of this production mechanism leads to a satisfactory, parameter-free description of the two-lepton mass distribution for 0-jet W+W- events and the four-lepton mass distribution for ZZ events. In such a scenario the Higgs-boson contribution is no longer necessary to describe the data. An experimental programme to prove or falsify such an assumption is proposed.


title: Fractional Hardy-Lieb-Thirring and related inequalities for interacting systems

abstract: We prove analogues of the Lieb-Thirring and Hardy-Lieb-Thirring inequalities for many-body quantum systems with fractional kinetic operators and homogeneous interaction potentials, where no anti-symmetry on the wave functions is assumed. These many-body inequalities imply interesting one-body interpolation inequalities, and we show that the corresponding one- and many-body inequalities are actually equivalent in certain cases.


title: Local Perturbations Perturb -Exponentially- Locally

abstract: We elaborate on the principle that for gapped quantum spin systems with local interaction "local perturbations [in the Hamiltonian] perturb locally [the ground state]". This principle was established in [Bachmann et al. 2012], relying on the `spectral flow technique' or `quasi-adiabatic continuation' [Hastings 2004] to obtain locality estimates with sub-exponential decay in the distance to the spatial support of the perturbation. We use ideas of [Hamza et Al. 2009] to obtain similarly a transformation between gapped eigenvectors and their perturbations that is local with exponential decay. This allows to improve locality bounds on the effect of perturbations on the low lying states in certain gapped models with a unique `bulk ground state' or `topological quantum order'. We also give some estimate on the exponential decay of correlations in models with impurities where some relevant correlations decay faster than one would naively infer from the global gap of the system, as one also expects in disordered systems with a localized ground state.

categories: math-ph cond-mat.stat-mech math.MP quant-ph


title: k variables are needed to define k-Clique in first-order logic

abstract: In an early paper, Immerman raised a proposal on developing model-theoretic techniques to prove lower bounds on ordered structures, which represents a long-standing challenge in finite model theory. An iconic question standing for such a challenge is how many variables are needed to define $k$-Clique in first-order logic on the class of finite ordered graphs? If $k$ variables are necessary, as widely believed, it would imply that the bounded (or finite) variable hierarchy in first-order logic is strict on the class of finite ordered graphs. In 2008, Rossman made a breakthrough by establishing an optimal average-case lower bound on the size of constant-depth unbounded fan-in circuits computing $k$-Clique. In terms of logic, this means that it needs greater than $\lfloor\frac{k} {4}\rfloor$ variables to describe the $k$-Clique problem in first-order logic on the class of finite ordered graphs, even in the presence of arbitrary arithmetic predicates. It follows, with an unpublished result of Immerman, that the bounded variable hierarchy in first-order logic is indeed strict. However, Rossman's methods come from circuit complexity and a novel notion of sensitivity by himself. And the challenge before finite model theory remains there. In this paper, we give an alternative proof for the strictness of bounded variable hierarchy in $\fo$ using pure model-theoretic toolkit, and anwser the question completely for first-order logic, i.e. $k$-variables are indeed needed to describe $k$-Clique in this logic. In contrast to Rossman's proof, our proof is purely constructive. Then we embed the main structures into a pure arithmetic structure to show a similar result where arbitrary arithmetic predicates are presented. Finally, we discuss its application in circuit complexity.


title: On the stability of cycles by delayed feedback control

abstract: We present a delayed feedback control (DFC) mechanism for stabilizing cycles of one dimensional discrete time systems. In particular, we consider a delayed feedback control for stabilizing $T$-cycles of a differentiable function $f: \mathbb{R}\rightarrow\mathbb{R}$ of the form $$x(k+1) = f(x(k)) + u(k)$$ where $$u(k) = (a_1 - 1)f(x(k)) + a_2 f(x(k-T)) + ... + a_N f(x(k-(N-1)T))\;,$$ with $a_1 + ... + a_N = 1$. Following an approach of Morg\"ul, we construct a map $F: \mathbb{R}^{T+1} \rightarrow \mathbb{R}^{T+1}$ whose fixed points correspond to $T$-cycles of $f$. We then analyze the local stability of the above DFC mechanism by evaluating the stability of the corresponding equilibrum points of $F$. We associate to each periodic orbit of $f$ an explicit polynomial whose Schur stability corresponds to the stability of the DFC on that orbit. An example indicating the efficacy of this method is provided.


title: Probing relevant ingredients in mean-field approaches for the athermal rheology of yield stress materials

abstract: Although the notion of mechanical noise is expected to play a key role in the non-linear rheology of athermally sheared amorphous systems, its characterization has so far remained elusive. Here, we show using molecular dynamic simulations that in spite of the presence of strong spatio-temporal correlations in the system, the local stress exhibits normal diffusion under the effect of the mechanical noise in the finite driving regime. The diffusion constant appears to be proportional to the mean plastic activity. Our data suggests that the corresponding proportionality constant is density independent, and can be directly related to the specific form of the rheological flow curve, pointing the way to a generic way of modeling mechanical noise in mean-field equations.


title: An optimal trading problem in intraday electricity markets

abstract: We consider the problem of optimal trading for a power producer in the context of intraday electricity markets. The aim is to minimize the imbalance cost induced by the random residual demand in electricity, i.e. the consumption from the clients minus the production from renewable energy. For a simple linear price impact model and a quadratic criterion, we explicitly obtain approximate optimal strategies in the intraday market and thermal power generation, and exhibit some remarkable properties of the trading rate. Furthermore, we study the case when there are jumps on the demand forecast and on the intraday price, typically due to error in the prediction of wind power generation. Finally, we solve the problem when taking into account delay constraints in thermal power production.


title: Rotationally symmetric biharmonic maps between models

abstract: The main aim of this paper is to study existence and stability properties of rotationally symmetric proper biharmonic maps between two $m$-dimensional models (in the sense of Greene and Wu). We obtain a complete classification of rotationally symmetric, proper biharmonic conformal diffeomorphisms in the special case that $m=4$ and the models have constant sectional curvature. Then, by introducing the Hamiltonian associated to this problem, we also obtain a complete description of conformal proper biharmonic solutions in the case that the domain model is ${\mathbb R}^4$. In the second part of the paper we carry out a stability study with respect to equivariant variations (equivariant stability). In particular, we prove that: (i) the inverse of the stereographic projection from the open $4$-dimensional Euclidean ball to the hyperbolic space is equivariant stable; (ii) the inverse of the stereographic projection from the closed $4$-dimensional Euclidean ball to the sphere is equivariant stable with respect to variations which preserve the boundary data.


title: Theory of Dipole Induced Electromagnetic Transparency

abstract: A detailed theory describing linear optics of vapors comprised of interacting multi-level quantum emitters is proposed. It is shown both by direct integration of Maxwell-Bloch equations and using a simple analytical model that at large densities narrow transparency windows appear in otherwise completely opaque spectra. The existence of such windows is attributed to overlapping resonances. This effect, first introduced for three-level systems in [R. Puthumpally-Joseph, M. Sukharev, O. Atabek and E. Charron, Phys. Rev. Lett. 113, 163603 (2014)], is due to strongly enhanced dipole-dipole interactions at high emitters' densities. The presented theory extends this effect to the case of multilevel systems. The theory is applied to the D1 transitions of interacting Rb-85 atoms. It is shown that at high atomic densities, Rb-85 atoms can behave as three-level emitters exhibiting all the properties of dipole induced electromagnetic transparency. Applications including slow light and laser pulse shaping are also proposed.


title: Stabilization of polynomial dynamical systems using linear programming based on Bernstein polynomials

abstract: In this paper, we deal with the problem of synthesizing static output feedback controllers for stabilizing polynomial systems. Our approach jointly synthesizes a Lyapunov function and a static output feedback controller that stabilizes the system over a given subset of the state-space. Specifically, our approach is simultaneously targeted towards two goals: (a) asymptotic Lyapunov stability of the system, and (b) invariance of a box containing the equilibrium. Our approach uses Bernstein polynomials to build a linear relaxation of polynomial optimization problems, and the use of a so-called "policy iteration" approach to deal with bilinear optimization problems. Our approach can be naturally extended to synthesizing hybrid feedback control laws through a combination of state-space decomposition and Bernstein polynomials. We demonstrate the effectiveness of our approach on a series of numerical benchmark examples.


title: Stability of Influence Maximization

abstract: The present article serves as an erratum to our paper of the same title, which was presented and published in the KDD 2014 conference. In that article, we claimed falsely that the objective function defined in Section 1.4 is non-monotone submodular. We are deeply indebted to Debmalya Mandal, Jean Pouget-Abadie and Yaron Singer for bringing to our attention a counter-example to that claim. Subsequent to becoming aware of the counter-example, we have shown that the objective function is in fact NP-hard to approximate to within a factor of $O(n^{1-\epsilon})$ for any $\epsilon > 0$. In an attempt to fix the record, the present article combines the problem motivation, models, and experimental results sections from the original incorrect article with the new hardness result. We would like readers to only cite and use this version (which will remain an unpublished note) instead of the incorrect conference version.


title: Thermalization with chemical potentials, and higher spin black holes

abstract: We study the long time behaviour of local observables following a quantum quench in 1+1 dimensional conformal field theories possessing additional conserved charges besides the energy. We show that the expectation value of an arbitrary string of {\it local} observables supported on a finite interval exponentially approaches an equilibrium value. The equilibrium is characterized by a temperature and chemical potentials defined in terms of the quenched state. For an infinite number of commuting conserved charges, the equilibrium ensemble is a generalized Gibbs ensemble (GGE). We compute the thermalization rate in a systematic perturbation in the chemical potentials, using a new technique to sum over an infinite number of Feynman diagrams. The above technique also allows us to compute relaxation times for thermal Green's functions in the presence of an arbitrary number of chemical potentials. In the context of a higher spin (hs[\lambda]) holography, the partition function of the final equilibrium GGE is known to agree with that of a higher spin black hole. The thermalization rate from the CFT computed in our paper agrees with the quasinormal frequency of a scalar field in this black hole.


title: Localization of charge carriers in monolayer graphene gradually disordered by ion irradiation

abstract: Gradual localization of charge carriers was studied in a series of micro-size samples of monolayer graphene fabricated on the common large scale film and irradiated by different doses of C$^+$ ions with energy 35 keV. Measurements of the temperature dependence of conductivity and magnetoresistance in fields up to 4 T showed that at low disorder, the samples are in the regime of weak localization and antilocalization. Further increase of disorder leads to strong localization regime, when conductivity is described by the variable-range-hopping (VRH) mechanism. A crossover from the Mott regime to the Efros-Shklovskii regime of VRH is observed with decreasing temperature. Theoretical analysis of conductivity in both regimes showed a remarkably good agreement with experimental data.


title: Exquisitely sensitive seal whisker-like sensors detect wakes at large distances

abstract: Blindfolded harbor seals are able to use their uniquely shaped whiskers to track vortex wakes left by moving animals and objects that passed by up to 30 seconds earlier; this is an impressive feat as the flow features they detect may have velocity as low as 1 mm/s, and the seals have some capacity to identify the shape of the object as well. They do so while swimming forward at high speed, hence their whiskers are sensitive enough to detect small-scale changes in the external flow field, while rejecting self-generated flow noise. Here we identify and illustrate a novel flow mechanism that allows artificial whiskers with the identical unique geometry as those of the harbor seal to detect the features of minute flow fluctuations in wakes produced by objects far away. This is shown through the study of a model problem, consisting of a harbor seal whisker model interacting with the wake of an upstream circular cylinder. We show that whereas in open water the whisker geometry results in very low vibration, once it enters a wake it oscillates with large amplitude and, remarkably, its response frequency coincides with the Strouhal frequency of the upstream cylinder, thus making the detection of an upstream wake as well as an estimation of the size and shape of the wake-generating body possible. An energy flow extraction mechanism causes the large amplitude whisker oscillations to lock in to the frequency of the oncoming wake, characterized by a slaloming motion among the oncoming wake vortices. This passive mechanism has some similarities with the flow mechanisms observed in actively controlled propulsive foils within upstream wakes and trout swimming behind bluff cylinders in a stream, but also differences due to the remarkable whisker morphology which causes it to operate passively and within a much wider parametric range.


title: An imbedding of spacetimes

abstract: It is shown that any two-dimensional spacetimes with compact Cauchy surfaces can be causally isomorphically imbedded into the two-dimensional Einstein's static universe. Also, it is shown that any two-dimensional globally hyperbolic spacetimes are conformally equivalent to a subset of the two-dimensional Einstein's static universe.


title: Entropy and Complexity of Polygonal Billiards with Spy Mirrors

abstract: We prove that a polygonal billiard with one-sided mirrors has zero topological entropy. In certain cases we show sub exponential and for other polynomial estimates on the complexity.


title: Multiplicative functions in short intervals

abstract: We introduce a general result relating "short averages" of a multiplicative function to "long averages" which are well understood. This result has several consequences. First, for the M\"obius function we show that there are cancellations in the sum of $\mu(n)$ in almost all intervals of the form $[x, x + \psi(x)]$ with $\psi(x) \rightarrow \infty$ arbitrarily slowly. This goes beyond what was previously known conditionally on the Density Hypothesis or the stronger Riemann Hypothesis. Second, we settle the long-standing conjecture on the existence of $x^{\epsilon}$-smooth numbers in intervals of the form $[x, x + c(\varepsilon) \sqrt{x}]$, recovering unconditionally a conditional (on the Riemann Hypothesis) result of Soundararajan. Third, we show that the mean-value of $\lambda(n)\lambda(n+1)$, with $\lambda(n)$ Liouville's function, is non-trivially bounded in absolute value by $1 - \delta$ for some $\delta > 0$. This settles an old folklore conjecture and constitutes progress towards Chowla's conjecture. Fourth, we show that a (general) real-valued multiplicative function $f$ has a positive proportion of sign changes if and only if $f$ is negative on at least one integer and non-zero on a positive proportion of the integers. This improves on many previous works, and is new already in the case of the M\"obius function. We also obtain some additional results on smooth numbers in almost all intervals, and sign changes of multiplicative functions in all intervals of square-root length.


title: Wigner translations and the observer-dependence of the position of masslesss spinning particles

abstract: The Wigner little group for massless particles is isomorphic to the Euclidean group ${\rm SE}(2)$. Applied to momentum eigenstates, or to infinite plane waves, the Euclidean "Wigner translations" act as the identity. We show that when applied to finite wavepackets the translation generators move the packet trajectory parallel to itself through a distance proportional to the particle's helicity. We relate this effect to the Hall effect of light and to the Lorentz-frame dependence of the position of a massless spinning particle.


title: Transferring Rich Feature Hierarchies for Robust Visual Tracking

abstract: Convolutional neural network (CNN) models have demonstrated great success in various computer vision tasks including image classification and object detection. However, some equally important tasks such as visual tracking remain relatively unexplored. We believe that a major hurdle that hinders the application of CNN to visual tracking is the lack of properly labeled training data. While existing applications that liberate the power of CNN often need an enormous amount of training data in the order of millions, visual tracking applications typically have only one labeled example in the first frame of each video. We address this research issue here by pre-training a CNN offline and then transferring the rich feature hierarchies learned to online tracking. The CNN is also fine-tuned during online tracking to adapt to the appearance of the tracked target specified in the first video frame. To fit the characteristics of object tracking, we first pre-train the CNN to recognize what is an object, and then propose to generate a probability map instead of producing a simple class label. Using two challenging open benchmarks for performance evaluation, our proposed tracker has demonstrated substantial improvement over other state-of-the-art trackers.


title: Correlation Trends in the Hyperfine Structures of $^{210,212}$Fr

abstract: We demonstrate the importance of electron correlation effects in the hyperfine structure constants of many low-lying states in $^{210}$Fr and $^{212}$Fr. This is achieved by calculating the magnetic dipole and electric quadrupole hyperfine structure constants using the Dirac-Fock approximation, second order many-body perturbation theory and the coupled-cluster method in the singles and doubles approximation in the relativistic framework. By combining our recommended theoretical results with the corresponding experimental values, improved nuclear magnetic dipole and electric quadrupole moments of the above isotopes are determined. In the present work, it is observed that there are large discrepancies between the hyperfine structure constants of the $7D_{5/2}$ state obtained from the experimental and theoretical studies whereas good agreements are found for the other $D_{5/2}$ states. Our estimated hyperfine constants for the $8P$, $6D$, $10S$ and $11S$ states could be very useful as benchmarks for the measurement of these quantities.


title: Zeeman deceleration of electron-impact-excited metastable helium atoms

abstract: We present experimental results that demonstrate - for the first time - the Zeeman deceleration of helium atoms in the metastable 2^3S_1state. A more than 40% decrease of the kinetic energy of the beam is achieved for deceleration from 490 m/s to a final velocity of 370 m/s. Metastable atom generation is achieved with an electron-impact-excitation source whose performance is enhanced through an additional discharge-type process which we characterize in detail. Comparison of deceleration data at different electron beam pulse durations confirms that a matching between the initial particle distribution and the phase-space acceptance of the decelerator is crucial for the production of a decelerated packet with a well-defined velocity distribution. The experimental findings are in good agreement with three-dimensional numerical particle trajectory simulations.


title: Nonlinear Model Reduction for Complex Systems using Sparse Optimal Sensor Locations from Learned Nonlinear Libraries

abstract: We demonstrate the synthesis of sparse sampling and machine learning to characterize and model complex, nonlinear dynamical systems over a range of bifurcation parameters. First, we construct modal libraries using the classical proper orthogonal decomposition to uncover dominant low-rank coherent structures. Here, nonlinear libraries are also constructed in order to take advantage of the discrete empirical interpolation method and projection that allows for the approximation of nonlinear terms in a low-dimensional way. The selected sampling points are shown to be nearly optimal sensing locations for characterizing the underlying dynamics, stability, and bifurcations of complex systems. The use of empirical interpolation points and sparse representation facilitate a family of local reduced-order models for each physical regime, rather than a higher-order global model, which has the benefit of physical interpretability of energy transfer between coherent structures. In particular, the discrete interpolation points and nonlinear modal libraries are used for sparse representation to classify the dynamic bifurcation regime in the complex Ginzburg-Landau equation. It is shown that nonlinear point measurements are more effective than linear measurements when sensor noise is present.


title: Rapidly rotating neutron stars in $R$-squared gravity

abstract: $f(R)$ theories of gravity are one of the most popular alternative explanations for dark energy and therefore studying the possible astrophysical implications of these theories is an important task. In the present paper we make a substantial advance in this direction by considering rapidly rotating neutron stars in $R^2$ gravity. The results are obtained numerically and the method we use is non-perturbative and self-consistent. The neutron star properties, such as mass, radius and moment of inertia, are studied in detail and the results show that rotation magnifies the deviations from general relativity and the maximum mass and moment of inertia can reach very high values. This observation is similar to previous studies of rapidly rotating neutron stars in other alternative theories of gravity, such as the scalar-tensor theories, and it can potentially lead to strong astrophysical manifestations.


title: Near-linear constructions of exact unitary 2-designs

abstract: A unitary 2-design can be viewed as a quantum analogue of a 2-universal hash function: it is indistinguishable from a truly random unitary by any procedure that queries it twice. We show that exact unitary 2-designs on n qubits can be implemented by quantum circuits consisting of ~O(n) elementary gates in logarithmic depth. This is essentially a quadratic improvement in size (and in width times depth) over all previous implementations that are exact or approximate (for sufficiently strong approximations).


title: Black hole spectra in holography: consequences for equilibration of dual gauge theories

abstract: For a closed system to equilibrate from a given initial condition there must exist an equilibrium state with the energy equal to the initial one. Equilibrium states of a strongly coupled gauge theory with a gravitational holographic dual are represented by black holes. We study the spectrum of black holes in Pilch-Warner geometry. These black holes are holographically dual to equilibrium states of strongly coupled $SU(N)$ ${\cal N}=2^*$ gauge theory plasma on $S^3$ in the planar limit. We find that there is no energy gap in the black hole spectrum. Thus, there is a priory no obstruction for equilibration of arbitrary low-energy states in the theory via a small black hole gravitational collapse. The latter is contrasted with phenomenological examples of holography with dual four-dimensional CFTs having non-equal central charges in the stress-energy tensor trace anomaly.


title: Optical observations of PSR J2021+3651 in the Dragonfly Nebula with the GTC

abstract: PSR J2021+3651 is a 17 kyr old rotation powered pulsar detected in the radio, X-rays, and $\gamma$-rays. It powers a torus-like pulsar wind nebula with jets, dubbed the Dragonfly, which is very similar to that of the Vela pulsar. The Dragonfly is likely associated with the extended TeV source VER J2019+368 and extended radio emission. We conducted first deep optical observations with the GTC in the Sloan $r'$ band to search for optical counterparts of the pulsar and its nebula. No counterparts were detected down to $r'\gtrsim27.2$ and $\gtrsim24.8$ for the point-like pulsar and the compact X-ray nebula, respectively. We also reanalyzed Chandra archival X-ray data taking into account an interstellar extinction--distance relation, constructed by us for the Dragonfly line of sight using the red-clump stars as standard candles. This allowed us to constrain the distance to the pulsar, $D=1.8^{+1.7}_{-1.4}$ kpc at 90% confidence. It is much smaller than the dispersion measure distance of $\sim$12 kpc but compatible with a $\gamma$-ray "pseudo-distance" of 1 kpc. Based on that and the optical upper limits, we conclude that PSR J2021+3651, similar to the Vela pulsar, is a very inefficient nonthermal emitter in the optical and X-rays, while its $\gamma$-ray efficiency is consistent with an average efficiency for $\gamma$-pulsars of similar age. Our optical flux upper limit for the pulsar is consistent with the long-wavelength extrapolation of its X-ray spectrum while the nebula flux upper limit does not constrain the respective extrapolation.


title: Asymptotics for the heat kernel in multicone domains

abstract: A multi cone domain $\Omega \subseteq \mathbb{R}^n$ is an open, connected set that resembles a finite collection of cones far away from the origin. We study the rate of decay in time of the heat kernel $p(t,x,y)$ of a Brownian motion killed upon exiting $\Omega$, using both probabilistic and analytical techniques. We find that the decay is polynomial and we characterize $\lim_{t\to\infty} t^{1+\alpha}p(t,x,y)$ in terms of the Martin boundary of $\Omega$ at infinity, where $\alpha>0$ depends on the geometry of $\Omega$. We next derive an analogous result for $t^{\kappa/2}\mathbb{P}_x(T >t)$, with $\kappa = 1+\alpha - n/2$, where $T$ is the exit time form $\Omega$. Lastly, we deduce the renormalized Yaglom limit for the process conditioned on survival.


title: Non-Global Logarithms, Factorization, and the Soft Substructure of Jets

abstract: An outstanding problem in QCD and jet physics is the factorization and resummation of logarithms that arise due to phase space constraints, so-called non-global logarithms (NGLs). In this paper, we show that NGLs can be factorized and resummed down to an unresolved infrared scale by making sufficiently many measurements on a jet or other restricted phase space region. Resummation is accomplished by renormalization group evolution of the objects in the factorization theorem and anomalous dimensions can be calculated to any perturbative accuracy and with any number of colors. To connect with the NGLs of more inclusive measurements, we present a novel perturbative expansion which is controlled by the volume of the allowed phase space for unresolved emissions. Arbitrary accuracy can be obtained by making more and more measurements so to resolve lower and lower scales. We find that even a minimal number of measurements produces agreement with Monte Carlo methods for leading-logarithmic resummation of NGLs at the sub-percent level over the full dynamical range relevant for the Large Hadron Collider. We also discuss other applications of our factorization theorem to soft jet dynamics and how to extend to higher-order accuracy.


title: Key (critical) relations preserved by a weak near-unanimity function

abstract: In the paper we introduce a notion of a key relation, which is similar to the notion of a critical relation introduced by Keith A.Kearnes and \'Agnes Szendrei. All clones on finite sets can be defined by only key relations. In addition there is a nice description of all key relations on 2 elements. These are exactly the relations that can be defined as a disjunction of linear equations. In the paper we show that, in general key relations do not have such a nice description. Nevertheless, we obtain a nice characterization of all key relations preserved by a weak near-unanimity function. This characterization is presented in the paper.


title: Proof of linear instability of the Reissner-Nordstr\"om Cauchy horizon under scalar perturbations

abstract: It has long been suggested that solutions to linear scalar wave equation $$\Box_g\phi=0$$ on a fixed subextremal Reissner-Nordstr\"om spacetime with non-vanishing charge are generically singular at the Cauchy horizon. We prove that generic smooth and compactly supported initial data on a Cauchy hypersurface indeed give rise to solutions with infinite nondegenerate energy near the Cauchy horizon in the interior of the black hole. In particular, the solution generically does not belong to $W^{1,2}_{loc}$. This instability is related to the celebrated blue shift effect in the interior of the black hole. The problem is motivated by the strong cosmic censorship conjecture and it is expected that for the full nonlinear Einstein-Maxwell system, this instability leads to a singular Cauchy horizon for generic small perturbations of Reissner-Nordstr\"om spacetime. Moreover, in addition to the instability result, we also show as a consequence of the proof that Price's law decay is generically sharp along the event horizon.


title: Deep Imaging of M51: a New View of the Whirlpool's Extended Tidal Debris

abstract: We present deep, wide-field imaging of the M51 system using CWRU's Burrell Schmidt telescope at KPNO to study the faint tidal features that constrain its interaction history. Our images trace M51's tidal morphology down to a limiting surface brightness of $\mu_{B,lim}\sim $30 mag arcsec$^{-2}$, and provide accurate colors ($\sigma_{B-V} < 0.1$) down to $\mu_B\sim 28$. We identify two new tidal streams in the system (the South and Northeast Plumes) with surface brightnesses of $\mu_B =29$ and luminosities of $\sim 10^6 L_{\odot,B}$. While the Northeast Plume may be a faint outer extension of the tidal "crown" north of NGC 5195 (M51b), the South Plume has no analogue in any existing M51 simulation and may represent a distinct tidal stream or disrupted dwarf galaxy. We also trace the extremely diffuse Northwest Plume out to a total extent of 20' (43 kpc) from NGC 5194(M51a), and show it to be physically distinct from the overlapping bright tidal streams from M51b. The Northwest Plume's morphology and red color ($B-V=0.8$) instead argue that it originated from tidal stripping of M51a's extreme outer disk. Finally, we confirm the strong segregation of gas and stars in the Southeast Tail, and do not detect any diffuse stellar component in the HI portion of the tail. Extant simulations of M51 have difficulty matching both the wealth of tidal structure in the system and the lack of stars in the HI tail, motivating new modeling campaigns to study the dynamical evolution of this classic interacting system.


title: An explicit open image theorem for products of elliptic curves

abstract: Let $K$ be a number field and $E_1, \ldots, E_n$ be elliptic curves over $K$, pairwise non-isogenous over $\overline{K}$ and without complex multiplication over $\overline{K}$. We study the image of the adelic representation of the absolute Galois group of $K$ naturally attached to $E_1 \times \cdots \times E_n$. The main result is an explicit bound for the index of this image in $\left\{ (x_1,\ldots,x_n) \in \operatorname{GL}_2(\hat{\mathbb{Z}})^n \bigm\vert \det x_i = \det x_j \;\; \forall i,j \right\}$.


title: Hilbert space inner products for PT-symmetric Su-Schrieffer-Heeger models

abstract: A Su-Schrieffer-Heeger model with added PT-symmetric boundary term is studied in the framework of pseudo-hermitian quantum mechanics. For two special cases, a complete set of pseudometrics is constructed in closed form. When complemented with a condition of positivity, the pseudometrics determine all the physical inner products of the considered model.


title: Triangular Spin-Orbit-Coupled Lattice with Strong Coulomb Correlations: Sn Atoms on a SiC(0001) Substrate

abstract: Two-dimensional (2D) atom lattices provide model setups for Coulomb correlations inducing competing ground states, partly with topological character. Hexagonal SiC(0001) is an intriguing wide-gap substrate, spectroscopically separated from the overlayer and hence reduced screening. We report the first study of an artificial high-Z atom lattice on SiC(0001) by Sn adatoms, based on combined experimental realization and theoretical modeling. Density-functional theory of our $\sqrt{3}$-structure model closely reproduces the scanning tunneling microscopy. Instead of metallic behavior, photoemission data show a deeply gapped state (~2 eV gap). Based on our calculations including dynamic mean-field theory, we argue that this reflects a pronounced Mott insulating scenario. We also find indications that the system is susceptible to antiferromagnetic superstructures. Such spin-orbit-coupled correlated heavy atom lattices on SiC(0001) thus form a novel testbed for peculiar quantum states of matter, with potential bearing for spin liquids and topological Mott insulators.


title: Single-stage reconstruction algorithm for quantitative photoacoustic tomography

abstract: The development of efficient and accurate image reconstruction algorithms is one of the cornerstones of computed tomography. Existing algorithms for quantitative photoacoustic tomography currently operate in a two-stage procedure: First an inverse source problem for the acoustic wave propagation is solved, whereas in a second step the optical parameters are estimated from the result of the first step. Such an approach has several drawbacks. In this paper we therefore propose the use of single-stage reconstruction algorithms for quantitative photoacoustic tomography, where the optical parameters are directly reconstructed from the observed acoustical data. In that context we formulate the image reconstruction problem of quantitative photoacoustic tomography as a single nonlinear inverse problem by coupling the radiative transfer equation with the acoustic wave equation. The inverse problem is approached by Tikhonov regularization with a convex penalty in combination with the proximal gradient iteration for minimizing the Tikhonov functional. We present numerical results, where the proposed single-stage algorithm shows an improved reconstruction quality at a similar computational cost.


title: Classification of Arnold-Beltrami Flows and their Hidden Symmetries

abstract: In the context of mathematical hydrodynamics, we consider the group theory structure which underlies the ABC-flow introduced by Beltrami, Arnold and Childress. Beltrami equation is the eigenstate equation for the first order Laplace-Beltrami operator *d, which we solve by using harmonic analysis. Taking torus T^3 constructed as R^3/L, where L is a crystallographic lattice, we present a general algorithm to construct solutions of Beltrami equation which utilizes as main ingredient the orbits under the action of the point group P_L of three-vectors in the momentum lattice L*. We introduce the new notion of a Universal Classifying Group GU_L which contains all crystallographic space groups as proper subgroups. We show that the *d-eigenfunctions are naturally arranged into irreducible representations of GU_L and by means of a systematic use of the branching rules with respect to various possible subgroups H of GU_L we search and find Beltrami fields with non trivial hidden symmetries. In the case of the cubic lattice the point group P_L is the proper octahedral group O_24 and the Universal Classifying Group is finite group G_1536 of order 1536 which we study in full detail deriving all of its 37 irreducible representations and the associated character table. We show that the O_24 orbits in the cubic lattice are arranged into 48 equivalence classes, the parameters of the corresponding Beltrami vector fields filling all the 37 irreducible representations of G_1536. In this way we obtain an exhaustive classification of all generalized ABC-flows and of their hidden symmetries. We make several conceptual comments about the possible relation of Arnold-Beltrami flows with (supersymmetric) Chern-Simons gauge theories. We also suggest linear generalizations of Beltrami equation to higher odd-dimensions that possibly make contact with M-theory and the geometry of flux-compactifications.

categories: math-ph hep-th math.DS math.GR math.MP nlin.CD


title: Statistical Problems in a Paper on Variation In Cancer Risk Among Tissues, and New Discoveries

abstract: Tomasetti and Vogelstein (2015) collected data on 31 different tissue types and found a correlation of 0.8 between the logarithms of the incidence of cancer (LCI), and the estimated number of stem cell divisions in those tissues (LSCD). Some of their conclusions however are statistically erroneous. Their excess risk score, "ERS" (log10 LCI x log10 LSCD), is non-monotonic under a change of time units for the rates, which renders meaningless the results derived from it, including a cluster of 22 "R-tumor" types for which they conclude, "primary prevention measures are not likely to be very effective". Further, r = 0.8 is consistent with the three orders of magnitude variation in other unmeasured factors, leaving room for the possibility of primary prevention if such factors can be intervened upon. Further exploration of the data reveals additional findings: (1) that LCI grows at approximately the square root of LSCD, which may provide a clue to the biology; (2) among different possible combinations of the primary data, the one maximizing the correlations with LCI is almost precisely the formula used by Tomasetti and Vogelstein to estimate LSCD, giving support to stem cell divisions as an independent factor in carcinogenesis, while not excluding other such factors.


title: No Open or Flat Bouncing Cosmologies in Einstein Gravity

abstract: We show that bouncing open or flat Friedmann-Robertson-Walker cosmologies are inconsistent with worldsheet string theory to first approximation. Specifically, the Virasoro constraint translates to the null energy condition in spacetime at leading order in the alpha-prime expansion. Thus one must go beyond minimally-coupled Einstein gravity in order to find bounce solutions.


title: Electronic correlations in Hund metals

abstract: To clarify the nature of correlations in Hund metals and its relationship with Mott physics we analyze the electronic correlations in multiorbital systems as a function of intraorbital interaction U, Hund's coupling JH and electronic filling n. We show that the main process behind the enhancement of correlations in Hund metals is the suppression of the double-occupancy of a given orbital, as it also happens in the Mott-insulator at half-filling. However, contrary to what happens in Mott correlated states the reduction of the quasiparticle weight Z with JH can happen on spite of increasing charge fluctuations. Therefore, in Hund metals the quasiparticle weight and the mass enhancement are not good measurements of the charge localization. Using simple energetic arguments we explain why the spin polarization induced by Hund's coupling produces orbital decoupling. We also discuss how the behavior at moderate interactions, with correlations controlled by the atomic spin polarization, changes at large $U$ and $J_H$ due to the proximity to a Mott insulating state.


title: Non-equilibrium thermodynamical framework for rate- and state-dependent friction

abstract: Rate- and state-dependent friction law for velocity-step and healing are analysed from a thermodynamic point of view. Assuming a logarithmic deviation from steady-state a unification of the classical Dieterich and Ruina models of rock friction is proposed.


title: New index transforms with the product of Bessel functions

abstract: New index transforms are investigated, which contain as the kernel products of the Bessel and modified Bessel functions. Mapping properties and invertibility in Lebesgue spaces are studied for these operators. Relationships with the Kontorovich-Lebedev and Fourier cosine transforms are established. Inversion theorems are proved. As an application, a solution of the initial value problem for the fourth order partial differential equation, involving the Laplacian is presented.


title: BiLipschitz decomposition of Lipschitz maps between Carnot groups

abstract: Let $f : G \to H$ be a Lipschitz map between two Carnot groups. We show that if $B$ is ball of $G$, then there exists a subset $Z \subset B$, whose image in $H$ under $f$ has small Hausdorff content, such that $B \backslash Z$ can be decomposed into a controlled number of pieces, the restriction of $f$ on each of which is quantitatively biLipschitz. This extends a result of \cite{meyerson}, which proved the same result, but with the restriction that $G$ has an appropriate discretization. We provide an example of a Carnot group not admitting such a discretization.


title: The crystalline period of a height one $p$-adic dynamical system over $\mathbf{Z}_p$

abstract: Let $f$ be a continuous ring endomorphism of $\mathbf{Z}_p[[x]]/\mathbf{Z}_p$ of degree $p.$ We prove that if $f$ acts on the tangent space at $0$ by a uniformizer and commutes with an automorphism of infinite order, then it is necessarily an endomorphism of a formal group over $\mathbf{Z}_p.$ The proof relies on finding a stable embedding of $\mathbf{Z}_p[[x]]$ in Fontaine's crystalline period ring with the property that $f$ appears in the monoid of endomorphisms generated by the Galois group of $\mathbf{Q}_p$ and crystalline Frobenius. Our result verifies, over $\mathbf{Z}_p,$ the height one case of a conjecture by Lubin.


title: Average-Value Tverberg Partitions via Finite Fourier Analysis

abstract: The long-standing topological Tverberg conjecture claimed, for any continuous map from the boundary of an $N(q,d):=(q-1)(d+1)$-simplex to $d$-dimensional Euclidian space, the existence of $q$ pairwise disjoint subfaces whose images have non-empty $q$-fold intersection. The affine cases, true for all $q$, constitute Tverberg's famous 1966 generalization of the classical Radon's Theorem. Although established for all prime powers in 1987 by \"Ozaydin, counterexamples to the conjecture, relying on 2014 work of Mabillard and Wagner, were first shown to exist for all non-prime-powers in 2015 by Frick. Starting with a reformulation of the topological Tverberg conjecture in terms of harmonic analysis on finite groups, we show that despite the failure of the conjecture, continuous maps \textit{below} the tight dimension $N(q,d)$ are nonetheless guaranteed $q$ pairwise disjoint subfaces -- including when $q$ is not a prime power -- which satisfy a variety of "average value" coincidences, the latter obtained as the vanishing of prescribed Fourier transforms.


title: Fra\"iss\'e Structures with Universal Automorphism Groups

abstract: We prove that the automorphism group of a Fra\"iss\'e structure M equipped with a notion of stationary independence is universal for the class of automorphism groups of substructures of M. Furthermore, we show that this applies to certain homogeneous n-gons.


title: Knot Cabling and the Degree of the Colored Jones Polynomial II

abstract: We continue our study of the degree of the colored Jones polynomial under knot cabling started in "Knot Cabling and the Degree of the Colored Jones Polynomial" (arXiv:1501.01574). Under certain hypothesis on this degree, we determine how the Jones slopes and the linear term behave under cabling. As an application we verify Garoufalidis' Slope Conjecture and a conjecture of the authors for cables of a two-parameter family of closed 3-braids called 2-fusion knots.


title: Singular Values Distribution of Squares of Elliptic Random Matrices and Type B Narayana Polynomials

abstract: We consider Gaussian elliptic random matrices $X$ of a size $N \times N$ with parameter $\rho$, i.e., matrices whose pairs of entries $(X_{ij}, X_{ji})$ are mutually independent Gaussian vectors, $E X_{ij} = 0$, $E X^2_{ij} = 1$ and $E X_{ij} X_{ji} = \rho$. We are interested in the asymptotic distribution of eigenvalues of the matrix $W =\frac{1}{N^2} X^2 X^{*2}$. We have shown that this distribution is defined by its moments and we provide a recurrent relation for these moments. We have proven that the (symmetrized) asymptotic distribution is determined by its free cumulants, which are Narayana polynomials of type B: $$c_{2n} = \sum_{k=0}^n \binom{n}{k}^2 \rho^{2k}.$$


title: Discretization of div-curl Systems by Weak Galerkin Finite Element Methods on Polyhedral Partitions

abstract: In this paper, the authors devise a new discretization scheme for div-curl systems defined in connected domains with heterogeneous media by using the weak Galerkin finite element method. Two types of boundary value problems are considered in the algorithm development: (1) normal boundary condition, and (2) tangential boundary condition. A new variational formulation is developed for the normal boundary value problem by using the Helmholtz decomposition which avoids the computation of functions in the harmonic fields. Both boundary value problems are reduced to a general saddle-point problem involving the curl and divergence operators, for which the weak Galerkin finite element method is devised and analyzed. The novelty of the technique lies in the discretization of the divergence operator applied to vector fields with heterogeneous media. Error estimates of optimal order are established for the corresponding finite element approximations in various discrete Sobolev norms.


title: Microscopic Model for the Scintillation-Light Generation and Light-Quenching in CaWO$_4$ single crystals

abstract: Scintillators are employed for particle detection and identification using light-pulse shapes and light quenching factors. We developed a comprehensive model describing the light generation and quenching in CaWO$_4$ single crystals used for direct dark matter search. All observed particle-dependent light-emission characteristics can be explained quantitatively, light-quenching factors and light-pulse shapes are calculated on a microscopic basis. This model can be extended to other scintillators such as inorganic crystal scintillators, liquid noble gases or organic liquid scintillators.


title: Quantum non-locality - It ain't necessarily so...

abstract: Bell's theorem is 50 years old. Still there is a controversy about its implications. Much of it has its roots in confusion regarding the premises from which the theorem can be derived. Some claim that a derivation of Bell's inequalities requires just locality assumption, and nothing more. Violations of the inequalities are then interpreted as ``nonlocality'' or ``quantum nonlocality''. We show that such claims are unfounded and that every derivation of Bell's inequalities requires a premise -- in addition to locality and freedom of choice -- which is either assumed tacitly, or unconsciously, or is embedded in a single compound condition (like Bell's ``local causality''). The premise is equivalent to the assumption of existence of additional variables which do not appear in the quantum formalism (in form of determinism, or joint probability for outcomes of all conceivable measurements, or ``additional causes`'', or ``hidden variables'', ``complete description of the state'' or counterfactual definiteness, etc.). A certain irony is that perhaps the main message of violation of Bell's inequalities is that our notion of locality should be based on an operationally well-defined no-signalling condition, rather than on local causality.


title: Maximum mass of magnetic white dwarfs

abstract: We revisit in this work the problem of the maximum masses of magnetized White Dwarfs (WD). The impact of a strong magnetic field onto the structure equations is addressed. The pressures become anisotropic due to the presence of the magnetic field and split into a parallel and perpendicular components. We first construct stable solutions of TOV equations for the parallel pressures, and found that physical solutions vanish for the perpendicular pressure when $B \gtrsim 10^{13}$ G. This fact establishes an upper bound for a magnetic field and the stability of the configurations in the (quasi) spherical approximation. Our findings also indicate that it is not possible to obtain stable magnetized WD with super Chandrasekhar masses because the values of the magnetic field needed for them are higher than this bound. To proceed into the anisotropic regime, we derived structure equations appropriated for a cylindrical metric with anisotropic pressures. From the solutions of the structure equations in cylindrical symmetry we have confirmed the same bound for $B \sim 10^{13} $ G, since beyond this value no physical solutions are possible. Our tentative conclusion is that massive WD, with masses well beyond the Chandrasekhar limit do not constitute stable solutions and should not exist.


title: Holomorph of generalized Bol loops

abstract: The notion of the holomorph of a generalized Bol loop and generalized flexible-Bol loop are characterized. With the aid of two self-mappings on the holomorph of a loop, it is shown that: the loop is a generalized Bol loop if and only if its holomorph is a generalized Bol loop; the loop is a generalized flexible-Bol loop if and only if its holomorph is a generalized flexible-Bol loop. Furthermore, elements of the Bryant Schneider group of a generalized Bol loop are characterized in terms of pseudo-automorphism, and the automorphisms gotten are used to build the holomorph of the generalized Bol loop.


title: Sparse Bayesian Learning for EEG Source Localization

abstract: Purpose: Localizing the sources of electrical activity from electroencephalographic (EEG) data has gained considerable attention over the last few years. In this paper, we propose an innovative source localization method for EEG, based on Sparse Bayesian Learning (SBL). Methods: To better specify the sparsity profile and to ensure efficient source localization, the proposed approach considers grouping of the electrical current dipoles inside human brain. SBL is used to solve the localization problem in addition with imposed constraint that the electric current dipoles associated with the brain activity are isotropic. Results: Numerical experiments are conducted on a realistic head model that is obtained by segmentation of MRI images of the head and includes four major components, namely the scalp, the skull, the cerebrospinal fluid (CSF) and the brain, with appropriate relative conductivity values. The results demonstrate that the isotropy constraint significantly improves the performance of SBL. In a noiseless environment, the proposed method was 1 found to accurately (with accuracy of >75%) locate up to 6 simultaneously active sources, whereas for SBL without the isotropy constraint, the accuracy of finding just 3 simultaneously active sources was <75%. Conclusions: Compared to the state-of-the-art algorithms, the proposed method is potentially more consistent in specifying the sparsity profile of human brain activity and is able to produce better source localization for EEG.


title: Color confinement and screening in the $\mathbf{\theta}$ -vacuum

abstract: QCD perturbation theory ignores the compact nature of $SU(3)$ gauge group that gives rise to the periodic $\theta$-vacuum of the theory. We propose to modify the gluon propagator to reconcile perturbation theory with the anomalous Ward identities for the topological current in the $\theta$-vacuum. As a result, the gluon couples to the Veneziano ghost describing the tunneling transitions between different Chern-Simons sectors of the vacuum; we call the emerging gluon dressed by ghost loops a "glost". We evaluate the glost propagator and find that it has the form $G(p) = (p^2 + \chi_{top}/p^2)^{-1}$ where $\chi_{top}$ is the Yang-Mills topological susceptibility related to the $\eta'$ mass by Witten-Veneziano relation; this propagator describes confinement of gluons at distances $\sim \chi_{top}^{-1/4} \simeq 1$ fm. The same functional form of the propagator was originally proposed by Gribov as a solution to the gauge copies problem that plagues perturbation theory. The resulting running coupling coincides with the perturbative one at $p^2 \gg \sqrt{\chi_{top}}$, but in the infrared region either freezes (in pure Yang-Mills theory) or vanishes (in full QCD with light quarks), in accord with experimental evidence. Our scenario makes explicit the connection between confinement and topology of the QCD vacuum; we discuss the implications for spin physics, high energy scattering, and the physics of quark-gluon plasma.


title: Misaligned gas discs around eccentric black-hole binaries and implications for the final-parsec problem

abstract: We investigate the evolution of low mass (Md /Mb = 0.005) misaligned gaseous discs around eccentric supermassive black hole (SMBH) binaries. These are expected to form from randomly oriented accretion events onto a SMBH binary formed in a galaxy merger. When expanding the interaction terms between the binary and a circular ring to quadrupole order and averaging over the binary orbit, we expect four non-precessing disc orientations: aligned or counter-aligned with the binary, or polar orbits around the binary eccentricity vector with either sense of rotation. All other orientations precess around either of these, with the polar precession dominating for high eccentricity. These expectations are borne out by smoothed particle hydrodynamics simulations of initially misaligned viscous circumbinary discs, resulting in the formation of polar rings around highly eccentric binaries in contrast to the co-planar discs around circular binaries. Moreover, we observe disc tearing and violent interactions between differentially precessing rings in the disc significantly disrupting the disc structure and causing gas to fall onto the binary with little angular momentum. While accretion from a polar disc may not promote SMBH binary coalescence (solving the `final-parsec problem'), ejection of this infalling low-angular momentum material via gravitational slingshot is a possible mechanism to reduce the binary separation. Moreover, this process acts on dynamical rather than viscous time scales, and so is much faster.


title: Tracing cool molecular gas and star formation on $\sim 100$pc scales within a $z=2.3$ galaxy

abstract: We present new, high-angular resolution interferometric observations with the Karl G. Jansky Very Large Array of $^{12}$CO $J=1-0$ line emission and 4-8 GHz continuum emission in the strongly lensed, $z=2.3$ submillimetre galaxy, SMM J21352-0102. Using these data, we identify and probe the conditions in $\sim 100$pc clumps within this galaxy, which we consider to be potential giant molecular cloud complexes, containing up to half of the total molecular gas in this system. In combination with far-infrared and submillimetre data, we investigate the far-infrared/radio correlation, measuring $q_{IR} = 2.39 \pm 0.17$ across SMM J21352. We search for variations in the properties of the interstellar medium throughout the galaxy by measuring the spatially-resolved $q_{IR}$ and radio spectral index, ${\alpha}_{\rm radio}$, finding ranges $q_{IR} = [2.1, 2.6]$ and ${\alpha}_{\rm radio} = [-1.5, -0.7]$. We argue that these ranges in ${\alpha}_{\rm radio}$ and $q_{IR}$ may reflect variations in the age of the ISM material. Using multi-$J$ $^{12}$CO data, we quantitatively test a recent theoretical model relating the star-formation rate surface density to the excitation of $^{12}$CO, finding good agreement between the model and the data. Lastly, we study the Schmidt-Kennicutt relation, both integrated across the system and within the individual clumps. We find small offsets between SMM J21352 and its clumps relative to other star-forming galaxy populations on the Schmidt-Kennicutt plot - such offsets have previously been interpreted as evidence for a bi-modal star-formation law, but we argue that they can be equally-well explained as arising due to a combination of observational uncertainties and systematic biases in the choice of model used to interpret the data.


title: Linear mode stability of Kerr-Newman and its quasinormal modes

abstract: We provide strong evidence that, up to $99.999\%$ of extremality, Kerr-Newman black holes (KN BHs) are linear mode stable within Einstein-Maxwell theory. We derive and solve, numerically, a coupled system of two PDEs for two gauge invariant fields that describe the most general linear perturbations of a KN BH (except for trivial modes that shift the parameters of the solution). We determine the quasinormal mode (QNM) spectrum of the KN BH as a function of its three parameters and find no unstable modes. In addition, we find that the QNMs that are connected continuously to the gravitational $\ell=m=2$ Schwarzschild QNM dominate the spectrum for all values of the parameter space ($m$ is the azimuthal number of the wave function and $\ell$ measures the number of nodes along the polar direction). Furthermore, all QNMs with $\ell=m$ approach Re$\,\omega = m \Omega_H^{ext}$ and Im$\,\omega=0$ at extremality; this is a universal property for any field of arbitrary spin $|s|\leq 2$ propagating on a KN BH background ($\omega$ is the wave frequency and $\Omega_H^{ext}$ the BH angular velocity at extremality). We compare our results with available perturbative results in the small charge or small rotation regimes and find good agreement. We also present a simple proof that the Regge-Wheeler (odd) and Zerilli (even) sectors of Schwarzschild perturbations must be isospectral.


title: Integrated Light Chemical Tagging Analyses of Seven M31 Outer Halo Globular Clusters from the Pan-Andromeda Archaeological Survey

abstract: Detailed chemical abundances are presented for seven M31 outer halo globular clusters (with projected distances from M31 greater than 30 kpc), as derived from high resolution integrated light spectra taken with the Hobby Eberly Telescope. Five of these clusters were recently discovered in the Pan-Andromeda Archaeological Survey (PAndAS)---this paper presents the first determinations of integrated Fe, Na, Mg, Ca, Ti, Ni, Ba, and Eu abundances for these clusters. Four of the target clusters (PA06, PA53, PA54, and PA56) are metal-poor ([Fe/H] < -1.5), alpha-enhanced (though they are possibly less alpha-enhanced than Milky Way stars at the 1 sigma level), and show signs of star-to-star Na and Mg variations. The other three globular clusters (H10, H23, and PA17) are more metal rich, with metallicities ranging from [Fe/H] = -1.4 to -0.9. While H23 is chemically similar to Milky Way field stars, Milky Way globular clusters, and other M31 clusters, H10 and PA17 have moderately low [Ca/Fe], compared to Milky Way field stars and clusters. Additionally, PA17's high [Mg/Ca] and [Ba/Eu] ratios are distinct from Milky Way stars, and are in better agreement with the stars and clusters in the Large Magellanic Cloud (LMC). None of the clusters studied here can be conclusively linked to any of the identified streams from PAndAS; however, based on their locations, kinematics, metallicities, and detailed abundances, the most metal-rich PAndAS clusters H23 and PA17 may be associated with the progenitor of the Giant Stellar Stream, H10 may be associated with the SW Cloud, and PA53 and PA56 may be associated with the Eastern Cloud.


title: Time domain studies of Active Galactic Nuclei with the Square Kilometre Array

abstract: Variability of radio-emitting active galactic nuclei can be used to probe both intrinsic variations arising from shocks, flares, and other changes in emission from regions surrounding the central supermassive black hole, as well as extrinsic variations due to scattering by structures in our own Galaxy. Such interstellar scattering also probes the structure of the emitting regions, with microarcsecond resolution. Current studies have necessarily been limited to either small numbers of objects monitored over long periods of time, or large numbers of objects but with poor time sampling. The dramatic increase in survey speed engendered by the Square Kilometre Array will enable precision synoptic monitoring studies of hundreds of thousands of sources with a cadence of days or less. Statistics of variability, in particular concurrent observations at multiple radio frequencies and in other bands of the electromagnetic spectrum, will probe accretion physics over a wide range of AGN classes, luminosities, and orientations, as well as enabling a detailed understanding of the structures responsible for radio wave scattering in the Galactic interstellar medium.


title: HST+COS spectra of the double white dwarf CSS 41177 place the secondary inside the pulsational instability strip

abstract: We present Hubble Space Telescope + Cosmic Origins Spectrograph (HST+COS) data of the eclipsing double white dwarf binary CSS 41177. Due to the temperature difference between the two white dwarfs, the HST+COS far-ultraviolet data are dominated by the hot, primary white dwarf and allow us to precisely measure its temperature (T1). Using eclipse observations, we also tightly constrain the temperature of the cooler secondary white dwarf (T2). Our results, where T1 = 22439 +/- 59 K and T2 = 10876 +/- 32 K, with the uncertainties being purely statistical, place the secondary inside and close to the blue edge of the empirical instability strip for low temperature hydrogen-atmosphere white dwarfs. Dedicated high-speed photometry is encouraged to probe for the presence of pulsations, which will constrain the border of the instability strip as well as probe a new region of low gravity within the strip.


title: The SKA view of Gamma-ray Bursts

abstract: Gamma-ray bursts (GRBs) are some of the most extreme events in the Universe. As well as providing a natural laboratory for investigating fundamental physical processes, they might trace the cosmic star formation rate up to extreme redshifts and probe the composition of the intergalactic medium over most of the Universe's history. Radio observations of GRBs play a key part in determining their physical properties, but currently they are largely limited to follow-up observations of $\gamma$-ray-detected objects. The SKA will significantly increase our ability to study GRB afterglows, following up several hundred objects in the high frequency bands already in the "early science" implementation of the telescope. SKA1-MID Bands 4 (4 GHz) and 5 (9.2 GHz) will be particularly suited to the detection of these transient phenomena. The SKA will trace the peak of the emission, sampling the thick-to-thin transition of the evolving spectrum, and follow-up the afterglow down to the time the ejecta slow down to non-relativistic speeds. The full SKA will be able to observe the afterglows across the non-relativistic transition, for ~25% of the whole GRB population. This will allow us to get a significant insight into the true energy budget of GRBs, probe their surrounding density profile, and the shock microphysics. The SKA will also be able to routinely detect the elusive "orphan afterglow" emission, from the population of GRBs whose jets are not pointed towards the Earth. We expect that a deep all-sky survey such as SKA1-SUR will see around 300 orphan afterglows every week. We predict these detection to be >1000 when the full SKA telescope will be operational.


title: Formation of In Situ Stellar Haloes in Milky Way-Mass Galaxies

abstract: We study the formation of stellar haloes in three Milky Way-mass galaxies using cosmological smoothed particle hydrodynamics simulations, focusing on the subset of halo stars that form in situ, as opposed to those accreted from satellites. In situ stars in our simulations dominate the stellar halo out to 20 kpc and account for 30-40 per cent of its total mass. We separate in situ halo stars into three straightforward, physically distinct categories according to their origin: stars scattered from the disc of the main galaxy ('heated disc'), stars formed from gas smoothly accreted on to the halo ('smooth' gas) and stars formed in streams of gas stripped from infalling satellites ('stripped' gas). We find that most belong to the stripped gas category. Those originating in smooth gas outside the disc tend to form at the same time and place as the stripped-gas population, suggesting that their formation is associated with the same gas-rich accretion events. The scattered disc star contribution is negligible overall but significant in the Solar neighbourhood, where ~90 per cent of stars on eccentric orbits once belonged to the disc. However, the distinction between halo and thick disc in this region is highly ambiguous. The chemical and kinematic properties of the different components are very similar at the present day, but the global properties of the in situ halo differ substantially between the three galaxies in our study. In our simulations, the hierarchical buildup of structure is the driving force behind not only the accreted stellar halo, but also those halo stars formed in situ.


title: Using machine learning to classify the diffuse interstellar bands

abstract: Using over a million and a half extragalactic spectra we study the correlations of the Diffuse Interstellar Bands (DIBs) in the Milky Way. We measure the correlation between DIB strength and dust extinction for 142 DIBs using 24 stacked spectra in the reddening range E(B-V) < 0.2, many more lines than ever studied before. Most of the DIBs do not correlate with dust extinction. However, we find 10 weak and barely studied DIBs with correlations that are higher than 0.7 with dust extinction and confirm the high correlation of additional 5 strong DIBs. Furthermore, we find a pair of DIBs, 5925.9A and 5927.5A which exhibits significant negative correlation with dust extinction, indicating that their carrier may be depleted on dust. We use Machine Learning algorithms to divide the DIBs to spectroscopic families based on 250 stacked spectra. By removing the dust dependency we study how DIBs follow their local environment. We thus obtain 6 groups of weak DIBs, 4 of which are tightly associated with C2 or CN absorption lines.


title: Galaxies that Shine: radiation-hydrodynamical simulations of disk galaxies

abstract: Radiation feedback is typically implemented using subgrid recipes in hydrodynamical simulations of galaxies. Very little work has so far been performed using radiation-hydrodynamics (RHD), and there is no consensus on the importance of radiation feedback in galaxy evolution. We present RHD simulations of isolated galaxy disks of different masses with a resolution of 18 pc. Besides accounting for supernova feedback, our simulations are the first galaxy-scale simulations to include RHD treatments of photo-ionisation heating and radiation pressure, from both direct optical/UV radiation and multi-scattered, re-processed infrared (IR) radiation. Photo-heating smooths and thickens the disks and suppresses star formation about as much as the inclusion of ("thermal dump") supernova feedback does. These effects decrease with galaxy mass and are mainly due to the prevention of the formation of dense clouds, as opposed to their destruction. Radiation pressure, whether from direct or IR radiation, has little effect, but for the IR radiation we show that its impact is limited by our inability to resolve the high optical depths for which multi-scattering becomes important. While artificially boosting the IR optical depths does reduce the star formation, it does so by smoothing the gas rather than by generating stronger outflows. We conclude that although higher-resolution simulations, and potentially also different supernova implementations, are needed for confirmation, our findings suggest that radiation feedback is more gentle and less effective than is often assumed in subgrid prescriptions.