Multi-mass solvers for lattice QCD on GPUs

Graphical Processing Units (GPUs) are more and more frequently used for lattice QCD calculations. Lattice studies often require computing the quark propagators for several masses. These systems can be solved using multi-shift inverters but these algorithms are memory intensive which limits the size of the problem that can be solved using GPUs. In this paper, we show how to efficiently use a memory-lean single-mass inverter to solve multi-mass problems. We focus on the BiCGstab algorithm for Wilson fermions and show that the single-mass inverter not only requires less memory but also outperforms the multi-shift variant by a factor of two.Comment: 27 pages, 6 figures, 3 Table

Two-band superconductors: Extended Ginzburg-Landau formalism by a systematic expansion in small deviation from the critical temperature

We derive the extended Ginzburg-Landau (GL) formalism for a clean s-wave two-band superconductor by employing a systematic expansion of the free-energy functional and the corresponding matrix gap equation in powers of the small deviation from the critical temperature tau = 1-T/T_c. The two lowest orders of this expansion produce the equation for T_c and the GL theory. It is shown that in agreement with previous studies, the two-band GL theory maps onto the single-band GL model and thus fails to describe the difference in the spatial profiles of the two band condensates. We prove that except for some very special cases, this difference appears already in the leading correction to the GL theory, which constitutes the extended GL formalism. We derive linear differential equations that determine the leading corrections to the band order parameters and magnetic field, discuss the validity of these equations, and consider examples of an important interplay between the band condensates. Finally, we present numerical results for the thermodynamic critical magnetic field and temperature-dependent band gaps (at zero field), which are in a very good agreement with those obtained from the full BCS approach in a wide temperature range. To this end, we emphasize the advantages of our extended GL theory in comparison with the often used two-component GL-like model based on an unreconstructed two-band generalization of the Gor'kov derivation

Limits on non-Gaussianities from WMAP data

We develop a method to constrain the level of non-Gaussianity of density perturbations when the 3-point function is of the "equilateral" type. Departures from Gaussianity of this form are produced by single field models such as ghost or DBI inflation and in general by the presence of higher order derivative operators in the effective Lagrangian of the inflaton. We show that the induced shape of the 3-point function can be very well approximated by a factorizable form, making the analysis practical. We also show that, unless one has a full sky map with uniform noise, in order to saturate the Cramer-Rao bound for the error on the amplitude of the 3-point function, the estimator must contain a piece that is linear in the data. We apply our technique to the WMAP data obtaining a constraint on the amplitude f_NL^equil of "equilateral" non-Gaussianity: -366 < f_NL^equil < 238 at 95% C.L. We also apply our technique to constrain the so-called "local" shape, which is predicted for example by the curvaton and variable decay width models. We show that the inclusion of the linear piece in the estimator improves the constraint over those obtained by the WMAP team, to -27 < f_NL^local < 121 at 95% C.L.Comment: 20 pages, 12 eps figure

Exact Solution for the Critical State in Thin Superconductor Strips with Field Dependent or Anisotropic Pinning

An exact analytical solution is given for the critical state problem in long thin superconductor strips in a perpendicular magnetic field, when the critical current density j_c(B) depends on the local induction B according to a simple three-parameter model. This model describes both isotropic superconductors with this j_c(B) dependence, but also superconductors with anisotropic pinning described by a dependence j_c(theta) where theta is the tilt angle of the flux lines away from the normal to the specimen plane

Probing local non-Gaussianities within a Bayesian framework

Aims: We outline the Bayesian approach to inferring f_NL, the level of non-Gaussianity of local type. Phrasing f_NL inference in a Bayesian framework takes advantage of existing techniques to account for instrumental effects and foreground contamination in CMB data and takes into account uncertainties in the cosmological parameters in an unambiguous way. Methods: We derive closed form expressions for the joint posterior of f_NL and the reconstructed underlying curvature perturbation, Phi, and deduce the conditional probability densities for f_NL and Phi. Completing the inference problem amounts to finding the marginal density for f_NL. For realistic data sets the necessary integrations are intractable. We propose an exact Hamiltonian sampling algorithm to generate correlated samples from the f_NL posterior. For sufficiently high signal-to-noise ratios, we can exploit the assumption of weak non-Gaussianity to find a direct Monte Carlo technique to generate independent samples from the posterior distribution for f_NL. We illustrate our approach using a simplified toy model of CMB data for the simple case of a 1-D sky. Results: When applied to our toy problem, we find that, in the limit of high signal-to-noise, the sampling efficiency of the approximate algorithm outperforms that of Hamiltonian sampling by two orders of magnitude. When f_NL is not significantly constrained by the data, the more efficient, approximate algorithm biases the posterior density towards f_NL = 0.Comment: 11 pages, 7 figures. Accepted for publication in Astronomy and Astrophysic

A smooth bouncing cosmology with scale invariant spectrum

We present a bouncing cosmology which evolves from the contracting to the expanding phase in a smooth way, without developing instabilities or pathologies and remaining in the regime of validity of 4d effective field theory. A nearly scale invariant spectrum of perturbations is generated during the contracting phase by an isocurvature scalar with a negative exponential potential and then converted to adiabatic. The model predicts a slightly blue spectrum, n_S >~ 1, no observable gravitational waves and a high (but model dependent) level of non-Gaussianities with local shape. The model represents an explicit and predictive alternative to inflation, although, at present, it is clearly less compelling.Comment: 20 pages, 1 fig. v2: references added, JCAP published versio

An Improved Calculation of the Non-Gaussian Halo Mass Function

The abundance of collapsed objects in the universe, or halo mass function, is an important theoretical tool in studying the effects of primordially generated non-Gaussianities on the large scale structure. The non-Gaussian mass function has been calculated by several authors in different ways, typically by exploiting the smallness of certain parameters which naturally appear in the calculation, to set up a perturbative expansion. We improve upon the existing results for the mass function by combining path integral methods and saddle point techniques (which have been separately applied in previous approaches). Additionally, we carefully account for the various scale dependent combinations of small parameters which appear. Some of these combinations in fact become of order unity for large mass scales and at high redshifts, and must therefore be treated non-perturbatively. Our approach allows us to do this, and to also account for multi-scale density correlations which appear in the calculation. We thus derive an accurate expression for the mass function which is based on approximations that are valid over a larger range of mass scales and redshifts than those of other authors. By tracking the terms ignored in the analysis, we estimate theoretical errors for our result and also for the results of others. We also discuss the complications introduced by the choice of smoothing filter function, which we take to be a top-hat in real space, and which leads to the dominant errors in our expression. Finally, we present a detailed comparison between the various expressions for the mass functions, exploring the accuracy and range of validity of each.Comment: 28 pages, 13 figures; v2: text reorganized and some figured modified for clarity, results unchanged, references added. Matches version published in JCA