title: A Universal Machine for Biform Theory Graphs

abstract: Broadly speaking, there are two kinds of semantics-aware assistant systems for mathematics: proof assistants express the semantic in logic and emphasize deduction, and computer algebra systems express the semantics in programming languages and emphasize computation. Combining the complementary strengths of both approaches while mending their complementary weaknesses has been an important goal of the mechanized mathematics community for some time. We pick up on the idea of biform theories and interpret it in the MMTt/OMDoc framework which introduced the foundations-as-theories approach, and can thus represent both logics and programming languages as theories. This yields a formal, modular framework of biform theory graphs which mixes specifications and implementations sharing the module system and typing information. We present automated knowledge management work flows that interface to existing specification/programming tools and enable an OpenMath Machine, that operationalizes biform theories, evaluating expressions by exhaustively applying the implementations of the respective operators. We evaluate the new biform framework by adding implementations to the OpenMath standard content dictionaries.


title: The MMT API: A Generic MKM System

abstract: The MMT language has been developed as a scalable representation and interchange language for formal mathematical knowledge. It permits natural representations of the syntax and semantics of virtually all declarative languages while making MMT-based MKM services easy to implement. It is foundationally unconstrained and can be instantiated with specific formal languages. The MMT API implements the MMT language along with multiple backends for persistent storage and frontends for machine and user access. Moreover, it implements a wide variety of MMT-based knowledge management services. The API and all services are generic and can be applied to any language represented in MMT. A plugin interface permits injecting syntactic and semantic idiosyncrasies of individual formal languages.


title: Synthesis of unitaries with Clifford+T circuits

abstract: We describe a new method for approximating an arbitrary $n$ qubit unitary with precision $\varepsilon$ using a Clifford and T circuit with $O(4^{n}n(\log(1/\varepsilon)+n))$ gates. The method is based on rounding off a unitary to a unitary over the ring $\mathbb{Z}[i,1/\sqrt{2}]$ and employing exact synthesis. We also show that any $n$ qubit unitary over the ring $\mathbb{Z}[i,1/\sqrt{2}]$ with entries of the form $(a+b\sqrt{2}+ic+id\sqrt{2})/2^{k}$ can be exactly synthesized using $O(4^{n}nk)$ Clifford and T gates using two ancillary qubits. This new exact synthesis algorithm is an improvement over the best known exact synthesis method by B. Giles and P. Selinger requiring $O(3^{2^{n}}nk)$ elementary gates.


title: Spatiospectral concentration of vector fields on a sphere

abstract: We construct spherical vector bases that are bandlimited and spatially concentrated, or, alternatively, spacelimited and spectrally concentrated, suitable for the analysis and representation of real-valued vector fields on the surface of the unit sphere, as arises in the natural and biomedical sciences, and engineering. Building on the original approach of Slepian, Landau, and Pollak we concentrate the energy of our function bases into arbitrarily shaped regions of interest on the sphere, and within certain bandlimits in the vector spherical-harmonic domain. As with the concentration problem for scalar functions on the sphere, which has been treated in detail elsewhere, a Slepian vector basis can be constructed by solving a finite-dimensional algebraic eigenvalue problem. The eigenvalue problem decouples into separate problems for the radial and tangential components. For regions with advanced symmetry such as polar caps, the spectral concentration kernel matrix is very easily calculated and block-diagonal, lending itself to efficient diagonalization. The number of spatiospectrally well-concentrated vector fields is well estimated by a Shannon number that only depends on the area of the target region and the maximal spherical-harmonic degree or bandwidth. The spherical Slepian vector basis is doubly orthogonal, both over the entire sphere and over the geographic target region. Like its scalar counterparts it should be a powerful tool in the inversion, approximation and extension of bandlimited fields on the sphere: vector fields such as gravity and magnetism in the earth and planetary sciences, or electromagnetic fields in optics, antenna theory and medical imaging.


title: Nuclear superfluidity in isospin asymmetric matter within the Skyrme model

abstract: The phase diagram of the superfluid phase coupled to spin singlet (S=0) and isospin triplet (T=1) states in infinite nuclear matter is analyzed within the nonrelativistic Skyrme model. We use an approach that allows a unified and consistent treatment of the particle-hole and particle-particle channels. The gap equation is solved for the full range of accessible densities, isospin asymmetries, and temperatures. The characteristic features of each of the components Tz=0, +1, -1 are emphasized. The search for signatures of a BCS-BEC crossover yields negative results within the scheme proposed.


title: Bregman Alternating Direction Method of Multipliers

abstract: The mirror descent algorithm (MDA) generalizes gradient descent by using a Bregman divergence to replace squared Euclidean distance. In this paper, we similarly generalize the alternating direction method of multipliers (ADMM) to Bregman ADMM (BADMM), which allows the choice of different Bregman divergences to exploit the structure of problems. BADMM provides a unified framework for ADMM and its variants, including generalized ADMM, inexact ADMM and Bethe ADMM. We establish the global convergence and the $O(1/T)$ iteration complexity for BADMM. In some cases, BADMM can be faster than ADMM by a factor of $O(n/\log(n))$. In solving the linear program of mass transportation problem, BADMM leads to massive parallelism and can easily run on GPU. BADMM is several times faster than highly optimized commercial software Gurobi.


title: The Heliometer of Rio de Janeiro in Operation - 2010 to 2013

abstract: Out of the three quantities that characterize the state of an isolated gaseous body: pressure, temperature and volume the radius is the only one directly measurable for the Sun, what is specially true in the optical window and for ground base measurements. The Heliometer of Observatorio Nacional, in Rio de Janeiro, measures the distance between two opposite limbs of the Sun in the same field of view, through the reflection on a 10 cm parabolic mirror split on its half and forming an appropriate angle. This configuration is free from optical aberrations and focal variations along the measurement direction. The mirrors are made on CCZ vitro-ceramic and the telescope structure is of carbon steel, resulting that there is no flexion or temperature deformation. The instrument is compact, and can perform hundreds of measurements per duty day, around all heliolatitudes. It attains an accuracy on the solar radius of 0.01 arcsec, becoming the ideal instrument to monitor from the ground the solar diameter, and to bridge satellites and astrolabes historical series of data. We discuss the first years of regular observation, with emphasis on the instrumental calibrations and on the statistic study of the derived time series, attitude series, and solar geometry series. On basis of these series we obtain how well the Heliometer and Solar Astrolabe results are matched.


title: Lifshitz Tails for Anderson Models with Sign-Indefinite Single-Site Potentials

abstract: We study the spectral minimum and Lifshitz tails for continuum random Schr\"{o}dinger operators of the form \begin{equation*} H_{\om}=-\De+V_{0}+\sum_{i\in\Z^{d}}\om_{i}u(\cdot-i), \end{equation*} where $V_{0}$ is the periodic potential, $\{\om_{i}\}_{i\in\Z^{d}}$ are i.i.d random variables and $u$ is the sign-indefinite impurity potential. Recently, this model has been proven to exhibit Lifshitz tails near the bottom of the spectrum under the small support assuption of $u$ and the reflection symmetry assumption of $V_{0}$ and $u$. We here drop the reflection symmetry assumption of $V_{0}$ and $u$. We first give characterizations of the bottom of the spectrum. Then, we show the existence of Lifshitz tails in the regime where the characterization of the bottom of the spectrum is explicit. In particular, this regime covers the reflection symmetry case.


title: A detailed census of variable stars in the globular cluster NGC 6333 (M9) from CCD differential photometry

abstract: We report CCD $V$ and $I$ time-series photometry of the globular cluster NGC 6333 (M9). The technique of difference image analysis has been used, which enables photometric precision better than 0.05 mag for stars brighter than $V \sim 19.0$ mag, even in the crowded central regions of the cluster. The high photometric precision has resulted in the discovery of two new RRc stars, three eclipsing binaries, seven long-term variables and one field RRab star behind the cluster. A detailed identification chart and equatorial coordinates are given for all the variable stars in the field of our images of the cluster. Our data together with literature $V$-data obtained in 1994 and 1995 allowed us to refine considerably the periods for all RR Lyrae stars. The nature of the new variables is discussed. We argue that variable V12 is a cluster member and an Anomalous Cepheid. Secular period variations, double mode pulsations and/or the Blazhko-like modulations in some RRc variables are addressed. Through the light curve Fourier decomposition of 12 RR Lyrae stars we have calculated a mean metallicity of [Fe/H]$_{\rm ZW}$=$-1.70 \pm 0.01{\rm(statistical)} \pm 0.14{\rm(systematic)}$ or [Fe/H]$_{UVES}=-1.67 \pm 0.01{\rm(statistical)} \pm 0.19{\rm(systematic)}$.Absolute magnitudes, radii and masses are also estimated for the RR Lyrae stars. A detailed search for SX Phe stars in the Blue Straggler region was conducted but none were discovered. If SX Phe exist in the cluster then their amplitudes must be smaller than the detection limit of our photometry. The CMD has been corrected for heavy differential reddening using the detailed extinction map of the cluster of Alonso-Garc\'ia et al. (2012). This has allowed us to set the mean cluster distance from two independent estimates; from the RRab and RRc absolute magnitudes, we find $8.04\pm 0.19$ kpc and $7.88\pm0.30$ kpc respectively.


title: Hermite spectral method with hyperbolic cross approximations to high-dimensional parabolic PDEs

abstract: It is well-known that sparse grid algorithm has been widely accepted as an efficient tool to overcome the "curse of dimensionality" in some degree. In this note, we first give the error estimate of hyperbolic cross (HC) approximations with generalized Hermite functions. The exponential convergence in both regular and optimized hyperbolic cross approximations has been shown. Moreover, the error estimate of Hermite spectral method to high-dimensional linear parabolic PDEs with HC approximations has been investigated in the properly weighted Korobov spaces. The numerical result verifies the exponential convergence of this approach.


title: Timescape cosmology with radiation fluid

