title: Riemann problems with non--local point constraints and capacity drop

abstract: In the present note we discuss in details the Riemann problem for a one--dimensional hyperbolic conservation law subject to a point constraint. We investigate how the regularity of the constraint operator impacts the well--posedness of the problem, namely in the case, relevant for numerical applications, of a discretized exit capacity. We devote particular attention to the case in which the constraint is given by a non--local operator depending on the solution itself. We provide several explicit examples. We also give the detailed proof of some results announced in the paper [Andreainov, Donadello, Rosini, "Crowd dynamics and conservation laws with non--local point constraints and capacity drop", which is devoted to existence and stability for a more general class of Cauchy problems subject to Lipschitz continuous non--local point constraints.


title: PT-symmetry breaking with divergent potentials: lattice and continuum cases

abstract: We investigate the parity- and time-reversal ($\mathcal{PT}$)-symmetry breaking in lattice models in the presence of long-ranged, non-hermitian, $\mathcal{PT}$-symmetric potentials that remain finite or become divergent in the continuum limit. By scaling analysis of the fragile $\mathcal{PT}$ threshold for an open finite lattice, we show that continuum loss-gain potentials $V_\alpha(x)\propto i |x|^\alpha \mathrm{sign}(x)$ have a positive $\mathcal{PT}$-breaking threshold for $\alpha>-2$, and a zero threshold for $\alpha\leq -2$. When $\alpha<0$ localized states with complex (conjugate) energies in the continuum energy-band occur at higher loss-gain strengths. We investigate the signatures of $\mathcal{PT}$-symmetry breaking in coupled waveguides, and show that the emergence of localized states dramatically shortens the relevant time-scale in the $\mathcal{PT}$-symmetry broken region.


title: G\'en\'eralisations Quantitatives du Crit\'ere D'ind\'ependance Lin\'eaire De Nesterenko

abstract: In this paper we extend Fischler's quantitative generalization of Nesterenko's linear independence criterion, by weakening the hypotheses on the divisors of the coe cients of the linear forms and allowing (to some extent) the linear forms not to tend to 0. Another version of this result is proved, in the spirit of Siegel's criterion, with a recurrence relation veri ed by the linear forms. Finally, the results are restated in a more general setting in terms of convex bodies and lattices of $\mathbb{R}^n$.


title: A reversible infinite HMM using normalised random measures

abstract: We present a nonparametric prior over reversible Markov chains. We use completely random measures, specifically gamma processes, to construct a countably infinite graph with weighted edges. By enforcing symmetry to make the edges undirected we define a prior over random walks on graphs that results in a reversible Markov chain. The resulting prior over infinite transition matrices is closely related to the hierarchical Dirichlet process but enforces reversibility. A reinforcement scheme has recently been proposed with similar properties, but the de Finetti measure is not well characterised. We take the alternative approach of explicitly constructing the mixing measure, which allows more straightforward and efficient inference at the cost of no longer having a closed form predictive distribution. We use our process to construct a reversible infinite HMM which we apply to two real datasets, one from epigenomics and one ion channel recording.


title: DISCO: a Spatio-Spectral Recombiner for Pupil Remapping Interferometry

abstract: Pupil-remapping is a new high-dynamic range imaging technique that has recently demonstrated feasibility on sky. The current prototypes present however deceiving limiting magnitude, restricting the current use to the brightest stars in the sky. We propose to combine pupil-remapping with spatio-spectral encoding, a technique first applied to the VEGA/CHARA interferometer. The result is an instrument proposal, called "Dividing Interferometer for Stars Characterizations and Observations" (DISCO). The idea is to take profit of wavelength multiplexing when using a spectrograph in order to pack as much as possible the available information, yet providing a potential boost of 1.5 magnitude if used in existing prototypes. We detail in this paper the potential of such a concept.


title: Macroscopic quantum tunneling and quantum-classical phase transitions of the escape rate in large spin systems

abstract: This article presents a review on the theoretical and the experimental developments on macroscopic quantum tunneling and phase transition of the escape rate in spin systems. We present the basic ideas with simplified calculations so that it is readable to both specialists and nonspecialists in this area of research. A brief derivation of the path integral formulation of quantum mechanics in its original form using the orthonormal position and momentum basis is reviewed. For spin systems such as single molecule magnets, the formulation of path integral requires the use of non-orthonormal spin coherent state in $(2s+1)$ dimensional Hilbert space, the coordinate independent and the coordinate dependent form of the spin coherent state path integral is derived. These two forms of spin coherent state path integral are applied to the tunneling of single molecule magnets through its magnetic anisotropy barrier. Most experimental and numerical results are presented. The suppression of tunneling for half-odd integer spin (spin-parity effect) at zero magnetic field is derived from both forms, which shows that this result (spin-parity effect) is independent of the coordinate. At nonzero magnetic field we present both the experimental and the theoretical results of the oscillation of tunneling splitting as a function of the applied magnetic field applied along the spin hard anisotropy axis direction. The experimental and the theoretical results of the tunneling in antiferromagnetic exchange coupled dimer model are also reviewed. As the spin coherent state path integral formalism is a semi-classical method, an alternative exact mapping of a spin system to a particle in a potential field (effective potential method) is derived. This effective potential method allows for the investigation of phase transition of the escape rate in spin systems.


title: High-Performance Image Synthesis for Radio Interferometry

abstract: A radio interferometer indirectly measures the intensity distribution of the sky over the celestial sphere. Since measurements are made over an irregularly sampled Fourier plane, synthesising an intensity image from interferometric measurements requires substantial processing. Furthermore there are distortions that have to be corrected. In this thesis, a new high-performance image synthesis tool (imaging tool) for radio interferometry is developed. Implemented in C++ and CUDA, the imaging tool achieves unprecedented performance by means of Graphics Processing Units (GPUs). The imaging tool is divided into several components, and the back-end handling numerical calculations is generalised in a new framework. A new feature termed compression arbitrarily increases the performance of an already highly efficient GPU-based implementation of the w-projection algorithm. Compression takes advantage of the behaviour of oversampled convolution functions and the baseline trajectories. A CPU-based component prepares data for the GPU which is multi-threaded to ensure maximum use of modern multi-core CPUs. Best performance can only be achieved if all hardware components in a system do work in parallel. The imaging tool is designed such that disk I/O and work on CPU and GPUs is done concurrently. Test cases show that the imaging tool performs nearly 100$\times$ faster than another general CPU-based imaging tool. Unfortunately, the tool is limited in use since deconvolution and A-projection are not yet supported. It is also limited by GPU memory. Future work will implement deconvolution and A-projection, whilst finding ways of overcoming the memory limitation.


title: K\"ahler manifolds of semi-negative holomorphic sectional curvature

abstract: In an earlier work, we investigated some consequences of the existence of a K\"ahler metric of negative holomorphic sectional curvature on a projective manifold. In the present work, we extend our results to the case of semi-negative (i.e., non-positive) holomorphic sectional curvature. In doing so, we define a new invariant that records the largest codimension of maximal subspaces in the tangent spaces on which the holomorphic sectional curvature vanishes. Using this invariant, we establish lower bounds for the nef dimension and, under certain additional assumptions, for the Kodaira dimension of the manifold. In dimension two, a precise structure theorem is obtained.


title: Evolution of Water Reservoirs on Mars: Constraints from Hydrogen Isotopes in Martian Meteorites

abstract: Martian surface morphology implies that Mars was once warm enough to maintain persistent liquid water on its surface. While the high D/H ratios (~6 times the Earth's ocean water) of the current martian atmosphere suggest that significant water has been lost from the surface during martian history, the timing, processes, and the amount of the water loss have been poorly constrained. Recent technical developments of ion-microprobe analysis of martian meteorites have provided accurate estimation of hydrogen isotope compositions (D/H) of martian water reservoirs at the time when the meteorites formed. Based on the D/H data from the meteorites, this study demonstrates that the water loss during the pre-Noachian (>41-99 m global equivalent layers, GEL) was more significant than in the rest of martian history (>10-53 m GEL). Combining our results with geological and geomorphological evidence for ancient oceans, we propose that undetected subsurface water/ice (~100-1000 m GEL) should have existed, and it exceeds the observable present water inventory (~20-30 m GEL) on Mars.


title: Managing gamma_5 in Dimensional Regularization and ABJ Anomaly

abstract: An integral representation is proposed for the trace involving gamma_5 in dimensional regularization. Lorentz covariance is preserved. ABJ anomaly naturally follows. The Local Functional Equation associated to x-dependent chiral transformations is verified.


title: Conformal Carroll groups

abstract: Conformal extensions of Levy-Leblond's Carroll group, based on geometric properties analogous to those of Newton-Cartan space-time are proposed. The extensions are labelled by an integer $k$. This framework includes and extends our recent study of the Bondi-Metzner-Sachs (BMS) and Newman-Unti (NU) groups. The relation to Conformal Galilei groups is clarified. Conformal Carroll symmetry is illustrated by "Carrollian photons". Motion both in the Newton-Cartan and Carroll spaces may be related to that of strings in the Bargmann space.


title: Destruction of wide binary stars in low mass elliptical galaxies: implications for initial mass function estimates

abstract: We discuss the effects of destruction of wide binaries in the nuclei of the lower mass giant elliptical galaxies. We show that the numbers of barium stars and extrinsic S stars should be dramatically reduced in these galaxies compared to what is seen in the largest elliptical galaxies. Given that the extrinsic S stars show strong Wing-Ford band and Na I D absorption, we argue that the recent claims of different initial mass functions from the most massive elliptical galaxies versus lower mass ellipticals may be the result of extrinsic S stars, rather than bottom-heavy initial mass function.


title: Galaxy Shapes and Intrinsic Alignments in The MassiveBlack-II Simulation

abstract: The intrinsic alignment of galaxy shapes with the large-scale density field is a contaminant to weak lensing measurements, as well as being an interesting signature of galaxy formation and evolution (albeit one that is difficult to predict theoretically). Here we investigate the shapes and relative orientations of the stars and dark matter of halos and subhalos (central and satellite) extracted from the MassiveBlack-II simulation, a state-of-the-art high resolution hydrodynamical cosmological simulation which includes stellar and AGN feedback in a volume of $(100{h^{-1}\mathrm{Mpc}})^3$. We consider redshift evolution from $z=1$ to $0.06$ and mass evolution within the range of subhalo masses, $10^{10} -6.0 \times 10^{14.0}{h^{-1}M_{\odot}}$. The shapes of the dark matter distributions are generally more round than the shapes defined by stellar matter. The projected root-mean-square (RMS) ellipticity per component for stellar matter is measured to be $e_\text{rms} = 0.28$ at $z=0.3$ for $M_{subhalo}> 10^{12.0}{h^{-1}M_{\odot}}$, which compares favourably with observational measurements. We find that the shapes of stellar and dark matter are more round for less massive subhalos and at lower redshifts. By directly measuring the relative orientation of the stellar matter and dark matter of subgroups, we find that, on average, the misalignment between the two components is larger for less massive subhalos. The mean misalignment angle varies from $\sim 30^{\circ}-10^{\circ}$ for $M \sim 10^{10} - 10^{14} {h^{-1}M_{\odot}}$ and shows a weak dependence on redshift. We also compare the misalignment angles in central and satellite subhalos at fixed subhalo mass, and find that centrals are more misaligned than satellites. We present fitting formulae for the shapes of dark and stellar matter in subhalos and also the probability distributions of misalignment angles.


title: Tensor Detection Severely Constrains Axion Dark Matter

abstract: The recent detection of B-modes by BICEP2 has non-trivial implications for axion dark matter implied by combining the tensor interpretation with isocurvature constraints from Planck. In this paper the measurement is taken as fact, and its implications considered, though further experimental verification is required. In the simplest inflation models $r=0.2$ implies $H_I=1.1\times 10^{14}\text{ GeV}$. If the axion decay constant $f_a<H_I/2\pi$ constraints on the dark matter (DM) abundance alone rule out the QCD axion as DM for $m_a \lesssim 52\chi^{6/7}\,\mu\text{eV}$ (where $\chi>1$ accounts for theoretical uncertainty). If $f_a>H_I/2\pi$ then vacuum fluctuations of the axion field place conflicting demands on axion DM: isocurvature constraints require a DM abundance which is too small to be reached when the back reaction of fluctuations is included. High $f_a$ QCD axions are thus ruled out. Constraints on axion-like particles, as a function of their mass and DM fraction, are also considered. For heavy axions with $m_a\gtrsim 10^{-22}\text{ eV}$ we find $\Omega_a/\Omega_d\lesssim 10^{-3}$, with stronger constraints on heavier axions. Lighter axions, however, are allowed and (inflationary) model-independent constraints from the CMB temperature power spectrum and large scale structure are stronger than those implied by tensor modes.


title: Socio-economic applications of finite state mean field games

abstract: In this paper we present different applications of finite state mean field games to socio-economic sciences. Examples include paradigm shifts in the scientific community or the consumer choice behaviour in the free market. The corresponding finite state mean field game models are hyperbolic systems of partial differential equations, for which we present and validate different numerical methods. We illustrate the behaviour of solutions with various numerical experiments, which show interesting phenomena like shock formation. Hence we conclude with an investigation of the shock structure in the case of two-state problems.


title: Review of COMPASS results on transverse-spin effects in SIDIS

abstract: The transversity parton distribution remains a poorly known cornerstone in the nucleon spin structure. While the Collins effect in spin asymmetries in Semi-Inclusive DIS (SIDIS) is one crucial tool to address the transversity function, the most promising alternative is the azimuthal asymmetry in SIDIS when a hadron pair is detected in the final state. In this case, the chiral-odd transversity function is coupled to another chiral-odd function, i.e. the hadron-pair interference fragmentation function (IFF). The measurement of azimuthal asymmetries in hadron-pair production on a transversely polarised nucleon target has been performed at COMPASS using a 160 GeV/c muon beam of CERN's M2 beam line. Results from the 2007 and 2010 recent measurements will be presented and compared to model predictions.


title: Electrodynamic modeling of strong coupling between a metasurface and intersubband transitions in quantum wells

abstract: Strong light-matter coupling has recently been demonstrated in sub-wavelength volumes by coupling engineered optical transitions in semiconductor heterostructures (e.g., quantum wells) to metasurface resonances via near fields. It has also been shown that different resonator shapes may lead to different Rabi splittings, though this has not yet been well explained. In this paper, our aim is to understand the correlation between resonator shape and Rabi splitting, and in particular determine and quantify the physical parameters that affect strong coupling by developing an equivalent circuit network model whose elements describe energy and dissipation. Because of the subwavelength dimension of each metasurface element, we resort to the quasi-static (electrostatic) description of the near-field and hence define an equivalent capacitance associated to each dipolar element of a flat metasurface, and we show that this is also able to accurately model the phenomenology involved in strong coupling between the metasurface and the intersubband transitions in quantum wells. We show that the spectral properties and stored energy of a metasurface/quantum-well system obtained using our model are in good agreement with both full-wave simulation and experimental results. We then analyze metasurfaces made of three different resonator geometries and observe that the magnitude of the Rabi splitting increases with the resonator capacitance in agreement with our theory, providing a phenomenological explanation for the resonator shape dependence of the strong coupling process.


title: The Jenkins-Serrin problem for constant mean curvature graphs in the Heisenberg space $\mathbf{Nil_3}(\tau)$

abstract: In this paper we find functions over bounded domains in the 2-dimensional Euclidean space, whose graphs (in the Heisenberg space) has constant mean curvature different from zero and taking on (possibly) infinite boundary values over the boundary of the domain.


title: On bundles that admit fiberwise hyperbolic dynamics

abstract: This paper is devoted to rigidity of smooth bundles which are equipped with fiberwise geometric or dynamical structure. We show that the fiberwise associated sphere bundle to a bundle whose leaves are equipped with (continuously varying) metrics of negative curvature is a topologically trivial bundle when either the base space is simply connected or, more generally, when the bundle is fiber homotopically trivial. We present two very different proofs of this result: a geometric proof and a dynamical proof. We also establish a number of rigidity results for bundles which are equipped with fiberwise Anosov dynamical systems. Finally, we present a number of examples which show that our results are sharp in certain ways or illustrate necessity of various assumptions.


title: General parity-odd CMB bispectrum estimation

abstract: We develop a methodology for estimating parity-odd bispectra in the cosmic microwave background (CMB). This is achieved through the extension of the original separable modal methodology to parity-odd bispectrum domains ($\ell_1 + \ell_2 + \ell_3 = {\rm odd}$). Through numerical tests of the parity-odd modal decomposition with some theoretical bispectrum templates, we verify that the parity-odd modal methodology can successfully reproduce the CMB bispectrum, without numerical instabilities. We also present simulated non-Gaussian maps produced by modal-decomposed parity-odd bispectra, and show the consistency with the exact results. Our new methodology is applicable to all types of parity-odd temperature and polarization bispectra.


title: On the distribution of interspecies correlation for Markov models of character evolution on Yule trees

abstract: Efforts to reconstruct phylogenetic trees and understand evolutionary processes depend fundamentally on stochastic models of speciation and mutation. The simplest continuous-time model for speciation in phylogenetic trees is the Yule process, in which new species are "born" from existing lineages at a constant rate. Recent work has illuminated some of the structural properties of Yule trees, but it remains mostly unknown how these properties affect sequence and trait patterns observed at the tips of the phylogenetic tree. Understanding the interplay between speciation and mutation under simple models of evolution is essential for deriving valid phylogenetic inference methods and gives insight into the optimal design of phylogenetic studies. In this work, we derive the probability distribution of interspecies covariance under Brownian motion and Ornstein-Uhlenbeck models of phenotypic change on a Yule tree. We compute the probability distribution of the number of mutations shared between two randomly chosen taxa in a Yule tree under discrete Markov mutation models. Our results suggest summary measures of phylogenetic information content, illuminate the correlation between site patterns in sequences or traits of related organisms, and provide heuristics for experimental design and reconstruction of phylogenetic trees.


title: Learning Negative Mixture Models by Tensor Decompositions

abstract: This work considers the problem of estimating the parameters of negative mixture models, i.e. mixture models that possibly involve negative weights. The contributions of this paper are as follows. (i) We show that every rational probability distributions on strings, a representation which occurs naturally in spectral learning, can be computed by a negative mixture of at most two probabilistic automata (or HMMs). (ii) We propose a method to estimate the parameters of negative mixture models having a specific tensor structure in their low order observable moments. Building upon a recent paper on tensor decompositions for learning latent variable models, we extend this work to the broader setting of tensors having a symmetric decomposition with positive and negative weights. We introduce a generalization of the tensor power method for complex valued tensors, and establish theoretical convergence guarantees. (iii) We show how our approach applies to negative Gaussian mixture models, for which we provide some experiments.


title: Thermal photon $v_2$ with slow quark chemical equilibration

abstract: Elliptic flow of direct photons in high-energy heavy ion collisions has been a topic of great interest as it is experimentally found larger than most hydrodynamic expectations. I discuss the implication of possible late formation of quark component in the hot QCD medium on the photon elliptic flow as quarks are the source of thermal photons in the deconfined phase. Hydrodynamic equations are numerically solved with the evolution equations for quark and gluon number densities. The numerical results suggest that thermal photon $v_2$ is visibly enhanced by the slow chemical equilibration of quarks and gluons, reducing the aforementioned problem.


title: Agravity

abstract: We explore the possibility that the fundamental theory of nature does not contain any scale. This implies a renormalizable quantum gravity theory where the graviton kinetic term has 4 derivatives, and can be reinterpreted as gravity minus an anti-graviton. We compute the super-Planckian RGE of adimensional gravity coupled to a generic matter sector. The Planck scale and a flat space can arise dynamically at quantum level provided that a quartic scalar coupling and its $\beta$ function vanish at the Planck scale. This is how the Higgs boson behaves for $M_h\approx 125$ GeV at $M_t\approx 171$ GeV. Within agravity, inflation is a generic phenomenon: the slow-roll parameters are given by the $\beta$-functions of the theory, and are small if couplings are perturbative. The predictions $n_s\approx 0.967$ and $r\approx 0.13$ arise if the inflaton is identified with the Higgs of gravity. Furthermore, quadratically divergent corrections to the Higgs mass vanish: a small weak scale is natural and can be generated by agravity quantum corrections.


title: A Scanning, All-Fiber Sagnac Interferometer for High Resolution Magneto-Optic Measurements at 820 nm

abstract: The Sagnac Interferometer has historically been used for detecting non-reciprocal phenomena, such as rotation. We demonstrate an apparatus in which this technique is employed for high resolution measurements of the Magneto-Optical Polar Kerr effect--a direct indicator of magnetism. Previous designs have incorporated free-space components which are bulky and difficult to align. We improve upon this technique by using all fiber-optic coupled components and demonstrate operation at a new wavelength, 820 nm, with which we can achieve better than 1 $\mu$rad resolution. Mounting the system on a piezo-electric scanner allows us to acquire diffraction limited images with 1.5 $\mu$m spatial resolution. We also provide extensive discussion on the details and of the Sagnac Interferometer's construction.


title: Consistency Tests of Classical and Quantum Models for a Quantum Annealer

abstract: Recently the question of whether the D-Wave processors exhibit large-scale quantum behavior or can be described by a classical model has attracted significant interest. In this work we address this question by studying a 503 qubit D-Wave Two device in the "black box" model, i.e., by studying its input-output behavior. Our work generalizes an approach introduced in Boixo et al. [Nat. Commun. 4, 2067 (2013)], and uses groups of up to 20 qubits to realize a transverse Ising model evolution with a ground state degeneracy whose distribution acts as a sensitive probe that distinguishes classical and quantum models for the D-Wave device. Our findings rule out all classical models proposed to date for the device and provide evidence that an open system quantum dynamical description of the device that starts from a quantized energy level structure is well justified, even in the presence of relevant thermal excitations and a small value of the ratio of the single-qubit decoherence time to the annealing time.


title: Infinite Systems of Competing Brownian Particles

abstract: Consider a system of infinitely many Brownian particles on the real line. At any moment, these particles can be ranked from the bottom upward. Each particle moves as a Brownian motion with drift and diffusion coefficients depending on its current rank. The gaps between consecutive particles form the (infinite-dimensional) gap process. We find a stationary distribution for the gap process. We also show that if the initial value of the gap process is stochastically larger than this stationary distribution, this process converges back to this distribution as time goes to infinity. This continues the work by Pal and Pitman (2008). Also, this includes infinite systems with asymmetric collisions, similar to the finite ones from Karatzas, Pal and Shkolnikov (2016).


title: VIMOS Ultra-Deep Survey (VUDS): Witnessing the Assembly of a Massive Cluster at z~3.3

abstract: Using new spectroscopic observations obtained as part of the VIMOS Ultra-Deep Survey (VUDS), we perform a systematic search for overdense environments in the early universe ($z>2$) and report here on the discovery of Cl J0227-0421, a massive protocluster at $z=3.29$. This protocluster is characterized by both the large overdensity of spectroscopically confirmed members, $\delta_{gal}=10.5\pm2.8$, and a significant overdensity in photometric redshift members. The halo mass of this protocluster is estimated, by a variety of methods, to be roughly $3\times10^{14}$ $\mathcal{M}_{\odot}$ at $z\sim3.3$, which, evolved to $z=0$ results in a halo mass rivaling or exceeding that of the Coma cluster. The properties of 19 spectroscopically confirmed member galaxies are compared with a large sample of VUDS/VVDS galaxies in lower density field environments at similar redshifts. We find tentative evidence for an excess of redder, brighter, and more massive galaxies within the confines of the protocluster relative to the field population, which suggests that we may be observing the beginning of environmentally-induced quenching. The properties of these galaxies are investigated, including a discussion of the brightest protocluster galaxy which appears to be undergoing vigorous coeval nuclear and starburst activity. The remaining member galaxies appear to have characteristics which are largely similar to the field population. Though we find weaker evidence of the suppression of the median star formation rates amongst and differences in stacked spectra of member galaxies with respect to the field, we defer any conclusions of these trends to future work with the ensemble of protostructures that are found in the full VUDS sample.


title: Semiquantisation Functor and Poisson-Riemannian Geometry, I

abstract: We study noncommutative bundles and Riemannian geometry at the semiclassical level of first order in a deformation parameter $\lambda$, using a functorial approach. The data for quantisation of the cotangent bundle is known to be a Poisson structure and Poisson preconnection and we now show that this data defines to a functor $Q$ from the monoidal category of classical vector bundles equipped with connections to the monodial category of bimodules equipped with bimodule connections over the quantised algebra. We adapt this functor to quantise the wedge product of the exterior algebra and in the Riemannian case, the metric and the Levi-Civita connection. Full metric compatibility requires vanishing of an obstruction in the classical data, expressed in terms of a generalised Ricci 2-form, without which our quantum Levi-Civita connection is still the best possible. We apply the theory to the Schwarzschild black-hole and to Riemann surfaces as examples, as well as verifying our results on the 2D bicrossproduct model quantum spacetime. The quantized Schwarzschild black-hole in particular has features similar to those encountered in $q$-deformed models, notably the necessity of nonassociativity of any rotationally invariant quantum differential calculus of classical dimensions.


title: Automatic Image Registration in Infrared-Visible Videos using Polygon Vertices

abstract: In this paper, an automatic method is proposed to perform image registration in visible and infrared pair of video sequences for multiple targets. In multimodal image analysis like image fusion systems, color and IR sensors are placed close to each other and capture a same scene simultaneously, but the videos are not properly aligned by default because of different fields of view, image capturing information, working principle and other camera specifications. Because the scenes are usually not planar, alignment needs to be performed continuously by extracting relevant common information. In this paper, we approximate the shape of the targets by polygons and use affine transformation for aligning the two video sequences. After background subtraction, keypoints on the contour of the foreground blobs are detected using DCE (Discrete Curve Evolution)technique. These keypoints are then described by the local shape at each point of the obtained polygon. The keypoints are matched based on the convexity of polygon's vertices and Euclidean distance between them. Only good matches for each local shape polygon in a frame, are kept. To achieve a global affine transformation that maximises the overlapping of infrared and visible foreground pixels, the matched keypoints of each local shape polygon are stored temporally in a buffer for a few number of frames. The matrix is evaluated at each frame using the temporal buffer and the best matrix is selected, based on an overlapping ratio criterion. Our experimental results demonstrate that this method can provide highly accurate registered images and that we outperform a previous related method.


title: Noise correlations of two-dimensional Bose gases

abstract: We analyze density-density correlations of expanding clouds of weakly interacting two-dimensional Bose gases below and above the Berezinskii-Kosterlitz-Thouless transition, with particular focus on short-time expansions. During time-of-flight expansion, phase fluctuations of the trapped system translate into density fluctuations, in addition to the density fluctuations that exist in in-situ. We calculate the correlations of these fluctuations both in real space and in momentum space, and derive analytic expressions in momentum space. Below the transition, the correlation functions show an oscillatory behavior, controlled by the scaling exponent of the quasi-condensed phase, due to constructive interference. We argue that this can be used to extract the scaling exponent of the quasi-condensate experimentally. Above the transition, the interference is rapidly suppressed when the atoms travel an average distance beyond the correlation length. This can be used to distinguish the two phases qualitatively.


title: Dust aerosol, clouds, and the atmospheric optical depth record over 5 Mars years of the Mars Exploration Rover mission

abstract: Dust aerosol plays a fundamental role in the behavior and evolution of the Martian atmosphere. The first five Mars years of Mars Exploration Rover data provide an unprecedented record of the dust load at two sites. This record is useful for characterization of the atmosphere at the sites and as ground truth for orbital observations. Atmospheric extinction optical depths have been derived from solar images after calibration and correction for time-varying dust that has accumulated on the camera windows. The record includes local, regional, and globally extensive dust storms. Comparison with contemporaneous thermal infrared data suggests significant variation in the size of the dust aerosols, with a 1 {\mu}m effective radius during northern summer and a 2 {\mu}m effective radius at the onset of a dust lifting event. The solar longitude (LS) 20-136{\deg} period is also characterized by the presence of cirriform clouds at the Opportunity site, especially near LS=50 and 115{\deg}. In addition to water ice clouds, a water ice haze may also be present, and carbon dioxide clouds may be present early in the season. Variations in dust opacity are important to the energy balance of each site, and work with seasonal variations in insolation to control dust devil frequency at the Spirit site.


title: Entanglement and reduction in two-photon correlations

abstract: For a number of basic experiments on two-photon intensity correlations it is pointed out that the results, which are usually explained in terms of the formalism of canonical field quantization, can also be explained in terms of quantum wave packets, where reduction and entanglement are taken into account.


title: On $p$-adic Hard Sphere model with three states on a Cayley tree

abstract: We study $p$-adic model of hard spheres with three states on the Cayley tree. We show that there exist three translation-invariant $p$-adic Gibbs measures and two periodic measures on a Cayley tree of oreder two.


title: Movable high Q nanoresonators realized by semiconductor nanowires on a Si photonic crystal platform

abstract: Subwavelength semiconductor nanowires have recently attracted interest for photonic applications because they possess various unique optical properties and offer great potential for miniaturizing devices. However, realizing tight light confinement or efficient coupling with photonic circuits is not straightforward and remains a challenge. Here we show that a high Q nanocavity can be created by placing a single III/V semiconductor nanowire with a diameter of under 100 nm in a grooved waveguide in a Si photonic crystal, by means of nanoprobe manipulation. We observe very fast spontaneous emission (91 ps) from nanowires accelerated by the strong Purcell enhancement in nanocavities, which proves that very strong light confinement can be achieved. Furthermore, this system enables us to move the nanocavity anywhere along the waveguide. This configuration provides a significant degree of flexibility in integrated photonics and permits the addition and displacement of various functionalities of III/V nanocavity devices in Si photonic circuits.


title: Computer Vision Accelerators for Mobile Systems based on OpenCL GPGPU Co-Processing

abstract: In this paper, we present an OpenCL-based heterogeneous implementation of a computer vision algorithm -- image inpainting-based object removal algorithm -- on mobile devices. To take advantage of the computation power of the mobile processor, the algorithm workflow is partitioned between the CPU and the GPU based on the profiling results on mobile devices, so that the computationally-intensive kernels are accelerated by the mobile GPGPU (general-purpose computing using graphics processing units). By exploring the implementation trade-offs and utilizing the proposed optimization strategies at different levels including algorithm optimization, parallelism optimization, and memory access optimization, we significantly speed up the algorithm with the CPU-GPU heterogeneous implementation, while preserving the quality of the output images. Experimental results show that heterogeneous computing based on GPGPU co-processing can significantly speed up the computer vision algorithms and makes them practical on real-world mobile devices.


title: On the eigenvalues of some non-Hermitian Hamiltonians with space-time symmetry

abstract: We calculate the eigenvalues of some two-dimensional non-Hermitian Hamiltonians by means of a pseudospectral method and straightforward diagonalization of the Hamiltonian matrix in a suitable basis set. Both sets of results agree remarkably well but differ considerably from the eigenvalues obtained some time ago by other authors. In particular, we do not observe the multiple phase transitions claimed to occur in one of the anharmonic oscillators.


title: Multi-epoch spectropolarimetry of SN 2009ip: direct evidence for aspherical circumstellar material

abstract: We present spectropolarimetry of SN 2009ip throughout the evolution of its 2012 explosion. During the initial 2012a phase, when the source spectrum exhibits broad P-Cygni lines, we measure a V-band polarization of P~0.9% at a position angle of ~166 deg, indicating substantial asphericity for the 2012a outflow. Near the subsequent peak of the 2012b phase, when the spectrum shows signs of intense interaction with circumstellar material (CSM), we measure P~1.7% at a position angle of 72 deg, indicating a separate physical component of polarization, which has a higher degree of asphericity than the 2012a outflow and an orthogonal axis of symmetry on the sky. Around 30 days past peak, which is coincident with a bump in the declining light curve, we measure P~0.7% and significant variations in P across some broad lines, particularly HeI/NaI. By 60 days past peak the source appears to be approaching a low value of interstellar polarization (P<0.2%). The results are consistent with a scenario in which a potentially bipolar explosion during 2012a impacts a toroidal distribution of CSM, thereby causing the 2012b brightening. Orthogonality between the 2012a and 2012b geometries is inconsistent with the hypothesis that the 2012a event launched the CSM that was hit 40 days later by material from the 2012b event. Rather, the aspherical CSM probably has a separate origin associated with the outbursts of the stellar progenitor during prior years. Previous calculations that assumed spherical symmetry for the CSM have underestimated the required explosion energy, as the results suggest that <10% of the fast SN ejecta participated in strong CSM interaction during the 2012b phase. In light of the spectropolarimetric results, a kinetic energy of 1e51 erg for the ejecta is difficult to avoid, supporting the interpretation that the 2012 outburst was the result of a core-collapse supernova explosion (abridged)


title: WFIRST Ultra-Precise Astrometry I: Kuiper Belt Objects

abstract: I show that the WFIRST microlensing survey will enable detection and precision orbit determination of Kuiper Belt Objects (KBOs) down to H_=28.2 over an effective area of about 17 sq.deg. Typical fractional period errors will be about 1.5% x 10^(0.4(H-28.2)) with similar errors in other parameters for roughly 5000 KBOs. Binary companions to detected KBOs can be detected to even fainter limits, H_vega=29, corresponding to $R\sim 30.5$ and effective diameters D~7 km. For KBOs H~23, binary companions can be found with separations down to 10 mas. This will provide an unprecedented probe of orbital resonance and KBO mass measurements. More than a thousand stellar occultations by KBOs can be combined to determine the mean size as a function of KBO magnitude down to H~25. Current ground-based microlensing surveys can make a significant start on finding and characterizing KBOs using existing and soon-to-be-acquired data.


title: Stabilizing the false vacuum: Mott skyrmions

abstract: Topological excitations keep fascinating physicists since many decades. While individual vortices and solitons emerge and have been observed in many areas of physics, their most intriguing higher dimensional topological relatives, skyrmions (smooth, topologically stable textures) and magnetic monopoles -- emerging almost necessarily in any grand unified theory and responsible for charge quantization -- remained mostly elusive. Here we propose that loading a three-component nematic superfluid such as $^{23}$Na into a deep optical lattice and thereby creating an insulating core, one can create topologically stable skyrmion textures and investigate their properties in detail. We show furthermore that the spectrum of the excitations of the superfluid and their quantum numbers change dramatically in the presence of the skyrmion, and they reflect the presence of a trapped monopole, as imposed by the skyrmion's topology.


title: Emergent Cosmology Revisited

abstract: We explore the possibility of emergent cosmology using the effective potential formalism. We discover new models of emergent cosmology which satisfy the constraints posed by the cosmic microwave background (CMB). We demonstrate that, within the framework of modified gravity, the emergent scenario can arise in a universe which is spatially open/closed. By contrast, in general relativity (GR) emergent cosmology arises from a spatially closed past-eternal Einstein Static Universe (ESU). In GR the ESU is unstable, which creates fine tuning problems for emergent cosmology. However, modified gravity models including Braneworld models, Loop Quantum Cosmology (LQC) and Asymptotically Free Gravity result in a stable ESU. Consequently, in these models emergent cosmology arises from a larger class of initial conditions including those in which the universe eternally oscillates about the ESU fixed point. We demonstrate that such an oscillating universe is necessarily accompanied by graviton production. For a large region in parameter space graviton production is enhanced through a parametric resonance, casting serious doubts as to whether this emergent scenario can be past-eternal.


title: Irradiation Instability at the Inner Edges of Accretion Disks

abstract: An instability can potentially operate in highly irradiated disks where the disk sharply transitions from being radially transparent to opaque (the 'transition region'). Such conditions may exist at the inner edges of transitional disks around T Tauri stars and accretion disks around AGNs. We derive the criterion for this instability, which we term the 'irradiation instability', or IRI. We also present the linear growth rate as a function of beta, the ratio between radiation force and gravity, and c_s, the sound speed of the disk, obtained using two methods: a semi-analytic analysis of the linearized equations and a numerical simulation using the GPU-accelerated hydrodynamical code PEnGUIn. In particular, we find that IRI occurs at beta~0.1 if the transition region extends as wide as ~0.05r, and at higher beta values if it is wider. This threshold value applies to c_s ranging from 3% of the Keplerian orbital speed to 5%, and becomes higher if c_s is lower. Furthermore, in the nonlinear evolution of the instability, disks with a large beta and small c_s exhibit 'clumping', extreme local surface density enhancements that can reach over ten times the initial disk surface density.


title: Cosmic Reionization On Computers I. Design and Calibration of Simulations

abstract: Cosmic Reionization On Computers (CROC) is a long-term program of numerical simulations of cosmic reionization. Its goal is to model fully self-consistently (albeit not necessarily from the first principles) all relevant physics, from radiative transfer to gas dynamics and star formation, in simulation volumes of up to 100 comoving Mpc, and with spatial resolution approaching 100 pc in physical units. In this method paper we describe our numerical method, the design of simulations, and the calibration of numerical parameters. Using several sets (ensembles) of simulations in 20 Mpc/h and 40 Mpc/h boxes with spatial resolution reaching 125 pc at z=6, we are able to match the observed galaxy UV luminosity functions at all redshifts between 6 and 10, as well as obtain reasonable agreement with the observational measurements of the Gunn-Peterson optical depth at z<6.


title: X-ray constraints on the local supermassive black hole occupation fraction

abstract: Distinct seed formation mechanisms are imprinted upon the fraction of dwarf galaxies currently containing a central supermassive black hole. Seeding by Pop III remnants is expected to produce a higher occupation fraction than is generated with direct gas collapse precursors. Chandra observations of nearby early-type galaxies can directly detect even low-level supermassive black hole activity, and the active fraction immediately provides a firm lower limit to the occupation fraction. Here, we use the volume-limited AMUSE surveys of ~200 optically-selected early-type galaxies to characterize simultaneously, for the first time, the occupation fraction and the scaling of nuclear X-ray luminosity with stellar mass, accounting for intrinsic scatter, measurement uncertainties, and X-ray limits. For early-type galaxies with log(M_star/M_sun)<10, we obtain a lower limit to the occupation fraction of >20% (at 95% confidence), but full occupation cannot be excluded. The preferred dependence of log(L_X) upon log(M_star) has a slope of about 0.7-0.8, consistent with the "downsizing" trend previously identified from the AMUSE dataset, and a uniform Eddington efficiency is disfavored at ~2 sigma. We provide guidelines for the future precision with which these parameters may be refined with larger or more sensitive samples.


title: The scale height of gas traced by [CII] in the Galactic plane

abstract: The distribution of various interstellar gas components and the pressure in the interstellar medium (ISM) is a result of the interplay of different dynamical mechanisms and energy sources on the gas in the Milky Way. The scale heights of the different gas tracers, such as HI and CO, are a measure of these processes. The scale height of [CII] emission in the Galactic plane is important for understanding those ISM components not traced by CO or HI. We determine the average distribution of [CII] perpendicular to the plane in the inner Galactic disk and compare it to the distributions of other key gas tracers, such as CO and HI. We calculated the vertical, z, distribution of [CII] in the inner Galactic disk by adopting a model for the emission that combines the latitudinal, b, spectrally unresolved BICE survey, with the spectrally resolved $Herschel$ Galactic plane survey of [CII] at b = 0 deg. Our model assumed a Gaussian emissivity distribution vertical to the plane, and related the distribution in z to that of the latitude b using the spectrally resolved [CII] Herschel survey as the boundary solution for the emissivity at b=0 deg. We find that the distribution of [CII] perpendicular to the plane has a full-width half-maximum of 172 pc, larger than that of CO, which averages ~110 pc in the inner Galaxy, but smaller than that of HI, ~230 pc, and is offset by -28 pc. We explain the difference in distributions of [CII], CO, and HI as due to [CII] tracing a mix of ISM components. Models of hydrostatic equilibrium of clouds in the disk predict different scale heights, for the same interstellar pressure. The diffuse molecular clouds with [CII] but no CO emission likely have a scale height intermediate between the low density atomic hydrogen HI clouds and the dense CO molecular clouds.


title: Divergent Thermopower without a Quantum Phase Transition

abstract: A general principle of modern statistical physics is that divergences of either thermodynamic or transport properties are only possible if the correlation length diverges. We show by explicit calculation that the thermopower in the quantum XY model $d=1+1$ and the Kitaev model in $d=2+1$ can 1) diverge even when the correlation length is finite and 2) remain finite even when the correlation length diverges, thereby providing a counterexample to the standard paradigm.Two conditions are necessary: 1) the sign of the charge carriers and that of the group velocity must be uncorrelated and 2) the current operator defined formally as the derivative of the Hamiltonian with respect to the gauge field does not describe a set of excitations that have a particle interpretation, as in strongly correlated electron matter. The recent experimental\cite{2dtp} and theoretical\cite{kirkpatrick} findings on the divergent thermopower of a 2D electron gas are discussed in this context.


