title: Positive and negative electrocaloric effect in BaTiO$_3$ in the presence of defect dipoles

abstract: The influence of defect dipoles on the electrocaloric effect (ECE) in acceptor doped BaTiO$_3$ is studied by means of lattice-based Monte-Carlo simulations. A Ginzburg-Landau type effective Hamiltonian is used. Oxygen vacancy-acceptor associates are described by fixed defect dipoles with orientation parallel or anti-parallel to the external field. By a combination of canonical and microcanoncial simulations the ECE is directly evaluated. Our results show that in the case of anti-parallel defect dipoles the ECE can be positive or negative depending on the density of defect dipoles. Moreover, a transition from a negative to positive ECE can be observed from a certain density of anti-parallel dipoles on when the external field increases. These transitions are due to the delicate interplay of internal and external fields, and are explained by the domain structure evolution and related field-induced entropy changes. The results are compared to those obtained by MD simulations employing an {\it{ab initio}} based effective Hamiltonian, and a good qualitative agreement is found. In addition, a novel electrocaloric cycle, which makes use of the negative ECE and defect dipoles, is proposed to enhance the cooling effect.


title: Time evolution of the anisotropies of the hydrodynamically expanding sQGP

abstract: In high energy heavy ion collisions of RHIC and LHC, a strongly interacting quark gluon plasma (sQGP) is created. This medium undergoes a hydrodynamic evolution, before it freezes out to form a hadronic matter. The initial state of the sQGP is determined by the initial distribution of the participating nucleons and their interactions. Due to the finite number of nucleons, the initial distribution fluctuates on an event-by-event basis. The transverse plane anisotropy of the initial state can be translated into a series of anisotropy coefficients or eccentricities: second, third, fourth-order anisotropy etc. These anisotropies then evolve in time, and result in measurable momentum-space anisotropies, to be measured with respect to their respective symmetry planes. In this paper we investigate the time evolution of the anisotropies. With a numerical hydrodynamic code, we analyze how the speed of sound and viscosity influence this evolution.


title: Local interstellar cosmic-ray spectra derived from gamma-ray emissivities

abstract: Precise gamma-ray emissivities from cosmic-ray interactions with interstellar gas have been recently derived using Fermi-LAT data, and used to constrain the local interstellar spectra of protons and leptons. We report on a continuing effort to exploit these emissivities combined with the latest hadronic gamma-ray production cross-sections and other constraints such as synchrotron emission for the leptonic component. The interstellar spectra provide important information for heliospheric modulation, and cosmic-ray origin and propagation.


title: Non-equilibrium steady-states in a driven-dissipative superfluid

abstract: We experimentally study a driven-dissipative Josephson junction array, realized with a weakly interacting Bose Einstein condensate residing in a one-dimensional optical lattice. Engineered losses on one site act as a local dissipative process, while tunneling from the neighboring sites constitutes the driving force. We characterize the emerging steady-states of this atomtronic device. With increasing dissipation strength $\gamma$ the system crosses from a superfluid state, characterized by a coherent Josephson current into the lossy site to a resistive state, characterized by an incoherent hopping transport. For intermediate values of $\gamma$, the system exhibits bistability, where a superfluid and a resistive branch coexist. We also study the relaxation dynamics towards the steady-state, where we find a critical slowing down, indicating the presence of a non-equilibrium phase transition.


title: Upper bounds for number of removed edges in the Erased Configuration Model

abstract: Models for generating simple graphs are important in the study of real-world complex networks. A well established example of such a model is the erased configuration model, where each node receives a number of half-edges that are connected to half-edges of other nodes at random, and then self-loops are removed and multiple edges are concatenated to make the graph simple. Although asymptotic results for many properties of this model, such as the limiting degree distribution, are known, the exact speed of convergence in terms of the graph sizes remains an open question. We provide a first answer by analyzing the size dependence of the average number of removed edges in the erased configuration model. By combining known upper bounds with a Tauberian Theorem we obtain upper bounds for the number of removed edges, in terms of the size of the graph. Remarkably, when the degree distribution follows a power-law, we observe three scaling regimes, depending on the power law exponent. Our results provide a strong theoretical basis for evaluating finite-size effects in networks.


title: Invariance of the Noether charge

abstract: Surprisingly, an interesting property of the Noether charge that it is by itself invariant under the corresponding symmetry transformation is never discussed in quantum field theory or classical mechanics textbooks we have checked. This property is also almost never mentioned in articles devoted to Noether's theorem. Nevertheless, to prove this property in the context of Lagrangian formalism is not quite trivial and the proof, outlined in this article, can constitute an useful and interesting exercise for students.


title: Precision Estimation of Source Dimensions from Higher-Order Intensity Correlations

abstract: An important topic of interest in imaging is the construction of protocols that are not diffraction limited. This can be achieved in a variety of ways, including classical superresolution techniques or quantum entanglement-based protocols. Here, we consider superresolving imaging in the far field using higher-order intensity correlations. We show that third and fourth order correlations can improve upon the first and second order correlations that are traditionally used in classical optics and Hanbury Brown and Twiss type experiments. The improvement is achieved entirely by post-processing of the data. As a demonstrator, we simulate the far field intensity distribution of a circular aperture that emits thermal light and use maximum likelihood estimation to determine the radius of the aperture. We compare the achieved precision to the Cram\'er-Rao lower bound and find that the variance of measurements for the third and fourth order correlation functions are indeed closer to the Cram\'er-Rao bound than that of the second order correlation function. The method presented here is general, and can be used for all kinds of incoherent emitters, geometries, and types of noise.


title: Multi-Stage Pricing for Coordination of Thermostatically Controlled Loads: A Dynamic Stackelberg Game Approach

abstract: This paper focuses on multi-stage coordination for a population of thermostatically controlled loads (TCL). Each load maximizes the individual utility in response to an energy price, while the coordinator determines the price to maximize the social welfare subject to a peak energy constraint. The coordination problem is formulated as a dynamic Stackelberg game. The Stackelberg solution is derived using an indirect approach: we first obtain an upper bound of the Stackelberg game, then we prove that this upper bound is attainable. Sufficient conditions for the optimal solution are characterized. Based on these conditions, a pricing algorithm with guaranteed convergence is developed to compute the Stackelberg solution employing the monotonicity of the user responses. Our control strategy is validated using real market data and weather information, and realistic simulation results are presented to demonstrate the effectiveness of the proposed control framework.


title: Crystal excitations features in the photon emission spectrum of the quantum channeled particle

abstract: The processes of generation of phonons and electron excitations with quantum oriented particle in crystal are considered. Expressions are obtained for the probability of emission of phonons and plasmons with channeled particle. The theory of the processes of the crystal excitation with the channeled particle, accompanied by the simultaneous emission of a photon is developed. It is proved that all the specific features of the electron and phonon crystal excitation structure appear as components of radiation of the oriented fast charged particle. The photon emission of a quantum channeled particle, accompanied by a plasmon excitation in a crystal target is considered. It is taken into account that the plasmon energy in the crystal is of the same order with the depth of the potential well in which the channeled particle moves. A weak dispersion of the plasmons in the crystal is taken into account, which leads to the selection of the fixed plasmon energy. It is shown that fast charged particle in the rest frame emits photons with an energy equal to the energy difference between two quantized levels of transverse motion with the deduction of the net plasmon energy. The characteristics of photon-plasmon radiation are investigated.


title: Viscosity solutions of systems of variational inequalities with interconnected bilateral obstacles of non-local type

abstract: In this paper, we study systems of nonlinear second-order variational inequalities with interconnected bilateral obstacles with non-local terms. They are of min-max and max-min types and related to a multiple modes zero-sum switching game in the jump-diffusion model. Using systems of penalized reflected backward SDEs with jumps and unilateral interconnected obstacles, and their associated deterministic functions, we construct for each system a continuous viscosity solution which is unique in the class of functions with polynomial growth.


title: Compositional Construction of Approximate Abstractions of Interconnected Control Systems

abstract: We consider a compositional construction of approximate abstractions of interconnected control systems. In our framework, an abstraction acts as a substitute in the controller design process and is itself a continuous control system. The abstraction is related to the concrete control system via a so-called simulation function: a Lyapunov-like function, which is used to establish a quantitative bound between the behavior of the approximate abstraction and the concrete system. In the first part of the paper, we provide a small gain type condition that facilitates the compositional construction of an abstraction of an interconnected control system together with a simulation function from the abstractions and simulation functions of the individual subsystems. In the second part of the paper, we restrict our attention to linear control system and characterize simulation functions in terms of controlled invariant, externally stabilizable subspaces. Based on those characterizations, we propose a particular scheme to construct abstractions for linear control systems. We illustrate the compositional construction of an abstraction on an interconnected system consisting of four linear subsystems. We use the abstraction as a substitute to synthesize a controller to enforce a certain linear temporal logic specification.


title: Injective tests of low complexity in the plane

abstract: We study injective versions of the characterization of sets potentially in a Wadge class of Borel sets, for the first Borel and Lavrentieff classes. We also study the case of oriented graphs in terms of continuous homomorphisms, injective or not.


title: Incremental Variational Inference for Latent Dirichlet Allocation

abstract: We introduce incremental variational inference and apply it to latent Dirichlet allocation (LDA). Incremental variational inference is inspired by incremental EM and provides an alternative to stochastic variational inference. Incremental LDA can process massive document collections, does not require to set a learning rate, converges faster to a local optimum of the variational bound and enjoys the attractive property of monotonically increasing it. We study the performance of incremental LDA on large benchmark data sets. We further introduce a stochastic approximation of incremental variational inference which extends to the asynchronous distributed setting. The resulting distributed algorithm achieves comparable performance as single host incremental variational inference, but with a significant speed-up.


title: Signatures of the Current Blockade Instability in Suspended Carbon Nanotubes