abstract: The timescape cosmology represents a potentially viable alternative to the standard homogeneous cosmology, without the need for dark energy. Although average cosmic evolution in the timescape scenario only differs substantially from that of Friedmann-Lemaitre model at relatively late epochs when the contribution from the energy density of radiation is negligible, a full solution of the Buchert equations to incorporate radiation is necessary to smoothly match parameters to the epoch of photon decoupling and to obtain constraints from cosmic microwave background data. Here we extend the matter-dominated solution found in earlier work to include radiation, providing series solutions at early times and an efficient numerical integration strategy for generating the complete solution. The numerical solution is used to directly calculate the scale of the sound horizon at decoupling, and at the baryon drag epoch. The constraints on these scales from the Planck satellite data yield bounds on the timescape cosmological parameters, which are found to also agree with the best-fit values from a recent analysis of SDSS-II supernova data, while avoiding the problem of a primordial lithium-7 abundance anomaly.


title: Anisotropic Cosmic Ray Diffusion and its Implications for Gamma-Ray Astronomy

abstract: Analyses of TeV-PeV cosmic ray (CR) diffusion around their sources usually assume either isotropic diffusion or anisotropic diffusion due to the regular Galactic magnetic field. We show that none of them are adequate on distances smaller than the maximal scale Lmax ~ 100 pc of fluctuations in the turbulent interstellar magnetic field. As a result, we predict anisotropic gamma-ray emissions around CR proton and electron sources, even for uniform densities of target gas. The centers of extended emission regions may have non-negligible offsets from their sources, leading to risks of misidentification. Gamma-rays from CR filaments have steeper energy spectra than those from surrounding regions. We point out that gamma-ray telescopes can be used in the future as a new way to probe and deduce the parameters of the interstellar magnetic field.


title: Dark matter halo formation in the multicomponent dark matter models

abstract: This work investigates a set of cosmological collisionless N-body simulations with featured power spectra of initial perturbations in the context of the core-cusp and satellites problems. On the studied power spectra some scales of fluctuations were suppressed. Such spectral features can be caused by multicomponent dark matter. The density profiles innermost resolved slopes $\alpha\equiv d \log(\rho) /d \log(r) $ of the five largest haloes were measured and its dependence on the parameters of the suppression was traced. In a certain range of the parameters the slopes flatten from initial value of about -1.2 to -0.6 or even to -0.2 in one of the cases. This qualitatively demonstrates that (i) profiles shape depends on initial power spectrum and (ii) this effect may be responsible for the solution of the core-cusp problem. The suppression of some part of the initial power spectrum also leads to the decrease of the number of massive subhaloes.


title: Superconformal Generalization of the Chaotic Inflation Model \lambda \phi^4/4 - \xi/ 2 \phi^2 R

abstract: A model of chaotic inflation based on the theory of a scalar field with potential \lambda\phi^4 perfectly matches the observational data if one adds to it a tiny non-minimal coupling to gravity -\xi/2 \phi^2 R with \xi > 0.002. We describe embedding of this model into the superconformal theory with spontaneous breaking of superconformal symmetry, and into supergravity. A model with small \xi is technically natural: setting the small parameter \xi to zero leads to a point of enhanced symmetry in the underlying superconformal theory.


title: Sparse Inverse Covariance Matrix Estimation Using Quadratic Approximation

abstract: The L1-regularized Gaussian maximum likelihood estimator (MLE) has been shown to have strong statistical guarantees in recovering a sparse inverse covariance matrix, or alternatively the underlying graph structure of a Gaussian Markov Random Field, from very limited samples. We propose a novel algorithm for solving the resulting optimization problem which is a regularized log-determinant program. In contrast to recent state-of-the-art methods that largely use first order gradient information, our algorithm is based on Newton's method and employs a quadratic approximation, but with some modifications that leverage the structure of the sparse Gaussian MLE problem. We show that our method is superlinearly convergent, and present experimental results using synthetic and real-world application data that demonstrate the considerable improvements in performance of our method when compared to other state-of-the-art methods.


title: Sets of Complex Unit Vectors with Two Angles and Distance-Regular Graphs

abstract: We study {0,\alpha}-sets, which are sets of unit vectors of $\mathbb{C}^m$ in which any two distinct vectors have angle 0 or \alpha. We investigate some distance-regular graphs that provide new constructions of {0,\alpha}-sets using a method by Godsil and Roy. We prove bounds for the sizes of {0,\alpha}-sets of flat vectors, and characterize all the distance-regular graphs that yield {0,\alpha}-sets meeting the bounds at equality.


title: Superconformal generalizations of the Starobinsky model

abstract: We find a way to represent the Starobinsky model in terms of a simple conformally invariant theory with spontaneous symmetry breaking. We also present a superconformal theory, which, upon spontaneous breaking of the superconformal symmetry, provides a consistent supergravity generalization of the Starobinsky model.


title: Weights for Monoids and Actions of Monoids

abstract: The main objective of the paper is to define the category of monoids as a weighted limit. We also define the category of actions of monoids along the action of a monoidal category as a weighted limit.


title: On the 'Reality' of Observable Properties