title: A $10^{10}$ Solar Mass Flow of Molecular Gas in the Abell 1835 Brightest Cluster Galaxy

abstract: We report ALMA Early Science observations of the Abell 1835 brightest cluster galaxy (BCG) in the CO (3-2) and CO (1-0) emission lines. We detect $5\times 10^{10}~\rm M_\odot$ of molecular gas within 10 kpc of the BCG. Its ensemble velocity profile width of $\sim 130 ~\rm km~s^{-1}$ FWHM is too narrow for the molecular cloud sto be supported in the galaxy by dynamic pressure. The gas may instead be supported in a rotating, turbulent disk oriented nearly face-on. Roughly $10^{10}~\rm M_\odot$ of molecular gas is projected $3-10 ~\rm kpc$ to the north-west and to the east of the nucleus with line of sight velocities lying between $-250 ~\rm km~s^{-1}$ to $+480 ~\rm km~s^{-1}$ with respect to the systemic velocity. The high velocity gas may be either inflowing or outflowing. However, the absence of high velocity gas toward the nucleus that would be expected in a steady inflow, and its bipolar distribution on either side of the nucleus, are more naturally explained as outflow. Star formation and radiation from the AGN are both incapable of driving an outflow of this magnitude. If so, the molecular outflow may be associated a hot outflow on larger scales reported by Kirkpatrick and colleagues. The molecular gas flow rate of approximately $200~\rm M_\odot ~yr^{-1}$ is comparable to the star formation rate of $100-180~\rm M_\odot ~yr^{-1}$ in the central disk. How radio bubbles would lift dense molecular gas in their updrafts, how much gas will be lost to the BCG, and how much will return to fuel future star formation and AGN activity are poorly understood. Our results imply that radio-mechanical (radio mode) feedback not only heats hot atmospheres surrounding elliptical galaxies and BCGs, it is able to sweep higher density molecular gas away from their centers.