abstract: Transport measurements allow sensitive detection of nanomechanical motion of suspended carbon nanotubes. It has been predicted that when the electro-mechanical coupling is sufficiently large a bistability with a current blockade appears. Unambiguous observation of this transition by current measurements may be difficult. Instead, we investigate the mechanical response of the system, namely the displacement spectral function; the linear response to a driving; and the ring-down behavior. We find that by increasing the electro-mechanical coupling the peak in the spectral function broadens and shifts at low frequencies while the oscillator dephasing time shortens. These effects are maximum at the transition where non-linearities dominate the dynamics. These strong signatures open the way to detect the blockade transition in devices currently studied by several groups.


title: Spectral Synchronicity in Brain Signals

abstract: Brain activity following stimulus presentation and during resting state are often the result of highly coordinated responses of large numbers of neurons both locally and globally. Coordinated activity of neurons can give rise to oscillations which are captured by electroencephalograms (EEG). In this paper, we examine EEGs as this is the primary data being used by our collaborators who are studying coordination of neuronal response during the execution of tasks such as learning, and memory formation, retention and retrieval. In this paper, we develop the spectral merger clustering (SMC) method that identifies synchronized brain regions during resting state in a sense that these regions share similar oscillations or waveforms. The SMC method, produces clusters of EEGs which serve as a proxy for segmenting the brain cortical surface since the EEGs capture neuronal activity over a locally distributed region on the cortical surface. The extent of desynchronicity between a pair of EEGs is measured using the total variation distance (TVD) which gives the largest possible difference between the spectral densities of the pair of EEGs. We considered the spectral merger algorithm for clustering EEGs, which updates the spectral estimate of the cluster from a weighted average of the spectral estimate obtained from each EEG in the cluster. Numerical experiments suggest that the SMC method performs very well in producing the correct clusters. When applied to resting state EEG data, the method showed how some regions, though not contiguously connected on the cortical surface, are spectrally synchronized during resting state. Moreover, the method demonstrates that brain organization, as expressed in cluster formation, evolves over resting state.


title: Support Vector Machine in Prediction of Building Energy Demand Using Pseudo Dynamic Approach

abstract: Building's energy consumption prediction is a major concern in the recent years and many efforts have been achieved in order to improve the energy management of buildings. In particular, the prediction of energy consumption in building is essential for the energy operator to build an optimal operating strategy, which could be integrated to building's energy management system (BEMS). This paper proposes a prediction model for building energy consumption using support vector machine (SVM). Data-driven model, for instance, SVM is very sensitive to the selection of training data. Thus the relevant days data selection method based on Dynamic Time Warping is used to train SVM model. In addition, to encompass thermal inertia of building, pseudo dynamic model is applied since it takes into account information of transition of energy consumption effects and occupancy profile. Relevant days data selection and whole training data model is applied to the case studies of Ecole des Mines de Nantes, France Office building. The results showed that support vector machine based on relevant data selection method is able to predict the energy consumption of building with a high accuracy in compare to whole data training. In addition, relevant data selection method is computationally cheaper (around 8 minute training time) in contrast to whole data training (around 31 hour for weekend and 116 hour for working days) and reveals realistic control implementation for online system as well.


title: Recent extensions to GALPROP

abstract: Some recent extensions to the GALPROP cosmic-ray propagation package are described. The enhancements include: an accurate solution option, improved convection formulation, alternative spatial boundary conditions, polarized synchrotron emission, new magnetic field models, updated gamma-ray production cross-sections, free-free radio emission and absorption, primary positrons, additional injection spectral breaks, deuterium production by pp fusion, hadronic energy losses, improved HEALPix skymap format, compatibility with latest HEALPix release, and various bug fixes. The Explanatory Supplement has been extensively updated, including details of these enhancements. A compatible plot package GALPLOT for GALPROP output is also provided, as well as other related software.


title: Non-asymptotic convergence analysis for the Unadjusted Langevin Algorithm

abstract: In this paper, we study a method to sample from a target distribution $\pi$ over $\mathbb{R}^d$ having a positive density with respect to the Lebesgue measure, known up to a normalisation factor. This method is based on the Euler discretization of the overdamped Langevin stochastic differential equation associated with $\pi$. For both constant and decreasing step sizes in the Euler discretization, we obtain non-asymptotic bounds for the convergence to the target distribution $\pi$ in total variation distance. A particular attention is paid to the dependency on the dimension $d$, to demonstrate the applicability of this method in the high dimensional setting. These bounds improve and extend the results of (Dalalyan 2014).


title: Tetraquarks from the Bethe-Salpeter equation

abstract: We present a numerical solution of the four-quark Bethe-Salpeter equation for a scalar tetraquark. We find that the four-body equation dynamically generates pseudoscalar poles in the Bethe-Salpeter amplitude. The sensitivity to the pion poles leads to a light isoscalar tetraquark mass M ~ 400 MeV, which is comparable to that of the sigma/f0(500). The masses of its multiplet partners kappa and a0/f0 follow a similar pattern, thereby providing support for the tetraquark interpretation of the light scalar nonet.


title: A study of the magnetotransport properties of the graphene (II. Fractional Quantum Hall Effect)

abstract: We present an approach to the fractional quantum Hall effect observed in grapheme (GFQHE), basing us on the model developed previously for the fractional quantum Hall effect in a two-dimensional electron system embedded in a quantum well (FQHE) [Hidalgo, 2013(*)]. The main idea in the view proposed f0r the FQHE is the breaking of the symmetry of the Hamiltonian of every electron in a two-dimensional electron gas (2DEG) under the application of a magnetic field and in the presence of an electrostatic potential due to the ionized impurities. As the magnetic field increases the effect of that electrostatic potential evolves; changing in turn the spatial symmetry of the Hamiltonian: from continuous to discrete one. The model provides the odd fractional states, and corresponding plateaux, p/3, p being any integer, observed in graphene.


title: Long term monitoring of the optical background in the Capo Passero deep-sea site with the NEMO tower prototype

abstract: The NEMO Phase-2 tower is the first detector which was operated underwater for more than one year at the "record" depth of 3500 m. It was designed and built within the framework of the NEMO (NEutrino Mediterranean Observatory) project. The 380 m high tower was successfully installed in March 2013 80 km offshore Capo Passero (Italy). This is the first prototype operated on the site where the italian node of the KM3NeT neutrino telescope will be built. The installation and operation of the NEMO Phase-2 tower has proven the functionality of the infrastructure and the operability at 3500 m depth. A more than one year long monitoring of the deep water characteristics of the site has been also provided. In this paper the infrastructure and the tower structure and instrumentation are described. The results of long term optical background measurements are presented. The rates show stable and low baseline values, compatible with the contribution of 40K light emission, with a small percentage of light bursts due to bioluminescence. All these features confirm the stability and good optical properties of the site.


title: High-pressure, temperature elasticity of Fe- and Al-bearing MgSiO3: implications for the Earth's lower mantle

abstract: Fe and Al are two of the most important rock-forming elements other than Mg, Si, and O. Their presence in the lower mantle's most abundant minerals, MgSiO_3 bridgmanite, MgSiO_3 post-perovskite and MgO periclase, alters their elastic properties. However, knowledge on the thermoelasticity of Fe- and Al-bearing MgSiO_3 bridgmanite, and post-perovskite is scarce. In this study, we perform ab initio molecular dynamics to calculate the elastic and seismic properties of pure, Fe^{3+}- and Fe^{2+}-, and Al^{3+}-bearing MgSiO_3 perovskite and post-perovskite, over a wide range of pressures, temperatures, and Fe/Al compositions. Our results show that a mineral assemblage resembling pyrolite fits a 1D seismological model well, down to, at least, a few hundred kilometers above the core-mantle boundary, i.e. the top of the D'' region. In D'', a similar composition is still an excellent fit to the average velocities and fairly approximate to the density. We also implement polycrystal plasticity with a geodynamic model to predict resulting seismic anisotropy, and find post-perovskite with predominant (001) slip across all compositions agrees best with seismic observations in the D''.


title: Vacuum polarization on the brane

abstract: We compute the renormalized expectation value of the square of a massless, conformally coupled, quantum scalar field on the brane of a higher-dimensional black hole. Working in the AADD brane-world scenario, the extra dimensions are flat and we assume that the compactification radius is large compared with the size of the black hole. The four-dimensional on-brane metric corresponds to a slice through a higher-dimensional Schwarzschild-Tangherlini black hole geometry and depends on the number of bulk space-time dimensions. The quantum scalar field is in a thermal state at the Hawking temperature. An exact, closed-form expression is derived for the renormalized expectation value of the square of the quantum scalar field on the event horizon of the black hole. Outside the event horizon, this renormalized expectation value is computed numerically. The answer depends on the number of bulk space-time dimensions, with a magnitude which increases rapidly as the number of bulk space-time dimensions increases.


title: Blocks for general linear supergroup $GL(m|n)$

abstract: We prove the linkage principle and describe blocks of the general linear supergroups $GL(m|n)$ over the ground field $K$ of characteristic $p\neq 2$.


title: Survey of Higgs interpretations of the diboson excesses

abstract: We investigate diboson signals in the Standard Model (SM) with an extended Higgs sector, motivated by the excesses in the diboson channels at the LHC. We begin with the unitarity sum-rules of the weak gauge boson scattering assuming the Higgs sector is extended. According to the sum-rules, we discuss the Higgs interpretations of the diboson signals and the consistency with the ATLAS diboson anomaly and other experimental constraints. As a concrete example, we propose a two-Higgs-doublet model where the Yukawa coupling of an extra CP-even scalar with up-type quark is relatively large. The diboson (WW and ZZ) signals can be explained by 2 TeV CP-even Higgs boson, while the partners, the CP-odd and the charged Higgs bosons in the extra doublet, are severely constrained by both the LHC direct search and the indirect search via flavor observables. Especially, in order to avoid the diboson resonance search in the Vh channel, further extensions of the model are required. The diboson excess is correlated with the SM Higgs signals in our framework, so that the precise measurement of the SM Higgs boson is also important to test the Higgs interpretation.