abstract: This note contains some initial work on attempting to bring recent developments in the foundations of quantum mechanics concerning the nature of the wavefunction within the scope of more logical and structural methods. A first step involves generalising and reformulating a criterion for the reality of the wavefunction proposed by Harrigan & Spekkens, which was central to the PBR theorem. The resulting criterion has several advantages, including the avoidance of certain technical difficulties relating to sets of measure zero. By considering the 'reality' not of the wavefunction but of the observable properties of any ontological physical theory a novel characterisation of non-locality and contextuality is found. Secondly, a careful analysis of preparation independence, one of the key assumptions of the PBR theorem, leads to an analogy with Bell locality, and thence to a proposal to weaken it to an assumption of `no-preparation-signalling' in analogy with no-signalling. This amounts to introducing non-local correlations in the joint ontic state, which is, at least, consistent with the Bell and Kochen-Specker theorems. The question of whether the PBR result can be strengthened to hold under this relaxed assumption is therefore posed.


title: LHC constraints on dark matter with (130 GeV) gamma ray lines

abstract: Dark matter annihilation into photons in our galaxy would constitute an exciting indirect signal of its existence, as underscored by tentative evidence for 130 GeV dark matter in Fermi/LAT data. Models that give a large annihilation cross section into photons typically require the dark matter to couple to, or be composed of, new charged particles, that can be produced in colliders. We consider the LHC constraints on some representative models of these types, including the signals of same-sign dileptons, opposite-sign dileptons, events mimicking the production and decay of excited leptons, four-photon events, resonant production of composite vectors decaying into two photons, and monophoton events.


title: Black hole variability and the star formation-AGN connection: Do all star-forming galaxies host an AGN?

abstract: We investigate the effect of active galactic nucleus (AGN) variability on the observed connection between star formation and black hole accretion in extragalactic surveys. Recent studies have reported relatively weak correlations between observed AGN luminosities and the properties of AGN hosts, which has been interpreted to imply that there is no direct connection between AGN activity and star formation. However, AGNs may be expected to vary significantly on a wide range of timescales (from hours to Myr) that are far shorter than the typical timescale for star formation (>~100 Myr). This variability can have important consequences for observed correlations. We present a simple model in which all star-forming galaxies host an AGN when averaged over ~100 Myr timescales, with long-term average AGN accretion rates that are perfectly correlated with the star formation rate (SFR). We show that reasonable prescriptions for AGN variability reproduce the observed weak correlations between SFR and L_AGN in typical AGN host galaxies, as well as the general trends in the observed AGN luminosity functions, merger fractions, and measurements of the average AGN luminosity as a function of SFR. These results imply there may be a tight connection between AGN activity and SFR over galaxy evolution timescales, and that the apparent similarities in rest-frame colors, merger rates, and clustering of AGNs compared to "inactive" galaxies may be due primarily to AGN variability. The results provide motivation for future deep, wide extragalactic surveys that can measure the distribution of AGN accretion rates as a function of SFR.


title: Nonlinear structure formation in the Cubic Galileon gravity model

abstract: We model the linear and nonlinear growth of large scale structure in the Cubic Galileon gravity model, by running a suite of N-body cosmological simulations using the {\tt ECOSMOG} code. Our simulations include the Vainshtein screening effect, which reconciles the Cubic Galileon model with local tests of gravity. In the linear regime, the amplitude of the matter power spectrum increases by $\sim 25%$ with respect to the standard $\Lambda$CDM model today. The modified expansion rate accounts for $\sim 20%$ of this enhancement, while the fifth force is responsible for only $\sim 5%$. This is because the effective unscreened gravitational strength deviates from standard gravity only at late times, even though it can be twice as large today. In the nonlinear regime ($k \gtrsim 0.1 h\rm{Mpc}^{-1}$), the fifth force leads to only a modest increase ($\lesssim 8%$) in the clustering power on all scales due to the very efficient operation of the Vainshtein mechanism. Such a strong effect is typically not seen in other models with the same screening mechanism. The screening also results in the fifth force increasing the number density of halos by less than 10%, on all mass scales. Our results show that the screening does not ruin the validity of linear theory on large scales which anticipates very strong constraints from galaxy clustering data. We also show that, whilst the model gives an excellent match to CMB data on small angular scales ($l \gtrsim 50$), the predicted integrated Sachs-Wolf effect is in tension with Planck/WMAP results.


title: Classical 3-loop 2-body diagrams

abstract: As part of the study of the two-body problem in Einstein's gravity, the fourth post-Newtonian order (4PN) of the two-body effective action is being computed presently by both effective field theory (EFT) methods and others. Diagrams with 3 (or 4) classical loops appear to be a significant obstacle. In this paper we develop a method to compute such 3-loop diagrams and demonstrate it through a specific diagram. We reduce the classical diagrams through shrinking the body worldlines to a form more familiar in Quantum Field Theory. A key ingredient in the evaluation is the Integration By Parts method for Feynman integrals.


title: A simple, quantitative method to infer the minimum atmospheric height of small exoplanets

abstract: Amongst the many hundreds of transiting planet candidates discovered by the Kepler Mission, one finds a large number of candidates with sizes between that of the Earth and Neptune. The composition of these worlds is not immediately obvious with no Solar System analog to draw upon and there exists some ambiguity as to whether a given candidate is a rocky Super-Earth or a gas-enveloped Mini-Neptune. The potential scientific value and observability of the atmospheres of these two classes of worlds varies significantly and given the sheer number of candidates in this size-range, there is evidently a need for a quick, simple metric to rank whether the planets have an extended atmosphere or not. In this work, we propose a way to calculate the 'minimum atmospheric height' (R_MAH) using only a planet's radius and mass as inputs. We assume and exploit the boundary condition that the bulk composition of a solid/liquid Super-Earth cannot be composed of a material lighter than that of water. Mass-radius loci above a pure-water composition planet correspond to R_MAH>0. The statistical confidence of a planet maintaining an extended atmosphere can be therefore easily calculated to provide a simple ranking of target planets for follow-up observations. We also discuss how this metric can be useful in the interpretation of the spectra of observed planetary atmospheres.


title: Turbulence In the Outer Regions of Protoplanetary Disks. II. Strong Accretion Driven by a Vertical Magnetic Field

abstract: We carry out a series of local, vertically stratified shearing box simulations of protoplanetary disks that include ambipolar diffusion and a net vertical magnetic field. The ambipolar diffusion profiles we employ correspond to 30AU and 100AU in a minimum mass solar nebula (MMSN) disk model, which consists of a far-UV-ionized surface layer and low-ionization disk interior. These simulations serve as a follow up to Simon et al. (2013), in which we found that without a net vertical field, the turbulent stresses that result from the magnetorotational instability (MRI) are too weak to account for observed accretion rates. The simulations in this work show a very strong dependence of the accretion stresses on the strength of the background vertical field; as the field strength increases, the stress amplitude increases. For gas to magnetic pressure ratios of 1e4 and 1e5, we find accretion rates between 1e-8 and 1e-7 solar masses per year. These accretion rates agree with observational constraints, suggesting a vertical magnetic field strength between 60 and 200 microgauss at 30AU and 10 and 30 microgauss at 100AU in a MMSN disk. Furthermore, the stress has a non-negligible component due to a magnetic wind. For sufficiently strong vertical field strengths, MRI turbulence is quenched, and the flow becomes largely laminar, with accretion proceeding through large scale correlations in the radial and toroidal field components as well as through the magnetic wind. In all simulations, the presence of a low ionization region near the disk mid-plane, which we call the ambipolar damping zone, results in reduced stresses there.


title: Gapped Symmetry Preserving Surface-State for the Electron Topological Insulator

abstract: It is well known that the 3D electronic topological insulator (TI) with charge-conservation and time-reversal symmetry cannot have a trivial insulating surface that preserves symmetry. It is often implicitly assumed that if the TI surface preserves both symmetries then it must be gapless. Here we show that it is possible for the TI surface to be both gapped and symmetry-preserving, at the expense of having surface-topological order. In contrast to analogous bosonic topological insulators, this symmetric surface topological order is intrinsically non-Abelian. We show that the surface-topological order provides a complete non-perturbative definition of the electron TI that transcends a free-particle band-structure picture, and could provide a useful perspective for studying strongly correlated topological Mott insulators.


title: Landau Levels, Anisotropy and Holography

abstract: We analyze properties of field theories dual to extremal black branes in (4+1) dimensions with anisotropic near-horizon geometries. Such gravity solutions were recently shown to fall into nine classes which align with the Bianchi classification of real three-dimensional Lie algebras. As a warmup we compute constraints on critical exponents from energy conditions in the bulk and scalar two point functions for a general type I metric, which has translation invariance but broken rotations. We also comment on the divergent nature of tidal forces in general Bianchi-type metrics. Then we come to our main focus: extremal branes whose near-horizon isometries are those of the Heisenberg algebra (type II). We find hyperscaling-violating solutions with type II isometries in (4+1)-dimensions. We show that scale invariant (4+1)-dimensional type II metrics are related by Kaluza-Klein reduction to more symmetric $AdS_2\times R^2$ and (3+1)-dimensional hyperscaling-violating spacetimes. These solutions generically have $\theta \leq 0$. We discuss how one can obtain flows from UV CFTs to Bianchi-type spacetimes in the IR via the Higgs mechanism, as well as potential inhomogeneous instabilities of type II. Finally, we compute two-point functions of massive and massless scalar operators in the dual field theory and find that they exhibit the behavior of a theory with Landau levels.


title: Astrometric orbit of a low-mass companion to an ultracool dwarf

abstract: Little is known about the existence of extrasolar planets around ultracool dwarfs. Furthermore, binary stars with Sun-like primaries and very low-mass binaries composed of ultracool dwarfs show differences in the distributions of mass ratio and orbital separation that can be indicative of distinct formation mechanisms. Using FORS2/VLT optical imaging for high precision astrometry we are searching for planets and substellar objects around ultracool dwarfs to investigate their multiplicity properties for very low companion masses. Here we report astrometric measurements with an accuracy of two tenths of a milli-arcsecond over two years that reveal orbital motion of the nearby L1.5 dwarf DENIS-P J082303.1-491201 located at 20.77 +/- 0.08 pc caused by an unseen companion that revolves about its host on an eccentric orbit in 246.4 +/- 1.4 days. We estimate the L1.5 dwarf to have 7.5 +/- 0.7 % of the Sun's mass that implies a companion mass of 28 +/- 2 Jupiter masses. This new system has the smallest mass ratio (0.36 +/- 0.02) of known very low-mass binaries with characterised orbits. With this discovery we demonstrate 200 micro-arcsecond astrometry over an arc-minute field and over several years that is sufficient to discover sub-Jupiter mass planets around ultracool dwarfs. We also show that the achieved parallax accuracy of < 0.4 % makes it possible to remove distance as a dominant source of uncertainty in the modelling of ultracool dwarfs.


title: A geometric crescent model for black hole images

abstract: The Event Horizon Telescope (EHT), a global very long baseline interferometry array operating at millimetre wavelengths, is spatially resolving the immediate environments of black holes for the first time. The current observations of the Galactic center black hole, Sagittarius A* (Sgr A*), and M87 have been interpreted in terms of either geometric models (e.g., a symmetric Gaussian) or detailed calculations of the appearance of black hole accretion flows. The former are not physically motivated, while the latter are subject to large systematic uncertainties. Motivated by the dominant relativistic effects of Doppler beaming and gravitational lensing in many calculations, we propose a geometric crescent model for black hole images. We show that this simple model provides an excellent statistical description of the existing EHT data of Sgr A* and M87, superior to other geometric models for Sgr A*. It also qualitatively matches physically predicted models, bridging accretion theory and observation. Based on our results, we make predictions for the detectability of the black hole shadow, a signature of strong gravity, in future observations.


title: KAT-7 Science Verification: Using HI Observations of NGC 3109 to Understand its Kinematics and Mass Distribution

abstract: HI observations of the Magellanic-type spiral NGC 3109, obtained with the seven dish Karoo Array Telescope (KAT-7), are used to analyze its mass distribution. Our results are compared to what is obtained using VLA data. KAT-7 is the precursor of the SKA pathfinder MeerKAT, which is under construction. The short baselines and low system temperature of the telescope make it sensitive to large scale low surface brightness emission. The new observations with KAT-7 allow the measurement of the rotation curve of NGC 3109 out to 32', doubling the angular extent of existing measurements. A total HI mass of 4.6 x 10^8 Msol is derived, 40% more than what was detected by the VLA observations. The observationally motivated pseudo-isothermal dark matter (DM) halo model can reproduce very well the observed rotation curve but the cosmologically motivated NFW DM model gives a much poorer fit to the data. While having a more accurate gas distribution has reduced the discrepancy between the observed RC and the MOdified Newtonian Dynamics (MOND) models, this is done at the expense of having to use unrealistic mass-to-light ratios for the stellar disk and/or very large values for the MOND universal constant a0. Different distances or HI contents cannot reconcile MOND with the observed kinematics, in view of the small errors on those two quantities. As for many slowly rotating gas-rich galaxies studied recently, the present result for NGC 3109 continues to pose a serious challenge to the MOND theory.


title: Resummation and S-duality in N=4 SYM

abstract: We consider the problem of resumming the perturbative expansions for anomalous dimensions of low twist, non-BPS operators in four dimensional N=4 supersymmetric Yang-Mills theories. The requirement of S-duality invariance imposes considerable restrictions on any such resummation. We introduce several prescriptions that produce interpolating functions on the upper half plane that are compatible with a subgroup of the full duality group. These lead to predictions for the anomalous dimensions at all points in the fundamental domain of the complex gauge coupling, and in particular at the duality-invariant values \tau=i and \tau=exp(i\pi/3). For low-rank gauge groups, the predictions are compatible with the bounds derived by conformal bootstrap methods for these anomalous dimensions; within numerical errors, they are in good agreement with the conjecture that said bounds are saturated at a duality-invariant point. We also find that the anomalous dimensions of the lowest twist operators lie within an extremely narrow window around a straight line as we vary the moduli of the theory over the two dimensional fundamental domain.


title: Non-Decoupling MSSM Higgs Sector and Light Superpartners

abstract: In the ``non-decoupling'' region of the Higgs sector in MSSM, the heavier CP-even Higgs boson (H^0) is Standard-Model-like and close to the charged Higgs bosons (H^\pm) in mass, while other neutral Higgs bosons (h^0,A^0) are lighter and near the Z mass. This scenario is consistent with the current Higgs search limits, although the improved sensitivity for a light charged Higgs boson search t -> H^+ b may result in certain degree of tension. We demonstrate that it can pass the stringent flavor constraints, provided there are other light SUSY particles to contribute in the loop induced processes. In turn, the non-decoupling Higgs sector implies the existence of light (left-handed) stop, sbottom and Wino-like gauginos, with mass all below 250 GeV. These light super-partners can still escape the current SUSY searches at the LHC. Dedicated searches for soft decay products should be devised for the LHC experiments to improve the searching sensitivity. The ILC would be able to cover the full spectrum region. The solutions for the viable SUSY parameters result from subtle cancellations and are often missed by the generic multiple dimensional scans, highlighting the importance of theoretical guidance in search for such special cases.


title: A Time-Reversal Invariant Topological Phase at the Surface of a 3D Topological Insulator

abstract: A 3D fermionic topological insulator has a gapless Dirac surface state protected by time-reversal symmetry and charge conservation symmetry. The surface state can be gapped by introducing ferromagnetism to break time-reversal symmetry, introducing superconductivity to break charge conservation, or entering a topological phase. In this paper, we construct a minimal gapped topological phase that preserves both time-reversal and charge conservation symmetries and supports Ising-type non-Abelian anyons. This phase can be understood heuristically as emerging from a surface s-wave superconducting state via the condensation of eight-vortex composites. The topological phase inherits vortices supporting Majorana zero modes from the surface superconducting state. However, since it is time-reversal invariant, the surface topological phase is a distinct phase from the Ising topological phase, which can be viewed as a quantum-disordered spin-polarized p_x + i p_y superconductor. We discuss the anyon model of this topological phase and the manner in which time-reversal symmetry is realized in it. We also study the interfaces between the topological state and other surface gapped phases.


title: A preferred ground state for the scalar field in de Sitter space

abstract: We investigate a recent proposal for a distinguished vacuum state of a free scalar quantum field in an arbitrarily curved spacetime, known as the Sorkin-Johnston (SJ) vacuum, by applying it to de Sitter space. We derive the associated two-point functions on both the global and Poincar\'e (cosmological) patches in general d+1 dimensions. In all cases where it is defined, the SJ vacuum belongs to the family of de Sitter invariant alpha-vacua. We obtain different states depending on the spacetime dimension, mass of the scalar field, and whether the state is evaluated on the global or Poincar\'e patch. We find that the SJ vacuum agrees with the Euclidean/Bunch-Davies state for heavy ("principal series") fields on the global patch in even spacetime dimensions. We also compute the SJ vacuum on a causal set corresponding to a causal diamond in 1+1 dimensional de Sitter space. Our simulations show that the mean of the SJ two-point function on the causal set agrees well with its expected continuum counterpart.


title: Quantum Instability of the Emergent Universe

abstract: We perform a semi-classical analysis of the Emergent Universe scenario for inflation. Fixing the background, and taking the inflaton to be homogenous, we cast the inflaton's evolution as a one-dimensional quantum mechanics problem. We find that the tuning required over a long time scale for this inflationary scenario is unstable quantum mechanically. Considering the inflaton field value as a wavepacket, the spreading of the wavepacket destroys any chance of both starting and ending with a well-formed state. Thus, one cannot have an Einstein static universe to begin with that evolves into a well-defined beginning to inflation a long time later.


title: Spherical collapse and halo mass function in f(R) theories

abstract: We compute the critical density of collapse for spherically symmetric overdensities in a class of f(R) modified gravity models. For the first time we evolve the Einstein, scalar field and non-linear fluid equations, making the minimal simplifying assumptions that the metric potentials and scalar field remain quasi-static throughout the collapse. Initially evolving a top hat profile, we find that the density threshold for collapse depends significantly on the initial conditions imposed, specifically the choice of size and shape. By imposing `natural' initial conditions, we obtain a fitting function for the spherical collapse delta_c as a function of collapse redshift, mass of the overdensity and f_{R0}, the background scalar field value at z=0. By extending delta_c into drifting and diffusing barrier within the context of excursion set theory, we obtain a realistic mass function that might be used to confront this class of scalar-tensor models with observations of dark matter halos. The proposed analytic formula for the halo mass function was tested against Monte Carlo random walks for a wide class of moving barriers and can therefore be applied to other modified gravity theories.