title: Constraints on the Progenitor System of the Type Ia Supernova 2014J from Pre-Explosion Hubble Space Telescope Imaging

abstract: We constrain the properties of the progenitor system of the highly reddened Type Ia supernova (SN) 2014J in Messier 82 (M82; d ~ 3.5 Mpc). We determine the SN location using Keck-II K-band adaptive optics images, and we find no evidence for flux from a progenitor system in pre-explosion near-ultraviolet through near-infrared Hubble Space Telescope (HST) images. Our upper limits exclude systems having a bright red giant companion, including symbiotic novae with luminosities comparable to that of RS Ophiuchi. While the flux constraints are also inconsistent with predictions for comparatively cool He-donor systems (T < ~35,000 K), we cannot preclude a system similar to V445 Puppis. The progenitor constraints are robust across a wide range of R_V and A_V values, but significantly greater values than those inferred from the SN light curve and spectrum would yield proportionally brighter luminosity limits. The comparatively faint flux expected from a binary progenitor system consisting of white dwarf stars would not have been detected in the pre-explosion HST imaging. Infrared HST exposures yield more stringent constraints on the luminosities of very cool (T < 3000 K) companion stars than was possible in the case of SN Ia 2011fe.


title: Cosmic Reionization On Computers II. Reionization History and Its Back Reaction on Early Galaxies