title: Flavour, Electroweak Symmetry Breaking and Dark Matter: state of the art and future prospects

abstract: With the discovery of the Higgs boson the Standard Model has become a complete and comprehensive theory, which has been verified with unparalleled precision and in principle might be valid at all scales. However, several reasons remain why we firmly believe that there should be physics beyond the Standard Model. Experiments such as the LHC, new $B$ factories, and earth- and space-based astro-particle experiments provide us with unique opportunities to discover a coherent framework for many of the long-standing puzzles of our field. Here we explore several significant interconnections between the physics of the Higgs boson, the physics of flavour, and the experimental clues we have about dark matter.


title: Maximum principles for the relativistic heat equation

abstract: The classical heat equation is incompatible with relativity, since the strong maximum principle allows for disturbances to propagate instantaneously. Some authors have proposed limiting the propagation speed by adding a linear hyperbolic correction term, but then even a weak maximum principle fails to hold. We study a more recently introduced relativistic heat equation, which replaces the Laplace operator by a quasilinear elliptic operator, and show that strong and weak maximum principles hold for stationary and time-varying solutions, respectively, as well as for sub- and supersolutions. Moreover, by transforming the equation into an equivalent form, related to the mean curvature operator, we prove even stronger tangency and comparison principles.


title: First- and second-order error estimates in Monte Carlo integration

abstract: In Monte Carlo integration an accurate and reliable determination of the numerical intregration error is essential. We point out the need for an independent estimate of the error on this error, for which we present an unbiased estimator. In contrast to the usual (first-order) error estimator, this second-order estimator can be shown to be not necessarily positive in an actual Monte Carlo computation. We propose an alternative and indicate how this can be computed in linear time without risk of large rounding errors. In addition, we comment on the relatively very slow convergence of the second-order error estimate.


title: Strata Hasse invariants, Hecke algebras and Galois representations

abstract: We construct group-theoretical generalizations of the Hasse invariant on strata closures of the stacks $G$-Zip$^{\mu}$. Restricting to zip data of Hodge type, we obtain a group-theoretical Hasse invariant on every Ekedahl-Oort stratum closure of a general Hodge-type Shimura variety. A key tool is the construction of a stack of zip flags $G$-ZipFlag$^\mu$, fibered in flag varieties over $G$-Zip$^{\mu}$. It provides a simultaneous generalization of the "classical case" homogeneous complex manifolds studied by Griffiths-Schmid and the "flag space" for Siegel varieties studied by Ekedahl-van der Geer. Four applications are obtained: (1) Pseudo-representations are attached to the coherent cohomology of Hodge-type Shimura varieties modulo a prime power. (2) Galois representations are associated to many automorphic representations with non-degenerate limit of discrete series archimedean component. (3) It is shown that all Ekedahl-Oort strata in the minimal compactification of a Hodge-type Shimura variety are affine, thereby proving a conjecture of Oort. (4) Part of Serre's letter to Tate on mod $p$ modular forms is generalized to general Hodge-type Shimura varieties.


title: Classification of Complex Wishart Matrices with a Diffusion-Reaction System guided by Stochastic Distances

abstract: We propose a new method for PolSAR (Polarimetric Synthetic Aperture Radar) imagery classification based on stochastic distances in the space of random matrices obeying complex Wishart distributions. Given a collection of prototypes $\{Z_m\}_{m=1}^M$ and a stochastic distance $d(.,.)$, we classify any random matrix $X$ using two criteria in an iterative setup. Firstly, we associate $X$ to the class which minimizes the weighted stochastic distance $w_md(X,Z_m)$, where the positive weights $w_m$ are computed to maximize the class discrimination power. Secondly, we improve the result by embedding the classification problem into a diffusion-reaction partial differential system where the diffusion term smooths the patches within the image, and the reaction term tends to move the pixel values towards the closest class prototype. In particular, the method inherits the benefits of speckle reduction by diffusion-like methods. Results on synthetic and real PolSAR data show the performance of the method.


title: Bootstrap tuning in ordered model selection

abstract: In the problem of model selection for a given family of linear estimators, ordered by their variance, we offer a new "smallest accepted" approach motivated by Lepski's method and multiple testing theory. The procedure selects the smallest model which satisfies an acceptance rule based on comparison with all larger models. The method is completely data-driven and does not use any prior information about the variance structure of the noise: its parameters are adjusted to the underlying possibly heterogeneous noise by the so-called "propagation condition" using a wild bootstrap method. The validity of the bootstrap calibration is proved for finite samples with an explicit error bound. We provide a comprehensive theoretical study of the method and describe in detail the set of possible values of the selector \( \hat{m} \). We also establish some precise oracle error bounds for the corresponding estimator \( \hat{\theta} = \tilde{\theta}_{\hat{m}} \) which equally applies to estimation of the whole parameter vectors, some subvector or linear mapping, as well as the estimation of a linear functional.


title: Fractional Riesz-Hilbert transforms and fractional monogenic signals

abstract: The fractional Hilbert transforms plays an important role in optics and signal processing. In particular the analytic signal proposed by Gabor has as a key component the Hilbert transform. The higher dimensional Hilbert transform is the Riesz-Hilbert transform which was used by Felsberg and Sommer to construct the monogenic signal. We will construct fractional and quaternionic fractional Riesz-Hilbert transforms based on a eigenvalue decomposition. We will prove properties of these transformations such as shift and scale invariance, orthogonality and the semigroup property. Based on the fractional/quaternionic fractional Riesz-Hilbert transform we construct (quaternionic) fractional monogenic signals. These signals are rotated and modulated monogenic signals.


title: Testing the Kerr Paradigm with X-ray Observations

abstract: Astrophysical black hole candidates are thought to be the Kerr black holes of general relativity, but the actual nature of these objects has still to be confirmed. The continuum-fitting and the iron line methods are currently the only available techniques to probe the spacetime geometry around these bodies and test the Kerr black hole paradigm. The continuum-fitting method is a robust approach, but the shape of the disk's thermal spectrum is in general too simple to measure the spin and to constrain possible deviations from the Kerr solution at the same time. The iron line analysis is potentially a powerful technique, but at the moment we do not have high quality data and a robust astrophysical model.


title: Set Theory for The (Smart) Masses

abstract: Proof Designer is a computer software program designed to help Mathematics students learn to write mathematical proofs. Under the guidance of the user, Proof Designer assists in writing outlines of proofs in elementary set theory. Proof Designer was designed by Daniel Velleman in association with his book "How To Prove It: A Structured Approach" to help students apply the methods discussed in the book, making classes based on the book more interactive. This article is an early report on the progress of our effort to "bring set theory to the masses" by developing Proof Maker, a new Proof Designer-inspired software that ports Proof Designer to hand-held devices such as smart-phones and tablets. Proof Maker, when completed, will allow students to use Proof Designer with the ease of a touch, literally, on their smart devices. Our goal behind developing Proof Maker is to enable any one who is interested enough to develop elementary set theory proofs anywhere he or she might be (think of doing proofs while waiting at a bus stop!) and at any time he or she wishes (think of writing proofs before going to bed, or even in bed!). In this paper we report on the improvements we made to Proof Designer so far, and on the (many) steps remaining for us to have a fully-functioning Proof Maker "in our hands".


title: Exponential Convergence through Linear Finite Element Discretization of Stratified Subdomains

abstract: Motivated by problems where the response is needed at select localized regions in a large computational domain, we devise a novel finite element discretization that results in exponential convergence at pre-selected points. The two key features of the discretization are (a) use of midpoint integration to evaluate the contribution matrices, and (b) an unconventional bending of the mesh into complex space. Named complex-length finite element method (CFEM), the technique is linked to Pade approximants that provide exponential convergence of the Dirichlet-to-Neumann maps and thus the solution at specified points in the domain. Exponential convergence facilitates drastic reduction in the number of elements. This, combined with sparse computation associated with linear finite elements, results in significant reduction in the computational cost. The paper presents the basic ideas of the method as well as illustration of its effectiveness for a variety of problems involving Laplace, Helmholtz and elastodynamic equations.


title: Trial for a proof of the Syracuse conjecture

abstract: The infamous 3x+1 conjecture spread by Lothar Collatz in 1952, despite its elementary formulation, remained unproved for over 60 years. From the heuristical probabilistic approach to the complex mapping of the algorithm, the scientific community has fetched for many methods to try to prove it formally, and thus, mathematicians like Erdos tend to believe that "mathematics are not yet ready for such problems". In this research report, covering domains like algebra and graph theory, it is shown a trial of proof of the conjecture by disproval of its two antitheses: the existence of an ever-growing Syracuse suite and the existence of a cycle different from the cycle 4,2,1.


title: Modeling the self-energy and wavefunction relaxation in the orbitals

abstract: The strong boundary normalized condition of wavefunction for fully occupied semicore 3d orbitals leads the linear response DFT+U on such metal oxide to have an insurmountable obstacle in Hubbard U determination. We treated the orbital self-energy and orbital relaxation as components of eigenvalues with respective orbital occupation number that follows the Fermi-Dirac distribution. By self-consistently solving the second partial deviation of total energy based on the most simple local density formalism with Hubbard U correction, we found the local density exchange-correlation potential functional can only give a minimized residue of the self-energy and orbital relaxation on the focus orbital if the Janak theorem maintained. Such residue turns to well counteracted in the fully occupied orbitals and non-zero the partially occupied orbitals. With keeping the validation of Janak theorem on localized orbitals, the self-consistent cycle by local density functional with Hubbard U correction cannot find out a set of orbital occupation that simultaneously offsets the orbital self-energy and relaxations in the empty and partially filled shell, but returns a unique set of the occupation for fully occupied shell. The band gap calculations on fully occupied orbital based compounds are thus improved and the relaxed lattices are also shown based on minimization of the self-energy error, which shows a possible route for accurate excited state studies.