title: Standard Model Vacuum Stability and Weyl Consistency Conditions

abstract: At high energy the standard model possesses conformal symmetry at the classical level. This is reflected at the quantum level by relations between the different beta functions of the model. These relations are known as the Weyl consistency conditions. We show that it is possible to satisfy them order by order in perturbation theory, provided that a suitable coupling constant counting scheme is used. As a direct phenomenological application, we study the stability of the standard model vacuum at high energies and compare with previous computations violating the Weyl consistency conditions.


title: Lagrangian structures on mapping stacks and semi-classical TFTs

abstract: We extend a recent result of Pantev-Toen-Vaquie-Vezzosi, who constructed shifted symplectic structures on derived mapping stacks having a Calabi-Yau source and a shifted symplectic target. Their construction gives a clear conceptual framework for the so-called AKSZ formalism. We extend the PTVV construction to derived mapping stacks with boundary conditions, which is required in most applications to quantum field theories (see e.g. the work of Cattaneo-Felder on the Poisson sigma model, and the recent work of Cattaneo-Mnev-Reshetikhin). We provide many examples of Lagrangian and symplectic structures that can be recovered in this way. We finally give an application to topological field theories (TFTs). We expect that our approach will help to rigorously constuct a 2 dimensional TFT introduced by Moore and Tachikawa. A subsequent paper will be devoted to the construction of fully extended TFTs (in the sense of Baez-Dolan and Lurie) from mapping stacks.


title: Loops in reheating and cosmological perturbations

abstract: We show that in scalar field inflationary models, the loop corrections in reheating corresponding to the decay of the inflaton can cause nontrivial superhorizon evolution of the curvature perturbation. The effect turns out to be prominent when the decay occurs via parametric resonance, even indicating the breakdown of the perturbation theory, as we demonstrate in a specific model.


title: Pre-main-sequence isochrones -- II. Revising star and planet formation timescales

abstract: We have derived ages for 13 young (<30 Myr) star-forming regions and find they are up to a factor two older than the ages typically adopted in the literature. This result has wide-ranging implications, including that circumstellar discs survive longer (~10-12 Myr) and that the average Class I lifetime is greater (~1 Myr) than currently believed. For each star-forming region we derived two ages from colour-magnitude diagrams. First we fitted models of the evolution between the zero-age main-sequence and terminal-age main-sequence to derive a homogeneous set of main-sequence ages, distances and reddenings with statistically meaningful uncertainties. Our second age for each star-forming region was derived by fitting pre-main-sequence stars to new semi-empirical model isochrones. For the first time (for a set of clusters younger than 50 Myr) we find broad agreement between these two ages, and since these are derived from two distinct mass regimes that rely on different aspects of stellar physics, it gives us confidence in the new age scale. This agreement is largely due to our adoption of empirical colour-Teff relations and bolometric corrections for pre-main-sequence stars cooler than 4000 K. The revised ages for the star-forming regions in our sample are: ~2 Myr for NGC 6611 (Eagle Nebula; M 16), IC 5146 (Cocoon Nebula), NGC 6530 (Lagoon Nebula; M 8), and NGC 2244 (Rosette Nebula); ~6 Myr for {\sigma} Ori, Cep OB3b, and IC 348; ~10 Myr for {\lambda} Ori (Collinder 69); ~11 Myr for NGC 2169; ~12 Myr for NGC 2362; ~13 Myr for NGC 7160; ~14 Myr for {\chi} Per (NGC 884); and ~20 Myr for NGC 1960 (M 36).


title: Classification of interacting electronic topological insulators in three dimensions