abstract: We compare the results from several sets of cosmological simulations of cosmic reionization, produced under Cosmic Reionization On Computers (CROC) project, with existing observational data on the high-redshift Ly-alpha forest and the abundance of Ly-alpha emitters. We find good consistency with the observational measurements and the previous simulation work. By virtue of having several independent realizations for each set of numerical parameters, we are able to explore the effect of cosmic variance on observable quantities. One unexpected conclusion we are forced into is that cosmic variance is unusually large at z>6, with both our simulations and, most likely, observational measurements are still not fully converged for even such basic quantities as the average Gunn-Peterson optical depth or the volume-weighted neutral fraction. We also find that reionization has little effect on the early galaxies or on global cosmic star formation history, because galaxies whose gas content is affected by photoionization contain no molecular (i.e. star-forming) gas in the first place. In particular, measurements of the faint end of the galaxy luminosity function by JWST are unlikely to provide a useful constraint on reionization.


title: The Spatial Structure of Young Stellar Clusters. I. Subclusters

abstract: The clusters of young stars in massive star-forming regions show a wide range of sizes, morphologies, and numbers of stars. Their highly subclustered structures are revealed by the MYStIX project's sample of 31,754 young stars in nearby sites of star formation (regions at distances <3.6 kpc that contain at least one O-type star.) In 17 of the regions surveyed by MYStIX, we identify subclusters of young stars using finite mixture models -- collections of isothermal ellipsoids that model individual subclusters. Maximum likelihood estimation is used to estimate the model parameters, and the Akaike Information Criterion is used to determine the number of subclusters. This procedure often successfully finds famous subclusters, such as the BN/KL complex behind the Orion Nebula Cluster and the KW-object complex in M17. A catalog of 142 subclusters is presented, with 1 to 20 subclusters per region. The subcluster core radius distribution for this sample is peaked at 0.17 pc with a standard deviation of 0.43 dex, and subcluster core radius is negatively correlated with gas/dust absorption of the stars -- a possible age effect. Based on the morphological arrangements of subclusters, we identify four classes of spatial structure: long chains of subclusters, clumpy structures, isolated clusters with a core-halo structure, and isolated clusters well fit by a single isothermal ellipsoid.