title: Tight Span of Path Connected Subsets of the Manhattan Plane

abstract: We show that the tight span of a path connected subset of the Manhattan plane can be constructed in a very simple way by "hatching" the subset in horizontal and vertical directions and then taking the closure of the resulting set.


title: Molecular geometry of alkaloids present in seeds of mexican prickly poppy

abstract: The work is a study of the geometry of the molecules via molecular mechanics of the main alkaloids found in the seeds of Argemone Mexicana Linn, a prickly poppy, which is considered one of the most important species of plants in traditional Mexican and Indian medicine system. The seeds have toxic properties as well as bactericide, hallucinogenic, fungicide, insecticide, in isoquinolines and sanguinarine alkaloids such as berberine. A computational study of the molecular geometry of the molecules through molecular mechanics of the main alkaloids compounds present in plant seeds is described in a computer simulation. The plant has active ingredients compounds: allocryptopine, berberine, chelerythrine, copsitine, dihydrosanguinarine, protopine and sanguinarine. The studied alkaloids form two groups having similar charge distribution among themselves, which have dipole moments of these two times higher than in the other group.


title: Replies to Comments on "Tetraquarks as Diquark-Antidiquark Bound Systems"

abstract: The present paper is written in response to the critical comments of M.R. Hadizadeh on our original paper Tetraquarks as diquark-antidiquark bound systems [Phys. Lett. B 741, 124 (2015), arXiv:1505.07510]. We present our clarifications on his arguments against the accuracy of the procedure and results of the study. The arguments turn out to stem mainly from the critic's misunderstanding of the issues discussed in the original paper.


title: Expansive maps are isometries

abstract: We show that expansive maps from a dense subset of a compact metric space into the metric space itself are isometries


title: On the ranking of a Swiss system chess team tournament

abstract: The paper suggests a family of paired comparison-based scoring procedures for ranking the participants of a Swiss system chess team tournament. We present the challenges of ranking in Swiss system, the features of individual and team competitions as well as the failures of the official rankings based on lexicographical order. The tournament is represented as a ranking problem such that the linearly-solvable row sum (score), generalized row sum, and least squares methods have favourable axiomatic properties. Two chess team European championships are analysed as case studies. Final rankings are compared by their distances and visualized with multidimensional scaling (MDS). Differences to the official ranking are revealed by the decomposition of the least squares method. Rankings are evaluated by prediction power, retrodictive performance, and stability. The paper argues for the use of least squares method with a results matrix favouring match points on the basis of its relative insensitivity to the choice between match and board points, retrodictive accuracy, and robustness.


title: On $n$-norm preservers and the Aleksandrov conservative $n$-distance problem

abstract: The goal of this paper is to point out that the results obtained in the recent papers [7,8,10,11] can be seriously strengthened in the sense that we can significantly relax the assumptions of the main results so that we still get the same conclusions. In order to do this first, we prove that for $n \geq 3$ any transformation which preserves the $n$-norm of any $n$ vectors is automatically plus-minus linear. This will give a re-proof of the well-known Mazur--Ulam-type result that every $n$-isometry is automatically affine ($n \geq 2$) which was proven in several papers, e.g. in [9]. Second, following the work of Rassias and \v{S}emrl [23], we provide the solution of a natural Aleksandrov-type problem in $n$-normed spaces, namely, we show that every surjective transformation which preserves the unit $n$-distance in both directions ($n\geq 2$) is automatically an $n$-isometry.


title: Piercing the Vainshtein screen with anomalous gravitational wave speed: Constraints on modified gravity from binary pulsars

abstract: By using observations of the Hulse-Taylor pulsar we constrain the gravitational wave (GW) speed to the level of $10^{-2}$. We apply this result to scalar-tensor theories that generalize Galileon 4 and 5 models, which display anomalous propagation speed and coupling to matter for GWs. We argue that this effect survives conventional screening due to the persistence of a scalar field gradient inside virialized overdensities, which effectively "pierces" the Vainshtein screening. In specific branches of solutions, our result allows to directly constrain the cosmological couplings in the effective field theory of dark energy formalism.


title: Black Hole - Moving Mirror II: Particle Creation

abstract: There is an exact correspondence between the simplest solution to Einstein's equations describing the formation of a black hole and a particular moving mirror trajectory. In both cases the Bogolubov coefficients in 1+1 dimensions are identical and can be computed analytically. Particle creation is investigated by using wave packets. The entire particle creation history is computed, incorporating the early-time non-thermal emission due to the formation of the black hole (or the early-time acceleration of the moving mirror) and the evolution to a Planckian spectrum.


title: Computer-assisted independent study in mutivariate calculus

abstract: Learning mathematics requires students to work in an independent way which is particularly challenging for such an abstract subject. Advancements in technology and, taking the student as the focus of his own learning, led to a change of paradigm in the design and development of educational contents. In this paper we describe the first experience with an interactive feedback and assessment tool (Siacua), based on parameterized math exercises, and explain how we use it to motivate student independent study in a multivariate calculus environment. We have defined an index about the subject, trying to make it consensual enough for being used in other courses about multivariate calculus. Then we have created a concept map, selected some existing parameterized true/false questions from PmatE project and classified them using our concept map, for being reused in our system. For complementing the course we have created about one hundred parameterized multiple choice question templates in system Megua and generated about one thousand instances for using in Siacua. Results based on data collected by this tool and also based on an informal survey are presented. This first experience allows us to conclude our approach has an important impact on student motivation and contributes to the success on learning multivariate calculus.


title: Six General Ecosystem Properties are more Intense in Biogeochemical Cycling Networks than Food Webs

abstract: Network analysis has revealed whole-network properties hypothesized to be general characteristics of ecosystems including pathway proliferation, and network non-locality, homogenization, amplification, mutualism, and synergism. Collectively these properties characterize the impact of indirect interactions among ecosystem elements. While ecosystem networks generally trace a thermodynamically conserved unit through the system, there appear to be several model classes. For example, trophic (TRO) networks are built upon a food web, usually follow energy or carbon, and are the most abundant in the literature. Biogeochemical cycling (BGC) networks trace nutrients like nitrogen or phosphorus and tend to have more recycling than TRO. We tested (1) the hypothesized generality of the properties in BGC networks and (2) that they tend to be more strongly expressed in BGC networks than in the TRO networks due to increased recycling. We compared the properties in 22 BGC and 57 TRO ecosystem networks from the literature using enaR. The results generally support the hypotheses. First, five of the properties occurred in all 22 BGC models, while network mutualism occurred in 86% of the models. Further, these results were generally robust to a $\pm$50% uncertainty in the model parameters. Second, the mean network statistics for the six properties were statistically significantly greater in the BGC models than the TRO models. These results (1) confirm the general presence of these properties in ecosystem networks, (2) highlight the significance of model types in determining property intensities, (3) reinforce the importance of recycling, and (4) provide a set of indicator benchmarks for future systems comparisons. Further, this work highlights how indirect effects distributed by network connectivity can transform whole-ecosystem functioning, and adds to the growing domain of network ecology.


title: Using a biased qubit to probe complex systems

