Caching and Interpolated Likelihoods: Accelerating Cosmological Monte Carlo Markov Chains
We describe a novel approach to accelerating Monte Carlo Markov Chains. Our focus is cosmological parameter estimation, but the algorithm is applicable to any problem for which the likelihood surface is a smooth function of the free parameters and computationally expensive to evaluate. We generate a high-order interpolating polynomial for the log-likelihood using the first points gathered by the Markov chains as a training set. This polynomial then accurately computes the majority of the likelihoods needed in the latter parts of the chains. We implement a simple version of this algorithm as a patch (InterpMC) to CosmoMC and show that it accelerates parameter estimatation by a factor of between two and four for well-converged chains. The current code is primarily intended as a "proof of concept", and we argue that there is considerable room for further performance gains. Unlike other approaches to accelerating parameter fits, we make no use of precomputed training sets or special choices of variables, and InterpMC is almost entirely transparent to the user.
Nonlinear matter spectra in growing neutrino quintessence
We investigate the nonlinear power spectra of density perturbations and acoustic oscillations in growing neutrino quintessence. In this scenario, the neutrino mass has a strong dependence on the quintessence field. The induced coupling stops the evolution of the field when the neutrinos become nonrelativistic, and triggers the transition to the accelerating phase of the cosmological expansion. At redshifts around five, the neutrino fluctuations are still linear and acoustic oscillations are present in the neutrino power spectrum, induced by the acoustic oscillations in the baryonic and dark-matter sectors. The neutrino perturbations become nonlinear at redshifts around three. The mode coupling generated by the nonlinearities erases the oscillations in the neutrino spectrum at some redshift above two. There is a potential danger that at later times the influence of the gravitational potentials induced by the neutrino inhomogeneities could erase the oscillations from the baryonic and dark-matter spectra, making the scenario incompatible with observations. For the scenario to be viable, the neutrino-induced gravitational potentials in the range of baryonic acoustic oscillations should not grow to average values much larger than 10^{-4}. The magnitude of the expected potentials is still not known reliably, as the process of structure formation is poorly understood in growing neutrino quintessence.
Physics of the Cosmic Microwave Background and the Planck Mission
This lecture is a sketch of the physics of the cosmic microwave background. The observed anisotropy can be divided into four main contributions: variations in the temperature and gravitational potential of the primordial plasma, Doppler effect from its motion, and a net red/blueshift the photons accumulate from traveling through evolving gravitational potentials on their way from the primordial plasma to here. These variations are due to primordial perturbations, probably caused by quantum fluctuations in the very early universe. The ongoing Planck satellite mission to observe the cosmic microwave background is also described.
Towards a UV Completion for Chameleon Scalar Theories
Chameleons are scalar fields that couple directly to ordinary matter with gravitational strength, but which nevertheless evade the stringent constraints on tests of gravity because of properties they acquire in the presence of high ambient matter density. Chameleon theories were originally constructed in a bottom-up, phenomenological fashion, with potentials and matter couplings designed to hide the scalar from experiments. In this paper, we attempt to embed the chameleon scenario within string compactifications, thus UV completing the scenario. We look for stabilized potentials that can realize a screening mechanism, and we find that the volume modulus rather generically works as a chameleon, and in fact the supersymmetric potential used by Kachru, Kallosh, Linde and Trivedi (KKLT) is an example of this type. We consider all constraints from tests of gravity, allowing us to put experimental constraints on the KKLT parameters.
Theory of planet formation
We review the current theoretical understanding how growth from micro-meter sized dust to massive giant planets occurs in disks around young stars. After introducing a number of observational constraints from the solar system, from observed protoplanetary disks, and from the extrasolar planets, we simplify the problem by dividing it into a number of discrete stages which are assumed to occur in a sequential way. In the first stage - the growth from dust to kilometer sized planetesimals - the aerodynamics of the bodies are of central importance. We discuss both a purely coagulative growth mode, as well as a gravoturbulent mode involving a gravitational instability of the dust. In the next stage, planetesimals grow to protoplanets of roughly 1000 km in size. Gravity is now the dominant force. The mass accretion can be strongly non-linear, leading to the detachment of a few big bodies from the remaining planetesimals. In the outer planetary system (outside a few AU), some of these bodies can become so massive that they eventually accrete a large gaseous envelope. This is the stage of giant planet formation, as understood within the core accretion-gas capture paradigm. We also discuss the direct gravitational collapse model where giant planets are thought to form directly via a gravitational fragmentation of the gas disk. In the inner system, protoplanets collide in the last stage - probably after the dispersal of the gaseous disk - in giant impacts until the separations between the remaining terrestrial planets become large enough to allow long term stability. We finish the review with some selected questions.
A much lower density for the transiting extrasolar planet WASP-7
John Southworth, M. Dominik, U. G. Jorgensen, S. Rahvar, C. Snodgrass, K. Alsubai, V. Bozza, P. Browne, M. Burgdorf, S. Calchi Novati, P. Dodds, S. Dreizler, F. Finet, T. Gerner, S. Hardis, K. Harpsoe,C. Hellier, T. C. Hinse, M. Hundertmark, N. Kains, E. Kerins, C. Liebig, L. Mancini, M. Mathiasen, M. T. Penny, S. Proft, D. Ricci, K. Sahu, G. Scarpetta, S. Schafer, F. Schonebeck, J. Surdej
We present the first high-precision photometry of the transiting extrasolar planetary system WASP-7, obtained using telescope defocussing techniques and reaching a scatter of 0.68 mmag per point. We find that the transit depth is greater and that the host star is more evolved than previously thought. The planet has a significantly larger radius (1.330 +/- 0.093 Rjup versus 0.915 +0.046 -0.040 Rjup) and much lower density (0.41 +/- 0.10 rhojup versus 1.26 +0.25 -0.21 rhojup) and surface gravity (13.4 +/- 2.6 m/s2 versus 26.4 +4.4 -4.0 m/s2) than previous measurements showed. Based on the revised properties it is no longer an outlier in planetary mass--radius and period--gravity diagrams. We also obtain a more precise transit ephemeris for the WASP-7 system.
No comments:
Post a Comment