title: Long Term Wind-Driven X-Ray Spectral Variability of NGC 1365 with Swift

abstract: We present long-term (months-years) X-ray spectral variability of the Seyfert 1.8 galaxy NGC 1365 as observed by Swift, which provides well sampled observations over a much longer timescale (6 years) and a much larger flux range than is afforded by other observatories. At very low luminosities the spectrum is very soft, becoming rapidly harder as the luminosity increases and then, above a particular luminosity, softening again. At a given flux level, the scatter in hardness ratio is not very large, meaning that the spectral shape is largely determined by the luminosity. The spectra were therefore summed in luminosity bins and fitted with a variety of models. The best fitting model consists of two power laws, one unabsorbed and another, more luminous, which is absorbed. In this model, we find a range of intrinsic 0.5-10.0 keV luminosities of approximately 1.1-3.5 ergs/s, and a very large range of absorbing columns, of approximately 10^22 - 10^24 cm^-2. Interestingly, we find that the absorbing column decreases with increasing luminosity, but that this result is not due to changes in ionisation. We suggest that these observations might be interpreted in terms of a wind model in which the launch radius varies as a function of ionising flux and disc temperature and therefore moves out with increasing accretion rate, i.e. increasing X-ray luminosity. Thus, depending on the inclination angle of the disc relative to the observer, the absorbing column may decrease as the accretion rate goes up. The weaker, unabsorbed, component may be a scattered component from the wind.


title: Uncovering the Intrinsic Variability of Gamma-ray Bursts

abstract: We develop a robust technique to determine the minimum variability timescale for Gamma-ray Burst light curves, utilizing Haar wavelets. Our approach averages over the data for a given GRB, providing an aggregate measure of signal variation while also retaining sensitivity to narrow pulses within complicated time-series. In contrast to previous studies using wavelets, which simply define the minimum timescale in reference to the measurement noise floor, our approach identifies the signature of temporally-smooth features in the wavelet scaleogram and then additionally identifies a break in the scaleogram on longer timescales as signature of a true, temporally-unsmooth light curve feature or features. We apply our technique to the large sample of Swift GRB Gamma-ray light curves and for the first time -- due to the presence of a large number of GRBs with measured redshift -- determine the distribution of minimum variability timescales in the source frame. We find a median minimum timescale for long-duration GRBs in the source frame of $\Delta t_{\rm min}=0.5$ s, with the shortest timescale found being on the order of 10 ms. This short timescale suggests a compact central engine ($3 \times 10^3$ km). We discuss further implications for the GRB fireball model and present a tantalizing correlation between minimum timescale and redshift, which may in part be due to cosmological time-dilation.