abstract: Complex mesoscopic systems play increasingly important roles in modern science -- from understanding biological functions at the molecular level, to designing solid-state information processing devices. The operation of these systems typically depends on their energetic structure, yet probing their energy-landscape can be extremely challenging; they have many degrees of freedom, which may be hard to isolate and measure independently. Here we show that a qubit (a two-level quantum system) with a biased energy-splitting can directly probe the spectral properties of a complex system, without knowledge of how they couple. Our work is based on the completely-positive and trace-preserving map formalism, which treats any unknown dynamics as a `black-box' process. This black box contains information about the system with which the probe interacts, which we access by measuring the survival probability of the initial state of the probe as function of the energy-splitting and the process time. Fourier transforming the results yields the energy spectrum of the complex system. Without making assumptions about the strength or form of its coupling, our probe could determine aspects of a complex molecule's energy landscape as well as, in many cases, test for coherent superposition of its energy eigenstates.


title: On the power graphs which are Cayley graphs of some groups

abstract: In 2013, Jemal Abawajy, Andrei Kelarev and Morshed Chowdhury [1] proposed a problem to characterize the finite groups whose power graphs are Cayley graphs of some groups. Here we give a complete answer to this question.


title: Massively Deep Artificial Neural Networks for Handwritten Digit Recognition

abstract: Greedy Restrictive Boltzmann Machines yield an fairly low 0.72% error rate on the famous MNIST database of handwritten digits. All that was required to achieve this result was a high number of hidden layers consisting of many neurons, and a graphics card to greatly speed up the rate of learning.


title: Equivariant Chow classes of matrix orbit closures

abstract: Let $G$ be the product $GL_r(C) \times (C^\times)^n$. We show that the $G$-equivariant Chow class of a $G$ orbit closure in the space of $r$-by-$n$ matrices is determined by a matroid. To do this, we split the natural surjective map from the $G$ equvariant Chow ring of the space of matrices to the torus equivariant Chow ring of the Grassmannian. The splitting takes the class of a Schubert variety to the corresponding factorial Schur polynomial, and also has the property that the class of a subvariety of the Grassmannian is mapped to the class of the closure of those matrices whose row span is in the variety.


title: Pricing American and Asian Options

abstract: An analytic method for pricing American call options is provided; followed by an empirical method for pricing Asian call options. The methodology is the pricing theory presented in "A Modern Theory of Random Variation", by Patrick Muldowney, 2012.


title: On some higher degree sign-definite multivariate polynomials associated with definite quadratic forms

abstract: Positive and negative quadratic forms are well known and widely used. They are multivariate homogeneous polynomials of degree two taking positive or negative values respectively for any values of their arguments not all zero. In the present paper a certain higher degree polynomial is associated with each quadratic form such that the form is definite if and only if this polynomial is sign-definite.


title: Visualizing Probabilistic Proof

abstract: The author revisits the Blue Bus Problem, a famous thought-experiment in law involving probabilistic proof, and presents simple Bayesian solutions to different versions of the blue bus hypothetical. In addition, the author expresses his solutions in standard and visual formats, i.e. in terms of probabilities and natural frequencies.


title: A new entropy based on a group-theoretical structure

abstract: A multi-parametric version of the nonadditive entropy $S_{q}$ is introduced. This new entropic form, denoted by $S_{a,b,r}$, possesses many interesting statistical properties, and it reduces to the entropy $S_q$ for $b=0$, $a=r:=1-q$ (hence Boltzmann-Gibbs entropy $S_{BG}$ for $b=0$, $a=r \to 0$). The construction of the entropy $S_{a,b,r}$ is based on a general group-theoretical approach recently proposed by one of us \cite{Tempesta2}. Indeed, essentially all the properties of this new entropy are obtained as a consequence of the existence of a rational group law, which expresses the structure of $S_{a,b,r}$ with respect to the composition of statistically independent subsystems. Depending on the choice of the parameters, the entropy $S_{a,b,r}$ can be used to cover a wide range of physical situations, in which the measure of the accessible phase space increases say exponentially with the number of particles $N$ of the system, or even stabilizes, by increasing $N$, to a limiting value. This paves the way to the use of this entropy in contexts where a system "freezes" some or many of its degrees of freedom by increasing the number of its constituting particles or subsystems.


title: Corrigendum to: Essential normality, essential norms and hyperrigidity

abstract: In our paper "Essential normality, essential norms and hyperrigidity" we claimed that the restriction of the identity representation of a certain operator system (constructed from a polynomial ideal) has the unique extension property, however the justification we gave was insufficient. In this note we provide the required justification under some additional assumptions. Fortunately, homogeneous ideals that are "sufficiently non-trivial" are covered by these assumptions. This affects the section of our paper relating essential normality and hyperrigidity. We show here that Proposition 4.11 and Theorem 4.12 hold under the additional assumptions. We do not know if they hold in the generality considered in our paper.


title: The origin of net electric currents in solar active regions

abstract: There is a recurring question in solar physics about whether or not electric currents are neutralized in active regions (ARs). This question was recently revisited using three-dimensional (3D) magnetohydrodynamic (MHD) numerical simulations of magnetic flux emergence into the solar atmosphere. Such simulations showed that flux emergence can generate a substantial net current in ARs. Another source of AR currents are photospheric horizontal flows. Our aim is to determine the conditions for the occurrence of net vs. neutralized currents with this second mechanism. Using 3D MHD simulations, we systematically impose line-tied, quasi-static, photospheric twisting and shearing motions to a bipolar potential magnetic field. We find that such flows: (1) produce both {\it direct} and {\it return} currents, (2) induce very weak compression currents - not observed in 2.5D - in the ambient field present in the close vicinity of the current-carrying field, and (3) can generate force-free magnetic fields with a net current. We demonstrate that neutralized currents are in general produced only in the absence of magnetic shear at the photospheric polarity inversion line - a special condition rarely observed. We conclude that, as magnetic flux emergence, photospheric flows can build up net currents in the solar atmosphere, in agreement with recent observations. These results thus provide support for eruption models based on pre-eruption magnetic fields possessing a net coronal current.


title: A Polynomial Time Bounded-error Quantum Algorithm for Boolean Satisfiability

abstract: The aim of the paper is to answer a long-standing open problem on the relationship between NP and BQP. The paper shows that BQP contains NP by proposing a BQP quantum algorithm for the MAX-E3-SAT problem which is a fundamental NP-hard problem. Given an E3-CNF Boolean formula, the aim of the MAX-E3-SAT problem is to find the variable assignment that maximizes the number of satisfied clauses. The proposed algorithm runs in $O(m^2)$ for an E3-CNF Boolean formula with $m$ clauses and in the worst case runs in $O(n^6)$ for an E3-CNF Boolean formula with $n$ inputs. The proposed algorithm maximizes the set of satisfied clauses using a novel iterative partial negation and partial measurement technique. The algorithm is shown to achieve an arbitrary high probability of success of $1-\epsilon$ for small $\epsilon>0$ using a polynomial resources. In addition to solving the MAX-E3-SAT problem, the proposed algorithm can also be used to decide if an E3-CNF Boolean formula is satisfiable or not, which is an NP-complete problem, based on the maximum number of satisfied clauses.


title: Destruction and Observational Signatures of Sun-Impacting Comets

abstract: Motivated by recent data on comets in the low corona, we discuss destruction of sun impacting comets in the dense lower solar atmosphere. Perihelion distances q less than the solar radius and incident masses Mo much greater than 1E12 g are required to reach such depths. Extending earlier work on planetary atmosphere impacts to solar conditions, we evaluate the mechanisms and spatial distribution of nucleus mass and energy loss as functions of Mo and q, and of parameter X = 2Q/CHvovo. Q is the total specific energy for ablative mass loss, CH the bow shock heat transfer efficiency, and vo the solar escape speed (619 km/s). We discuss factors affecting Q and CH and conclude that, for solar vo, X is most likely less than 1 so that solar impactors are mostly ablated before decelerating. Sun impacting comets have kinetic energies 2E30 erg x(Mo/1E15 g), comparable with the energies of magnetic flares. This is released as a localised explosive airburst within a few scale heights H around 200 km of the photosphere, depending weakly on Mo , q and X. For X = 0.01, Mo around 1E15 g, and a typical Kreutz Group shallow incidence angle, comet the airburst occurs around atmospheric density n around 3E15 per ml and this would be 1000 times larger (700 km deeper) for vertical entry. Such airbursts drive flare like phenomena including prompt radiation, hot rising plumes and photospheric ripples, the observability and diagnostic value of which we discuss.


title: An extrapolation cascadic multigrid method combined with a fourth order compact scheme for 3D poisson equation

abstract: In this paper, we develop an EXCMG method to solve the three-dimensional Poisson equation on rectangular domains by using the compact finite difference (FD) method with unequal meshsizes in different coordinate directions. The resulting linear system from compact FD discretization is solved by the conjugate gradient (CG) method with a relative residual stopping criterion. By combining the Richardson extrapolation and tri-quartic Lagrange interpolation for the numerical solutions from two-level of grids (current and previous grids), we are able to produce an extremely accurate approximation of the actual numerical solution on the next finer grid, which can greatly reduce the number of relaxation sweeps needed. Additionally, a simple method based on the midpoint extrapolation formula is used for the fourth-order FD solutions on two-level of grids to achieve sixth-order accuracy on the entire fine grid cheaply and directly. The gradient of the numerical solution can also be easily obtained through solving a series of tridiagonal linear systems resulting from the fourth-order compact FD discretizations. Numerical results show that our EXCMG method is much more efficient than the classical V-cycle and W-cycle multigrid methods. Moreover, only few CG iterations are required on the finest grid to achieve full fourth-order accuracy in both the $L^2$-norm and $L^{\infty}$-norm for the solution and its gradient when the exact solution belongs to $C^6$. Finally, numerical result shows that our EXCMG method is still effective when the exact solution has a lower regularity, which widens the scope of applicability of our EXCMG method.


title: Analytical mass formula and nuclear surface properties in the ETF approximation

abstract: The problem of the determination of the nuclear surface and surface symmetry energy is addressed in the framework of the Extended Thomas Fermi (ETF) approximation using Skyrme functionals. We propose an analytical model for the density profiles with variationally determined diffuseness parameters. For the case of symmetric nuclei, the resulting ETF functional can be exactly integrated, leading to an analytical formula expressing the surface energy as a function of the couplings of the energy functional. The importance of non-local terms is stressed, which cannot be simply deduced from the local part of the functional. In the case of asymmetric nuclei, we propose an approximate expression for the diffuseness and the surface energy. These quantities are analytically related to the parameters of the energy functional. In particular, the influence of the different equation of state parameters can be explicitly quantified. Detailed analyses of the different energy components (local/non-local, isoscalar/isovector, surface/curvature and higher order) are also performed. Our analytical solution of the ETF integral improves over previous models and leads to a precision better than 200 keV per nucleon in the determination of the nuclear binding energy for dripline nuclei.


title: Non-backtracking loop soups and statistical mechanics on spin networks

abstract: We introduce and study a Markov field on the edges of a graph in dimension $d\geq2$ whose configurations are spin networks. The field arises naturally as the edge-occupation field of a Poissonian model (a soup) of non-backtracking loops and walks characterized by a spatial Markov property such that, conditionally on the value of the edge-occupation field on a boundary set, the distributions of the loops and arcs on either side of the boundary are independent of each other. The field has a Gibbs distribution with a Hamiltonian given by a sum of terms which involve only edges incident on the same vertex. Its free energy density and other quantities can be computed exactly, and their critical behavior analyzed, in any dimension.

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


title: Spatially adaptive, Bayesian estimation for probabilistic temperature forecasts

abstract: Uncertainty in the prediction of future weather is commonly assessed through the use of forecast ensembles that employ a numerical weather prediction model in distinct variants. Statistical postprocessing can correct for biases in the numerical model and improves calibration. We propose a Bayesian version of the standard ensemble model output statistics (EMOS) postprocessing method, in which spatially varying bias coefficients are interpreted as realizations of Gaussian Markov random fields. Our Markovian EMOS (MEMOS) technique utilizes the recently developed stochastic partial differential equation (SPDE) and integrated nested Laplace approximation (INLA) methods for computationally efficient inference. The MEMOS approach shows good predictive performance in a comparative study of 24-hour ahead temperature forecasts over Germany based on the 50-member ensemble of the European Centre for Medium-Range Weather Forecasting (ECMWF).


title: High Temperature Asymptotics of Orthogonal Mean-Field Spin Glasses

abstract: We evaluate the high temperature limit of the free energy of spin glasses on the hypercube with Hamiltonian $H_N(\sigma) = \sigma^T J \sigma$, where the coupling matrix $J$ is drawn from certain symmetric orthogonally invariant ensembles. Our derivation relates the annealed free energy of these models to a spherical integral, and expresses the limit of the free energy in terms of the limiting spectral measure of the coupling matrix $J$. As an application, we derive the limiting free energy of the Random Orthogonal Model (ROM) at high temperatures, which confirms non-rigorous calculations of Marinari et al. (1994). Our methods also apply to other well-known models of disordered systems, including the SK and Gaussian Hopfield models.


title: First Physics Results at the Physical Pion Mass from $N_f = 2$ Wilson Twisted Mass Fermions at Maximal Twist

abstract: We present physics results from simulations of QCD using $N_f = 2$ dynamical Wilson twisted mass fermions at the physical value of the pion mass. These simulations were enabled by the addition of the clover term to the twisted mass quark action. We show evidence that compared to previous simulations without this term, the pion mass splitting due to isospin breaking is almost completely eliminated. Using this new action, we compute the masses and decay constants of pseudoscalar mesons involving the dynamical up and down as well as valence strange and charm quarks at one value of the lattice spacing, $a \approx 0.09$ fm. Further, we determine renormalized quark masses as well as their scale-independent ratios, in excellent agreement with other lattice determinations in the continuum limit. In the baryon sector, we show that the nucleon mass is compatible with its physical value and that the masses of the $\Delta$ baryons do not show any sign of isospin breaking. Finally, we compute the electron, muon and tau lepton anomalous magnetic moments and show the results to be consistent with extrapolations of older ETMC data to the continuum and physical pion mass limits. We mostly find remarkably good agreement with phenomenology, even though we cannot take the continuum and thermodynamic limits.


title: Automorphisms of $p$-local compact groups

abstract: Self equivalences of classifying spaces of $p$-local compact groups are well understood by means of the algebraic structure that gives rise to them, but explicit descriptions are lacking. In this paper we use a construction of Robinson of an amalgam $G$, realizing a given fusion system, to produce a split epimorphism from the outer automorphism group of $G$ to the group of homotopy classes of self homotopy equivalences of the classifying space of the corresponding $p$-local compact group.


title: One-loop triple collinear splitting amplitudes in QCD

abstract: We study the factorisation properties of one-loop scattering amplitudes in the triple collinear limit and extract the universal splitting amplitudes for processes initiated by a gluon. The splitting amplitudes are derived from the analytic Higgs plus four partons amplitudes. We present compact results for primitive helicity splitting amplitudes making use of super-symmetric decompositions. The universality of the collinear factorisation is checked numerically against the full colour six parton squared matrix elements.


title: Crowding induced clustering under confinement

abstract: We present Langevin dynamics simulations that study the collective behavior of driven particles embedded in a densely packed background consisting of passive particles. Depending on the driving force, the densities of driven and passive particles, and the temperature we observe a dynamical phase separation of the driven particles, which cluster together in tight bands. We explore the mechanisms that drive this cluster formation, and determine the critical conditions for such phase separation. A simple physical picture explains the formation and subsequent growth of a jammed zone developing in front of the driven cluster. The model correctly captures the observed scaling with time. We analyze the implications of this clustering transition for the driven transport in dense particulate flows, which due to a non-monotonic dependence on the applied driving force is not straightforwardly optimized. We provide proof-of-concept for a direct application of the clustering effect, and propose a 'colloidal chromatograph'; a setup that permits the separation of colloids by mass or size.


title: The Random-Diluted Triangular Plaquette Model: study of phase transitions in a Kinetically Constrained Model

abstract: We study how the thermodynamic properties of the Triangular Plaquette Model (TPM) are influenced by the addition of extra interactions. The thermodynamics of the original TPM is trivial, while its dynamics is glassy, as usual in Kinetically Constrained Models. As soon as we generalize the model to include additional interactions, a thermodynamic phase transition appears in the system. The additional interactions we consider are either short ranged, forming a regular lattice in the plane, or long ranged of the small-world kind. In the case of long-range interactions we call the new model Random-Diluted TPM. We provide arguments that the model so modified should undergo a thermodynamic phase transition, and that in the long-range case this is a glass transition of the "Random First-Order" kind. Finally, we give support to our conjectures studying the finite temperature phase diagram of the Random-Diluted TPM in the Bethe approximation. This corresponds to the exact calculation on the random regular graph, where free-energy and configurational entropy can be computed by means of the cavity equations.


title: Sequential Quantiles via Hermite Series Density Estimation

abstract: Sequential quantile estimation refers to incorporating observations into quantile estimates in an incremental fashion thus furnishing an online estimate of one or more quantiles at any given point in time. Sequential quantile estimation is also known as online quantile estimation. This area is relevant to the analysis of data streams and to the one-pass analysis of massive data sets. Applications include network traffic and latency analysis, real time fraud detection and high frequency trading. We introduce new techniques for online quantile estimation based on Hermite series estimators in the settings of static quantile estimation and dynamic quantile estimation. In the static quantile estimation setting we apply the existing Gauss-Hermite expansion in a novel manner. In particular, we exploit the fact that Gauss-Hermite coefficients can be updated in a sequential manner. To treat dynamic quantile estimation we introduce a novel expansion with an exponentially weighted estimator for the Gauss-Hermite coefficients which we term the Exponentially Weighted Gauss-Hermite (EWGH) expansion. These algorithms go beyond existing sequential quantile estimation algorithms in that they allow arbitrary quantiles (as opposed to pre-specified quantiles) to be estimated at any point in time. In doing so we provide a solution to online distribution function and online quantile function estimation on data streams. In particular we derive an analytical expression for the CDF and prove consistency results for the CDF under certain conditions. In addition we analyse the associated quantile estimator. Simulation studies and tests on real data reveal the Gauss-Hermite based algorithms to be competitive with a leading existing algorithm.


title: Stochastic Loewner Evolution Relates Anomalous Diffusion and Anisotropic Percolation

abstract: We disclose the origin of anisotropic percolation perimeters in terms of the Stochastic Loewner Evolution (SLE) process. Precisely, our results from extensive numerical simulations indicate that the perimeters of multi-layered and directed percolation clusters at criticality are the scaling limits of the Loewner evolution of an anomalous Brownian motion, being subdiffusive and superdiffusive, respectively. The connection between anomalous diffusion and fractal anisotropy is further tested by using long-range power-law correlated time series (fractional Brownian motion) as driving functions in the evolution process. The fact that the resulting traces are distinctively anisotropic corroborates our hypothesis. Under the conceptual framework of SLE, our study therefore reveals new perspectives for mathematical and physical interpretations of non-Markovian processes in terms of anisotropic paths at criticality and vice-versa.


title: The External Magnetic Field Created by the Superposition of Identical Parallel Finite Solenoids

abstract: Using superposition and numerical approximations of a published analytical expression for the magnetic field generated by a finite solenoid, we show that the magnetic field external to parallel identical solenoids can be nearly uniform and substantial, even when the solenoids have lengths that are large compared to their radii. We study two arrangements of solenoids---a ring of parallel solenoids whose surfaces are tangent to a common cylindrical surface and to nearest neighbours, and a large finite hexagonal array of parallel solenoids---and summarize how the magnitude and uniformity of the resultant external field depend on the solenoid length and distances between solenoids. We also report some novel results about single solenoids, e.g., that the energy stored in the internal magnetic field exceeds the energy stored in the spatially infinite external magnetic field for even short solenoids. These results should be broadly interesting to undergraduates learning about electricity and magnetism as novel examples of superposition based on a familiar source of magnetic fields.


title: Precision age indicators that exploit chemically peculiar stars

abstract: We would like to find a way to improve the determination of galaxy star formation history from integrated light spectroscopy. To this end, several classes of chemically peculiar (CP) stars arise during the course of normal evolution in single stars and noninteracting binary stars. An aging stellar population has periods of time in which CP stars contribute to the integrated light, and others in which the contributions fade. The HgMn stars, for example, occupy a narrow temperature range of 10500 to 16000 K, which maps to a narrow range of ages. Wolf-Rayet stars, He-poor stars, Bp-Ap stars, Am-Fm stars, and C stars all become very common in a normal stellar population at various ages between zero and several Gyr, fading in and out in a way that is analogous to features used in stellar spectral classification. We examine population fractions and light fractions in order to assess the feasibility of using CP stars as age tracers. We find that, even though CP stars do not usually dominate in number, there are enough of them so that the CP spectral features are detectable in high-quality integrated spectra of young and intermediate age stellar populations. The new technique should be calibratable and useful. Furthermore, using CP signatures as age dating tools sidesteps reliance on photometry that is susceptible to dust and Balmer features that are susceptible to nebular fill-in.


title: Magnetotransport in Dirac metals: chiral magnetic effect and quantum oscillations

abstract: Dirac metals are characterized by the linear dispersion of fermionic quasi-particles, with the Dirac point hidden inside a Fermi surface. We study the magnetotransport in these materials using chiral kinetic theory to describe within the same framework both the negative magnetoresistance caused by chiral magnetic effect and quantum oscillations in the magnetoresistance due to the existence of the Fermi surface. We discuss the relevance of obtained results to recent measurements on ${\rm Cd_3As_2}$.


title: A Non-Equilibrium Kinetic Theory for Trapped Binary Condensates

abstract: We derive a non-equilibrium finite-temperature kinetic theory for a binary mixture of two interacting atomic Bose-Einstein condensates and use it to explore the degree of hydrodynamicity attainable in realistic experimental geometries. Based on the standard separation of timescale argument of kinetic theory, the dynamics of the condensates of the multi-component system are shown to be described by dissipative Gross-Pitaevskii equations, self-consistently coupled to corresponding Quantum Boltzmann equations for the non-condensate atoms: on top of the usual mean field contributions, our scheme identifies a total of eight distinct collisional processes, whose dynamical interplay is expected to be responsible for the systems equilibration. In order to provide their first characterization, we perform a detailed numerical analysis of the role of trap frequency and geometry on collisional rates for experimentally accessible mixtures of $^{87}$Rb-$^{41}$K and $^{87}$Rb-$^{85}$Rb, discussing the extent to which the system may approach the hydrodynamic regime with regard to some of those processes, as a guide for future experimental investigations of ultracold Bose gas mixtures.


title: Riemann Manifold Langevin Methods on Stochastic Volatility Estimation

abstract: In this paper we perform Bayesian estimation of stochastic volatility models with heavy tail distributions using Metropolis adjusted Langevin (MALA) and Riemman manifold Langevin (MMALA) methods. We provide analytical expressions for the application of these methods, assess the performance of these methodologies in simulated data and illustrate their use on two financial time series data sets.


title: Primes represented by incomplete norm forms

abstract: Let $K=\mathbb{Q}(\omega)$ with $\omega$ the root of a degree $n$ monic irreducible polynomial $f\in\mathbb{Z}[X]$. We show the degree $n$ polynomial $N(\sum_{i=1}^{n-k}x_i\omega^{i-1})$ in $n-k$ variables formed by setting the final $k$ coefficients to 0 takes the expected asymptotic number of prime values if $n\ge 4k$. In the special case $K=\mathbb{Q}(\sqrt[n]{\theta})$, we show $N(\sum_{i=1}^{n-k}x_i\sqrt[n]{\theta^{i-1}})$ takes infinitely many prime values provided $n\ge 22k/7$. Our proof relies on using suitable `Type I' and `Type II' estimates in Harman's sieve, which are established in a similar overall manner to the previous work of Friedlander and Iwaniec on prime values of $X^2+Y^4$ and of Heath-Brown on $X^3+2Y^3$. Our proof ultimately relies on employing explicit elementary estimates from the geometry of numbers and algebraic geometry to control the number of highly skewed lattices appearing in our final estimates.