abstract: A fundamental open problem in condensed matter physics is how the dichotomy between conventional and topological band insulators is modified in the presence of strong electron interactions. We show that there are 6 new electronic topological insulators that have no non-interacting counterpart. Combined with the previously known band-insulators, these produce a total of 8 topologically distinct phases. Two of the new topological insulators have a simple physical description as Mott insulators in which the electron spins form spin analogs of the familiar topological band-insulator. The remaining are obtained as combinations of these two `topological paramagnets' and the topological band insulator. We prove that these 8 phases form a complete list of all possible interacting topological insulators, and are classified by a Z_2^3 group-structure. Experimental signatures are also discussed for these phases.


title: A Higgs--Saw Mechanism as a Source for Dark Energy

abstract: Motivated by the see-saw mechanism for neutrinos which naturally generates small neutrino masses, we explore how a small GUT-scale mixing between the Standard Model Higgs and an otherwise massless hidden sector scalar can naturally generate a small mass and VEV for the new scalar which produces a false vacuum energy density contribution comparable to that of the observed dark energy dominating the current expansion of the Universe. The new scalar produces no discernible signatures in existing terrestrial experiments so that one may have to rely on other cosmological tests of this idea.


title: The impact of edges and dopants on the work function of graphene nanostructures. The way to high electronic emission from pure carbon medium

abstract: The impact of the edges and presence of dopants to the work function (WF) of graphene nanoribbons (GNR) and nanoflakes was studied by an ab initio approach. The strong dependence of the WF upon the GNR structure was found and a promising character for the field emission by the donor type impurities was observed. Basing on the predominant impact of the nanostructure edges to the emission properties, the small graphene flakes were investigated as a possible source for the electron emission. The obtained weak dependence of the low WF values of the graphene flakes on their size and shape allows to suggest that the pure carbon medium with high and uniform emission properties can be fabricated by today technology.


title: Testing the role of the Barbero-Immirzi parameter and the choice of connection in Loop Quantum Gravity

abstract: We study the role of the Barbero-Immirzi parameter $\gamma$ and the choice of connection in the construction of (a symmetry-reduced version of) loop quantum gravity. We start with the four-dimensional Lorentzian Holst action that we reduce to three dimensions in a way that preserves the presence of $\gamma$. In the time gauge, the phase space of the resulting three-dimensional theory mimics exactly that of the four-dimensional one. Its quantization can be performed, and on the kinematical Hilbert space spanned by SU(2) spin network states the spectra of geometric operators are discrete and $\gamma$-dependent. However, because of the three-dimensional nature of the theory, its SU(2) Ashtekar-Barbero Hamiltonian constraint can be traded for the flatness constraint of an sl(2,C) connection, and we show that this latter has to satisfy a linear simplicity-like condition analogous to the one used in the construction of spin foam models. The physically relevant solution to this constraint singles out the non-compact subgroup SU(1,1), which in turn leads to the disappearance of the Barbero-Immirzi parameter and to a continuous length spectrum, in agreement with what is expected from Lorentzian three-dimensional gravity.


title: Is F-term hybrid inflation natural within minimal supersymmetric SO(10)?

abstract: We examine whether F-term supersymmetric hybrid inflation can be embedded within the minimal SO(10) model, in a natural way. We show that none of the singlets of the Standard Model symmetries in the minimal set of SO(10) representations can satisfy the conditions which are necessary for a scalar field to play the role of the inflaton. As a consequence, one has to introduce an extra scalar field, which however may spoil the naturalness of inflation within the context of SO(10). Nevertheless, if we add an extra scalar field, we are then able to construct a model that can accommodate flat directions while it preserves the stability of the inflationary valley.


title: ABJM amplitudes and WL at finite N

abstract: We evaluate ABJM observables at two loops, for any value of the rank N of the gauge group. We compute the color subleading contributions to the four-point scattering amplitude in ABJM at two loops. Contrary to the four dimensional case, IR divergent N-subleading contributions are proportional to leading poles in the regularization parameter. We then exploit the non-planar calculation for the amplitude to derive an expression for the two-loop Sudakov form factor at any N. In the planar limit the result coincides with the one recently obtained in literature by using Feynman diagrams and unitarity. Finally, we analyze the subleading contributions to the light-like four-cusps Wilson loop and interpret the result in terms of the non-abelian exponentiation theorem. All these perturbative results satisfy the uniform transcendentality principle, hinting at its validity in ABJM beyond the planar limit.


title: The Scientific Reach of Multi-Ton Scale Dark Matter Direct Detection Experiments

abstract: The next generation of large scale WIMP direct detection experiments have the potential to go beyond the discovery phase and reveal detailed information about both the particle physics and astrophysics of dark matter. We report here on early results arising from the development of a detailed numerical code modeling the proposed DARWIN detector, involving both liquid argon and xenon targets. We incorporate realistic detector physics, particle physics and astrophysical uncertainties and demonstrate to what extent two targets with similar sensitivities can remove various degeneracies and allow a determination of dark matter cross sections and masses while also probing rough aspects of the dark matter phase space distribution. We find that, even assuming dominance of spin-independent scattering, multi-ton scale experiments still have degeneracies that depend sensitively on the dark matter mass, and on the possibility of isospin violation and inelasticity in interactions. We find that these experiments are best able to discriminate dark matter properties for dark matter masses less than around 200 GeV. In addition, and somewhat surprisingly, the use of two targets gives only a small improvement (aside from the advantage of different systematics associated with any claimed signal) in the ability to pin down dark matter parameters when compared with one target of larger exposure.


title: Chiral asymmetry in QED matter in a magnetic field

abstract: We calculate the electron self-energy in a magnetized QED plasma to the leading perturbative order in the coupling constant and to the linear order in an external magnetic field. We find that the chiral asymmetry of the normal ground state of the system is characterized by two new Dirac structures. One of them is the familiar chiral shift previously discussed in the Nambu--Jona-Lasinio model. The other structure is new. It formally looks like that of the chiral chemical potential but is an odd function of the longitudinal component of the momentum, directed along the magnetic field. The origin of this new parity-even chiral structure is directly connected with the long-range character of the QED interaction. The form of the Fermi surface in the weak magnetic field is determined.


title: Spectra of geometric operators in three-dimensional LQG: From discrete to continuous

abstract: We study and compare the spectra of geometric operators (length and area) in the quantum kinematics of two formulations of three-dimensional Lorentzian loop quantum gravity. In the SU(2) Ashtekar-Barbero framework, the spectra are discrete and depend on the Barbero-Immirzi parameter $\gamma$ exactly like in the four-dimensional case. However, we show that when working with the self-dual variables and imposing the reality conditions the spectra become continuous and $\gamma$-independent.


title: Program for quantum wave-packet dynamics with time-dependent potentials

abstract: We present a program to simulate the dynamics of a wave packet interacting with a time-dependent potential. The time-dependent Schr\"odinger equation is solved on a one-, two-, or three-dimensional spatial grid using the split operator method. The program can be compiled for execution either on a single processor or on a distributed-memory parallel computer.


title: Detecting initial system-environment correlations: Performance of various distance measures for quantum states

abstract: We study the time evolution of four distance measures in the presence of initial systemenvironment correlations. It is well-known that the trace distance between two quantum states of an open system may increase due to initial correlations which leads to a breakdown of the contractivity of the reduced dynamics. Here we compare and analyze, for two different models, the time evolution of the trace distance, the Bures metric, the Hellinger distance and the Jensen-Shannon divergence regarding an increase above their initial values, witnessing initial correlations. This work generalizes, deepens and corrects the study performed by Dajka et al. [Phys. Rev. A 84 032120 (2011)] and thereby reveals generic features of the considered distance measures with respect to the capability of detecting initial system-environment correlations.


title: Coisotropic submanifolds and dual pairs

abstract: The Poisson sigma model is a widely studied two-dimensional topological field theory. This note shows that boundary conditions for the Poisson sigma model are related to coisotropic submanifolds (a result announced in [math.QA/0309180]) and that the corresponding reduced phase space is a (possibly singular) dual pair between the reduced spaces of the given two coisotropic submanifolds. In addition the generalization to a more general tensor field is considered and it is shown that the theory produces Lagrangian evolution relations if and only if the tensor field is Poisson.


title: Symmetry Enforced Non-Abelian Topological Order at the Surface of a Topological Insulator

abstract: The surfaces of three dimensional topological insulators (3D TIs) are generally described as Dirac metals, with a single Dirac cone. It was previously believed that a gapped surface implied breaking of either time reversal $\mathcal T$ or U(1) charge conservation symmetry. Here we discuss a novel possibility in the presence of interactions, a surface phase that preserves all symmetries but is nevertheless gapped and insulating. Then the surface must develop topological order of a kind that cannot be realized in a 2D system with the same symmetries. We discuss candidate surface states - non-Abelian Quantum Hall states which, when realized in 2D, have $\sigma_{xy}=1/2$ and hence break $\mathcal T$ symmetry. However, by constructing an exactly soluble 3D lattice model, we show they can be realized as $\mathcal T$ symmetric surface states. The corresponding 3D phases are confined, and have $\theta=\pi$ magnetoelectric response. Two candidate states have the same 12 particle topological order, the (Read-Moore) Pfaffian state with the neutral sector reversed, which we term T-Pfaffian topological order, but differ in their $\mathcal T$ transformation. Although we are unable to connect either of these states directly to the superconducting TI surface, we argue that one of them describes the 3D TI surface, while the other differs from it by a bosonic topological phase. We also discuss the 24 particle Pfaffian-antisemion topological order (which can be connected to the superconducting TI surface) and demonstrate that it can be realized as a $\mathcal T$ symmetric surface state.


title: On Brauer groups of double covers of ruled surfaces

abstract: Let X be a smooth double cover of a geometrically ruled surface defined over a separably closed field of characteristic different from 2. The main result of this paper is a finite presentation of the 2-torsion in the Brauer group of X with generators given by central simple algebras over the function field of X and relations coming from the N\'eron-Severi group of X. In particular, the result gives a central simple algebra representative for the unique nontrivial Brauer class on any Enriques surface. An example demonstrating the applications to the study of rational points is given.


title: Global Stabilization of Nonlinear Delay Systems With a Compact Absorbing Set

abstract: Predictor-based stabilization results are provided for nonlinear systems with input delays and a compact absorbing set. The control scheme consists of an inter-sample predictor, a global observer, an approximate predictor, and a nominal controller for the delay-free case. The control scheme is applicable even to the case where the measurement is sampled and possibly delayed. The closed-loop system is shown to have the properties of global asymptotic stability and exponential convergence in the disturbance-free case, robustness with respect to perturbations of the sampling schedule, and robustness with respect to measurement errors. In contrast to existing predictor feedback laws, the proposed control scheme utilizes an approximate predictor of a dynamic type which is expressed by a system described by Integral Delay Equations. Additional results are provided for systems that can be transformed to systems with a compact absorbing set by means of a preliminary predictor feedback.


title: Pointing LISA-like gravitational wave detectors

abstract: Space-based gravitational wave detectors based on the Laser Interferometer Space Antenna (LISA) design operate by synthesizing one or more interferometers from fringe velocity measurements generated by changes in the light travel time between three spacecraft in a special set of drag-free heliocentric orbits. These orbits determine the inclination of the synthesized interferometer with respect to the ecliptic plane. Once these spacecraft are placed in their orbits, the orientation of the interferometers at any future time is fixed by Kepler's Laws based on the initial orientation of the spacecraft constellation, which may be freely chosen. Over the course of a full solar orbit, the initial orientation determines a set of locations on the sky were the detector has greatest sensitivity to gravitational waves as well as a set of locations where nulls in the detector response fall. By artful choice of the initial orientation, we can choose to optimize or suppress the antennas sensitivity to sources whose location may be known in advance (e.g., the Galactic Center or globular clusters).


title: Minimum-variance multitaper spectral estimation on the sphere

abstract: We develop a method to estimate the power spectrum of a stochastic process on the sphere from data of limited geographical coverage. Our approach can be interpreted either as estimating the global power spectrum of a stationary process when only a portion of the data are available for analysis, or estimating the power spectrum from local data under the assumption that the data are locally stationary in a specified region. Restricting a global function to a spatial subdomain -- whether by necessity or by design -- is a windowing operation, and an equation like a convolution in the spectral domain relates the expected value of the windowed power spectrum to the underlying global power spectrum and the known power spectrum of the localization window. The best windows for the purpose of localized spectral analysis have their energy concentrated in the region of interest while possessing the smallest effective bandwidth as possible. Solving an optimization problem in the sense of Slepian (1960) yields a family of orthogonal windows of diminishing spatiospectral localization, the best concentrated of which we propose to use to form a weighted multitaper spectrum estimate in the sense of Thomson (1982). Such an estimate is both more representative of the target region and reduces the estimation variance when compared to estimates formed by any single bandlimited window. We describe how the weights applied to the individual spectral estimates in forming the multitaper estimate can be chosen such that the variance of the estimate is minimized.


title: A spatiospectral localization approach to estimating potential fields on the surface of a sphere from noisy, incomplete data taken at satellite altitudes

abstract: Satellites mapping the spatial variations of the gravitational or magnetic fields of the Earth or other planets ideally fly on polar orbits, uniformly covering the entire globe. Thus, potential fields on the sphere are usually expressed in spherical harmonics, basis functions with global support. For various reasons, however, inclined orbits are favorable. These leave a "polar gap": an antipodal pair of axisymmetric polar caps without any data coverage, typically smaller than 10 degrees in diameter for terrestrial gravitational problems, but 20 degrees or more in some planetary magnetic configurations. The estimation of spherical harmonic field coefficients from an incompletely sampled sphere is prone to error, since the spherical harmonics are not orthogonal over the partial domain of the cut sphere. Although approaches based on wavelets have gained in popularity in the last decade, we present a method for localized spherical analysis that is firmly rooted in spherical harmonics. We construct a basis of bandlimited spherical functions that have the majority of their energy concentrated in a subdomain of the unit sphere by solving Slepian's (1960) concentration problem in spherical geometry, and use them for the geodetic problem at hand. Most of this work has been published by us elsewhere. Here, we highlight the connection of the "spherical Slepian basis" to wavelets by showing their asymptotic self-similarity, and focus on the computational considerations of calculating concentrated basis functions on irregularly shaped domains.


title: Black Hole Formation in Fuzzy Sphere Collapse

abstract: We study the collapse of a fuzzy sphere, that is a spherical membrane built out of D0-branes, in the BFSS model. At weak coupling, as the sphere shrinks, open strings are produced. If the initial radius is large then open string production is not important and the sphere behaves classically. At intermediate initial radius the back-reaction from open string production is important but the fuzzy sphere retains its identity. At small initial radius the sphere collapses to form a black hole. The crossover between the later two regimes is smooth and occurs at the correspondence point of Horowitz and Polchinski.


title: 3-color Bounded Patterned Self-assembly

abstract: Patterned self-assembly tile set synthesis PATS is the problem of finding a minimal tile set which uniquely self-assembles into a given pattern. Czeizler and Popa proved the NP-completeness of PATS and Seki showed that the PATS problem is already NP-complete for patterns with 60 colors. In search for the minimal number of colors such that PATS remains NP-complete, we introduce multiple bound PATS (mbPATS) where we allow bounds for the numbers of tile types of each color. We show that mbPATS is NP-complete for patterns with just three colors and, as a byproduct of this result, we also obtain a novel proof for the NP-completeness of PATS which is more concise than the previous proofs.


title: Dimerized Solids and Resonating Plaquette Order in SU(N)-Dirac Fermions

abstract: We study the quantum phases of fermions with an explicit SU(N)-symmetric, Heisenberg-like nearest-neighbor flavor exchange interaction on the honeycomb lattice at half-filling. Employing projective (zero temperature) quantum Monte Carlo simulations for even values of N, we explore the evolution from a weak-coupling semimetal into the strong-coupling, insulating regime. Furthermore, we compare our numerical results to a saddle-point approximation in the large-N limit. From the large-N regime down to the SU(6) case, the insulating state is found to be a columnar valence bond crystal, with a direct transition to the semimetal at weak, finite coupling, in agreement with the mean-field result in the large-N limit. At SU(4) however, the insulator exhibits a subtly different valence bond crystal structure, stabilized by resonating valence bond plaquettes. In the SU(2) limit, our results support a direct transition between the semimetal and an antiferromagnetic insulator.


title: The Submillimeter Polarization Spectrum of M17

abstract: We present 450 {\mu}m polarimetric observations of the M17 molecular cloud obtained with the SHARP polarimeter at the Caltech Submillimeter Observatory. Across the observed region, the magnetic field orientation is consistent with previous submillimeter and far-infrared polarization measurements. Our observations are centered on a region of the molecular cloud that has been compressed by stellar winds from a cluster of OB stars. We have compared these new data with previous 350 {\mu}m polarimetry and find an anti-correlation between the 450 to 350 {\mu}m polarization magnitude ratio and the ratio of 21 cm to 450 {\mu}m intensity. The polarization ratio is lower near the east end of the studied region where the cloud is exposed to stellar winds and radiation. At the west end of the region, the polarization ratio is higher. We interpret the varying polarization spectrum as evidence supporting the radiative alignment torque (RAT) model for grain alignment, implying higher alignment efficiency in the region that is exposed to a higher anisotropic radiation field.


title: Semilinearity and Context-Freeness of Languages Accepted by Valence Automata

abstract: Valence automata are a generalization of various models of automata with storage. Here, each edge carries, in addition to an input word, an element of a monoid. A computation is considered valid if multiplying the monoid elements on the visited edges yields the identity element. By choosing suitable monoids, a variety of automata models can be obtained as special valence automata. This work is concerned with the accepting power of valence automata. Specifically, we ask for which monoids valence automata can accept only context-free languages or only languages with semilinear Parikh image, respectively. First, we present a characterization of those graph products (of monoids) for which valence automata accept only context-free languages. Second, we provide a necessary and sufficient condition for a graph product of copies of the bicyclic monoid and the integers to yield only languages with semilinear Parikh image when used as a storage mechanism in valence automata. Third, we show that all languages accepted by valence automata over torsion groups have a semilinear Parikh image.


title: arXiv e-prints and the journal of record: An analysis of roles and relationships

abstract: Since its creation in 1991, arXiv has become central to the diffusion of research in a number of fields. Combining data from the entirety of arXiv and the Web of Science (WoS), this paper investigates (a) the proportion of papers across all disciplines that are on arXiv and the proportion of arXiv papers that are in the WoS, (b) elapsed time between arXiv submission and journal publication, and (c) the aging characteristics and scientific impact of arXiv e-prints and their published version. It shows that the proportion of WoS papers found on arXiv varies across the specialties of physics and mathematics, and that only a few specialties make extensive use of the repository. Elapsed time between arXiv submission and journal publication has shortened but remains longer in mathematics than in physics. In physics, mathematics, as well as in astronomy and astrophysics, arXiv versions are cited more promptly and decay faster than WoS papers. The arXiv versions of papers - both published and unpublished - have lower citation rates than published papers, although there is almost no difference in the impact of the arXiv versions of both published and unpublished papers.


title: Internal Gravity Waves in Massive Stars: Angular Momentum Transport

abstract: We present numerical simulations of internal gravity waves (IGW) in a star with a convective core and extended radiative envelope. We report on amplitudes, spectra, dissipation and consequent angular momentum transport by such waves. We find that these waves are generated efficiently and transport angular momentum on short timescales over large distances. We show that, as in the Earth's atmosphere, IGW drive equatorial flows which change magnitude and direction on short timescales. These results have profound consequences for the observational inferences of massive stars, as well as their long term angular momentum evolution. We suggest IGW angular momentum transport may explain many observational mysteries, such as: the misalignment of hot Jupiters around hot stars, the Be class of stars, Ni enrichment anomalies in massive stars and the non-synchronous orbits of interacting binaries.


title: Emery vs. Hubbard model for cuprate superconductors: a Composite Operator Method study

abstract: Within the Composite Operator Method (COM), we report the solution of the Emery model (also known as p-d or three band model), which is relevant for the cuprate high-Tc superconduc- tors. We also discuss the relevance of the often-neglected direct oxygen-oxygen hopping for a more accurate, sometimes unique, description of this class of materials. The benchmark of the solution is performed by comparing our results with the available quantum Monte Carlo ones. Both single- particle and thermodynamic properties of the model are studied in detail. Our solution features a metal-insulator transition at half filling. The resulting metal-insulator phase diagram agrees qual- itatively very well with the one obtained within Dynamical Mean-Field Theory. We discuss the type of transition (Mott-Hubbard (MH) or charge-transfer (CT)) for the microscopic (ab-initio) parameter range relevant for cuprates getting, as expected a CT type. The emerging single-particle scenario clearly suggests a very close relation between the relevant sub-bands of the three- (Emery) and the single- band (Hubbard) models, thus providing an independent and non-perturbative proof of the validity of the mapping between the two models for the model parameters optimal to describe cuprates. Such a result confirms the emergence of the Zhang-Rice scenario, which has been recently questioned. We also report the behavior of the specific heat and of the entropy as functions of the temperature on varying the model parameters as these quantities, more than any other, depend on and, consequently, reveal the most relevant energy scales of the system.


title: The Intrinsic Shape of Galaxies in SDSS/Galaxy Zoo

abstract: By modelling the axis ratio distribution of SDSS DR8 galaxies we find the intrinsic 3D shapes of spirals and ellipticals. We use morphological information from the Galaxy Zoo project and assume a non-parametric distribution intrinsic of shapes, while taking into account dust extinction. We measure the dust extinction of the full sample of spiral galaxies and find a smaller value than previous estimations, with an edge-on extinction of $E_0 = 0.284^{+0.015}_{-0.026}$ in the SDSS r band. We also find that the distribution of minor to major axis ratio has a mean value of $ 0.267 \pm 0.009 $, slightly larger than previous estimates mainly due to the lower extinction used; the same affects the circularity of galactic discs, which are found to be less round in shape than in previous studies, with a mean ellipticity of $0.215 \pm 0.013$. For elliptical galaxies, we find that the minor to major axis ratio, with a mean value of $0.584 \pm 0.006$, is larger than previous estimations due to the removal of spiral interlopers present in samples with morphological information from photometric profiles. These interlopers are removed when selecting ellipticals using Galaxy Zoo data. We find that the intrinsic shapes of galaxies and their dust extinction vary with absolute magnitude, colour and physical size. We find that bright elliptical galaxies are more spherical than faint ones, a trend that is also present with galaxy size, and that there is no dependence of elliptical galaxy shape with colour. For spiral galaxies we find that the reddest ones have higher dust extinction as expected, due to the fact that this reddening is mainly due to dust. We also find that the thickness of discs increases with luminosity and size, and that brighter, smaller and redder galaxies have less round discs.


title: Painleve Field Theory

abstract: We propose multidimensional versions of the Painlev\'e VI equation and its degenerations. These field theories are related to the isomonodromy problems of flat holomorphic infinite rank bundles over elliptic curves and take the form of non-autonomous Hamiltonian equations. The modular parameter of curves plays the role of "time". Reduction of the field equations to the zero modes leads to ${\rm SL}(N, {\mathbb C})$ monodromy preserving equations. The latter coincide with the Painlev\'e VI equation for $N\!=\!2$. We consider two types of the bundles. In the first one the group of automorphisms is the centrally and cocentrally extended loop group $L({\rm SL}(N, {\mathbb C}))$ or some multiloop group. In the case of the Painlev\'e VI field theory in D=1+1 four constants of the Painlev\'e VI equation become dynamical fields. The second type of bundles are defined by the group of automorphisms of the noncommutative torus. They lead to the equations in dimension 2+1. In both cases we consider trigonometric, rational and scaling limits of the theories. Generically (except some degenerate cases) the derived equations are nonlocal. We consider Whitham quasiclassical limit to integrable systems. In this way we derive two and three dimensional integrable nonlocal versions of the integrable Euler-Arnold tops.


title: Water loss from terrestrial planets with CO2-rich atmospheres

abstract: Water photolysis and hydrogen loss from the upper atmospheres of terrestrial planets is of fundamental importance to climate evolution but remains poorly understood in general. Here we present a range of calculations we performed to study the dependence of water loss rates from terrestrial planets on a range of atmospheric and external parameters. We show that CO2 can only cause significant water loss by increasing surface temperatures over a narrow range of conditions, with cooling of the middle and upper atmosphere acting as a bottleneck on escape in other circumstances. Around G-stars, efficient loss only occurs on planets with intermediate CO2 atmospheric partial pressures (0.1 to 1 bar) that receive a net flux close to the critical runaway greenhouse limit. Because G-star total luminosity increases with time but XUV/UV luminosity decreases, this places strong limits on water loss for planets like Earth. In contrast, for a CO2-rich early Venus, diffusion limits on water loss are only important if clouds caused strong cooling, implying that scenarios where the planet never had surface liquid water are indeed plausible. Around M-stars, water loss is primarily a function of orbital distance, with planets that absorb less flux than ~270 W/m^2 (global mean) unlikely to lose more than one Earth ocean of H2O over their lifetimes unless they lose all their atmospheric N2/CO2 early on. Because of the variability of H2O delivery during accretion, our results suggest that many 'Earth-like' exoplanets in the habitable zone may have ocean-covered surfaces, stable CO2/H2O-rich atmospheres, and high mean surface temperatures.


title: Oscillatory instability of fully 3D flow in a cubic diagonally lid-driven cavity

abstract: A transition to unsteadiness of a flow inside a cubic diagonally lid-driven cavity with no-slip boundaries is numerically investigated by a series of direct numerical simulations (DNS) performed on 100^3 and 200^3 stretched grids. It is found that the observed oscillatory instability is setting in via a subcritical symmetry-breaking Hopf bifurcation. The instability evolves on two vortices in a coupled manner. Critical values of Reynolds number Recr=2320 and non-dimensional angular oscillating frequency omegacr=0.249 for transition from steady to oscillatory flow are accurately estimated. Characteristic patterns of the 3D oscillatory flow are presented.


title: Repulsive Interactions in Quantum Hall Systems as a Pairing Problem

abstract: A subtle relation between Quantum Hall physics and the phenomenon of pairing is unveiled. By use of second quantization, we establish a connection between (i) a broad class of rotationally symmetric two-body interactions within the lowest Landau level and (ii) integrable hyperbolic Richardson-Gaudin type Hamiltonians that arise in (p_{x}+ip_{y}) superconductivity. Specifically, we show that general Haldane pseudopotentials (and their sums) can be expressed as a sum of repulsive non-commuting (p_{x}+ip_{y})-type pairing Hamiltonians. For the Laughlin sequence, it is observed that this problem is frustration free and zero energy ground states lie in the common null space of all of these non-commuting Hamiltonians. This property allows for the use of a new truncated basis of pairing configurations in which to express Laughlin states at general filling factors. We prove separability of arbitrary Haldane pseudopotentials, providing explicit expressions for their second quantized forms, and further show by explicit construction how to exploit the topological equivalence between different geometries (disk, cylinder, and sphere) sharing the same topological genus number, in the second quantized formalism, through similarity transformations. As an application of the second quantized approach, we establish a "squeezing principle" that applies to the zero modes of a general class of Hamiltonians, which includes but is not limited to Haldane pseudopotentials. We also show how one may establish (bounds on) "incompressible filling factors" for those Hamiltonians. By invoking properties of symmetric polynomials, we provide explicit second quantized quasi-hole generators; the generators that we find directly relate to bosonic chiral edge modes and further make aspects of dimensional reduction in the Quantum Hall systems precise.

categories: cond-mat.str-el cond-mat.supr-con math-ph math.MP quant-ph


title: Nuclear Polarization Corrections to mu-d Atoms in Zero-Range Approximation

abstract: Nuclear polarization corrections to the 2P-2S Lamb shift in mu-d atoms are developed in order alpha^5 and are shown to agree with a recent calculation. The nuclear physics in the resulting corrections is then evaluated in zero-range approximation. The dominant part of the correction is very simple in form and differs from a recent potential model calculation by less than 1%. It is also demonstrated how the third-Zemach-moment contribution largely cancels against part of the polarization correction, as it did in e-d atoms and does so exactly for point-like nucleons. This suggests that it may be possible to reduce the uncertainty in the theory (of which nuclear polarization is the largest contributor) to less than 1%.


title: Interacting viscous dark fluids

abstract: We revise the conditions for the physical viability of a cosmological model in which dark matter has bulk viscosity and also interacts with dark energy. We have also included radiation and baryonic matter components; all matter components are represented by perfect fluids, except the dark matter, that is treated as an imperfect fluid. We impose upon the model the condition of a complete cosmological dynamics that results in an either null or negative bulk viscosity, but the latter also disagrees with the Local Second Law of Thermodynamics. The model is also compared with cosmological observations at different redshifts: type Ia supernova, the shift parameter of CMB, the acoustic peak of BAO, and the Hubble parameter H(z). In general, observations consistently point out to a negative value of the bulk viscous coefficient, and in overall the fitting procedure shows no preference for the model over the standard LCDM model.


title: The calculation of prompt fission neutron spectrum for 233U(n,f) reaction by the semi-empirical method

abstract: The prompt fission neutron spectra for neutron-induced fission of 233U for low energy neutrons (below 6 MeV) are calculated using the nuclear evaporation theory with a semi-empirical method, in which the partition of the total excitation energy between the fission fragments for the nth+233U fission reactions are determined with the available experimental and evaluation data. The calculated prompt fission neutron spectra agree well with the experimental data. The proportions of high- energy outgoing neutrons of prompt fission neutron spectrum versus incident neutron energies are investigated with the theoretical spectra, and the results are consistent with the systematics. We conclude that the semi-empirical method is sound, and it could be a useful tool for prompt fission neutron spectra evaluation.


title: A method for evaluating the expectation value of a power spectrum using the probability density function of phases

abstract: Here, we present a new method to evaluate the expectation value of the power spectrum of a time series. A statistical approach is adopted to define the method. After its demonstration, it is validated showing that it leads to the known properties of the power spectrum when the time series contains a periodic signal. The approach is also validated in general with numerical simulations. The method puts into evidence the importance that is played by the probability density function of the phases associated to each time stamp for a given frequency, and how this distribution can be perturbed by the uncertainties of the parameters in the pulsar ephemeris. We applied this method to solve the power spectrum in the case the first derivative of the pulsar frequency is unknown and not negligible. We also undertook the study of the most general case of a blind search, in which both the frequency and its first derivative are uncertain. We found the analytical solutions of the above cases invoking the sum of Fresnel's integrals squared.


title: Performance of commercial CMOS cameras for high-speed multicolor photometry

abstract: We present some results of testing of commercial color CMOS cameras for astronomical applications. CMOS sensors allow to perform photometry in three filters simultaneously that gives a great advantage compared with monochrome CCD detectors. The Bayer BGR colour system realized in CMOS sensors is close to the Johnson BVR system. We demonstrate transformation from the Bayer color system to the Johnson one. Our photometric measurements with color CMOS cameras coupled to small telescopes (11 - 30 inch) reveal that in video mode stars up to V $\sim$ 9 can be shot at 24 frames per second. Using a high-speed CMOS camera with short exposure times (10 - 20 ms) we can perform an imaging mode called "lucky imaging". We can pick out high quality frames and combine them into a single image using "shift-and-add" technique. This allows us obtain an image with much higher resolution than would be possible shooting a single image with long exposure. For image selection we use the Strehl-selection method. We demonstrates advantage of the lucky imaging technique in comparison with long exposure shooting. The FWHM of the blurred image caused by atmosphere turbulence can be decreased twice and more.


title: The exit time of planar Brownian motion and the Phragmen-Lindelof principle

abstract: In this note, we prove a version of the Phragmen-Lindelof principle using probabilistic techniques. In particular, we will show that if the p-th moment of the exit time of Brownian motion from a planar domain is finite, then an analytic function on that domain is either bounded by its supremum on the boundary or else goes to infinity along some sequence more rapidly than $e^{|z|^{2p}}$. We also present a method for constructing domains whose exit time has finite p-th moment, thereby giving a general Phragmen-Lindelof principle for spiral-like and star-like domains, and in the process giving a probabilistic proof of a theorem of Hansen. A number of auxiliary results are presented as well.


title: A slope conjecture for links

abstract: The slope conjecture gives a precise relation between the degree of the colored Jones polynomial of a knot and the boundary slopes of essential surfaces in the knot complement. In this note we propose a generalization of the slope conjecture to links. We prove the conjecture for all alternating and more generally adequate links. We also verify the conjecture for torus links.


title: Virtual neighborhood technique for pseudo-holomorphic spheres

abstract: This is the first part of a trilogy where we apply the theory of virtual manifold/orbifolds developed by the first named author and Tian to study the Gromov-Witten moduli spaces. In this paper, we resolve the main analytic issue arising from the lack of differentiability of $PSL(2, \C)$-action on spaces of $W^{1, p}$-maps from the Riemann sphere to a symplectic manifold $(X, \omega, J)$ with a non-zero homology class $A$. In particular, we establish the slice and tubular neighbourhood theorems for $PSL(2, \C)$-action along smooth maps, and construct a $PSL(2, \C)$-obstruction bundle along $PSL(2, \C)$-orbit of a pseudo-holomorphic map representing a point in the moduli space $\cM_{0, 0}(X, A)$. In Sections 2 and 3 of this paper, we explain an integration theory on virtual orbifolds using proper \'etale groupoids and establish the virtual neighborhood technique for a general orbifold Fredholm system. When the moduli space $\cM_{0, 0}(X, A)$ of pseudo-holomorphic spheres in $(X, \omega, J)$ is compact, applying the virtual neighborhood technique developed in Section 3, we obtain a virtual system for the moduli space $\cM_{0, 0}(X, A)$ of pseudo-holomorphic spheres in $(X, \omega, J)$ and show that the genus zero Gromov-Witten invariant is well-defined.


title: Bayesian State-Space Modelling on High-Performance Hardware Using LibBi

abstract: LibBi is a software package for state-space modelling and Bayesian inference on modern computer hardware, including multi-core central processing units (CPUs), many-core graphics processing units (GPUs) and distributed-memory clusters of such devices. The software parses a domain-specific language for model specification, then optimises, generates, compiles and runs code for the given model, inference method and hardware platform. In presenting the software, this work serves as an introduction to state-space models and the specialised methods developed for Bayesian inference with them. The focus is on sequential Monte Carlo (SMC) methods such as the particle filter for state estimation, and the particle Markov chain Monte Carlo (PMCMC) and SMC^2 methods for parameter estimation. All are well-suited to current computer hardware. Two examples are given and developed throughout, one a linear three-element windkessel model of the human arterial system, the other a nonlinear Lorenz '96 model. These are specified in the prescribed modelling language, and LibBi demonstrated by performing inference with them. Empirical results are presented, including a performance comparison of the software with different hardware configurations.


title: A decomposition theorem for immersions of product manifolds

abstract: We introduce polar metrics on a product manifold, which have product and warped product metrics as special cases. We prove a de Rham-type theorem characterizing Riemannian manifolds that can be locally decomposed as a product manifold endowed with a polar metric. For a product manifold endowed with a polar metric, our main result gives a complete description of all its isometric immersions into a space form whose second fundamental forms are adapetd to its product structure, in the sense that the tangent spaces to each factor are preserved by all shape operators. This is a far-reaching generalization of a basic decomposition theorem for isometric immersions of Riemannian products due to Moore as well as its extension by N\"olker to isometric immersions of warped products.


title: An analysis of B meson decays to two light pseudoscalar nonets with deltaS=1 in QCD factorization approach

abstract: We have performed a numerical analysis of branching ratios and direct CP-asymmetries of B meson decays to two light pseudoscalar nonets with deltaS=1 using QCD improved factorization. The parameters related to B meson including its semileptonic form factors, strange quark mass and the CKM phase angle have been varied over a limited range and phenomenological parameters used to parameterize divergences in hard spectator scattering and weak annihilation have been varied over a wide range simultaneously with an aim to achieve a fit of branching ratios and some direct CP-asymmetries. A good fit for the majority of the parameters has been achieved indicating that the standard model is quite satisfactory in accounting for these processes.


abstract: We define Eisenstein series on rank 2 hyperbolic Kac--Moody groups over R, induced from quasi--characters. We prove convergence of the constant term and hence the almost everywhere convergence of the Eisenstein series. We define and calculate the degenerate Fourier coefficients. We also consider Eisenstein series induced from cusp forms and show that these are entire functions.


title: The orbital periods of subdwarf B binaries produced by the first stable Roche overflow channel

abstract: Long-orbital-period subdwarf B (sdB) stars with main-sequence companions are believed to be the product of stable Roche Lobe overflow (RLOF), a scenario challenged by recent observations. Here we represent the results of a systematic study of the orbital-period distribution of sdB binaries in this channel using detailed binary evolution calculations. We show that the observed orbital-period distribution of long-period sdB binaries can be well explained by this scenario. Furthermore, we find that, if the progenitors of the sdB stars have initial masses below the helium flash mass, the sdB binaries produced from stable RLOF follow a unique mass -- orbital period relation for a given metallicity $Z$; increasing the orbital period from $\sim 400$ to $\sim 1100$\,d corresponds to increasing the mass of the sdB star from $\sim 0.40$ to $\sim 0.49\,M_\odot$ for $Z=0.02$. We suggest that the longest sdB binaries (with orbital period $> 1100$\,d) could be the result of atmospheric RLOF. The mass -- orbital period relation can be tested observationally if the mass of the sdB star can be determined precisely, e.g.\ from asteroseismology. Using this relation, we revise the orbital period distribution of sdB binaries produced by the first stable RLOF channel for the best fitting model of Han et al (2003), and show that the orbital period has a peak around 830\,d.


title: A Pseudo su(1,1)-Algebraic Deformation of the Cooper-Pair in the su(2)-Algebraic Many-Fermion Model

abstract: A pseudo su(1,1)-algebra is formulated as a possible deformation of the Cooper-pair in the su(2)-algebraic many-fermion system. With the aid of this algebra, it is possible to describe behavior of individual fermions which are generated as the result of interaction with the external environment. The form presented in this paper is a generalization of a certain simple case developed recently by the present authors. Basic idea follows the su(1,1)-algebra in the Schwinger boson representation for treating energy transfer between the harmonic oscillator and the external environment. Hamiltonian is given under the idea of the phase space doubling in the thermo-field dynamics formalism and the time-dependent variational method is applied to this Hamiltonian. Its trial state is constructed in the frame deformed from the BCS-Bogoliubov approach to the superconductivity. Several numerical results are shown.


title: Quark Propagation in the Instantons of Lattice QCD

abstract: We quantitatively examine the extent to which instanton degress of freedom, contained within standard Monte-carlo generated gauge-field configurations, can maintain the characteristic features of the mass and renormalisation functions of the non-perturbative quark propagator. We use over-improved stout-link smearing to isolate instanton effects on the lattice. Using a variety of measures, we illustrate how gauge fields consisting almost solely of instanton-like objects are produced after only 50 sweeps of smearing. We find a full vacuum, with a packing fraction more than three times larger than phenomenological models predict. We calculate the overlap quark propagator on these smeared configurations, and find that even at high levels of smearing the majority of the characteristic features of the propagator are reproduced. We thus conclude that instantons contained within standard Monte-carlo generated gauge-field configurations are the degrees of freedom responsible for the dynamical generation of mass observed in lattice QCD.


title: All-Distances Sketches, Revisited: HIP Estimators for Massive Graphs Analysis

abstract: Graph datasets with billions of edges, such as social and Web graphs, are prevalent, and scalable computation is critical. All-distances sketches (ADS) [Cohen 1997], are a powerful tool for scalable approximation of statistics. The sketch is a small size sample of the distance relation of a node which emphasizes closer nodes. Sketches for all nodes are computed using a nearly linear computation and estimators are applied to sketches of nodes to estimate their properties. We provide, for the first time, a unified exposition of ADS algorithms and applications. We present the Historic Inverse Probability (HIP) estimators which are applied to the ADS of a node to estimate a large natural class of statistics. For the important special cases of neighborhood cardinalities (the number of nodes within some query distance) and closeness centralities, HIP estimators have at most half the variance of previous estimators and we show that this is essentially optimal. Moreover, HIP obtains a polynomial improvement for more general statistics and the estimators are simple, flexible, unbiased, and elegant. For approximate distinct counting on data streams, HIP outperforms the original estimators for the HyperLogLog MinHash sketches (Flajolet et al. 2007), obtaining significantly improved estimation quality for this state-of-the-art practical algorithm.


title: Third-Order Edge Statistics: Contour Continuation, Curvature, and Cortical Connections

abstract: Association field models have attempted to explain human contour grouping performance, and to explain the mean frequency of long-range horizontal connections across cortical columns in V1. However, association fields only depend on the pairwise statistics of edges in natural scenes. We develop a spectral test of the sufficiency of pairwise statistics and show there is significant higher order structure. An analysis using a probabilistic spectral embedding reveals curvature-dependent components.


title: A symmetry-respecting topologically-ordered surface phase of 3d electron topological insulators

abstract: A 3d electron topological insulator (ETI) is a phase of matter protected by particle-number conservation and time-reversal symmetry. It was previously believed that the surface of an ETI must be gapless unless one of these symmetries is broken. A well-known symmetry-preserving, gapless surface termination of an ETI supports an odd number of Dirac cones. In this paper we deduce a symmetry-respecting, gapped surface termination of an ETI, which carries an intrinsic 2d topological order, Moore-Read x anti-semion. The Moore-Read sector supports non-Abelian charge 1/4 anyons, while the Abelian, anti-semion sector is electrically neutral. Time-reversal symmetry is implemented in this surface phase in a highly non-trivial way. Moreover, it is impossible to realize this phase strictly in 2d, simultaneously preserving its implementation of both the particle number and time-reversal symmetries. A 1d edge on the ETI surface between the topologically-ordered phase and the topologically trivial time-reversal-broken phase with a Hall conductivity sigma_H = 1/2, carries a right-moving neutral Majorana mode, a right-moving bosonic charge mode and a left-moving bosonic neutral mode. The topologically-ordered phase is separated from the surface superconductor by a direct second order phase transition in the XY* universality class, which is driven by the condensation of a charge 1/2 boson, when approached from the topologically-ordered side, and proliferation of a flux 4\pi-vortex, when approached from the superconducting side. In addition, we prove that time-reversal invariant (interacting) electron insulators with no intrinsic topological order and electromagnetic response characterized by a \theta-angle, \theta = \pi, do not exist if the electrons transform as Kramers singlets under time-reversal.