title: The highest-frequency detection of a radio relic: 16-GHz AMI observations of the `Sausage' cluster

abstract: We observed the cluster CIZA J2242.8+5301 with the Arcminute Microkelvin Imager at $16$ GHz and present the first high radio-frequency detection of diffuse, non-thermal cluster emission. This cluster hosts a variety of bright, extended, steep-spectrum synchrotron-emitting radio sources, associated with the intra-cluster medium, called radio relics. Most notably, the northern, Mpc-wide, narrow relic provides strong evidence for diffusive shock acceleration in clusters. We detect a puzzling, flat-spectrum, diffuse extension of the southern relic, which is not visible in the lower radio-frequency maps. The northern radio relic is unequivocally detected and measures an integrated flux of $1.2\pm0.3$ mJy. While the low-frequency ($<2$ GHz) spectrum of the northern relic is well represented by a power-law, it clearly steepens towards $16$ GHz. This result is inconsistent with diffusive shock acceleration predictions of ageing plasma behind a uniform shock front. The steepening could be caused by an inhomogeneous medium with temperature/density gradients or by lower acceleration efficiencies of high energy electrons. Further modelling is necessary to explain the observed spectrum.


title: Mobility versus quality in 2D semiconductor structures

abstract: We consider theoretically effects of random charged impurity disorder on the {\it quality} of high-mobility two dimensional (2D) semiconductor structures, explicitly demonstrating that the sample mobility is not necessarily a reliable or universal indicator of the sample quality in high-mobility modulation-doped 2D GaAs structures because, depending on the specific system property of interest, mobility and quality may be controlled by different aspects of the underlying disorder distribution, particularly since these systems are dominated by long-range Coulomb disorder from both near and far random quenched charged impurities. We show that in the presence of both channel and remote charged impurity scattering, which is a generic situation in modulation-doped high-mobility 2D carrier systems, it is quite possible for higher (lower) mobility structures to have lower (higher) quality as measured by the disorder-induced single-particle level broadening. In particular, we establish that there is no reason to expect a unique relationship between mobility and quality in 2D semiconductor structures as both are independent functionals of the disorder distribution, and are therefore, in principle, independent of each other. Using a simple, but reasonably realistic, ``2-impurity'' minimal model of the disorder distribution, we provide concrete examples of situations where higher (lower) mobilities correspond to lower (higher) sample qualities. We discuss experimental implications of our theoretical results and comment on possible strategies for future improvement of 2D sample quality.


title: CPV Phenomenology of Flavor Conserving Two Higgs Doublet Models

abstract: We analyze the constraints on CP-violating, flavor conserving Two Higgs Doublet Models (2HDMs) implied by measurements of Higgs boson properties at the Large Hadron Collider (LHC) and by the non-observation of permanent electric dipole moments (EDMs) of molecules, atoms and the neutron. We find that the LHC and EDM constraints are largely complementary, with the LHC studies constraining the mixing between the neutral CP-even states and EDMs probing the effect of mixing between the CP-even and CP-odd scalars. The presently most stringent constraints are implied by the non-observation of the ThO molecule EDM signal. Future improvements in the sensitivity of neutron and diamagnetic atom EDM searches could yield competitive or even more severe constraints. We analyze the quantitative impact of hadronic and nuclear theory uncertainties on the interpretation of the latter systems and conclude that these uncertainties cloud the impact of projected improvements in the corresponding experimental sensitivities.


title: A continuous Mott transition between a metal and a quantum spin liquid

abstract: More than half a century after first being proposed by Sir Nevill Mott, the deceptively simple question of whether the interaction-driven electronic metal-insulator transition may be continuous remains enigmatic. Recent experiments on two-dimensional materials suggest that when the insulator is a quantum spin liquid, lack of magnetic long-range order on the insulating side may cause the transition to be continuous, or only very weakly first order. Motivated by this, we study a half-filled extended Hubbard model on a triangular lattice strip geometry. We argue, through use of large-scale numerical simulations and analytical bosonization, that this model harbors a continuous (Kosterlitz-Thouless-like) quantum phase transition between a metal and a gapless spin liquid characterized by a spinon Fermi surface, i.e., a "spinon metal." These results may provide a rare insight into the development of Mott criticality in strongly interacting two-dimensional materials and represent one of the first numerical demonstrations of a Mott insulating quantum spin liquid phase in a genuinely electronic microscopic model.


title: A z=2.5 protocluster associated with the radio galaxy MRC 2104-242: star formation and differing mass functions in dense environments

abstract: We present results from a narrow-band survey of the field around the high redshift radio galaxy MRC 2104-242. We have selected Halpha emitters in a 7sq.arcmin field and compared the measured number density with that of a field sample at similar redshift. We find that MRC 2104-242 lies in an overdensity of galaxies that is 8.0 +/- 0.8 times the average density of a blank field, suggesting it resides in a large-scale structure that may eventually collapse to form a massive cluster. We find that there is more dust obscured star formation in the protocluster galaxies than in similarly selected control field galaxies and there is tentative evidence of a higher fraction of starbursting galaxies in the denser environment. However, on average we do not find a difference between the star formation rate (SFR)-mass relations of the protocluster and field galaxies and so conclude that the SFR of these galaxies at z~2.5 is governed predominantly by galaxy mass and not the host environment. We also find that the stellar mass distribution of the protocluster galaxies is skewed towards higher masses and there is a significant lack of galaxies at M < 10^10Msun within our small field of view. Based on the level of overdensity we expect to find ~22 star forming galaxies below 10^10Msun in the protocluster and do not detect any. This lack of low mass galaxies affects the level of overdensity which we detect. If we only consider high mass (M > 10^10.5Msun) galaxies, the density of the protocluster field increases to ~55 times the control field density.


title: Waving transport and propulsion in a generalized Newtonian fluid

abstract: Cilia and flagella are hair-like appendages that protrude from the surface of a variety of eukaryotic cells and deform in a wavelike fashion to transport fluids and propel cells. Motivated by the ubiquity of non-Newtonian fluids in biology, we address mathematically the role of shear-dependent viscosities on both the waving flagellar locomotion and ciliary transport by metachronal waves. Using a two-dimensional waving sheet as model for the kinematics of a flagellum or an array of cilia, and allowing for both normal and tangential deformation of the sheet, we calculate the flow field induced by a small-amplitude deformation of the sheet in a generalized Newtonian Carreau fluid up to order four in the dimensionless waving amplitude. Shear-thinning and shear-thickening fluids are seen to always induce opposite effects. When the fluid is shear-thinning, the rate of working of the sheet against the fluid is always smaller than in the Newtonian fluid, and the largest gain is obtained for antiplectic metachronal waves. Considering a variety of deformation kinematics for the sheet, we further show that in all cases transport by the sheet is more efficiency in a shear-thinning fluid, and in most cases the transport speed in the fluid is also increased. Comparing the order of magnitude of the shear-thinning contributions with past work on elastic effects as well as the magnitude of the Newtonian contributions, our theoretical results suggest that the impact of shear-dependent viscosities on transport could play a major biological role.


title: Locality and universality of quantum memory effects

abstract: Recently, a series of different measures quantifying memory effects in the quantum dynamics of open systems has been proposed. Here, we derive a mathematical representation for the non-Markovianity measure based on the exchange of information between the open system and its environment which substantially simplifies its numerical and experimental determination, and fully reveals the locality and universality of non-Markovianity in the quantum state space. We further illustrate the application of this representation by means of an all-optical experiment which allows the measurement of the degree of memory effects in a photonic quantum process with high accuracy.


title: Electroweak Symmetry Breaking via QCD

abstract: We propose a new mechanism to generate the electroweak scale within the framework of QCD, which is extended to include conformally invariant scalar degrees of freedom belonging to a larger irreducible representation of $SU(3)_c$. The electroweak symmetry breaking is triggered dynamically via the Higgs portal by the condensation of the colored scalar field around 1 TeV. The mass of the colored boson is restricted to be 350 GeV $\lesssim m_S\lesssim $ 3 TeV, with the upper bound obtained from perturbative renormalization group evolution. This implies that the colored boson can be produced at LHC. If the colored boson is electrically charged, the branching fraction of the Higgs decaying into two photons can slightly increase, and moreover, it can be produced at future linear colliders. Our idea of non-perturbative EW scale generation can serve as a new starting point for more realistic model building in solving the hierarchy problem.


title: Specific Features of Reflection of Radiation of a Compact Source Rotating about its Axis

abstract: Reprocessing of X-ray pulsar radiation by the atmosphere of a companion in a binary system may result in reflected pulse radiation with a period of the pulses equal to the period of the pulsar under suitable conditions. In this paper the influence of the rotation of the source - the pulsar - about its axis on the parameters of thus reflected pulses is investigated. The binary system is modeled by the spherical reflective screen and the compact source uniformly rotating about its axis; the beam pattern (BP) of the source periodically runs over the surface of the screen. Irradiation of the screen by the pulses infinitely narrow in time and by the rectangular pulses is considered. In this model parameters of the pulses reflected in some directions are calculated. The main conclusion provided by the consideration of this model is that the properties of reflected pulses - their profile, and the moments of reaching the observer - substantially depend on the correlation between the light speed and the speed of the BP passing over the companion surface. The possibility of applying of the obtained results to the known X-ray accretion-powered pulsars and rotation-powered pulsars in binary systems is examined.


title: Identifying recombination hotspots using population genetic data

abstract: Motivation: Recombination rates vary considerably at the fine scale within mammalian genomes, with the majority of recombination occurring within hotspots of ~2 kb in width. We present a method for inferring the location of recombination hotspots from patterns of linkage disequilibrium within samples of population genetic data. Results: Using simulations, we show that our method has hotspot detection power of approximately 50-60%, but depending on the magnitude of the hotspot. The false positive rate is between 0.24 and 0.56 false positives per Mb for data typical of humans. Availability: http://github.com/auton1/LDhot


title: Single-parameter spin-pumping in driven metallic rings with spin-orbit coupling

abstract: We consider the generation of a pure spin-current at zero bias voltage with a single time-dependent potential. To such end we study a device made of a mesoscopic ring connected to electrodes and clarify the interplay between a magnetic flux, spin-orbit coupling and non-adiabatic driving in the production of a spin and electrical current. By using Floquet theory, we show that the generated spin to charge current ratio can be controlled by tuning the spin-orbit coupling.


title: Large tunable image-charge effects in single-molecule junctions

abstract: The characteristics of molecular electronic devices are critically determined by metal-organic interfaces, which influence the arrangement of the orbital levels that participate in charge transport. Studies on self-assembled monolayers (SAMs) show (molecule-dependent) level shifts as well as transport-gap renormalization, suggesting that polarization effects in the metal substrate play a key role in the level alignment with respect to the metal's Fermi energy. Here, we provide direct evidence for an electrode-induced gap renormalization in single-molecule junctions. We study charge transport in single porphyrin-type molecules using electrically gateable break junctions. In this set-up, the position of the occupied and unoccupied levels can be followed in situ and with simultaneous mechanical control. When increasing the electrode separation, we observe a substantial increase in the transport gap with level shifts as high as several hundreds of meV for displacements of a few \aa ngstroms. Analysis of this large and tunable gap renormalization with image-charge calculations based on atomic charges obtained from density functional theory confirms and clarifies the dominant role of image-charge effects in single-molecule junctions.


title: Balancing Sparsity and Rank Constraints in Quadratic Basis Pursuit

abstract: We investigate the methods that simultaneously enforce sparsity and low-rank structure in a matrix as often employed for sparse phase retrieval problems or phase calibration problems in compressive sensing. We propose a new approach for analyzing the trade off between the sparsity and low rank constraints in these approaches which not only helps to provide guidelines to adjust the weights between the aforementioned constraints, but also enables new simulation strategies for evaluating performance. We then provide simulation results for phase retrieval and phase calibration cases both to demonstrate the consistency of the proposed method with other approaches and to evaluate the change of performance with different weights for the sparsity and low rank structure constraints.


title: Constraining chameleon models with cosmology

abstract: Chameleon fields may modify gravity on cluster scales while recovering general relativity locally. This article reviews signatures of chameleon modifications in the nonlinear cosmological structure, comparing different techniques to model them, summarising the current state of observational constraints, and concluding with an outlook on prospective constraints from future observations and applications of the analytic tools developed in the process to more general scalar-tensor theories. Particular focus is given to the Hu-Sawicki and designer models of f(R) gravity.


title: Photoelastic coupling in gallium arsenide optomechanical disk resonators

abstract: We analyze the magnitude of the radiation pressure and electrostrictive stresses exerted by light confined inside GaAs semiconductor WGM optomechanical disk resonators, through analytical and numerical means, and find the electrostrictive force to be of prime importance. We investigate the geometric and photoelastic optomechanical coupling resulting respectively from the deformation of the disk boundary and from the strain-induced refractive index changes in the material, for various mechanical modes of the disks. Photoelastic optomechanical coupling is shown to be a predominant coupling mechanism for certain disk dimensions and mechanical modes, leading to total coupling g$_{om}$ and g$_0$ reaching respectively 3 THz/nm and 4 MHz. Finally, we point towards ways to maximize the photoelastic coupling in GaAs disk resonators, and we provide some upper bounds for its value in various geometries.


title: Exciton Binding Energy and Nonhydrogenic Rydberg Series in Monolayer WS2

abstract: We have experimentally determined the energies of the ground and first four excited excitonic states of the fundamental optical transition in monolayer WS2, a model system for the growing class of atomically thin two-dimensional semiconductor crystals. From the spectra, we establish a large exciton binding energy of 0.32 eV and a pronounced deviation from the usual hydrogenic Rydberg series of energy levels of the excitonic states. We explain both of these results using a microscopic theory in which the non-local nature of the effective dielectric screening modifies the functional form of the Coulomb interaction. These strong but unconventional electron-hole interactions are expected to be ubiquitous in atomically thin materials.


title: Electronic correlations stabilizing time-reversal broken chiral superconductivity in single-trilayer TiSe$_2$

abstract: Bulk TiSe$_2$ is an intrinsically layered transition metal dichalcogenide (TMD) hosting both superconducting and charge density wave (CDW) ordering. Motivated by the recent progress in preparing two-dimensional TMDs, we study these frustrated orderings in {\it single} trilayer of TiSe$_2$ within a renormalization group approach. We establish that a novel state with time-reversal symmetry broken chiral superconductivity can emerge from the strong competition between CDW formation and superconductivity. Its stability depends on the precise strength and screening of the electron-electron interactions in two-dimensional TiSe$_2$.


title: The Herschel Exploitation of Local Galaxy Andromeda (HELGA): IV. Dust scaling relations at sub-kpc resolution

abstract: The imprints of dust-starlight interactions are visible in scaling relations between stellar characteristics, star formation parameters and dust properties. We aim to examine dust scaling relations on a sub-kpc resolution in the Andromeda galaxy (M31) by comparing the properties on a local and global scale to other galaxies of the local universe. New Herschel observations are combined with available data from GALEX, SDSS, WISE and Spitzer to construct a dataset covering UV to submm wavelengths. We work at the resolution of the SPIRE $500\; \mu$m beam, with pixels corresponding to physical regions of 137 x 608 pc in the galaxy's disk. A panchromatic spectral energy distribution was modelled for each pixel and several dust scaling relations are investigated. We find, on a sub-kpc scale, strong correlations between $M_d/M_\star$ and NUV-r, and between $M_d/M_\star$ and $\mu_\star$ (the stellar mass surface density). Striking similarities with corresponding relations based on integrated galaxies are found. We decompose M31 in four macro-regions based on their FIR morphology; the bulge, inner disk, star forming ring and the outer disk. All regions closely follow the galaxy-scale average trends. The specific star formation characteristics we derive for these macro-regions give strong hints of an inside-out formation of the bulge-disk geometry, as well as an internal downsizing process. However, within each macro-region, a great diversity in individual micro-regions is found. Furthermore, we confirm that dust in the bulge of M31 is heated only by the old stellar populations. In general, the local dust scaling relations indicate that the dust content in M31 is maintained by a subtle interplay of past and present star formation. The similarity with galaxy-based relations strongly suggests that they are in situ correlations, with underlying processes that must be local in nature. (Abriged)


title: Asymmetries in Silicon Microstrip Response Function and Lorentz Angle

abstract: An experimental set up, dedicated to isolate an error present in the $\eta$-algorithm, gave an unexpected result. The average of a center of gravity algorithm at orthogonal particle incidence turns out to be non zero. This non zero average signals an asymmetry in the response function of the strips, and introduces a further parameter in the corrections: the shift of the strip response center of gravity respect its geometrical position. A strategy to extract this parameter from a standard data set is discussed. Some simulations with various asymmetric response functions are explored for this test. The method is able to detect easily the asymmetry parameters introduced in the simulations. Its robustness is tested against angular rotations, and we see an almost linear variation with the angle. This simple property is used to simulate a determination of a Lorentz angle with and without the asymmetry of the response function.


title: Numerical integrators for motion under a strong constraining force

abstract: This paper deals with the numerical integration of Hamiltonian systems in which a stiff anharmonic potential causes highly oscillatory solution behavior with solution-dependent frequencies. The impulse method, which uses micro- and macro-steps for the integration of fast and slow parts, respectively, does not work satisfactorily on such problems. Here it is shown that variants of the impulse method with suitable projection preserve the actions as adiabatic invariants and yield accurate approximations, with macro-stepsizes that are not restricted by the stiffness parameter.


title: Deforming solutions of geometric variational problems with varying symmetry groups

abstract: We prove an equivariant implicit function theorem for variational problems that are invariant under a varying symmetry group (corresponding to a bundle of Lie groups). Motivated by applications to families of geometric variational problems lacking regularity, several non-smooth extensions of the result are discussed. Among such applications is the submanifold problem of deforming the ambient metric preserving a given variational property of a prescribed family of submanifolds, e.g., constant mean curvature, up to the action of the corresponding ambient isometry groups.


title: M-theoretic matrix models

abstract: Some matrix models admit, on top of the usual 't Hooft expansion, an M-theory-like expansion, i.e. an expansion at large N but where the rest of the parameters are fixed, instead of scaling with N. These models, which we call M-theoretic matrix models, appear in the localization of Chern-Simons-matter theories, and also in two-dimensional statistical physics. Generically, their partition function receives non-perturbative corrections which are not captured by the 't Hooft expansion. In this paper, we discuss general aspects of these type of matrix integrals and we analyze in detail two different examples. The first one is the matrix model computing the partition function of N=4 supersymmetric Yang-Mills theory in three dimensions with one adjoint hypermultiplet and N_f fundamentals, which has a conjectured M-theory dual, and which we call the N_f matrix model. The second one, which we call the polymer matrix model, computes form factors of the 2d Ising model and is related to the physics of 2d polymers. In both cases we determine their exact planar limit. In the N_f matrix model, the planar free energy reproduces the expected behavior of the M-theory dual. We also study their M-theory expansion by using Fermi gas techniques, and we find non-perturbative corrections to the 't Hooft expansion.


title: Competition between monomeric and dimeric crystals in schematic models for globular proteins

abstract: Advances in experimental techniques and in theoretical models have improved our understanding of protein crystallization. But they have also left open questions regarding the protein phase behavior and self-assembly kinetics, such as why (nearly) identical crystallization conditions sometimes result in the formation of different crystal forms. Here, we develop a patchy particle model with competing sets of patches that provides a microscopic explanation of this phenomenon. We identify different regimes in which one or two crystal forms can coexist with a low-density fluid. Using analytical approximations, we extend our findings to different crystal phases, providing a general framework for treating protein crystallization when multiple crystal forms compete. Our results also suggest different experimental routes for targeting a specific crystal form, and for reducing the dynamical competition between the two forms, thus facilitating protein crystal assembly.


title: Multilevel methods for nonuniformly elliptic operators

abstract: We develop and analyze multilevel methods for nonuniformly elliptic operators whose ellipticity holds in a weighted Sobolev space with an $A_2$--Muckenhoupt weight. Using the so-called Xu-Zikatanov (XZ) identity, we derive a nearly uniform convergence result, under the assumption that the underlying mesh is quasi-uniform. We also consider the so-called $\alpha$-harmonic extension to localize fractional powers of elliptic operators. Motivated by the scheme proposed in [R.H. Nochetto, E. Otarola, and A.J. Salgado. A PDE approach to fractional diffusion in general domains: a priori error analysis. arXiv:1302.0698, 2013] we present a multilevel method with line smoothers and obtain a nearly uniform convergence result on anisotropic meshes. Numerical experiments reveal a competitive performance of our method.


title: Effective Field Theory of Relativistic Quantum Hall Systems

abstract: Motivated by the observation of the fractional quantum Hall effect in graphene, we consider the effective field theory of relativistic quantum Hall states. We find that, beside the Chern-Simons term, the effective action also contains a term of topological nature, which couples the electromagnetic field with a topologically conserved current of $2+1$ dimensional relativistic fluid. In contrast to the Chern-Simons term, the new term involves the spacetime metric in a nontrivial way. We extract the predictions of the effective theory for linear electromagnetic and gravitational responses. For fractional quantum Hall states at the zeroth Landau level, additional holomorphic constraints allow one to express the results in terms of two dimensionless constants of topological nature.


title: Quantum bidding in Bridge

abstract: Quantum methods allow to reduce communication complexity of some computational tasks, with several separated partners, beyond classical constraints. Nevertheless, experimental demonstrations of this fact are thus far limited to some abstract problems, far away from real-life tasks. We show here, and demonstrate experimentally, that the power of reduction of communication complexity can be harnessed to gain advantage in famous, immensely popular, card game - Bridge. The essence of a winning strategy in Bridge is efficient communication between the partners. The rules of the game allow only specific form of communication, of a very low complexity (effectively one has a strong limitations on number of exchanged bits). Surprisingly, our quantum technique is not violating the existing rules of the game (as there is no increase in information flow). We show that our quantum Bridge auction corresponds to a biased nonlocal Clauser-Horne-Shimony-Holt (CHSH) game, which is equivalent to a $2\to 1$ quantum random access code. Thus our experiment is also a realization of such protocols. However, this correspondence is not full which enables the Bridge players to have efficient strategies regardless of the quality of their detectors.


title: A small, infinitely-ended 2-knot group

abstract: We show that a 2-knot group discovered in the course of a census of 4-manifolds with small triangulations is an HNN extension with finite base and proper associated subgroups, and has the smallest base among such knot groups.


title: Global phase and frequency comb structures in nonlinear Compton and Thomson scattering

abstract: The Compton and Thomson radiation spectra, generated in collisions of an electron beam with a powerful laser beam, are studied in the framework of quantum and classical electrodynamics, respectively. We show that there are frequency regimes where both radiation spectra are nearly identical, which for Compton scattering relates to the process which preserves the electron spin. Although the radiation spectra are nearly identical, the corresponding probability amplitudes exhibit different global phases. This has pronounced consequences, which we demonstrate by investigating temporal power distributions in both cases. We show that, contrary to Thomson scattering, it is not always possible to synthesize short laser pulses from Compton radiation. This happens when the global phase of the Compton amplitude varies in a nonlinear way with the frequency of emitted photons. We also demonstrate that while the Compton process driven by a non-chirped laser pulse can generate chirped bursts of radiation, this is not the case for the Thomson process. In principle, both processes can lead to a generation of coherent frequency combs when single or multiple driving laser pulses collide with electrons. Once we synthesize these combs into short bursts of radiation, we can control them, for instance, by changing the time delay between the driving pulses.

categories: hep-ph physics.acc-ph physics.atom-ph physics.optics physics.plasm-ph


title: A poset view of the major index

abstract: We introduce the Major MacMahon map and show how this map interacts with the pyramid and bipyramid operators. When the Major MacMahon map is applied to the ab-index of a simplicial poset, it yields the q-analogue of n! times the h-polynomial of the poset. Applying the map to the Boolean algebra gives the distribution of the major index on the symmetric group, a seminal result due to MacMahon. Similarly, when applied to the cross-polytope we obtain the distribution of one of the major indexes on signed permutations due to Reiner.


title: Anomalous magnetic moment of the positronium ion

abstract: We determine the gyromagnetic factor of the positronium ion, a three-body system consisting of two electrons and a positron, including first relativistic corrections. We find that the g-factor is modified by a term -0.51(1)alpha^2, exceeding 15 times the alpha^2 correction for a free electron. We compare this effect with analogous results found previously in atomic positronium and in hydrogen-like ions.


title: Loop measures and the Gaussian free field

abstract: Loop measures and their associated loop soups are generally viewed as arising from finite state Markov chains. We generalize several results to loop measures arising from potentially complex edge weights. We discuss two applications: Wilson's algorithm to produce uniform spanning trees and an isomorphism theorem due to Le Jan.


title: Absolutely $k$-convex domains and holomorphic foliations on homogeneous manifolds

abstract: We consider a holomorphic foliation $\mathcal{F}$ of codimension $k\geq 1$ on a homogeneous compact K\"ahler manifold $X$ of dimension $n>k$. Assuming that the singular set $Sing(\mathcal{F})$ of $\mathcal{F}$ is contained in an absolutely $k$-convex domain $U\subset X$, we prove that the determinant of normal bundle $\det(N_{\mathcal{F}})$ of $\mathcal{F}$ cannot be an ample line bundle, provided $[n/k]\geq 2k+3$. Here $[n/k]$ denotes the largest integer $\leq n/k.$


title: Factorization of Graded Traces on Nichols Algebras

abstract: We study the factorization of the Hilbert series and more general graded traces of a Nichols algebra into cyclotomic polynomials. Nichols algebras play the role of Borel subalgebras of finite-dimensional quantum groups, so this observation can be viewed as an analog to the factorization of the order of finite Lie groups into cyclotomic polynomials. We prove results on this factorization and give a table of many examples of rank 1 over nonabelian groups.


title: On characterizing nonlocality and anisotropy for the magnetorotational instability

abstract: The extent to which angular momentum transport in accretion discs is primarily local or non-local and what determines this is an important avenue of study for understanding accretion engines. Taking a step along this path, we analyze simulations of the magnetorotational instability (MRI) by calculating energy and stress power spectra in stratified isothermal shearing box simulations in several new ways. We divide our boxes in two regions, disc and corona where the disc is the MRI unstable region and corona is the magnetically dominated region. We calculate the fractional power in different quantities, including magnetic energy and Maxwell stresses and find that they are dominated by contributions from the lowest wave numbers. This is even more dramatic for the corona than the disc, suggesting that transport in the corona region is dominated by larger structures than the disc. By calculating averaged power spectra in one direction of $k$ space at a time, we also show that the MRI turbulence is strongly anisotropic on large scales when analyzed by this method, but isotropic on small scales. Although the shearing box itself is meant to represent a local section of an accretion disc, the fact that the stress and energy are dominated by the largest scales highlights that the locality is not captured within the box. This helps to quantify the intuitive importance of global simulations for addressing the question of locality of transport, for which similar analyses can be performed.

categories: astro-ph.SR astro-ph.GA astro-ph.HE physics.plasm-ph


title: Telling Breaking News Stories from Wikipedia with Social Multimedia: A Case Study of the 2014 Winter Olympics

abstract: With the ability to watch Wikipedia and Wikidata edits in realtime, the online encyclopedia and the knowledge base have become increasingly used targets of research for the detection of breaking news events. In this paper, we present a case study of the 2014 Winter Olympics, where we tell the story of breaking news events in the context of the Olympics with the help of social multimedia stemming from multiple social network sites. Therefore, we have extended the application Wikipedia Live Monitor-a tool for the detection of breaking news events-with the capability of automatically creating media galleries that illustrate events. Athletes winning an Olympic competition, a new country leading the medal table, or simply the Olympics themselves are all events newsworthy enough for people to concurrently edit Wikipedia and Wikidata-around the world in many languages. The Olympics being an event of common interest, an even bigger majority of people share the event in a multitude of languages on global social network sites, which makes the event an ideal subject of study. With this work, we connect the world of Wikipedia and Wikidata with the world of social network sites, in order to convey the spirit of the 2014 Winter Olympics, to tell the story of victory and defeat, and always following the Olympic motto Citius, Altius, Fortius. The proposed system-generalized for all sort of breaking news stories-has been put in production in form of the Twitter bot @mediagalleries, available and archived at https://twitter.com/mediagalleries.


title: Data-Driven Model Reduction for the Bayesian Solution of Inverse Problems

abstract: One of the major challenges in the Bayesian solution of inverse problems governed by partial differential equations (PDEs) is the computational cost of repeatedly evaluating numerical PDE models, as required by Markov chain Monte Carlo (MCMC) methods for posterior sampling. This paper proposes a data-driven projection-based model reduction technique to reduce this computational cost. The proposed technique has two distinctive features. First, the model reduction strategy is tailored to inverse problems: the snapshots used to construct the reduced-order model are computed adaptively from the posterior distribution. Posterior exploration and model reduction are thus pursued simultaneously. Second, to avoid repeated evaluations of the full-scale numerical model as in a standard MCMC method, we couple the full-scale model and the reduced-order model together in the MCMC algorithm. This maintains accurate inference while reducing its overall computational cost. In numerical experiments considering steady-state flow in a porous medium, the data-driven reduced-order model achieves better accuracy than a reduced-order model constructed using the classical approach. It also improves posterior sampling efficiency by several orders of magnitude compared to a standard MCMC method.


title: An importance sampling approach for copula models in insurance

abstract: An importance sampling approach for sampling copula models is introduced. We propose two algorithms that improve Monte Carlo estimators when the functional of interest depends mainly on the behaviour of the underlying random vector when at least one of the components is large. Such problems often arise from dependence models in finance and insurance. The importance sampling framework we propose is general and can be easily implemented for all classes of copula models from which sampling is feasible. We show how the proposal distribution of the two algorithms can be optimized to reduce the sampling error. In a case study inspired by a typical multivariate insurance application, we obtain variance reduction factors between 10 and 30 in comparison to standard Monte Carlo estimators.


title: Decay Oscillations in Electron Capture and the Neutrino Mass Difference

abstract: Quantum mechanical theory disallows the model that has been used to infer the neutrino mass difference from the reported "GSI oscillations" in the rates of decay of hydrogen-like ions by electron capture. It has not been proved that the existence of mass-difference-dependent oscillations conflicts with quantum mechanics but no consistent quantum mechanical model has been shown to predict them.


title: On a condition number of random polynomial systems

abstract: Condition numbers of random polynomial systems have been widely studied in the literature under certain coefficient ensembles of invariant type. In this note we introduce a method that allows us to study these numbers for a broad family of probability distributions. Our work also extends to perturbed systems.