title: Toward a new language of legal drafting

abstract: Lawyers should write in document markup language just like web developers, digital publishers, scientists, and almost everyone else


title: Gromov-Hausdorff limits of Kahler manifolds and algebraic geometry, II

abstract: In this paper we continue to study Gromov-Hausdorff limits of Kahler manifolds and algebraic geometry. Our main focus is on the algebro-geometric meaning of Riemannian tangent cones and rescaled limits.


title: Identifying Codes in the Complementary Prism of Cycles

abstract: We show that an identifying code of minimum order in the complementary prism of a cycle of order $n$ has order $7n/9+\Theta(1)$. Furthermore, we observe that the clique-width of the complementary prism of a graph of clique-width $k$ is at most $4k$, and discuss some algorithmic consequences.


title: Detection of magnetic fields in both B-type components of the $\epsilon$ Lupi system: a new constraint on the origin of fossil fields?

abstract: High-resolution circular spectropolarimetric observations, obtained with ESPaDOnS in the context of the BinaMIcS Large Program, have revealed a magnetic field in the B3V secondary component of the SB2 binary system $\epsilon$ Lupi (B2/B3). As the B2V primary is already known to be magnetic, this is the first detection of a magnetic field in both components of an early-type binary system. The longitudinal magnetic field of the primary is $\sim -200$ G; that of the secondary $\sim +100$ G. Observations can be approximately reproduced by a model assuming the magnetic axes of the two stars are anti-aligned, and roughly parallel to their respective rotation axes. Estimated magnetospheric radii indicate a high probability that their magnetospheres are interacting. As many of the arguments for the different proposed formation scenarios of fossil magnetic fields rely upon evidence drawn from investigations of close binaries, in particular the rarity of magnetic ABO stars in close binaries and the previous absence of any known close binary with two magnetic, massive stars, this discovery may be an important new constraint on the origin of fossil magnetic fields.


title: Secure management of logs in internet of things

abstract: Ever since the advent of computing, managing data has been of extreme importance. With innumerable devices getting added to network infrastructure, there has been a proportionate increase in the data which needs to be stored. With the advent of Internet of Things (IOT) it is anticipated that billions of devices will be a part of the internet in another decade. Since those devices will be communicating with each other on a regular basis with little or no human intervention, plethora of real time data will be generated in quick time which will result in large number of log files. Apart from complexity pertaining to storage, it will be mandatory to maintain confidentiality and integrity of these logs in IOT enabled devices. This paper will provide a brief overview about how logs can be efficiently and securely stored in IOT devices.


title: Parallel Correlation Clustering on Big Graphs

abstract: Given a similarity graph between items, correlation clustering (CC) groups similar items together and dissimilar ones apart. One of the most popular CC algorithms is KwikCluster: an algorithm that serially clusters neighborhoods of vertices, and obtains a 3-approximation ratio. Unfortunately, KwikCluster in practice requires a large number of clustering rounds, a potential bottleneck for large graphs. We present C4 and ClusterWild!, two algorithms for parallel correlation clustering that run in a polylogarithmic number of rounds and achieve nearly linear speedups, provably. C4 uses concurrency control to enforce serializability of a parallel clustering process, and guarantees a 3-approximation ratio. ClusterWild! is a coordination free algorithm that abandons consistency for the benefit of better scaling; this leads to a provably small loss in the 3-approximation ratio. We provide extensive experimental results for both algorithms, where we outperform the state of the art, both in terms of clustering accuracy and running time. We show that our algorithms can cluster billion-edge graphs in under 5 seconds on 32 cores, while achieving a 15x speedup.


title: Type I and Type II Bayesian Methods for Sparse Signal Recovery using Scale Mixtures

abstract: In this paper, we propose a generalized scale mixture family of distributions, namely the Power Exponential Scale Mixture (PESM) family, to model the sparsity inducing priors currently in use for sparse signal recovery (SSR). We show that the successful and popular methods such as LASSO, Reweighted $\ell_1$ and Reweighted $\ell_2$ methods can be formulated in an unified manner in a maximum a posteriori (MAP) or Type I Bayesian framework using an appropriate member of the PESM family as the sparsity inducing prior. In addition, exploiting the natural hierarchical framework induced by the PESM family, we utilize these priors in a Type II framework and develop the corresponding EM based estimation algorithms. Some insight into the differences between Type I and Type II methods is provided and of particular interest in the algorithmic development is the Type II variant of the popular and successful reweighted $\ell_1$ method. Extensive empirical results are provided and they show that the Type II methods exhibit better support recovery than the corresponding Type I methods.


title: 4f fine-structure levels as the dominant error in the electronic structures of binary lanthanide oxides

abstract: The ground-state 4f fine-structure levels in the intrinsic optical transition gaps between the 2p and 5d orbitals of lanthanide sesquioxides (Ln2O3, Ln=La...Lu) were calculated by a two-way crossover search for the U parameters for DFT+U calculations. The original 4f-shell potential perturbation in the linear response method were reformulated within the constraint volume of thegiven solids. The band structures were also calculated. This method yields nearly constant optical transition gaps between Ln-5d and O-2p orbitals, with magnitudes of 5.3~5.5 eV. This result verifies that the error in the band structure calculations for Ln2O3 is dominated by the inaccuracies in the predicted 4f levels in the 2p-5d transition gaps, which strongly and nonlinearly depend on the on-site Hubbard U. The relationship between the 4f occupancies and Hubbard U is non-monotonic and is entirely different from that for materials with 3d or 4d orbitals, such as transition metal oxides. This new linear response DFT+U method can provide a simpler understanding of the electronic structure of Ln2O3 and enables a quick examination of the electronic structures of lanthanide solids prior to hybrid functional or GW calculations.


title: Supercurrent Reversal in Two-Dimensional Topological Insulators

abstract: We calculate supercurrent across a two-dimensional topological insulator subjected to an external magnetic field. When the edge states of a narrow two-dimensional topological insulator are hybridized, an external magnetic field can close the hybridization gap, thus driving a quantum phase transition from insulator to semimetal states of the topological insulator. We find a sign reversal of the supercurrent at the quantum phase transition revealing intrinsic properties of topological insulators via Josephson effect.


title: Weak lensing by galaxy troughs in DES Science Verification data

abstract: We measure the weak lensing shear around galaxy troughs, i.e. the radial alignment of background galaxies relative to underdensities in projections of the foreground galaxy field over a wide range of redshift in Science Verification data from the Dark Energy Survey. Our detection of the shear signal is highly significant (10 to 15$\sigma$ for the smallest angular scales) for troughs with the redshift range z in [0.2,0.5] of the projected galaxy field and angular diameters of 10 arcmin...1{\deg}. These measurements probe the connection between the galaxy, matter density, and convergence fields. By assuming galaxies are biased tracers of the matter density with Poissonian noise, we find agreement of our measurements with predictions in a fiducial Lambda cold dark matter model. The prediction for the lensing signal on large trough scales is virtually independent of the details of the underlying model for the connection of galaxies and matter. Our comparison of the shear around troughs with that around cylinders with large galaxy counts is consistent with a symmetry between galaxy and matter over- and underdensities. In addition, we measure the two-point angular correlation of troughs with galaxies which, in contrast to the lensing signal, is sensitive to galaxy bias on all scales. The lensing signal of troughs and their clustering with galaxies is therefore a promising probe of the statistical properties of matter underdensities and their connection to the galaxy field.


title: Silicon vacancy center in 4H-SiC: Electronic structure and spin-photon interfaces

abstract: Defects in silicon carbide are of intense and increasing interest for quantum-based applications due to this material's properties and technological maturity. We calculate the multi-particle symmetry adapted wave functions of the negatively charged silicon vacancy defect in hexagonal silicon carbide via use of group theory and density functional theory and find the effects of spin-orbit and spin-spin interactions on these states. Although we focused on $\textrm{V}_{\textrm{Si}}^-$ in 4H-SiC, because of its unique fine structure due to odd number of active electrons, our methods can be easily applied to other defect centers of different polytpes, especially to the 6H-SiC. Based on these results we identify the mechanism that polarizes the spin under optical drive, obtain the ordering of its dark doublet states, point out a path for electric field or strain sensing, and find the theoretical value of its ground-state zero field splitting to be 68 MHz, in good agreement with experiment. Moreover, we present two distinct protocols of a spin-photon interface based on this defect. Our results pave the way toward novel quantum information and quantum metrology applications with silicon carbide.


title: Interloper bias in future large-scale structure surveys

abstract: Next-generation spectroscopic surveys will map the large-scale structure of the observable universe, using emission line galaxies as tracers. While each survey will map the sky with a specific emission line, interloping emission lines can masquerade as the survey's intended emission line at different redshifts. Interloping lines from galaxies that are not removed can contaminate the power spectrum measurement, mixing correlations from various redshifts and diluting the true signal. We assess the potential for power spectrum contamination, finding that an interloper fraction worse than 0.2% could bias power spectrum measurements for future surveys by more than 10% of statistical errors, while also biasing power spectrum inferences. We also construct a formalism for predicting cosmological parameter bias, demonstrating that a 0.15%-0.3% interloper fraction could bias the growth rate by more than 10% of the error, which can affect constraints on gravity upcoming surveys. We use the COSMOS Mock Catalog (CMC), with the emission lines re-scaled to better reproduce recent data, to predict potential interloper fractions for the Prime Focus Spectrograph (PFS) and the Wide-Field InfraRed Survey Telescope (WFIRST). We find that secondary line identification, or confirming galaxy redshifts by finding correlated emission lines, can remove interlopers for PFS. For WFIRST, we use the CMC to predict that the 0.2% target can be reached for the WFIRST H$\alpha$ survey, but sensitive optical and near-infrared photometry will be required. For the WFIRST [OIII] survey, the predicted interloper fractions reach several percent and their effects will have to be estimated and removed statistically (e.g. with deep training samples). (Abridged)


title: Higgs mass predictions of public NMSSM spectrum generators

abstract: The publicly available spectrum generators for the NMSSM often lead to different predictions for the mass of the standard model-like Higgs boson even if using the same renormalization scheme and two-loop accuracy. Depending on the parameter point, the differences can exceed 5 GeV, and even reach 8 GeV for moderate superparticle masses of up to 2 TeV. It is shown here that these differences can be traced back to the calculation of the running standard model parameters entering all calculations, to the approximations used in the two-loop corrections included in the different codes, and to different choices for the renormalization conditions and scales. In particular, the importance of the calculation of the top Yukawa coupling is pointed out.


title: Mid-IR Enhanced Galaxies in the Coma & Virgo Cluster: lenticulars with a high star formation rate

abstract: We explore the properties of early-type galaxies (ETGs), including ellipticals (E) and lenticulars (S0), in rich environments such as clusters of galaxies (Virgo and Coma). The L_24/L_K distribution of ETGs in both Virgo and Coma clusters shows that some S0s have a much larger L_24/L_K ratio (0.5 to ~2 dex) than the bulk of the ETG population. This could be interpreted as an enhanced star formation rate in these lenticulars. We compare the optical colors of galaxies in these two clusters and investigate the nature of these sources with a large L24/L_K ratio by looking at their spatial distribution within the cluster, by analyzing their optical spectra and by looking at their optical colors compared to late-types. We obtain 10 Coma and 3 Virgo early-type sources with larger L24/L_K ratios than the bulk of their population. We call these sources Mid-Infrared Enhanced Galaxies (MIEGs). In Coma, they are mostly located in the South-West part of the cluster where a substructure is falling onto the main cluster. MIEGs present lower g-r color than the rest of the ETG sample, because of a blue continuum. We interpret the excess L24/L_K ratio as evidence for an enhanced star-formation induced as a consequence of their infall into the main cluster.


title: The Three-mm Ultimate Mopra Milky Way Survey. I. Survey Overview, Initial Data Releases, and First Results

abstract: We describe a new mm-wave molecular-line mapping survey of the southern Galactic Plane and its first data releases. The Three-mm Ultimate Mopra Milky Way Survey (ThrUMMS) maps a 60{\deg}x2{\deg} sector of our Galaxy's fourth quadrant, using a combination of fast mapping techniques with the Mopra radio telescope, simultaneously in the J=1-0 lines of $^{12}$CO, $^{13}$CO, C$^{18}$O, and CN near 112 GHz at ~arcminute and ~0.3 km s$^{-1}$ resolution, with ~2 K channel$^{-1}$ sensitivity for $^{12}$CO and ~1 K channel$^{-1}$ for the other transitions. The calibrated data cubes from these observations are made available to the community after processing through our pipeline. Here, we describe the motivation for ThrUMMS, the development of new observing techniques for Mopra, and how these techniques were optimised to the objectives of the survey. We showcase some sample data products and describe the first science results on CO-isotopologue line ratios. These vary dramatically across the Galactic Plane, indicating a very wide range of optical depth and excitation conditions, from warm and translucent to cold and opaque. The population of cold clouds in particular have optical depths for $^{12}$CO easily exceeding 100. We derive a new, nonlinear conversion law from $^{12}$CO integrated intensity to column density, which suggests that the molecular mass traced by CO in the Galactic disk may have been substantially underestimated. This further suggests that some global relationships in disk galaxies, such as star formation laws, may need to be recalibrated. The large ThrUMMS team is proceeding with several other science investigations.


title: Muon g-2 and galactic centre $\gamma$-ray excess in a scalar extension of the 2HDM type-X

abstract: We consider an extension of the lepto-specific 2HDM with an extra singlet $S$ as a dark matter candidate. Taking into account theoretical and experimental constraints, we investigate the possibility to address both the $\gamma$-ray excess detected at the Galactic centre and the discrepancy between the Standard Model prediction and experimental results of the anomalous magnetic moment of the muon. Our analyses reveal that the $SS \to \tau^+ \tau^-$ and $SS \to b \bar b$ channels reproduce the Galactic centre excess, with an emerging dark matter candidate which complies with the bounds from direct detection experiments, measurements of the Higgs boson invisible decay width and observations of the dark matter relic abundance. Addressing the anomalous magnetic moment of the muon imposes further strong constraints on the model. Remarkably, under these conditions, the $SS \to b \bar b$ channel still allows for the fitting of the Galactic centre. We also comment on a scenario allowed by the model where the $SS \to \tau^+ \tau^-$ and $SS \to b \bar b$ channels have comparable branching ratios, which possibly yield an improved fitting of the Galactic centre excess.


title: Measuring Transit Signal Recovery in the Kepler Pipeline II: Detection Efficiency as Calculated in One Year of Data

abstract: The Kepler planet sample can only be used to reconstruct the underlying planet occurrence rate if the detection efficiency of the Kepler pipeline is known, here we present the results of a second experiment aimed at characterising this detection efficiency. We inject simulated transiting planet signals into the pixel data of ~10,000 targets, spanning one year of observations, and process the pixels as normal. We compare the set of detections made by the pipeline with the expectation from the set of simulated planets, and construct a sensitivity curve of signal recovery as a function of the signal-to-noise of the simulated transit signal train. The sensitivity curve does not meet the hypothetical maximum detection efficiency, however it is not as pessimistic as some of the published estimates of the detection efficiency. For the FGK stars in our sample, the sensitivity curve is well fit by a gamma function with the coefficients a = 4.35 and b = 1.05. We also find that the pipeline algorithms recover the depths and periods of the injected signals with very high fidelity, especially for periods longer than 10 days. We perform a simplified occurrence rate calculation using the measured detection efficiency compared to previous assumptions of the detection efficiency found in the literature to demonstrate the systematic error introduced into the resulting occurrence rates. The discrepancies in the calculated occurrence rates may go some way towards reconciling some of the inconsistencies found in the literature.


title: The Lexicocalorimeter: Gauging public health through caloric input and output on social media

abstract: We propose and develop a Lexicocalorimeter: an online, interactive instrument for measuring the "caloric content" of social media and other large-scale texts. We do so by constructing extensive yet improvable tables of food and activity related phrases, and respectively assigning them with sourced estimates of caloric intake and expenditure. We show that for Twitter, our naive measures of "caloric input", "caloric output", and the ratio of these measures are all strong correlates with health and well-being measures for the contiguous United States. Our caloric balance measure in many cases outperforms both its constituent quantities, is tunable to specific health and well-being measures such as diabetes rates, has the capability of providing a real-time signal reflecting a population's health, and has the potential to be used alongside traditional survey data in the development of public policy and collective self-awareness. Because our Lexicocalorimeter is a linear superposition of principled phrase scores, we also show we can move beyond correlations to explore what people talk about in collective detail, and assist in the understanding and explanation of how population-scale conditions vary, a capacity unavailable to black-box type methods.


title: A geometrical approach to gravitational lensing magnification

abstract: The standard definition of gravitational lensing magnification is generalized to Lorentzian spacetimes, and it is shown how it can be interpreted geometrically in terms of the van Vleck determinant and the exponential map. This is joint work with Amir B. Aazami (Kavli IPMU, University of Tokyo).


title: Tree-level metastability bounds for the most general two Higgs doublet model

abstract: Within two Higgs doublet models, it is possible that the current vacuum is not the global minimum, in which case it could possibly decay at a later stage. We discuss the tree-level conditions which must be obeyed by the most general scalar potential in order to preclude that possibility. We propose a new procedure which is not only more general but also easier to implement than the previously published one, including CP conserving as well as CP violating scalar sectors. We illustrate these conditions within the context of the Z2 model, softly broken by a complex, CP violating parameter.