Theoretical accuracy in cosmological growth estimation

NASA Astrophysics Data System (ADS)

Bose, Benjamin; Koyama, Kazuya; Hellwing, Wojciech A.; Zhao, Gong-Bo; Winther, Hans A.

2017-07-01

We elucidate the importance of the consistent treatment of gravity-model specific nonlinearities when estimating the growth of cosmological structures from redshift space distortions (RSD). Within the context of standard perturbation theory (SPT), we compare the predictions of two theoretical templates with redshift space data from COLA (comoving Lagrangian acceleration) simulations in the normal branch of DGP gravity (nDGP) and general relativity (GR). Using COLA for these comparisons is validated using a suite of full N-body simulations for the same theories. The two theoretical templates correspond to the standard general relativistic perturbation equations and those same equations modeled within nDGP. Gravitational clustering nonlinear effects are accounted for by modeling the power spectrum up to one-loop order and redshift space clustering anisotropy is modeled using the Taruya, Nishimichi and Saito (TNS) RSD model. Using this approach, we attempt to recover the simulation's fiducial logarithmic growth parameter f . By assigning the simulation data with errors representing an idealized survey with a volume of 10 Gpc3/h3 , we find the GR template is unable to recover fiducial f to within 1 σ at z =1 when we match the data up to kmax=0.195 h /Mpc . On the other hand, the DGP template recovers the fiducial value within 1 σ . Further, we conduct the same analysis for sets of mock data generated for generalized models of modified gravity using SPT, where again we analyze the GR template's ability to recover the fiducial value. We find that for models with enhanced gravitational nonlinearity, the theoretical bias of the GR template becomes significant for stage IV surveys. Thus, we show that for the future large data volume galaxy surveys, the self-consistent modeling of non-GR gravity scenarios will be crucial in constraining theory parameters.

The real and apparent convergence of N-body simulations of the dark matter structures: Is the Navarro-Frenk-White profile real?

NASA Astrophysics Data System (ADS)

Baushev, A. N.

2015-03-01

While N-body simulations suggest a cuspy profile in the centra of the dark matter halos of galaxies, the majority of astronomical observations favor a relatively soft cored density distribution of these regions. The routine method of testing the convergence of N-body simulations (in particular, the negligibility of two-body scattering effect) is to find the conditions under which formed structures is insensitive to numerical parameters. The results obtained with this approach suggest a surprisingly minor role of the particle collisions: the central density profile remains untouched and close to the Navarro-Frenk-White shape, even if the simulation time significantly exceeds the collisional relaxation time τr . In order to check the influence of the unphysical test body collisions we use the Fokker-Planck equation. It turns out that a profile ρ ∝r-β where β ≃ 1 is an attractor: the Fokker-Planck diffusion transforms any reasonable initial distribution into it in a time shorter than τr , and then the cuspy profile should survive much longer than τr , since the Fokker-Planck diffusion is self-compensated if β ≃ 1 . Thus the purely numerical effect of test body scattering may create a stable NFW-like pseudosolution. Moreover, its stability may be mistaken for the simulation convergence. We present analytical estimations for this potential bias effect and call for numerical tests. For that purpose, we suggest a simple test that can be performed as the simulation progresses and would indicate the magnitude of the collisional influence and the veracity of the simulation results.

GRAPE-5: A Special-Purpose Computer for N-Body Simulations

NASA Astrophysics Data System (ADS)

Kawai, Atsushi; Fukushige, Toshiyuki; Makino, Junichiro; Taiji, Makoto

2000-08-01

We have developed a special-purpose computer for gravitational many-body simulations, GRAPE-5. GRAPE-5 accelerates the force calculation which dominates the calculation cost of the simulation. All other calculations, such as the time integration of orbits, are performed on a general-purpose computer (host computer) connected to GRAPE-5. A GRAPE-5 board consists of eight custom pipeline chips (G5 chip) and its peak performance is 38.4 Gflops. GRAPE-5 is the successor of GRAPE-3. The differences between GRAPE-5 and GRAPE-3 are: (1) The newly developed G5 chip contains two pipelines operating at 80 MHz, while the GRAPE chip, which was used for GRAPE-3, had one at 20 MHz. The calculation speed of GRAPE-5 is 8-times faster than that of GRAPE-3. (2) The GRAPE-5 board adopted a PCI bus as the interface to the host computer instead of VME of GRAPE-3, resulting in a communication speed one order of magnitude faster. (3) In addition to the pure 1/r potential, the G5 chip can calculate forces with arbitrary cutoff functions, so that it can be applied to the Ewald or P3M methods. (4) The pairwise force calculated on GRAPE-5 is about 10-times more accurate than that on GRAPE-3. On one GRAPE-5 board, one timestep with a direct summation algorithm takes 14 (N/128 k)2 seconds. With the Barnes-Hut tree algorithm (theta = 0.75), one timestep can be done in 15 (N/106) seconds.

A small-scale dynamo in feedback-dominated galaxies - III. Cosmological simulations

NASA Astrophysics Data System (ADS)

Rieder, Michael; Teyssier, Romain

2017-12-01

Magnetic fields are widely observed in the Universe in virtually all astrophysical objects, from individual stars to entire galaxies, even in the intergalactic medium, but their specific genesis has long been debated. Due to the development of more realistic models of galaxy formation, viable scenarios are emerging to explain cosmic magnetism, thanks to both deeper observations and more efficient and accurate computer simulations. We present here a new cosmological high-resolution zoom-in magnetohydrodynamic (MHD) simulation, using the adaptive mesh refinement technique, of a dwarf galaxy with an initially weak and uniform magnetic seed field that is amplified by a small-scale dynamo (SSD) driven by supernova-induced turbulence. As first structures form from the gravitational collapse of small density fluctuations, the frozen-in magnetic field separates from the cosmic expansion and grows through compression. In a second step, star formation sets in and establishes a strong galactic fountain, self-regulated by supernova explosions. Inside the galaxy, the interstellar medium becomes highly turbulent, dominated by strong supersonic shocks, as demonstrated by the spectral analysis of the gas kinetic energy. In this turbulent environment, the magnetic field is quickly amplified via a SSD process and is finally carried out into the circumgalactic medium by a galactic wind. This realistic cosmological simulation explains how initially weak magnetic seed fields can be amplified quickly in early, feedback-dominated galaxies, and predicts, as a consequence of the SSD process, that high-redshift magnetic fields are likely to be dominated by their small-scale components.

On the effect of galactic outflows in cosmological simulations of disc galaxies

NASA Astrophysics Data System (ADS)

Valentini, Milena; Murante, Giuseppe; Borgani, Stefano; Monaco, Pierluigi; Bressan, Alessandro; Beck, Alexander M.

2017-09-01

We investigate the impact of galactic outflow modelling on the formation and evolution of a disc galaxy, by performing a suite of cosmological simulations with zoomed-in initial conditions (ICs) of a Milky Way-sized halo. We verify how sensitive the general properties of the simulated galaxy are to the way in which stellar feedback triggered outflows are implemented, keeping ICs, simulation code and star formation (SF) model all fixed. We present simulations that are based on a version of the gadget3 code where our sub-resolution model is coupled with an advanced implementation of smoothed particle hydrodynamics that ensures a more accurate fluid sampling and an improved description of gas mixing and hydrodynamical instabilities. We quantify the strong interplay between the adopted hydrodynamic scheme and the sub-resolution model describing SF and feedback. We consider four different galactic outflow models, including the one introduced by Dalla Vecchia & Schaye (2012) and a scheme that is inspired by the Springel & Hernquist (2003) model. We find that the sub-resolution prescriptions adopted to generate galactic outflows are the main shaping factor of the stellar disc component at low redshift. The key requirement that a feedback model must have to be successful in producing a disc-dominated galaxy is the ability to regulate the high-redshift SF (responsible for the formation of the bulge component), the cosmological infall of gas from the large-scale environment, and gas fall-back within the galactic radius at low redshift, in order to avoid a too high SF rate at z = 0.

The diversity of evolutionary pathways of compact elliptical galaxies in cosmological simulations

NASA Astrophysics Data System (ADS)

Wellons, Sarah

2017-01-01

Observations of the high-redshift universe have revealed a population of galaxies which are already very massive (~1e11 solar masses at z=2) and have typical sizes of < 2 kpc, much smaller than their counterparts in the local universe. How such dense, massive galaxies form, and why they appear to be less common at low redshift, have been questions of interest for both theorists and observers. I will discuss these questions in the context of the Illustris simulation, a hydrodynamical cosmological simulation in which tens of thousands of galaxies form, evolve, and interact with each other, situated within a cosmological context. I select a group of massive compact galaxies at z=2 in the simulation and trace them back and forth in time to discover both how they formed at high redshift, and what they evolve into at the present day. I find a variety of both progenitors (compact galaxies form in the simulation either via central starbursts generally brought on by mergers, or by racing out to the tip of the SF main sequence and forming very early) and descendants (many formerly-compact galaxies lurk at the core of a more massive galaxy today, others were consumed in mergers, and some evolve passively and undisturbed). I will also discuss the implications of these results for observational methods of connecting galaxy populations across redshifts - in particular, the assumption of a constant cumulative comoving number density - and suggest an improvement to this method which takes the complexity and variety of galaxies' evolutionary paths into account.

An Investigation of Intracluster Light Evolution Using Cosmological Hydrodynamical Simulations

NASA Astrophysics Data System (ADS)

Tang, Lin; Lin, Weipeng; Cui, Weiguang; Kang, Xi; Wang, Yang; Contini, E.; Yu, Yu

2018-06-01

Intracluster light (ICL) in observations is usually identified through the surface brightness limit (SBL) method. In this paper, for the first time we produce mock images of galaxy groups and clusters, using a cosmological hydrodynamical simulation to investigate the ICL fraction and focus on its dependence on observational parameters, e.g., the SBL, the effects of cosmological redshift-dimming, point-spread function (PSF), and CCD pixel size. Detailed analyses suggest that the width of the PSF has a significant effect on the measured ICL fraction, while the relatively small pixel size shows almost no influence. It is found that the measured ICL fraction depends strongly on the SBL. At a fixed SBL and redshift, the measured ICL fraction decreases with increasing halo mass, while with a much fainter SBL, it does not depend on halo mass at low redshifts. In our work, the measured ICL fraction shows a clear dependence on the cosmological redshift-dimming effect. It is found that there is more mass locked in the ICL component than light, suggesting that the use of a constant mass-to-light ratio at high surface brightness levels will lead to an underestimate of ICL mass. Furthermore, it is found that the radial profile of ICL shows a characteristic radius that is almost independent of halo mass. The current measurement of ICL from observations has a large dispersion due to different methods, and we emphasize the importance of using the same definition when observational results are compared with theoretical predictions.

Counts of galaxy clusters as cosmological probes: the impact of baryonic physics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Balaguera-Antolínez, Andrés; Porciani, Cristiano, E-mail: abalan@astro.uni-bonn.de, E-mail: porciani@astro.uni-bonn.de

2013-04-01

The halo mass function from N-body simulations of collisionless matter is generally used to retrieve cosmological parameters from observed counts of galaxy clusters. This neglects the observational fact that the baryonic mass fraction in clusters is a random variable that, on average, increases with the total mass (within an overdensity of 500). Considering a mock catalog that includes tens of thousands of galaxy clusters, as expected from the forthcoming generation of surveys, we show that the effect of a varying baryonic mass fraction will be observable with high statistical significance. The net effect is a change in the overall normalizationmore » of the cluster mass function and a milder modification of its shape. Our results indicate the necessity of taking into account baryonic corrections to the mass function if one wants to obtain unbiased estimates of the cosmological parameters from data of this quality. We introduce the formalism necessary to accomplish this goal. Our discussion is based on the conditional probability of finding a given value of the baryonic mass fraction for clusters of fixed total mass. Finally, we show that combining information from the cluster counts with measurements of the baryonic mass fraction in a small subsample of clusters (including only a few tens of objects) will nearly optimally constrain the cosmological parameters.« less

High-precision spectra for dynamical Dark Energy cosmologies from constant-w models

NASA Astrophysics Data System (ADS)

Casarini, Luciano

2010-08-01

Spanning the whole functional space of cosmologies with any admissible DE state equations w(a) seems a need, in view of forthcoming observations, namely those aiming to provide a tomography of cosmic shear. In this paper I show that this duty can be eased and that a suitable use of results for constant-w cosmologies can be sufficient. More in detail, I ``assign'' here six cosmologies, aiming to span the space of state equations w(a) = wo+wa(1-a), for wo and wa values consistent with WMAP5 and WMAP7 releases and run N-body simulations to work out their non-linear fluctuation spectra at various redshifts z. Such spectra are then compared with those of suitable auxiliary models, characterized by constant w. For each z a different auxiliary model is needed. Spectral discrepancies between the assigned and the auxiliary models, up to k simeq 2-3 h Mpc-1, are shown to keep within 1 %. Quite in general, discrepancies are smaller at greater z and exhibit a specific trend across the wo and wa plane. Besides of aiming at simplifying the evaluation of spectra for a wide range of models, this paper also outlines a specific danger for future studies of the DE state equation, as models fairly distant on the w0-wa plane can be easily confused.

Simulating the Growth of a Disk Galaxy and its Supermassive Black Hole in a Cosmological Simulating the Growth of a Disk Galaxy and its Supermassive Black Hole in a Cosmological Context

DOE Office of Scientific and Technical Information (OSTI.GOV)

Levine, Robyn Deborah

2008-01-01

Supermassive black holes (SMBHs) are ubiquitous in the centers of galaxies. Their formation and subsequent evolution is inextricably linked to that of their host galaxies, and the study of galaxy formation is incomplete without the inclusion of SMBHs. The present work seeks to understand the growth and evolution of SMBHs through their interaction with the host galaxy and its environment. In the first part of the thesis (Chap. 2 and 3), we combine a simple semi-analytic model of outflows from active galactic nuclei (AGN) with a simulated dark matter density distribution to study the impact of SMBH feedback on cosmologicalmore » scales. We find that constraints can be placed on the kinetic efficiency of such feedback using observations of the filling fraction of the Lyα forest. We also find that AGN feedback is energetic enough to redistribute baryons over cosmological distances, having potentially significant effects on the interpretation of cosmological data which are sensitive to the total matter density distribution (e.g. weak lensing). However, truly assessing the impact of AGN feedback in the universe necessitates large-dynamic range simulations with extensive treatment of baryonic physics to first model the fueling of SMBHs. In the second part of the thesis (Chap. 4-6) we use a hydrodynamic adaptive mesh refinement simulation to follow the growth and evolution of a typical disk galaxy hosting a SMBH, in a cosmological context. The simulation covers a dynamical range of 10 million allowing us to study the transport of matter and angular momentum from super-galactic scales all the way down to the outer edge of the accretion disk around the SMBH. Focusing our attention on the central few hundred parsecs of the galaxy, we find the presence of a cold, self-gravitating, molecular gas disk which is globally unstable. The global instabilities drive super-sonic turbulence, which maintains local stability and allows gas to fuel a SMBH without first fragmenting

N-body simulations for f(R) gravity using a self-adaptive particle-mesh code

NASA Astrophysics Data System (ADS)

Zhao, Gong-Bo; Li, Baojiu; Koyama, Kazuya

2011-02-01

We perform high-resolution N-body simulations for f(R) gravity based on a self-adaptive particle-mesh code MLAPM. The chameleon mechanism that recovers general relativity on small scales is fully taken into account by self-consistently solving the nonlinear equation for the scalar field. We independently confirm the previous simulation results, including the matter power spectrum, halo mass function, and density profiles, obtained by Oyaizu [Phys. Rev. DPRVDAQ1550-7998 78, 123524 (2008)10.1103/PhysRevD.78.123524] and Schmidt [Phys. Rev. DPRVDAQ1550-7998 79, 083518 (2009)10.1103/PhysRevD.79.083518], and extend the resolution up to k˜20h/Mpc for the measurement of the matter power spectrum. Based on our simulation results, we discuss how the chameleon mechanism affects the clustering of dark matter and halos on full nonlinear scales.

Antimatter cosmic rays from dark matter annihilation: First results from an N-body experiment

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lavalle, J.; Nezri, E.; Athanassoula, E.

2008-11-15

While the particle hypothesis for dark matter may be very soon investigated at the LHC, and as the PAMELA and GLAST satellites are currently taking new data on charged and gamma cosmic rays, the need of controlling the theoretical uncertainties affecting the possible indirect signatures of dark matter annihilation is of paramount importance. The uncertainties which originate from the dark matter distribution are difficult to estimate because current astrophysical observations provide rather weak dynamical constraints and because, according to cosmological N-body simulations, dark matter is neither smoothly nor spherically distributed in galactic halos. Some previous studies made use of N-bodymore » simulations to compute the {gamma}-ray flux from dark matter annihilation, but such a work has never been performed for the antimatter (positron and antiproton) primary fluxes, for which transport processes complicate the calculations. We take advantage of the galaxylike 3D dark matter map extracted from the Horizon Project results to calculate the positron and antiproton fluxes from dark matter annihilation, in a model-independent approach as well as for dark matter particle benchmarks relevant at the LHC scale (from supersymmetric and extradimensional theories). We find that the flux uncertainties arise mainly from fluctuations of the local dark matter density, and are of {approx}1 order of magnitude. We compare our results to analytic descriptions of the dark matter halo, showing how the latter can well reproduce the former. The overall antimatter predictions associated with our benchmark models are shown to lie far below the existing measurements and, in particular, that of the positron fraction recently reported by PAMELA, and far below the background predictions as well. Finally, we stress the limits of the use of an N-body framework in this context.« less

Simulations of dusty plasmas using a special-purpose computer system designed for gravitational N-body problems

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yamamoto, K.; Mizuno, Y.; Hibino, S.

2006-01-15

Simulations of dusty plasmas were performed using GRAPE-6, a special-purpose computer designed for gravitational N-body problems. The collective behavior of dust particles, which are injected into the plasma, was studied by means of three-dimensional computer simulations. As an example of a dusty plasma simulation, experiments on Coulomb crystals in plasmas are simulated. Formation of a quasi-two-dimensional Coulomb crystal has been observed under typical laboratory conditions. Another example was to simulate movement of dust particles in plasmas under microgravity conditions. Fully three-dimensional spherical structures of dust clouds have been observed. For the simulation of a dusty plasma in microgravity with 3x10{supmore » 4} particles, GRAPE-6 can perform the whole operation 1000 times faster than by using a Pentium 4 1.6 GHz processor.« less

Introduction to Cosmology, Proceedings of the Polish Astronomical Society volume 4

NASA Astrophysics Data System (ADS)

Biernacka, Monika; Bajan, Katarzyna; Stachowski, Grzegorz; Pollo, Agnieszka

2016-07-01

On 11-23 July 2016, Jan Kochanowski University in Kielce was the host of the Second Cosmological School "Introduction to Cosmology". The main purpose of the School was to provide an introduction to a selection of the most interesting topics in modern cosmology, both in theory and observations. The program included a series of mini-workshops on cosmological simulations, Virtual Observatory database and tools and Spectral Energy Distribution tting. The School was intended for undergraduate, MSc and PhD students, as well as young postdoctoral researchers. The School was co-organized by the Polish Astronomical Society, the Jan Kochanowski University in Kielce, the Jagiellonian University in Cracow, the Nuclear Centre for Nuclear Research and the N. Copernicus Astronomical Center in Warsaw. The Interdisciplinary Centre for Mathematical and Computational Modeling kindly provided us with the possibility to remotely use their computing facilities.

Introduction to Cosmology, Proceedings of the Polish Astronomical Society volume 4

NASA Astrophysics Data System (ADS)

Biernacka, Monika; Bajan, Katarzyna; Stachowski, Grzegorz; Pollo, Agnieszka

2017-08-01

On 11-23 July 2016, Jan Kochanowski University in Kielce was the host of the Second Cosmological School "Introduction to Cosmology". The main purpose of the School was to provide an introduction to a selection of the most interesting topics in modern cosmology, both in theory and observations. The program included a series of mini-workshops on cosmological simulations, Virtual Observatory database and tools and Spectral Energy Distribution tting. The School was intended for undergraduate, MSc and PhD students, as well as young postdoctoral researchers. The School was co-organized by the Polish Astronomical Society, the Jan Kochanowski University in Kielce, the Jagiellonian University in Cracow, the Nuclear Centre for Nuclear Research and the N. Copernicus Astronomical Center in Warsaw. The Interdisciplinary Centre for Mathematical and Computational Modeling kindly provided us with the possibility to remotely use their computing facilities.

Dwarf galaxy mass estimators versus cosmological simulations

NASA Astrophysics Data System (ADS)

González-Samaniego, Alejandro; Bullock, James S.; Boylan-Kolchin, Michael; Fitts, Alex; Elbert, Oliver D.; Hopkins, Philip F.; Kereš, Dušan; Faucher-Giguère, Claude-André

2017-12-01

We use a suite of high-resolution cosmological dwarf galaxy simulations to test the accuracy of commonly used mass estimators from Walker et al. (2009) and Wolf et al. (2010), both of which depend on the observed line-of-sight velocity dispersion and the 2D half-light radius of the galaxy, Re. The simulations are part of the Feedback in Realistic Environments (FIRE) project and include 12 systems with stellar masses spanning 105-107 M⊙ that have structural and kinematic properties similar to those of observed dispersion-supported dwarfs. Both estimators are found to be quite accurate: M_Wolf/M_true = 0.98^{+0.19}_{-0.12} and M_Walker/M_true =1.07^{+0.21}_{-0.15}, with errors reflecting the 68 per cent range over all simulations. The excellent performance of these estimators is remarkable given that they each assume spherical symmetry, a supposition that is broken in our simulated galaxies. Though our dwarfs have negligible rotation support, their 3D stellar distributions are flattened, with short-to-long axis ratios c/a ≃ 0.4-0.7. The median accuracy of the estimators shows no trend with asphericity. Our simulated galaxies have sphericalized stellar profiles in 3D that follow a nearly universal form, one that transitions from a core at small radius to a steep fall-off ∝r-4.2 at large r; they are well fit by Sérsic profiles in projection. We find that the most important empirical quantity affecting mass estimator accuracy is Re. Determining Re by an analytic fit to the surface density profile produces a better estimated mass than if the half-light radius is determined via direct summation.

An N-body Integrator for Planetary Rings

NASA Astrophysics Data System (ADS)

Hahn, Joseph M.

2011-04-01

A planetary ring that is disturbed by a satellite's resonant perturbation can respond in an organized way. When the resonance lies in the ring's interior, the ring responds via an m-armed spiral wave, while a ring whose edge is confined by the resonance exhibits an m-lobed scalloping along the ring-edge. The amplitude of these disturbances are sensitive to ring surface density and viscosity, so modelling these phenomena can provide estimates of the ring's properties. However a brute force attempt to simulate a ring's full azimuthal extent with an N-body code will likely fail because of the large number of particles needed to resolve the ring's behavior. Another impediment is the gravitational stirring that occurs among the simulated particles, which can wash out the ring's organized response. However it is possible to adapt an N-body integrator so that it can simulate a ring's collective response to resonant perturbations. The code developed here uses a few thousand massless particles to trace streamlines within the ring. Particles are close in a radial sense to these streamlines, which allows streamlines to be treated as straight wires of constant linear density. Consequently, gravity due to these streamline is a simple function of the particle's radial distance to all streamlines. And because particles are responding to smooth gravitating streamlines, rather than discrete particles, this method eliminates the stirring that ordinarily occurs in brute force N-body calculations. Note also that ring surface density is now a simple function of streamline separations, so effects due to ring pressure and viscosity are easily accounted for, too. A poster will describe this N-body method in greater detail. Simulations of spiral density waves and scalloped ring-edges are executed in typically ten minutes on a desktop PC, and results for Saturn's A and B rings will be presented at conference time.

Resolving the Small-Scale Structure of the Circumgalactic Medium in Cosmological Simulations

NASA Astrophysics Data System (ADS)

Corlies, Lauren

2017-08-01

We propose to resolve the circumgalactic medium (CGM) of L* galaxies down to 100 Msun (250 pc) in a full cosmological simulation to examine how mixing and cooling shape the physical nature of this gas on the scales expected from observations. COS has provided the best characterization of the low-z CGM to date, revealing the extent and amount of low- and high-ions and hinting at the kinematic relations between them. Yet cosmological galaxy simulations that can reproduce the stellar properties of galaxies have all struggled to reproduce these results even qualitatively. However, while the COS data imply that the low-ion absorption is occurring on sub-kpc scales, such scales can not be traced by simulations with resolution between 1-5 kpc in the CGM. Our proposed simulations will, for the first time, reach the resolution required to resolve these structures in the outer halo of L* galaxies. Using the adaptive mesh refinement code enzo, we will experiment with the size, shape, and resolution of an enforced high refinement region extending from the disk into the CGM to identify the best configuration for probing the flows of gas throughout the CGM. Our test case has found that increasing the resolution alone can have dramatic consequences for the density, temperature, and kinematics along a line of sight. Coupling this technique with an independent feedback study already underway will help disentangle the roles of global and small scale physics in setting the physical state of the CGM. Finally, we will use the MISTY pipeline to generate realistic mock spectra for direct comparison with COS data which will be made available through MAST.

SPH/N-body simulations of small (D = 10 km) monolithic asteroidal breakups and improved parametric relations for Monte-Carlo collisional models

NASA Astrophysics Data System (ADS)

Ševecek, Pavel; Broz, Miroslav; Nesvorny, David; Durda, Daniel D.; Asphaug, Erik; Walsh, Kevin J.; Richardson, Derek C.

2016-10-01

Detailed models of asteroid collisions can yield important constrains for the evolution of the Main Asteroid Belt, but the respective parameter space is large and often unexplored. We thus performed a new set of simulations of asteroidal breakups, i.e. fragmentations of intact targets, subsequent gravitational reaccumulation and formation of small asteroid families, focusing on parent bodies with diameters D = 10 km.Simulations were performed with a smoothed-particle hydrodynamics (SPH) code (Benz & Asphaug 1994), combined with an efficient N-body integrator (Richardson et al. 2000). We assumed a number of projectile sizes, impact velocities and impact angles. The rheology used in the physical model does not include friction nor crushing; this allows for a direct comparison to results of Durda et al. (2007). Resulting size-frequency distributions are significantly different from scaled-down simulations with D = 100 km monolithic targets, although they may be even more different for pre-shattered targets.We derive new parametric relations describing fragment distributions, suitable for Monte-Carlo collisional models. We also characterize velocity fields and angular distributions of fragments, which can be used as initial conditions in N-body simulations of small asteroid families. Finally, we discuss various uncertainties related to SPH simulations.

The large-scale environment from cosmological simulations - I. The baryonic cosmic web

NASA Astrophysics Data System (ADS)

Cui, Weiguang; Knebe, Alexander; Yepes, Gustavo; Yang, Xiaohu; Borgani, Stefano; Kang, Xi; Power, Chris; Staveley-Smith, Lister

2018-01-01

Using a series of cosmological simulations that includes one dark-matter-only (DM-only) run, one gas cooling-star formation-supernova feedback (CSF) run and one that additionally includes feedback from active galactic nuclei (AGNs), we classify the large-scale structures with both a velocity-shear-tensor code (VWEB) and a tidal-tensor code (PWEB). We find that the baryonic processes have almost no impact on large-scale structures - at least not when classified using aforementioned techniques. More importantly, our results confirm that the gas component alone can be used to infer the filamentary structure of the universe practically un-biased, which could be applied to cosmology constraints. In addition, the gas filaments are classified with its velocity (VWEB) and density (PWEB) fields, which can theoretically connect to the radio observations, such as H I surveys. This will help us to bias-freely link the radio observations with dark matter distributions at large scale.

The redshift distribution of cosmological samples: a forward modeling approach

NASA Astrophysics Data System (ADS)

Herbel, Jörg; Kacprzak, Tomasz; Amara, Adam; Refregier, Alexandre; Bruderer, Claudio; Nicola, Andrina

2017-08-01

Determining the redshift distribution n(z) of galaxy samples is essential for several cosmological probes including weak lensing. For imaging surveys, this is usually done using photometric redshifts estimated on an object-by-object basis. We present a new approach for directly measuring the global n(z) of cosmological galaxy samples, including uncertainties, using forward modeling. Our method relies on image simulations produced using \\textsc{UFig} (Ultra Fast Image Generator) and on ABC (Approximate Bayesian Computation) within the MCCL (Monte-Carlo Control Loops) framework. The galaxy population is modeled using parametric forms for the luminosity functions, spectral energy distributions, sizes and radial profiles of both blue and red galaxies. We apply exactly the same analysis to the real data and to the simulated images, which also include instrumental and observational effects. By adjusting the parameters of the simulations, we derive a set of acceptable models that are statistically consistent with the data. We then apply the same cuts to the simulations that were used to construct the target galaxy sample in the real data. The redshifts of the galaxies in the resulting simulated samples yield a set of n(z) distributions for the acceptable models. We demonstrate the method by determining n(z) for a cosmic shear like galaxy sample from the 4-band Subaru Suprime-Cam data in the COSMOS field. We also complement this imaging data with a spectroscopic calibration sample from the VVDS survey. We compare our resulting posterior n(z) distributions to the one derived from photometric redshifts estimated using 36 photometric bands in COSMOS and find good agreement. This offers good prospects for applying our approach to current and future large imaging surveys.

Forming Circumbinary Planets: N-body Simulations of Kepler-34

NASA Astrophysics Data System (ADS)

Lines, S.; Leinhardt, Z. M.; Paardekooper, S.; Baruteau, C.; Thebault, P.

2014-02-01

Observations of circumbinary planets orbiting very close to the central stars have shown that planet formation may occur in a very hostile environment, where the gravitational pull from the binary should be very strong on the primordial protoplanetary disk. Elevated impact velocities and orbit crossings from eccentricity oscillations are the primary contributors to high energy, potentially destructive collisions that inhibit the growth of aspiring planets. In this work, we conduct high-resolution, inter-particle gravity enabled N-body simulations to investigate the feasibility of planetesimal growth in the Kepler-34 system. We improve upon previous work by including planetesimal disk self-gravity and an extensive collision model to accurately handle inter-planetesimal interactions. We find that super-catastrophic erosion events are the dominant mechanism up to and including the orbital radius of Kepler-34(AB)b, making in situ growth unlikely. It is more plausible that Kepler-34(AB)b migrated from a region beyond 1.5 AU. Based on the conclusions that we have made for Kepler-34, it seems likely that all of the currently known circumbinary planets have also migrated significantly from their formation location with the possible exception of Kepler-47(AB)c.

Non-linear matter power spectrum covariance matrix errors and cosmological parameter uncertainties

NASA Astrophysics Data System (ADS)

Blot, L.; Corasaniti, P. S.; Amendola, L.; Kitching, T. D.

2016-06-01

The covariance of the matter power spectrum is a key element of the analysis of galaxy clustering data. Independent realizations of observational measurements can be used to sample the covariance, nevertheless statistical sampling errors will propagate into the cosmological parameter inference potentially limiting the capabilities of the upcoming generation of galaxy surveys. The impact of these errors as function of the number of realizations has been previously evaluated for Gaussian distributed data. However, non-linearities in the late-time clustering of matter cause departures from Gaussian statistics. Here, we address the impact of non-Gaussian errors on the sample covariance and precision matrix errors using a large ensemble of N-body simulations. In the range of modes where finite volume effects are negligible (0.1 ≲ k [h Mpc-1] ≲ 1.2), we find deviations of the variance of the sample covariance with respect to Gaussian predictions above ˜10 per cent at k > 0.3 h Mpc-1. Over the entire range these reduce to about ˜5 per cent for the precision matrix. Finally, we perform a Fisher analysis to estimate the effect of covariance errors on the cosmological parameter constraints. In particular, assuming Euclid-like survey characteristics we find that a number of independent realizations larger than 5000 is necessary to reduce the contribution of sampling errors to the cosmological parameter uncertainties at subpercent level. We also show that restricting the analysis to large scales k ≲ 0.2 h Mpc-1 results in a considerable loss in constraining power, while using the linear covariance to include smaller scales leads to an underestimation of the errors on the cosmological parameters.

An iterative reconstruction of cosmological initial density fields

NASA Astrophysics Data System (ADS)

Hada, Ryuichiro; Eisenstein, Daniel J.

2018-05-01

We present an iterative method to reconstruct the linear-theory initial conditions from the late-time cosmological matter density field, with the intent of improving the recovery of the cosmic distance scale from the baryon acoustic oscillations (BAOs). We present tests using the dark matter density field in both real and redshift space generated from an N-body simulation. In redshift space at z = 0.5, we find that the reconstructed displacement field using our iterative method are more than 80% correlated with the true displacement field of the dark matter particles on scales k < 0.10h Mpc-1. Furthermore, we show that the two-point correlation function of our reconstructed density field matches that of the initial density field substantially better, especially on small scales (<40h-1 Mpc). Our redshift-space results are improved if we use an anisotropic smoothing so as to account for the reduced small-scale information along the line of sight in redshift space.

Fast Multipole Methods for Three-Dimensional N-body Problems

NASA Technical Reports Server (NTRS)

Koumoutsakos, P.

1995-01-01

We are developing computational tools for the simulations of three-dimensional flows past bodies undergoing arbitrary motions. High resolution viscous vortex methods have been developed that allow for extended simulations of two-dimensional configurations such as vortex generators. Our objective is to extend this methodology to three dimensions and develop a robust computational scheme for the simulation of such flows. A fundamental issue in the use of vortex methods is the ability of employing efficiently large numbers of computational elements to resolve the large range of scales that exist in complex flows. The traditional cost of the method scales as Omicron (N(sup 2)) as the N computational elements/particles induce velocities at each other, making the method unacceptable for simulations involving more than a few tens of thousands of particles. In the last decade fast methods have been developed that have operation counts of Omicron (N log N) or Omicron (N) (referred to as BH and GR respectively) depending on the details of the algorithm. These methods are based on the observation that the effect of a cluster of particles at a certain distance may be approximated by a finite series expansion. In order to exploit this observation we need to decompose the element population spatially into clusters of particles and build a hierarchy of clusters (a tree data structure) - smaller neighboring clusters combine to form a cluster of the next size up in the hierarchy and so on. This hierarchy of clusters allows one to determine efficiently when the approximation is valid. This algorithm is an N-body solver that appears in many fields of engineering and science. Some examples of its diverse use are in astrophysics, molecular dynamics, micro-magnetics, boundary element simulations of electromagnetic problems, and computer animation. More recently these N-body solvers have been implemented and applied in simulations involving vortex methods. Koumoutsakos and Leonard (1995

Constructing high-quality bounding volume hierarchies for N-body computation using the acceptance volume heuristic

NASA Astrophysics Data System (ADS)

Olsson, O.

2018-01-01

We present a novel heuristic derived from a probabilistic cost model for approximate N-body simulations. We show that this new heuristic can be used to guide tree construction towards higher quality trees with improved performance over current N-body codes. This represents an important step beyond the current practice of using spatial partitioning for N-body simulations, and enables adoption of a range of state-of-the-art algorithms developed for computer graphics applications to yield further improvements in N-body simulation performance. We outline directions for further developments and review the most promising such algorithms.

A Novel Approach to Visualizing Dark Matter Simulations.

PubMed

Kaehler, R; Hahn, O; Abel, T

2012-12-01

In the last decades cosmological N-body dark matter simulations have enabled ab initio studies of the formation of structure in the Universe. Gravity amplified small density fluctuations generated shortly after the Big Bang, leading to the formation of galaxies in the cosmic web. These calculations have led to a growing demand for methods to analyze time-dependent particle based simulations. Rendering methods for such N-body simulation data usually employ some kind of splatting approach via point based rendering primitives and approximate the spatial distributions of physical quantities using kernel interpolation techniques, common in SPH (Smoothed Particle Hydrodynamics)-codes. This paper proposes three GPU-assisted rendering approaches, based on a new, more accurate method to compute the physical densities of dark matter simulation data. It uses full phase-space information to generate a tetrahedral tessellation of the computational domain, with mesh vertices defined by the simulation's dark matter particle positions. Over time the mesh is deformed by gravitational forces, causing the tetrahedral cells to warp and overlap. The new methods are well suited to visualize the cosmic web. In particular they preserve caustics, regions of high density that emerge, when several streams of dark matter particles share the same location in space, indicating the formation of structures like sheets, filaments and halos. We demonstrate the superior image quality of the new approaches in a comparison with three standard rendering techniques for N-body simulation data.

Sixth- and eighth-order Hermite integrator for N-body simulations

NASA Astrophysics Data System (ADS)

Nitadori, Keigo; Makino, Junichiro

2008-10-01

We present sixth- and eighth-order Hermite integrators for astrophysical N-body simulations, which use the derivatives of accelerations up to second-order ( snap) and third-order ( crackle). These schemes do not require previous values for the corrector, and require only one previous value to construct the predictor. Thus, they are fairly easy to implement. The additional cost of the calculation of the higher-order derivatives is not very high. Even for the eighth-order scheme, the number of floating-point operations for force calculation is only about two times larger than that for traditional fourth-order Hermite scheme. The sixth-order scheme is better than the traditional fourth-order scheme for most cases. When the required accuracy is very high, the eighth-order one is the best. These high-order schemes have several practical advantages. For example, they allow a larger number of particles to be integrated in parallel than the fourth-order scheme does, resulting in higher execution efficiency in both general-purpose parallel computers and GRAPE systems.

The redshift distribution of cosmological samples: a forward modeling approach

DOE Office of Scientific and Technical Information (OSTI.GOV)

Herbel, Jörg; Kacprzak, Tomasz; Amara, Adam

Determining the redshift distribution n ( z ) of galaxy samples is essential for several cosmological probes including weak lensing. For imaging surveys, this is usually done using photometric redshifts estimated on an object-by-object basis. We present a new approach for directly measuring the global n ( z ) of cosmological galaxy samples, including uncertainties, using forward modeling. Our method relies on image simulations produced using \\textsc(UFig) (Ultra Fast Image Generator) and on ABC (Approximate Bayesian Computation) within the MCCL (Monte-Carlo Control Loops) framework. The galaxy population is modeled using parametric forms for the luminosity functions, spectral energy distributions, sizesmore » and radial profiles of both blue and red galaxies. We apply exactly the same analysis to the real data and to the simulated images, which also include instrumental and observational effects. By adjusting the parameters of the simulations, we derive a set of acceptable models that are statistically consistent with the data. We then apply the same cuts to the simulations that were used to construct the target galaxy sample in the real data. The redshifts of the galaxies in the resulting simulated samples yield a set of n ( z ) distributions for the acceptable models. We demonstrate the method by determining n ( z ) for a cosmic shear like galaxy sample from the 4-band Subaru Suprime-Cam data in the COSMOS field. We also complement this imaging data with a spectroscopic calibration sample from the VVDS survey. We compare our resulting posterior n ( z ) distributions to the one derived from photometric redshifts estimated using 36 photometric bands in COSMOS and find good agreement. This offers good prospects for applying our approach to current and future large imaging surveys.« less

OBERON: OBliquity and Energy balance Run on N-body systems

NASA Astrophysics Data System (ADS)

Forgan, Duncan H.

2016-08-01

OBERON (OBliquity and Energy balance Run on N-body systems) models the climate of Earthlike planets under the effects of an arbitrary number and arrangement of other bodies, such as stars, planets and moons. The code, written in C++, simultaneously computes N body motions using a 4th order Hermite integrator, simulates climates using a 1D latitudinal energy balance model, and evolves the orbital spin of bodies using the equations of Laskar (1986a,b).

Observing the baryon cycle in hydrodynamic cosmological simulations

NASA Astrophysics Data System (ADS)

Vander Vliet, Jacob Richard

An understanding of galaxy evolution requires an understanding of the flow of baryons in and out of a galaxy. The accretion of baryons is required for galaxies to form stars, while stars eject baryons out of the galaxy through stellar feedback mechanisms such as supernovae, stellar winds, and radiation pressure. The interplay between outfiowing and infalling material form the circumgalactic medium (CGM). Hydrodynamic simulations provide understanding of the connection between stellar feedback and the distribution and kinematics of baryons in the CGM. To compare simulations and observations properly the simulated CGI must be observed in the same manner as the real CGM. I have developed the Mockspec code to generate synthetic quasar absorption line observations of the CGM in cosmological hydrodynamic simulations. Mockspec generates synthetic spectra based on the phase; lnetallicity, and kinematics of CGM gas and mimics instrumental effects. Mockspec includes automated analysis of the spectra and identifies the gas responsible for the absorption. Mockspec was applied to simulations of dwarf galaxies at low redshift to examine the observable effect different feedback models have on the CGM. While the different feedback models had strong effects on the galaxy, they all produced a similar CGM that failed match observations. Mockspec was applied to the VELA simulation suite of high redshift, high mass galaxies to examine the variance of the CGM across different galaxies in different environments. The observed CGM showed little variation between the different galaxies and almost no evolution from z=4 to z=1. The VELAs were not able to generate a CGM to match the observations. The properties of cells responsible for the absorption were compared to the derived properties from Voigt Profile decomposition. VP modeling was found to accurately describe the HI and MgII absorbing gas but failed for high ionization species such as CIV and OVI, which do not arise in the coherent

Zonal methods for the parallel execution of range-limited N-body simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bowers, Kevin J.; Dror, Ron O.; Shaw, David E.

2007-01-20

Particle simulations in fields ranging from biochemistry to astrophysics require the evaluation of interactions between all pairs of particles separated by less than some fixed interaction radius. The applicability of such simulations is often limited by the time required for calculation, but the use of massive parallelism to accelerate these computations is typically limited by inter-processor communication requirements. Recently, Snir [M. Snir, A note on N-body computations with cutoffs, Theor. Comput. Syst. 37 (2004) 295-318] and Shaw [D.E. Shaw, A fast, scalable method for the parallel evaluation of distance-limited pairwise particle interactions, J. Comput. Chem. 26 (2005) 1318-1328] independently introducedmore » two distinct methods that offer asymptotic reductions in the amount of data transferred between processors. In the present paper, we show that these schemes represent special cases of a more general class of methods, and introduce several new algorithms in this class that offer practical advantages over all previously described methods for a wide range of problem parameters. We also show that several of these algorithms approach an approximate lower bound on inter-processor data transfer.« less

Phantom-GRAPE: Numerical software library to accelerate collisionless N-body simulation with SIMD instruction set on x86 architecture

NASA Astrophysics Data System (ADS)

Tanikawa, Ataru; Yoshikawa, Kohji; Nitadori, Keigo; Okamoto, Takashi

2013-02-01

We have developed a numerical software library for collisionless N-body simulations named "Phantom-GRAPE" which highly accelerates force calculations among particles by use of a new SIMD instruction set extension to the x86 architecture, Advanced Vector eXtensions (AVX), an enhanced version of the Streaming SIMD Extensions (SSE). In our library, not only the Newton's forces, but also central forces with an arbitrary shape f(r), which has a finite cutoff radius rcut (i.e. f(r)=0 at r>rcut), can be quickly computed. In computing such central forces with an arbitrary force shape f(r), we refer to a pre-calculated look-up table. We also present a new scheme to create the look-up table whose binning is optimal to keep good accuracy in computing forces and whose size is small enough to avoid cache misses. Using an Intel Core i7-2600 processor, we measure the performance of our library for both of the Newton's forces and the arbitrarily shaped central forces. In the case of Newton's forces, we achieve 2×109 interactions per second with one processor core (or 75 GFLOPS if we count 38 operations per interaction), which is 20 times higher than the performance of an implementation without any explicit use of SIMD instructions, and 2 times than that with the SSE instructions. With four processor cores, we obtain the performance of 8×109 interactions per second (or 300 GFLOPS). In the case of the arbitrarily shaped central forces, we can calculate 1×109 and 4×109 interactions per second with one and four processor cores, respectively. The performance with one processor core is 6 times and 2 times higher than those of the implementations without any use of SIMD instructions and with the SSE instructions. These performances depend only weakly on the number of particles, irrespective of the force shape. It is good contrast with the fact that the performance of force calculations accelerated by graphics processing units (GPUs) depends strongly on the number of particles

The origin of kinematically distinct cores and misaligned gas discs in galaxies from cosmological simulations

NASA Astrophysics Data System (ADS)

Taylor, Philip; Federrath, Christoph; Kobayashi, Chiaki

2018-06-01

Integral field spectroscopy surveys provide spatially resolved gas and stellar kinematics of galaxies. They have unveiled a range of atypical kinematic phenomena, which require detailed modelling to understand. We present results from a cosmological simulation that includes stellar and AGN feedback. We find that the distribution of angles between the gas and stellar angular momenta of galaxies is not affected by projection effects. We examine five galaxies (≈6 per cent of well resolved galaxies) that display atypical kinematics; two of the galaxies have kinematically distinct cores (KDC), while the other three have counter-rotating gas and stars. All five form the majority of their stars in the field, subsequently falling into cosmological filaments where the relative orientation of the stellar angular momentum and the bulk gas flow leads to the formation of a counter-rotating gas disc. The accreted gas exchanges angular momentum with pre-existing co-rotating gas causing it to fall to the centre of the galaxy. This triggers low-level AGN feedback, which reduces star formation. Later, two of the galaxies experience a minor merger (stellar mass ratio ˜1/10) with a galaxy on a retrograde orbit compared to the spin of the stellar component of the primary. This produces the KDCs, and is a different mechanism than suggested by other works. The role of minor mergers in the kinematic evolution of galaxies may have been under-appreciated in the past, and large, high-resolution cosmological simulations will be necessary to gain a better understanding in this area.

The Atacama Cosmology Telescope: Cosmology from Galaxy Clusters Detected Via the Sunyaev-Zel'dovich Effect

NASA Technical Reports Server (NTRS)

Sehgal, Neelima; Trac, Hy; Acquaviva, Viviana; Ade, Peter A. R.; Aguirre, Paula; Amiri, Mandana; Appel, John W.; Barrientos, L. Felipe; Battistelli, Elia S.; Bond, J. Richard;

The Formation And Evolution Of Milky Way Sized Galaxies In High-Resolution Cosmological Zoom Simulations

NASA Astrophysics Data System (ADS)

Grand, Robert

2016-09-01

Simulations are playing an increasingly important role in probing the formation history of the Milky Way, including the formation of the thick/thin disc and origin of the metal distribution and chemo-dynamical relations. We introduce the Auriga project, a suite of high resolution cosmological-zoom simulations of Milky Way-sized galaxies simulated with the state-of-the-art cosmological magneto-hydrodynamical code AREPO, and present an analysis of the formation and evolution of the stellar disc(s) from early times to present day. In particular, we show that 'thickened discs' are mainly driven by a bar (if present) and interactions with satellites of masses log10 (M/ Mo ) >= 10, whereas other potential heating mechanisms such as spiral arms, radial migration, and adiabatic heating from mid-plane density growth are all sub-dominant. Interestingly, we find that even in cases of violent satellite interactions the disc reforms quickly (within a few giga years), producing a well-defined disc-bulge system. In nearly all simulations the overall structure of the disc becomes gradually more radially extended and vertically thinner with time, in support of the inside-out, upside-down formation scenario, and without the presence of a thin/thick disc dichotomy. In addition, we comment on the mass distribution of mono-abundance populations and their relation to the bulge and disc components, which are readily comparable to observations from surveys such as APOGEE and Gaia.

Astronomy and Cosmology of the Guarani of Southern Brazil

NASA Astrophysics Data System (ADS)

de Mello, Flávia Cristina

The Guarani Indians in South Brazil have a sophisticated system of thought about the cosmos. Presented here are some elements of their cosmology and cosmogony and the influences of the heavenly bodies in this people's everyday life, which have been collected in ethnographic research in the first decade of the 21st century. The main themes of cosmology and the origin of the Sun, Moon, and Earth are described, approximating this anthropological research to ethnoastronomy and cultural anthropology discussions. This research seeks to analyze comparatively this cosmological concept with other indigenous cosmological systems and to compare them with the Western cosmological system, thus including it in studies of cultural astronomy.

SPMHD simulations of structure formation

NASA Astrophysics Data System (ADS)

Barnes, David J.; On, Alvina Y. L.; Wu, Kinwah; Kawata, Daisuke

2018-05-01

The intracluster medium of galaxy clusters is permeated by μ {G} magnetic fields. Observations with current and future facilities have the potential to illuminate the role of these magnetic fields play in the astrophysical processes of galaxy clusters. To obtain a greater understanding of how the initial seed fields evolve to the magnetic fields in the intracluster medium requires magnetohydrodynamic simulations. We critically assess the current smoothed particle magnetohydrodynamic (SPMHD) schemes, especially highlighting the impact of a hyperbolic divergence cleaning scheme and artificial resistivity switch on the magnetic field evolution in cosmological simulations of the formation of a galaxy cluster using the N-body/SPMHD code GCMHD++. The impact and performance of the cleaning scheme and two different schemes for the artificial resistivity switch is demonstrated via idealized test cases and cosmological simulations. We demonstrate that the hyperbolic divergence cleaning scheme is effective at suppressing the growth of the numerical divergence error of the magnetic field and should be applied to any SPMHD simulation. Although the artificial resistivity is important in the strong field regime, it can suppress the growth of the magnetic field in the weak field regime, such as galaxy clusters. With sufficient resolution, simulations with divergence cleaning can reproduce observed magnetic fields. We conclude that the cleaning scheme alone is sufficient for galaxy cluster simulations, but our results indicate that the SPMHD scheme must be carefully chosen depending on the regime of the magnetic field.

Numerical Convergence in the Dark Matter Halos Properties Using Cosmological Simulations

NASA Astrophysics Data System (ADS)

Mosquera-Escobar, X. E.; Muñoz-Cuartas, J. C.

2017-07-01

Nowadays, the accepted cosmological model is the so called -Cold Dark Matter (CDM). In such model, the universe is considered to be homogeneous and isotropic, composed of diverse components as the dark matter and dark energy, where the latter is the most abundant one. Dark matter plays an important role because it is responsible for the generation of gravitational potential wells, commonly called dark matter halos. At the end, dark matter halos are characterized by a set of parameters (mass, radius, concentration, spin parameter), these parameters provide valuable information for different studies, such as galaxy formation, gravitational lensing, etc. In this work we use the publicly available code Gadget2 to perform cosmological simulations to find to what extent the numerical parameters of the simu- lations, such as gravitational softening, integration time step and force calculation accuracy affect the physical properties of the dark matter halos. We ran a suite of simulations where these parameters were varied in a systematic way in order to explore accurately their impact on the structural parameters of dark matter halos. We show that the variations on the numerical parameters affect the structural pa- rameters of dark matter halos, such as concentration, virial radius, and concentration. We show that these modifications emerged when structures become non- linear (at redshift 2) for the scale of our simulations, such that these variations affected the formation and evolution structure of halos mainly at later cosmic times. As a quantitative result, we propose which would be the most appropriate values for the numerical parameters of the simulations, such that they do not affect the halo properties that are formed. For force calculation accuracy we suggest values smaller or equal to 0.0001, integration time step smaller o equal to 0.005 and for gravitational softening we propose equal to 1/60th of the mean interparticle distance, these values, correspond to the

HNBody: A Simulation Package for Hierarchical N-Body Systems

NASA Astrophysics Data System (ADS)

Rauch, Kevin P.

2018-04-01

HNBody (http://www.hnbody.org/) is an extensible software package forintegrating the dynamics of N-body systems. Although general purpose, itincorporates several features and algorithms particularly well-suited tosystems containing a hierarchy (wide dynamic range) of masses. HNBodyversion 1 focused heavily on symplectic integration of nearly-Kepleriansystems. Here I describe the capabilities of the redesigned and expandedpackage version 2, which includes: symplectic integrators up to eighth order(both leap frog and Wisdom-Holman type methods), with symplectic corrector andclose encounter support; variable-order, variable-timestep Bulirsch-Stoer andStörmer integrators; post-Newtonian and multipole physics options; advancedround-off control for improved long-term stability; multi-threading and SIMDvectorization enhancements; seamless availability of extended precisionarithmetic for all calculations; extremely flexible configuration andoutput. Tests of the physical correctness of the algorithms are presentedusing JPL Horizons ephemerides (https://ssd.jpl.nasa.gov/?horizons) andpreviously published results for reference. The features and performanceof HNBody are also compared to several other freely available N-body codes,including MERCURY (Chambers), SWIFT (Levison & Duncan) and WHFAST (Rein &Tamayo).

Supercomputer simulations of structure formation in the Universe

NASA Astrophysics Data System (ADS)

Ishiyama, Tomoaki

2017-06-01

We describe the implementation and performance results of our massively parallel MPI†/OpenMP‡ hybrid TreePM code for large-scale cosmological N-body simulations. For domain decomposition, a recursive multi-section algorithm is used and the size of domains are automatically set so that the total calculation time is the same for all processes. We developed a highly-tuned gravity kernel for short-range forces, and a novel communication algorithm for long-range forces. For two trillion particles benchmark simulation, the average performance on the fullsystem of K computer (82,944 nodes, the total number of core is 663,552) is 5.8 Pflops, which corresponds to 55% of the peak speed.

Cosmology Constraints from the Weak Lensing Peak Counts and the Power Spectrum in CFHTLenS

DOE PAGES

Liu, Jia; May, Morgan; Petri, Andrea; ...

2015-03-04

Lensing peaks have been proposed as a useful statistic, containing cosmological information from non-Gaussianities that is inaccessible from traditional two-point statistics such as the power spectrum or two-point correlation functions. Here we examine constraints on cosmological parameters from weak lensing peak counts, using the publicly available data from the 154 deg2 CFHTLenS survey. We utilize a new suite of ray-tracing N-body simulations on a grid of 91 cosmological models, covering broad ranges of the three parameters Ω m, σ 8, and w, and replicating the galaxy sky positions, redshifts, and shape noise in the CFHTLenS observations. We then build anmore » emulator that interpolates the power spectrum and the peak counts to an accuracy of ≤ 5%, and compute the likelihood in the three-dimensional parameter space (Ω m, σ 8, w) from both observables. We find that constraints from peak counts are comparable to those from the power spectrum, and somewhat tighter when different smoothing scales are combined. Neither observable can constrain w without external data. When the power spectrum and peak counts are combined, the area of the error “banana” in the (Ω m, σ 8) plane reduces by a factor of ≈ two, compared to using the power spectrum alone. For a flat Λ cold dark matter model, combining both statistics, we obtain the constraint σ 8(Ω m/0.27)0.63 = 0.85 +0.03 -0.03.« less

Cosmology Constraints from the Weak Lensing Peak Counts and the Power Spectrum in CFHTLenS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Liu, Jia; May, Morgan; Petri, Andrea

Lensing peaks have been proposed as a useful statistic, containing cosmological information from non-Gaussianities that is inaccessible from traditional two-point statistics such as the power spectrum or two-point correlation functions. Here we examine constraints on cosmological parameters from weak lensing peak counts, using the publicly available data from the 154 deg2 CFHTLenS survey. We utilize a new suite of ray-tracing N-body simulations on a grid of 91 cosmological models, covering broad ranges of the three parameters Ω m, σ 8, and w, and replicating the galaxy sky positions, redshifts, and shape noise in the CFHTLenS observations. We then build anmore » emulator that interpolates the power spectrum and the peak counts to an accuracy of ≤ 5%, and compute the likelihood in the three-dimensional parameter space (Ω m, σ 8, w) from both observables. We find that constraints from peak counts are comparable to those from the power spectrum, and somewhat tighter when different smoothing scales are combined. Neither observable can constrain w without external data. When the power spectrum and peak counts are combined, the area of the error “banana” in the (Ω m, σ 8) plane reduces by a factor of ≈ two, compared to using the power spectrum alone. For a flat Λ cold dark matter model, combining both statistics, we obtain the constraint σ 8(Ω m/0.27)0.63 = 0.85 +0.03 -0.03.« less

BOOSTED TIDAL DISRUPTION BY MASSIVE BLACK HOLE BINARIES DURING GALAXY MERGERS FROM THE VIEW OF N -BODY SIMULATION

DOE Office of Scientific and Technical Information (OSTI.GOV)

Li, Shuo; Berczik, Peter; Spurzem, Rainer

Supermassive black hole binaries (SMBHBs) are productions of the hierarchical galaxy formation model. There are many close connections between a central SMBH and its host galaxy because the former plays very important roles on galaxy formation and evolution. For this reason, the evolution of SMBHBs in merging galaxies is a fundamental challenge. Since there are many discussions about SMBHB evolution in a gas-rich environment, we focus on the quiescent galaxy, using tidal disruption (TD) as a diagnostic tool. Our study is based on a series of numerical, large particle number, direct N -body simulations for dry major mergers. According tomore » the simulation results, the evolution can be divided into three phases. In phase I, the TD rate for two well separated SMBHs in a merging system is similar to that for a single SMBH in an isolated galaxy. After two SMBHs approach close enough to form a bound binary in phase II, the disruption rate can be enhanced by ∼2 orders of magnitude within a short time. This “boosted” disruption stage finishes after the SMBHB evolves to a compact binary system in phase III, corresponding to a reduction in disruption rate back to a level of a few times higher than in phase I. We also discuss how to correctly extrapolate our N -body simulation results to reality, and the implications of our results to observations.« less

Cosmological constraints from the CFHTLenS shear measurements using a new, accurate, and flexible way of predicting non-linear mass clustering

NASA Astrophysics Data System (ADS)

Angulo, Raul E.; Hilbert, Stefan

2015-03-01

We explore the cosmological constraints from cosmic shear using a new way of modelling the non-linear matter correlation functions. The new formalism extends the method of Angulo & White, which manipulates outputs of N-body simulations to represent the 3D non-linear mass distribution in different cosmological scenarios. We show that predictions from our approach for shear two-point correlations at 1-300 arcmin separations are accurate at the ˜10 per cent level, even for extreme changes in cosmology. For moderate changes, with target cosmologies similar to that preferred by analyses of recent Planck data, the accuracy is close to ˜5 per cent. We combine this approach with a Monte Carlo Markov chain sampler to explore constraints on a Λ cold dark matter model from the shear correlation functions measured in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). We obtain constraints on the parameter combination σ8(Ωm/0.27)0.6 = 0.801 ± 0.028. Combined with results from cosmic microwave background data, we obtain marginalized constraints on σ8 = 0.81 ± 0.01 and Ωm = 0.29 ± 0.01. These results are statistically compatible with previous analyses, which supports the validity of our approach. We discuss the advantages of our method and the potential it offers, including a path to model in detail (i) the effects of baryons, (ii) high-order shear correlation functions, and (iii) galaxy-galaxy lensing, among others, in future high-precision cosmological analyses.

N-body simulations of gravitational redshifts and other relativistic distortions of galaxy clustering

NASA Astrophysics Data System (ADS)

Zhu, Hongyu; Alam, Shadab; Croft, Rupert A. C.; Ho, Shirley; Giusarma, Elena

2017-10-01

Large redshift surveys of galaxies and clusters are providing the first opportunities to search for distortions in the observed pattern of large-scale structure due to such effects as gravitational redshift. We focus on non-linear scales and apply a quasi-Newtonian approach using N-body simulations to predict the small asymmetries in the cross-correlation function of two galaxy different populations. Following recent work by Bonvin et al., Zhao and Peacock and Kaiser on galaxy clusters, we include effects which enter at the same order as gravitational redshift: the transverse Doppler effect, light-cone effects, relativistic beaming, luminosity distance perturbation and wide-angle effects. We find that all these effects cause asymmetries in the cross-correlation functions. Quantifying these asymmetries, we find that the total effect is dominated by the gravitational redshift and luminosity distance perturbation at small and large scales, respectively. By adding additional subresolution modelling of galaxy structure to the large-scale structure information, we find that the signal is significantly increased, indicating that structure on the smallest scales is important and should be included. We report on comparison of our simulation results with measurements from the SDSS/BOSS galaxy redshift survey in a companion paper.

Hydrodynamic Simulation of the Cosmological X-Ray Background

NASA Astrophysics Data System (ADS)

Croft, Rupert A. C.; Di Matteo, Tiziana; Davé, Romeel; Hernquist, Lars; Katz, Neal; Fardal, Mark A.; Weinberg, David H.

2001-08-01

We use a hydrodynamic simulation of an inflationary cold dark matter model with a cosmological constant to predict properties of the extragalactic X-ray background (XRB). We focus on emission from the intergalactic medium (IGM), with particular attention to diffuse emission from warm-hot gas that lies in relatively smooth filamentary structures between galaxies and galaxy clusters. We also include X-rays from point sources associated with galaxies in the simulation, and we make maps of the angular distribution of the emission. Although much of the X-ray luminous gas has a filamentary structure, the filaments are not evident in the simulated maps because of projection effects. In the soft (0.5-2 keV) band, our calculated mean intensity of radiation from intergalactic and cluster gas is 2.3×10-12 ergs-1 cm-2 deg-2, 35% of the total softband emission. This intensity is compatible at the ~1 σ level with estimates of the unresolved soft background intensity from deep ROSAT and Chandra measurements. Only 4% of the hard (2-10 keV) emission is associated with intergalactic gas. Relative to active galactic nuclei flux, the IGM component of the XRB peaks at a lower redshift (median z~0.45) and spans a narrower redshift range, so its clustering makes an important contribution to the angular correlation function of the total emission. The clustering on the scales accessible to our simulation (0.1‧-10') is significant, with an amplitude roughly consistent with an extrapolation of recent ROSAT results to small scales. A cross-correlation analysis of the XRB against nearby galaxies taken from a simulated redshift survey also yields a strong signal from the IGM. Our conclusions about the soft background intensity differ from those of some recent papers that have argued that the expected emission from gas in galaxy, group, and cluster halos would exceed the observed background unless much of the gas is expelled by supernova feedback. We obtain reasonable compatibility with

Self-consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations

DOE PAGES

Oñorbe, Jose; Hennawi, Joseph F.; Lukić, Zarija

2017-03-08

The ultraviolet background (UVB) emitted by quasars and galaxies governs the ionization and thermal state of the intergalactic medium (IGM), regulates the formation of high-redshift galaxies, and is thus a key quantity for modeling cosmic reionization. The vast majority of cosmological hydrodynamical simulations implement the UVB via a set of spatially uniform photoionization and photoheating rates derived from UVB synthesis models. In this paper, we show that simulations using canonical UVB rates reionize and, perhaps more importantly, spuriously heat the IGM, much earlier (more » $$z\\sim 15$$) than they should. This problem arises because at $$z\\gt 6$$, where observational constraints are nonexistent, the UVB amplitude is far too high. We introduce a new methodology to remedy this issue, and we generate self-consistent photoionization and photoheating rates to model any chosen reionization history. Following this approach, we run a suite of hydrodynamical simulations of different reionization scenarios and explore the impact of the timing of reionization and its concomitant heat injection on the thermal state of the IGM. We present a comprehensive study of the pressure smoothing scale of IGM gas, illustrating its dependence on the details of both hydrogen and helium reionization, and argue that it plays a fundamental role in interpreting Lyα forest statistics and the thermal evolution of the IGM. The premature IGM heating we have uncovered implies that previous work has likely dramatically overestimated the impact of photoionization feedback on galaxy formation, which sets the minimum halo mass able to form stars at high redshifts. Finally, we make our new UVB photoionization and photoheating rates publicly available for use in future simulations.« less

Self-consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Oñorbe, Jose; Hennawi, Joseph F.; Lukić, Zarija

The ultraviolet background (UVB) emitted by quasars and galaxies governs the ionization and thermal state of the intergalactic medium (IGM), regulates the formation of high-redshift galaxies, and is thus a key quantity for modeling cosmic reionization. The vast majority of cosmological hydrodynamical simulations implement the UVB via a set of spatially uniform photoionization and photoheating rates derived from UVB synthesis models. In this paper, we show that simulations using canonical UVB rates reionize and, perhaps more importantly, spuriously heat the IGM, much earlier (more » $$z\\sim 15$$) than they should. This problem arises because at $$z\\gt 6$$, where observational constraints are nonexistent, the UVB amplitude is far too high. We introduce a new methodology to remedy this issue, and we generate self-consistent photoionization and photoheating rates to model any chosen reionization history. Following this approach, we run a suite of hydrodynamical simulations of different reionization scenarios and explore the impact of the timing of reionization and its concomitant heat injection on the thermal state of the IGM. We present a comprehensive study of the pressure smoothing scale of IGM gas, illustrating its dependence on the details of both hydrogen and helium reionization, and argue that it plays a fundamental role in interpreting Lyα forest statistics and the thermal evolution of the IGM. The premature IGM heating we have uncovered implies that previous work has likely dramatically overestimated the impact of photoionization feedback on galaxy formation, which sets the minimum halo mass able to form stars at high redshifts. Finally, we make our new UVB photoionization and photoheating rates publicly available for use in future simulations.« less

N-body simulations of planet formation: understanding exoplanet system architectures

NASA Astrophysics Data System (ADS)

Coleman, Gavin; Nelson, Richard

2015-12-01

Observations have demonstrated the existence of a significant population of compact systems comprised of super-Earths and Neptune-mass planets, and a population of gas giants that appear to occur primarily in either short-period (<10 days) or longer period (>100 days) orbits. The broad diversity of system architectures raises the question of whether or not the same formation processes operating in standard disc models can explain these planets, or if different scenarios are required instead to explain the widely differing architectures. To explore this issue, we present the results from a comprehensive suite of N-body simulations of planetary system formation that include the following physical processes: gravitational interactions and collisions between planetary embryos and planetesimals; type I and II migration; gas accretion onto planetary cores; self-consistent viscous disc evolution and disc removal through photo-evaporation. Our results indicate that the formation and survival of compact systems of super-Earths and Neptune-mass planets occur commonly in disc models where a simple prescription for the disc viscosity is assumed, but such models never lead to the formation and survival of gas giant planets due to migration into the star. Inspired in part by the ALMA observations of HL Tau, and by MHD simulations that display the formation of long-lived zonal flows, we have explored the consequences of assuming that the disc viscosity varies in both time and space. We find that the radial structuring of the disc leads to conditions in which systems of giant planets are able to form and survive. Furthermore, these giants generally occupy those regions of the mass-period diagram that are densely populated by the observed gas giants, suggesting that the planet traps generated by radial structuring of protoplanetary discs may be a necessary ingredient for forming giant planets.

Ultralight Axion Dark Matter and Its Impact on Dark Halo Structure in N-body Simulations

NASA Astrophysics Data System (ADS)

Zhang, Jiajun; Sming Tsai, Yue-Lin; Kuo, Jui-Lin; Cheung, Kingman; Chu, Ming-Chung

2018-01-01

Ultralight axion is a dark matter candidate with mass { O }({10}-22){eV} and de Broglie wavelength of order kiloparsec. Such an axion, also called fuzzy dark matter (FDM), thermalizes via gravitational force and forms a Bose–Einstein condensate. Recent studies suggested that the quantum pressure from FDM can significantly affect structure formation in small scales, thus alleviating the so-called “small-scale crisis.” In this paper, we develop a new technique to discretize the quantum pressure and illustrate the interactions among FDM particles in an N-body simulation that accurately simulates the formation of the dark matter halo and its inner structure in the region outside the softening length. In a self-gravitationally bound virialized halo, we find a constant density solitonic core, which is consistent with theoretical prediction. The existence of the solitonic core reveals the nonlinear effect of quantum pressure and impacts structure formation in the FDM model.

Beta Dips in the Gaia Era: Simulation Predictions of the Galactic Velocity Anisotropy Parameter (β)

NASA Astrophysics Data System (ADS)

Loebman, Sarah; Valluri, Monica; Hattori, Kohei; Debattista, Victor P.; Bell, Eric F.; Stinson, Greg; Christensen, Charlotte; Brooks, Alyson; Quinn, Thomas R.; Governato, Fabio

2017-01-01

Milky Way (MW) science has entered a new era with the advent of Gaia. Combined with spectroscopic survey data, we have newfound access to full 6D phase space information for halo stars. Such data provides an invaluable opportunity to assess kinematic trends as a function of radius and confront simulations with these observations to draw insight about our merger history. I will discuss predictions for the velocity anisotropy parameter, β, drawn from three suites of state-of-the-art cosmological N-body and N-body+SPH MW-like simulations. On average, all three suites predict a monotonically increasing value of β that is radially biased, and beyond 10 kpc, β > 0.5. I will also discuss β as a function of time for individual simulated galaxies. I will highlight when "dips" in β form, the severity (the rarity of β < 0), origin (in situ versus accreted halo), and persistence of these dips. Thereby, I present a cohesive set of predictions of β from simulations for comparison to forthcoming observations.

Current Issues in Cosmology

NASA Astrophysics Data System (ADS)

Pecker, Jean-Claude; Narlikar, Jayant

2011-09-01

Part I. Observational Facts Relating to Discrete Sources: 1. The state of cosmology G. Burbidge; 2. The redshifts of galaxies and QSOs E. M. Burbidge and G. Burbidge; 3. Accretion discs in quasars J. Sulentic; Part II. Observational Facts Relating to Background Radiation: 4. CMB observations and consequences F. Bouchet; 5. Abundances of light nuclei K. Olive; 6. Evidence for an accelerating universe or lack of A. Blanchard; Part III. Standard Cosmology: 7. Cosmology, an overview of the standard model F. Bernardeau; 8. What are the building blocks of our universe? K. C. Wali; Part IV. Large-Scale Structure: 9. Observations of large-scale structure V. de Lapparent; 10. Reconstruction of large-scale peculiar velocity fields R. Mohayaee, B. Tully and U. Frisch; Part V. Alternative Cosmologies: 11. The quasi-steady state cosmology J. V. Narlikar; 12. Evidence for iron whiskers in the universe N. C. Wickramasinghe; 13. Alternatives to dark matter: MOND + Mach D. Roscoe; 14. Anthropic principle in cosmology B. Carter; Part VI. Evidence for Anomalous Redshifts: 15. Anomalous redshifts H. C. Arp; 16. Redshifts of galaxies and QSOs: the problem of redshift periodicities G. Burbidge; 17. Statistics of redshift periodicities W. Napier; 18. Local abnormal redshifts J.-C. Pecker; 19. Gravitational lensing and anomalous redshifts J. Surdej, J.-F. Claeskens and D. Sluse; Panel discussion; General discussion; Concluding remarks.

Simple inflationary models in Gauss-Bonnet brane-world cosmology

NASA Astrophysics Data System (ADS)

Okada, Nobuchika; Okada, Satomi

2016-06-01

In light of the recent Planck 2015 results for the measurement of the cosmic microwave background (CMB) anisotropy, we study simple inflationary models in the context of the Gauss-Bonnet (GB) brane-world cosmology. The brane-world cosmological effect modifies the power spectra of scalar and tensor perturbations generated by inflation and causes a dramatic change for the inflationary predictions of the spectral index (n s) and the tensor-to-scalar ratio (r) from those obtained in the standard cosmology. In particular, the predicted r values in the inflationary models favored by the Planck 2015 results are suppressed due to the GB brane-world cosmological effect, which is in sharp contrast with inflationary scenario in the Randall-Sundrum brane-world cosmology, where the r values are enhanced. Hence, these two brane-world cosmological scenarios are distinguishable. With the dramatic change of the inflationary predictions, the inflationary scenario in the GB brane-world cosmology can be tested by more precise measurements of n s and future observations of the CMB B-mode polarization.

Cosmological evolution of the nitrogen abundance

NASA Astrophysics Data System (ADS)

Vangioni, Elisabeth; Dvorkin, Irina; Olive, Keith A.; Dubois, Yohan; Molaro, Paolo; Petitjean, Patrick; Silk, Joe; Kimm, Taysun

2018-06-01

The abundance of nitrogen in the interstellar medium is a powerful probe of star formation processes over cosmological time-scales. Since nitrogen can be produced both in massive and intermediate-mass stars with metallicity-dependent yields, its evolution is challenging to model, as evidenced by the differences between theoretical predictions and observations. In this work, we attempt to identify the sources of these discrepancies using a cosmic evolution model. To further complicate matters, there is considerable dispersion in the abundances from observations of damped Lyα absorbers (DLAs) at z ˜ 2-3. We study the evolution of nitrogen with a detailed cosmic chemical evolution model and find good agreement with these observations, including the relative abundances of (N/O) and (N/Si). We find that the principal contribution of nitrogen comes from intermediate-mass stars, with the exception of systems with the lowest N/H, where nitrogen production might possibly be dominated by massive stars. This last result could be strengthened if stellar rotation which is important at low metallicity can produce significant amounts of nitrogen. Moreover, these systems likely reside in host galaxies with stellar masses below 108.5 M⊙. We also study the origin of the observed dispersion in nitrogen abundances using the cosmological hydrodynamical simulations Horizon-AGN. We conclude that this dispersion can originate from two effects: difference in the masses of the DLA host galaxies, and difference in their position inside the galaxy.

Properties of H I discs in the Auriga cosmological simulations

NASA Astrophysics Data System (ADS)

Marinacci, Federico; Grand, Robert J. J.; Pakmor, Rüdiger; Springel, Volker; Gómez, Facundo A.; Frenk, Carlos S.; White, Simon D. M.

2017-04-01

We analyse the properties of the H I gas distribution in the Auriga project, a set of magnetohydrodynamic cosmological simulations performed with the moving-mesh code arepo and a physics model for galaxy formation that succeeds in forming realistic late-type galaxies in the 30 Milky Way-sized haloes simulated in this project. We use a simple approach to estimate the neutral hydrogen fraction in our simulation set, which treats low-density and star-forming gas separately, and we explore two different prescriptions to subtract the contribution of molecular hydrogen from the total H I content. The H I gas in the vast majority of the systems forms extended discs although more disturbed morphologies are present. Notwithstanding the general good agreement with observed H I properties - such as radial profiles and the mass-diameter relation - the Auriga galaxies are systematically larger and more gas-rich than typical nearby galaxies. Interestingly, the amount of H I gas outside the disc plane correlates with the star formation rate, consistent with a picture where most of this extra-planar H I gas originates from a fountain-like flow. Our findings are robust with respect to the different assumptions adopted for computing the molecular hydrogen fraction and do not vary significantly over a wide range of numerical resolution. The H I modelling introduced in this paper can be used in future work to build artificial interferometric H I data cubes, allowing an even closer comparison of the gas dynamics in simulated galaxies with observations.

The impact of non-Gaussianity upon cosmological forecasts

NASA Astrophysics Data System (ADS)

Repp, A.; Szapudi, I.; Carron, J.; Wolk, M.

2015-12-01

The primary science driver for 3D galaxy surveys is their potential to constrain cosmological parameters. Forecasts of these surveys' effectiveness typically assume Gaussian statistics for the underlying matter density, despite the fact that the actual distribution is decidedly non-Gaussian. To quantify the effect of this assumption, we employ an analytic expression for the power spectrum covariance matrix to calculate the Fisher information for Baryon Acoustic Oscillation (BAO)-type model surveys. We find that for typical number densities, at kmax = 0.5h Mpc-1, Gaussian assumptions significantly overestimate the information on all parameters considered, in some cases by up to an order of magnitude. However, after marginalizing over a six-parameter set, the form of the covariance matrix (dictated by N-body simulations) causes the majority of the effect to shift to the `amplitude-like' parameters, leaving the others virtually unaffected. We find that Gaussian assumptions at such wavenumbers can underestimate the dark energy parameter errors by well over 50 per cent, producing dark energy figures of merit almost three times too large. Thus, for 3D galaxy surveys probing the non-linear regime, proper consideration of non-Gaussian effects is essential.

Novel kinematic methods to trace Spiral Arms nature using Gaia data

NASA Astrophysics Data System (ADS)

Roca-Fàbrega, S.; Figueras, F.; Valenzuela, O.; Romero-Gómez, M.; Antoja, T.; Colín, P.; Pichardo, B.; Velázquez, H.

2014-07-01

In this work we shed new light in the nature of spiral arm structures in galaxies. We present a disk kinematic and dynamic study of MW like galaxies using complementary approaches: analytical models, test-particle simulations, pure N-body and cosmological N-body plus hydrodynamic simulations. Using collisionless N-body data we have found that models with strong bar present a flat rotation frequency, i.e. rigid body rotation, whereas in the opposite extreme case, i.e. in unbarred systems, spiral arms are disk corotant (Roca-Fàbrega et al. 2013). Complementary to this work, we discuss how the vertex deviation parameter is a good tracer of corotation (CR) and outer Lindblad resonance radius (OLR) (Roca-Fàbrega et al. 2014). We have succeeded to produce MW like models in fully cosmological N-body plus hydrodynamic simulations with a high resolution (Roca-Fàbrega et al., in preparation). First results concerning disk phase space properties in terms of spiral arm nature using these simulations are presented (http://www.am.ub.edu/ sroca/shared/PosterRocaFabrega.pdf).

The metallicity distribution of H I systems in the EAGLE cosmological simulations

NASA Astrophysics Data System (ADS)

Rahmati, Alireza; Oppenheimer, Benjamin D.

2018-06-01

The metallicity of strong H I systems, spanning from damped Lyman α absorbers (DLAs) to Lyman-limit systems (LLSs), is explored between z = 5 → 0 using the EAGLE high-resolution cosmological hydrodynamic simulation of galaxy formation. The metallicities of LLSs and DLAs steadily increase with time in agreement with observations. DLAs are more metal rich than LLSs, although the metallicities in the LLS column density range (N_{H I }≈ 10^{17}-10^{20} cm^{-2}) are relatively flat, evolving from a median H I-weighted metallicity of {Z}≲ 10^{-2} Z_{⊙} at z = 3 to ≈10-0.5 Z⊙ by z = 0. The metal content of H I systems tracks the increasing stellar content of the Universe, holding ≈ 5 {per cent} of the integrated total metals released from stars at z = 0. We also consider partial LLS (pLLS, N_{H I}≈ 10^{16}-10^{17} cm^{-2}) metallicities, and find good agreement with Wotta et al. for the fraction of systems above (37 per cent) and below (63 per cent) 0.1 Z⊙. We also find a large dispersion of pLLS metallicities, although we do not reproduce the observed metallicity bimodality and instead we make the prediction that a larger sample will yield more pLLSs around 0.1 Z⊙. We underpredict the median metallicity of strong LLSs, and predict a population of Z < 10-3 Z⊙ DLAs at z > 3 that are not observed, which may indicate more widespread early enrichment in the real Universe compared to EAGLE.

N-body simulations of terrestrial planet formation under the influence of a hot Jupiter

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ogihara, Masahiro; Kobayashi, Hiroshi; Inutsuka, Shu-ichiro, E-mail: omasahiro@oca.eu, E-mail: ogihara@nagoya-u.jp

We investigate the formation of multiple-planet systems in the presence of a hot Jupiter (HJ) using extended N-body simulations that are performed simultaneously with semianalytic calculations. Our primary aims are to describe the planet formation process starting from planetesimals using high-resolution simulations, and to examine the dependences of the architecture of planetary systems on input parameters (e.g., disk mass, disk viscosity). We observe that protoplanets that arise from oligarchic growth and undergo type I migration stop migrating when they join a chain of resonant planets outside the orbit of an HJ. The formation of a resonant chain is almost independentmore » of our model parameters, and is thus a robust process. At the end of our simulations, several terrestrial planets remain at around 0.1 AU. The formed planets are not equal mass; the largest planet constitutes more than 50% of the total mass in the close-in region, which is also less dependent on parameters. In the previous work of this paper, we have found a new physical mechanism of induced migration of the HJ, which is called a crowding-out. If the HJ opens up a wide gap in the disk (e.g., owing to low disk viscosity), crowding-out becomes less efficient and the HJ remains. We also discuss angular momentum transfer between the planets and disk.« less

The Multi-SAG project: filling the MultiDark simulations with semi-analytic galaxies

NASA Astrophysics Data System (ADS)

Vega-Martínez, C. A.; Cora, S. A.; Padilla, N. D.; Muñoz Arancibia, A. M.; Orsi, A. A.; Ruiz, A. N.

2016-08-01

The semi-analytical model sag is a code of galaxy formation and evolution which is applied to halo catalogs and merger trees extracted from cosmological -body simulations of dark matter. This contribution describes the project of constructing a catalog of simulated galaxies by adapting and applying the model sag over two dark matter simulations of the spanish MultiDark Project publicly available. Those simulations have particles, each, in boxes with sizes of 1000 Mpc and 400 Mpc respectively with Planck cosmological parameters. They cover a large range of masses and have halo mass resolutions of , therefore each simulation is able to produce more than 150 millions of simulated galaxies. A detailed description of the method is explained, and the first statistical results are shown.

The metal enrichment of passive galaxies in cosmological simulations of galaxy formation

NASA Astrophysics Data System (ADS)

Okamoto, Takashi; Nagashima, Masahiro; Lacey, Cedric G.; Frenk, Carlos S.

2017-02-01

Massive early-type galaxies have higher metallicities and higher ratios of α elements to iron than their less massive counterparts. Reproducing these correlations has long been a problem for hierarchical galaxy formation theory, both in semi-analytic models and cosmological hydrodynamic simulations. We show that a simulation in which gas cooling in massive dark haloes is quenched by radio-mode active galactic nuclei (AGNs) feedback naturally reproduces the observed trend between α/Fe and the velocity dispersion of galaxies, σ. The quenching occurs earlier for more massive galaxies. Consequently, these galaxies complete their star formation before α/Fe is diluted by the contribution from Type Ia supernovae. For galaxies more massive than ˜1011 M⊙, whose α/Fe correlates positively with stellar mass, we find an inversely correlated mass-metallicity relation. This is a common problem in simulations in which star formation in massive galaxies is quenched either by quasar- or radio-mode AGN feedback. The early suppression of gas cooling in progenitors of massive galaxies prevents them from recapturing enriched gas ejected as winds. Simultaneously reproducing the [α/Fe]-σ relation and the mass-metallicity relation is, thus, difficult in the current framework of galaxy formation.

Investigating the physics and environment of Lyman limit systems in cosmological simulations

NASA Astrophysics Data System (ADS)

Erkal, Denis

2015-07-01

In this work, I investigate the properties of Lyman limit systems (LLSs) using state-of-the-art zoom-in cosmological galaxy formation simulations with on the fly radiative transfer, which includes both the cosmic UV background (UVB) and local stellar sources. I compare the simulation results to observations of the incidence frequency of LLSs and the H I column density distribution function over the redshift range z = 2-5 and find good agreement. I explore the connection between LLSs and their host haloes and find that LLSs reside in haloes with a wide range of halo masses with a nearly constant covering fraction within a virial radius. Over the range z = 2-5, I find that more than half of the LLSs reside in haloes with M < 1010 h-1 M⊙, indicating that absorption line studies of LLSs can probe these low-mass galaxies which H2-based star formation models predict to have very little star formation. I study the physical state of individual LLSs and test a simple model which encapsulates many of their properties. I confirm that LLSs have a characteristic absorption length given by the Jeans length and that they are in photoionization equilibrium at low column densities. Finally, I investigate the self-shielding of LLSs to the UVB and explore how the non-sphericity of LLSs affects the photoionization rate at a given N_{H I}. I find that at z ≈ 3, LLSs have an optical depth of unity at a column density of ˜1018 cm-2 and that this is the column density which characterizes the onset of self-shielding.

Cosmological hydrodynamical simulations of galaxy clusters: X-ray scaling relations and their evolution

NASA Astrophysics Data System (ADS)

Truong, N.; Rasia, E.; Mazzotta, P.; Planelles, S.; Biffi, V.; Fabjan, D.; Beck, A. M.; Borgani, S.; Dolag, K.; Gaspari, M.; Granato, G. L.; Murante, G.; Ragone-Figueroa, C.; Steinborn, L. K.

2018-03-01

We analyse cosmological hydrodynamical simulations of galaxy clusters to study the X-ray scaling relations between total masses and observable quantities such as X-ray luminosity, gas mass, X-ray temperature, and YX. Three sets of simulations are performed with an improved version of the smoothed particle hydrodynamics GADGET-3 code. These consider the following: non-radiative gas, star formation and stellar feedback, and the addition of feedback by active galactic nuclei (AGN). We select clusters with M500 > 1014 M⊙E(z)-1, mimicking the typical selection of Sunyaev-Zeldovich samples. This permits to have a mass range large enough to enable robust fitting of the relations even at z ˜ 2. The results of the analysis show a general agreement with observations. The values of the slope of the mass-gas mass and mass-temperature relations at z = 2 are 10 per cent lower with respect to z = 0 due to the applied mass selection, in the former case, and to the effect of early merger in the latter. We investigate the impact of the slope variation on the study of the evolution of the normalization. We conclude that cosmological studies through scaling relations should be limited to the redshift range z = 0-1, where we find that the slope, the scatter, and the covariance matrix of the relations are stable. The scaling between mass and YX is confirmed to be the most robust relation, being almost independent of the gas physics. At higher redshifts, the scaling relations are sensitive to the inclusion of AGNs which influences low-mass systems. The detailed study of these objects will be crucial to evaluate the AGN effect on the ICM.

BOOK REVIEW: Cosmology

NASA Astrophysics Data System (ADS)

Silk, Joseph

2008-11-01

The field of cosmology has been transformed since the glorious decades of the 1920's and 1930's when theory and observation converged to develop the current model of the expanding universe. It was a triumph of the theory of general relativity and astronomy. The first revolution came when the nuclear physicists entered the fray. This marked the debut of the hot big bang, in which the light elements were synthesized in the first three minutes. It was soon realised that elements like carbon and iron were synthesized in exploding stars. However helium, as well as deuterium and lithium, remain as George Gamow envisaged, the detritus of the big bang. The climax arrived with one of the most remarkable discoveries of the twentieth century, the cosmic microwave background radiation, in 1964. The fossil glow turned out to have the spectrum of an ideal black body. One could not imagine a stronger confirmation of the hot and dense origin of the universe. This discovery set the scene for the next major advance. It was now the turn of the particle physicists, who realized that the energies attained near the beginning of the universe, and unachievable in any conceivable terrestrial accelerator, provided a unique testing ground for theories of grand unification of the fundamental forces. This led Alan Guth and Andrei Linde in 1980 to propose the theory of inflation, which solved outstanding puzzles of the big bang. One could now understand why the universe is so large and homogeneous, and the origin of the seed fluctuations that gave rise to large-scale structure. A key prediction was that the universe should have Euclidean geometry, now verified to a precision of a few percent. Modern cosmology is firmly embedded in particle physics. It merits a text written by a particle physicist who can however appreciate the contributions of astronomy that provide the foundation and infrastructure for the theory of the expanding universe. There are now several such texts available. The most

Cosmological perturbation theory and the spherical collapse model - II. Non-Gaussian initial conditions

NASA Astrophysics Data System (ADS)

Gaztanaga, Enrique; Fosalba, Pablo

1998-12-01

In Paper I of this series, we introduced the spherical collapse (SC) approximation in Lagrangian space as a way of estimating the cumulants xi_J of density fluctuations in cosmological perturbation theory (PT). Within this approximation, the dynamics is decoupled from the statistics of the initial conditions, so we are able to present here the cumulants for generic non-Gaussian initial conditions, which can be estimated to arbitrary order including the smoothing effects. The SC model turns out to recover the exact leading-order non-linear contributions up to terms involving non-local integrals of the J-point functions. We argue that for the hierarchical ratios S_J, these non-local terms are subdominant and tend to compensate each other. The resulting predictions show a non-trivial time evolution that can be used to discriminate between models of structure formation. We compare these analytic results with non-Gaussian N-body simulations, which turn out to be in very good agreement up to scales where sigma<~1.

Precision Cosmology

NASA Astrophysics Data System (ADS)

Jones, Bernard J. T.

2017-04-01

Preface; Notation and conventions; Part I. 100 Years of Cosmology: 1. Emerging cosmology; 2. The cosmic expansion; 3. The cosmic microwave background; 4. Recent cosmology; Part II. Newtonian Cosmology: 5. Newtonian cosmology; 6. Dark energy cosmological models; 7. The early universe; 8. The inhomogeneous universe; 9. The inflationary universe; Part III. Relativistic Cosmology: 10. Minkowski space; 11. The energy momentum tensor; 12. General relativity; 13. Space-time geometry and calculus; 14. The Einstein field equations; 15. Solutions of the Einstein equations; 16. The Robertson-Walker solution; 17. Congruences, curvature and Raychaudhuri; 18. Observing and measuring the universe; Part IV. The Physics of Matter and Radiation: 19. Physics of the CMB radiation; 20. Recombination of the primeval plasma; 21. CMB polarisation; 22. CMB anisotropy; Part V. Precision Tools for Precision Cosmology: 23. Likelihood; 24. Frequentist hypothesis testing; 25. Statistical inference: Bayesian; 26. CMB data processing; 27. Parametrising the universe; 28. Precision cosmology; 29. Epilogue; Appendix A. SI, CGS and Planck units; Appendix B. Magnitudes and distances; Appendix C. Representing vectors and tensors; Appendix D. The electromagnetic field; Appendix E. Statistical distributions; Appendix F. Functions on a sphere; Appendix G. Acknowledgements; References; Index.

HBT+: an improved code for finding subhaloes and building merger trees in cosmological simulations

NASA Astrophysics Data System (ADS)

Han, Jiaxin; Cole, Shaun; Frenk, Carlos S.; Benitez-Llambay, Alejandro; Helly, John

2018-02-01

Dark matter subhalos are the remnants of (incomplete) halo mergers. Identifying them and establishing their evolutionary links in the form of merger trees is one of the most important applications of cosmological simulations. The HBT (Hierachical Bound-Tracing) code identifies haloes as they form and tracks their evolution as they merge, simultaneously detecting subhaloes and building their merger trees. Here we present a new implementation of this approach, HBT+ , that is much faster, more user friendly, and more physically complete than the original code. Applying HBT+ to cosmological simulations, we show that both the subhalo mass function and the peak-mass function are well fitted by similar double-Schechter functions. The ratio between the two is highest at the high-mass end, reflecting the resilience of massive subhaloes that experience substantial dynamical friction but limited tidal stripping. The radial distribution of the most-massive subhaloes is more concentrated than the universal radial distribution of lower mass subhaloes. Subhalo finders that work in configuration space tend to underestimate the masses of massive subhaloes, an effect that is stronger in the host centre. This may explain, at least in part, the excess of massive subhaloes in galaxy cluster centres inferred from recent lensing observations. We demonstrate that the peak-mass function is a powerful diagnostic of merger tree defects, and the merger trees constructed using HBT+ do not suffer from the missing or switched links that tend to afflict merger trees constructed from more conventional halo finders. We make the HBT+ code publicly available.

Low-redshift Lyman limit systems as diagnostics of cosmological inflows and outflows

NASA Astrophysics Data System (ADS)

Hafen, Zachary; Faucher-Giguère, Claude-André; Anglés-Alcázar, Daniel; Kereš, Dušan; Feldmann, Robert; Chan, T. K.; Quataert, Eliot; Murray, Norman; Hopkins, Philip F.

2017-08-01

We use cosmological hydrodynamic simulations with stellar feedback from the FIRE (Feedback In Realistic Environments) project to study the physical nature of Lyman limit systems (LLSs) at z ≤ 1. At these low redshifts, LLSs are closely associated with dense gas structures surrounding galaxies, such as galactic winds, dwarf satellites and cool inflows from the intergalactic medium. Our analysis is based on 14 zoom-in simulations covering the halo mass range Mh ≈ 109-1013 M⊙ at z = 0, which we convolve with the dark matter halo mass function to produce cosmological statistics. We find that the majority of cosmologically selected LLSs are associated with haloes in the mass range 1010 ≲ Mh ≲ 1012 M⊙. The incidence and H I column density distribution of simulated absorbers with columns in the range 10^{16.2} ≤ N_{H I} ≤ 2× 10^{20} cm-2 are consistent with observations. High-velocity outflows (with radial velocity exceeding the halo circular velocity by a factor of ≳ 2) tend to have higher metallicities ([X/H] ˜ -0.5) while very low metallicity ([X/H] < -2) LLSs are typically associated with gas infalling from the intergalactic medium. However, most LLSs occupy an intermediate region in metallicity-radial velocity space, for which there is no clear trend between metallicity and radial kinematics. The overall simulated LLS metallicity distribution has a mean (standard deviation) [X/H] = -0.9 (0.4) and does not show significant evidence for bimodality, in contrast to recent observational studies, but consistent with LLSs arising from haloes with a broad range of masses and metallicities.

JANUS: a bit-wise reversible integrator for N-body dynamics

NASA Astrophysics Data System (ADS)

Rein, Hanno; Tamayo, Daniel

2018-01-01

Hamiltonian systems such as the gravitational N-body problem have time-reversal symmetry. However, all numerical N-body integration schemes, including symplectic ones, respect this property only approximately. In this paper, we present the new N-body integrator JANUS , for which we achieve exact time-reversal symmetry by combining integer and floating point arithmetic. JANUS is explicit, formally symplectic and satisfies Liouville's theorem exactly. Its order is even and can be adjusted between two and ten. We discuss the implementation of JANUS and present tests of its accuracy and speed by performing and analysing long-term integrations of the Solar system. We show that JANUS is fast and accurate enough to tackle a broad class of dynamical problems. We also discuss the practical and philosophical implications of running exactly time-reversible simulations.

Phenomenology with fluctuating quantum geometries in loop quantum cosmology

NASA Astrophysics Data System (ADS)

Agullo, Ivan; Ashtekar, Abhay; Gupt, Brajesh

2017-04-01

The goal of this paper is to probe phenomenological implications of large fluctuations of quantum geometry in the Planck era, using cosmology of the early universe. For the background (Friedmann, Lemaître, Robertson, Walker) quantum geometry, we allow ‘widely spread’ states in which the relative dispersions are as large as 168 % in the Planck regime. By introducing suitable methods to overcome the ensuing conceptual and computational issues, we calculate the power spectrum {{P}R}(k) and the spectral index n s (k) of primordial curvature perturbations. These results generalize the previous work in loop quantum cosmology which focused on those states which were known to remain sharply peaked throughout the Planck regime. Surprisingly, even though the fluctuations we now consider are large, their presence does not add new features to the final {{P}R}(k) and n s (k): within observational error bars, their effect is degenerate with a different freedom in the theory, namely the number of pre-inflationary e-folds {{N}\\text{B\\star}} between the bounce and the onset of inflation. Therefore, with regard to observational consequences, one can simulate the freedom in the choice of states with large fluctuations in the Planck era using the simpler, sharply peaked states, simply by allowing for different values of {{N}\\text{B \\star}} .

Retrieving cosmological signal using cosmic flows

NASA Astrophysics Data System (ADS)

Bouillot, V.; Alimi, J.-M.

2011-12-01

To understand the origin of the anomalously high bulk flow at large scales, we use very large simulations in various cosmological models. To disentangle between cosmological and environmental effects, we select samples with bulk flow profiles similar to the observational data Watkins et al. (2009) which exhibit a maximum in the bulk flow at 53 h^{-1} Mpc. The estimation of the cosmological parameters Ω_M and σ_8, done on those samples, is correct from the rms mass fluctuation whereas this estimation gives completely false values when done on bulk flow measurements, hence showing a dependance of velocity fields on larger scales. By drawing a clear link between velocity fields at 53 h^{-1} Mpc and asymmetric patterns of the density field at 85 h^{-1} Mpc, we show that the bulk flow can depend largely on the environment. The retrieving of the cosmological signal is achieved by studying the convergence of the bulk flow towards the linear prediction at very large scale (˜ 150 h^{-1} Mpc).

Properties of Local Group galaxies in hydrodynamical simulations of sterile neutrino dark matter cosmologies

NASA Astrophysics Data System (ADS)

Lovell, Mark R.; Bose, Sownak; Boyarsky, Alexey; Crain, Robert A.; Frenk, Carlos S.; Hellwing, Wojciech A.; Ludlow, Aaron D.; Navarro, Julio F.; Ruchayskiy, Oleg; Sawala, Till; Schaller, Matthieu; Schaye, Joop; Theuns, Tom

2017-07-01

We study galaxy formation in sterile neutrino dark matter models that differ significantly from both cold and from 'warm thermal relic' models. We use the eagle code to carry out hydrodynamic simulations of the evolution of pairs of galaxies chosen to resemble the Local Group, as part of the APOSTLE simulations project. We compare cold dark matter (CDM) with two sterile neutrino models with 7 keV mass: one, the warmest among all models of this mass (LA120) and the other, a relatively cold case (LA10). We show that the lower concentration of sterile neutrino subhaloes compared to their CDM counterparts makes the inferred inner dark matter content of galaxies like Fornax (or Magellanic Clouds) less of an outlier in the sterile neutrino cosmologies. In terms of the galaxy number counts, the LA10 simulations are indistinguishable from CDM when one takes into account halo-to-halo (or 'simulation-to-simulation') scatter. In order for the LA120 model to match the number of Local Group dwarf galaxies, a higher fraction of low-mass haloes is required to form galaxies than is predicted by the eagle simulations. As the census of the Local Group galaxies nears completion, this population may provide a strong discriminant between cold and warm dark matter models.

COSMOG: Cosmology Oriented Sub-mm Modeling of Galactic Foregrounds

NASA Technical Reports Server (NTRS)

Kashlinsky, A.; Leisawitz, D.

2004-01-01

With upcoming missions in mid- and far-Infrared there is a need for software packages to reliably simulate the planned observations. This would help in both planning the observation and scanning strategy and in developing the concepts of the far-off missions. As this workshop demonstrated, many of the new missions are to be in the far-IR range of the electromagnetic spectrum and at the same time will map the sky with a sub-arcsec angular resolution. We present here a computer package for simulating foreground maps for the planned sub-mm and far-IR missions. such as SPECS. The package allows to study confusion limits and simulate cosmological observations for specified sky location interactively and in real time. Most of the emission at wavelengths long-ward of approximately 50 microns is dominated by Galactic cirrus and Zodiacal dust emission. Stellar emission at these wavelengths is weak and is for now neglected. Cosmological sources (distant and not-so-distant) galaxies for specified cosmologies will be added. Briefly, the steps that the algorithm goes through is described.

N-BODY SIMULATION OF PLANETESIMAL FORMATION THROUGH GRAVITATIONAL INSTABILITY AND COAGULATION. II. ACCRETION MODEL

DOE Office of Scientific and Technical Information (OSTI.GOV)

Michikoshi, Shugo; Kokubo, Eiichiro; Inutsuka, Shu-ichiro, E-mail: michikoshi@cfca.j, E-mail: kokubo@th.nao.ac.j, E-mail: inutsuka@tap.scphys.kyoto-u.ac.j

2009-10-01

The gravitational instability of a dust layer is one of the scenarios for planetesimal formation. If the density of a dust layer becomes sufficiently high as a result of the sedimentation of dust grains toward the midplane of a protoplanetary disk, the layer becomes gravitationally unstable and spontaneously fragments into planetesimals. Using a shearing box method, we performed local N-body simulations of gravitational instability of a dust layer and subsequent coagulation without gas and investigated the basic formation process of planetesimals. In this paper, we adopted the accretion model as a collision model. A gravitationally bound pair of particles ismore » replaced by a single particle with the total mass of the pair. This accretion model enables us to perform long-term and large-scale calculations. We confirmed that the formation process of planetesimals is the same as that in the previous paper with the rubble pile models. The formation process is divided into three stages: the formation of nonaxisymmetric structures; the creation of planetesimal seeds; and their collisional growth. We investigated the dependence of the planetesimal mass on the simulation domain size. We found that the mean mass of planetesimals formed in simulations is proportional to L {sup 3/2} {sub y}, where L{sub y} is the size of the computational domain in the direction of rotation. However, the mean mass of planetesimals is independent of L{sub x} , where L{sub x} is the size of the computational domain in the radial direction if L{sub x} is sufficiently large. We presented the estimation formula of the planetesimal mass taking into account the simulation domain size.« less

Cosmological Constraints from Fourier Phase Statistics

NASA Astrophysics Data System (ADS)

Ali, Kamran; Obreschkow, Danail; Howlett, Cullan; Bonvin, Camille; Llinares, Claudio; Oliveira Franco, Felipe; Power, Chris

2018-06-01

Most statistical inference from cosmic large-scale structure relies on two-point statistics, i.e. on the galaxy-galaxy correlation function (2PCF) or the power spectrum. These statistics capture the full information encoded in the Fourier amplitudes of the galaxy density field but do not describe the Fourier phases of the field. Here, we quantify the information contained in the line correlation function (LCF), a three-point Fourier phase correlation function. Using cosmological simulations, we estimate the Fisher information (at redshift z = 0) of the 2PCF, LCF and their combination, regarding the cosmological parameters of the standard ΛCDM model, as well as a Warm Dark Matter (WDM) model and the f(R) and Symmetron modified gravity models. The galaxy bias is accounted for at the level of a linear bias. The relative information of the 2PCF and the LCF depends on the survey volume, sampling density (shot noise) and the bias uncertainty. For a volume of 1h^{-3}Gpc^3, sampled with points of mean density \\bar{n} = 2× 10^{-3} h3 Mpc^{-3} and a bias uncertainty of 13%, the LCF improves the parameter constraints by about 20% in the ΛCDM cosmology and potentially even more in alternative models. Finally, since a linear bias only affects the Fourier amplitudes (2PCF), but not the phases (LCF), the combination of the 2PCF and the LCF can be used to break the degeneracy between the linear bias and σ8, present in 2-point statistics.

Cosmological Views of Anania Shirakatsi

NASA Astrophysics Data System (ADS)

Farmanyan, Sona V.; Mickaelian, Areg M.

2017-12-01

Since the ancient times the usage of cosmological ideas in mythology and poetry has contributed to the formation and development of human's philosophical thought. It is believed that before the M. Mashtots's alphabet, ancient Armenians have expressed their astronomical knowledge through stone structures and rock art.In the Armenian reality, the cosmological views, the idea of the spherical shape of the Earth and information of other celestial bodies more vividly were manifested in the works of Movses Khorenatsi, David Anhaght (5th century) and Anania Shirakatsi (7th century).Anania Shirakatsi is an Armenian Astronomer, Mathematician, Philosopher, Geographer and Alchemist.The importance of his work is also noted by foreign authors and he was called 7th century Cosmologist, First Scientist of Armenia and Middival Astronomer. Shirakatsi's works are united in his comprehensive knowledge, his insight of the mind, the ability of combining and analyzing facts and his literature talent.His works have simultaneous historical, cosmic, geographical, religious, literary and mystical significance. In the present study we will show Anania Shirakatsi's cosmological ideas and observations.

Cosmology with CLASS

NASA Astrophysics Data System (ADS)

Watts, Duncan; CLASS Collaboration

2018-01-01

The Cosmology Large Angular Scale Surveyor (CLASS) will use large-scale measurements of the polarized cosmic microwave background (CMB) to constrain the physics of inflation, reionization, and massive neutrinos. The experiment is designed to characterize the largest scales, which are inaccessible to most ground-based experiments, and remove Galactic foregrounds from the CMB maps. In this dissertation talk, I present simulations of CLASS data and demonstrate their ability to constrain the simplest single-field models of inflation and to reduce the uncertainty of the optical depth to reionization, τ, to near the cosmic variance limit, significantly improving on current constraints. These constraints will bring a qualitative shift in our understanding of standard ΛCDM cosmology. In particular, CLASS's measurement of τ breaks cosmological parameter degeneracies. Probes of large scale structure (LSS) test the effect of neutrino free-streaming at small scales, which depends on the mass of the neutrinos. CLASS's τ measurement, when combined with next-generation LSS and BAO measurements, will enable a 4σ detection of neutrino mass, compared with 2σ without CLASS data.. I will also briefly discuss the CLASS experiment's measurements of circular polarization of the CMB and the implications of the first-such near-all-sky map.

Cosmological Hydrodynamics on a Moving Mesh

NASA Astrophysics Data System (ADS)

Hernquist, Lars

We propose to construct a model for the visible Universe using cosmological simulations of structure formation. Our simulations will include both dark matter and baryons, and will employ two entirely different schemes for evolving the gas: smoothed particle hydrodynamics (SPH) and a moving mesh approach as incorporated in the new code, AREPO. By performing simulations that are otherwise nearly identical, except for the hydrodynamics solver, we will isolate and understand differences in the properties of galaxies, galaxy groups and clusters, and the intergalactic medium caused by the computational approach that have plagued efforts to understand galaxy formation for nearly two decades. By performing simulations at different levels of resolution and with increasingly complex treatments of the gas physics, we will identify the results that are converged numerically and that are robust with respect to variations in unresolved physical processes, especially those related to star formation, black hole growth, and related feedback effects. In this manner, we aim to undertake a research program that will redefine the state of the art in cosmological hydrodynamics and galaxy formation. In particular, we will focus our scientific efforts on understanding: 1) the formation of galactic disks in a cosmological context; 2) the physical state of diffuse gas in galaxy clusters and groups so that they can be used as high-precision probes of cosmology; 3) the nature of gas inflows into galaxy halos and the subsequent accretion of gas by forming disks; 4) the co-evolution of galaxies and galaxy clusters with their central supermassive black holes and the implications of related feedback for galaxy evolution and the dichotomy between blue and red galaxies; 5) the physical state of the intergalactic medium (IGM) and the evolution of the metallicity of the IGM; and 6) the reaction of dark matter around galaxies to galaxy formation. Our proposed work will be of immediate significance for

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Cosmological implications of the Fourier space wedges of the final sample

NASA Astrophysics Data System (ADS)

Grieb, Jan Niklas; Sánchez, Ariel G.; Salazar-Albornoz, Salvador; Scoccimarro, Román; Crocce, Martín; Dalla Vecchia, Claudio; Montesano, Francesco; Gil-Marín, Héctor; Ross, Ashley J.; Beutler, Florian; Rodríguez-Torres, Sergio; Chuang, Chia-Hsun; Prada, Francisco; Kitaura, Francisco-Shu; Cuesta, Antonio J.; Eisenstein, Daniel J.; Percival, Will J.; Vargas-Magaña, Mariana; Tinker, Jeremy L.; Tojeiro, Rita; Brownstein, Joel R.; Maraston, Claudia; Nichol, Robert C.; Olmstead, Matthew D.; Samushia, Lado; Seo, Hee-Jong; Streblyanska, Alina; Zhao, Gong-bo

2017-05-01

We extract cosmological information from the anisotropic power-spectrum measurements from the recently completed Baryon Oscillation Spectroscopic Survey (BOSS), extending the concept of clustering wedges to Fourier space. Making use of new fast-Fourier-transform-based estimators, we measure the power-spectrum clustering wedges of the BOSS sample by filtering out the information of Legendre multipoles ℓ > 4. Our modelling of these measurements is based on novel approaches to describe non-linear evolution, bias and redshift-space distortions, which we test using synthetic catalogues based on large-volume N-body simulations. We are able to include smaller scales than in previous analyses, resulting in tighter cosmological constraints. Using three overlapping redshift bins, we measure the angular-diameter distance, the Hubble parameter and the cosmic growth rate, and explore the cosmological implications of our full-shape clustering measurements in combination with cosmic microwave background and Type Ia supernova data. Assuming a Λ cold dark matter (ΛCDM) cosmology, we constrain the matter density to Ω M= 0.311_{-0.010}^{+0.009} and the Hubble parameter to H_0 = 67.6_{-0.6}^{+0.7} km s^{-1 Mpc^{-1}}, at a confidence level of 68 per cent. We also allow for non-standard dark energy models and modifications of the growth rate, finding good agreement with the ΛCDM paradigm. For example, we constrain the equation-of-state parameter to w = -1.019_{-0.039}^{+0.048}. This paper is part of a set that analyses the final galaxy-clustering data set from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. to produce the final cosmological constraints from BOSS.

Cosmology with the cosmic web

NASA Astrophysics Data System (ADS)

Forero-Romero, J. E.

2017-07-01

This talk summarizes different algorithms that can be used to trace the cosmic web both in simulations and observations. We present different applications in galaxy formation and cosmology. To finalize, we show how the Dark Energy Spectroscopic Instrument (DESI) could be a good place to apply these techniques.

Holographic cosmology and phase transitions of SYM theory

NASA Astrophysics Data System (ADS)

Ghoroku, Kazuo; Meyer, René; Toyoda, Fumihiko

2017-10-01

We study the time development of strongly coupled N =4 supersymmetric Yang Mills (SYM) theory on cosmological Friedmann-Robertson-Walker (FRW) backgrounds via the AdS/CFT correspondence. We implement the cosmological background as a boundary metric fulfilling the Friedmann equation with a four-dimensional cosmological constant and a dark radiation term. We analyze the dual bulk solution of the type IIB supergravity and find that the time dependence of the FRW background strongly influences the dynamical properties of the SYM theory. We in particular find a phase transition between a confined and a deconfined phase. We also argue that some cosmological solutions could be related to the inflationary scenario.

HACC: Simulating sky surveys on state-of-the-art supercomputing architectures

NASA Astrophysics Data System (ADS)

Habib, Salman; Pope, Adrian; Finkel, Hal; Frontiere, Nicholas; Heitmann, Katrin; Daniel, David; Fasel, Patricia; Morozov, Vitali; Zagaris, George; Peterka, Tom; Vishwanath, Venkatram; Lukić, Zarija; Sehrish, Saba; Liao, Wei-keng

2016-01-01

Current and future surveys of large-scale cosmic structure are associated with a massive and complex datastream to study, characterize, and ultimately understand the physics behind the two major components of the 'Dark Universe', dark energy and dark matter. In addition, the surveys also probe primordial perturbations and carry out fundamental measurements, such as determining the sum of neutrino masses. Large-scale simulations of structure formation in the Universe play a critical role in the interpretation of the data and extraction of the physics of interest. Just as survey instruments continue to grow in size and complexity, so do the supercomputers that enable these simulations. Here we report on HACC (Hardware/Hybrid Accelerated Cosmology Code), a recently developed and evolving cosmology N-body code framework, designed to run efficiently on diverse computing architectures and to scale to millions of cores and beyond. HACC can run on all current supercomputer architectures and supports a variety of programming models and algorithms. It has been demonstrated at scale on Cell- and GPU-accelerated systems, standard multi-core node clusters, and Blue Gene systems. HACC's design allows for ease of portability, and at the same time, high levels of sustained performance on the fastest supercomputers available. We present a description of the design philosophy of HACC, the underlying algorithms and code structure, and outline implementation details for several specific architectures. We show selected accuracy and performance results from some of the largest high resolution cosmological simulations so far performed, including benchmarks evolving more than 3.6 trillion particles.

HACC: Simulating sky surveys on state-of-the-art supercomputing architectures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Habib, Salman; Pope, Adrian; Finkel, Hal

2016-01-01

Current and future surveys of large-scale cosmic structure are associated with a massive and complex datastream to study, characterize, and ultimately understand the physics behind the two major components of the ‘Dark Universe’, dark energy and dark matter. In addition, the surveys also probe primordial perturbations and carry out fundamental measurements, such as determining the sum of neutrino masses. Large-scale simulations of structure formation in the Universe play a critical role in the interpretation of the data and extraction of the physics of interest. Just as survey instruments continue to grow in size and complexity, so do the supercomputers thatmore » enable these simulations. Here we report on HACC (Hardware/Hybrid Accelerated Cosmology Code), a recently developed and evolving cosmology N-body code framework, designed to run efficiently on diverse computing architectures and to scale to millions of cores and beyond. HACC can run on all current supercomputer architectures and supports a variety of programming models and algorithms. It has been demonstrated at scale on Cell- and GPU-accelerated systems, standard multi-core node clusters, and Blue Gene systems. HACC’s design allows for ease of portability, and at the same time, high levels of sustained performance on the fastest supercomputers available. We present a description of the design philosophy of HACC, the underlying algorithms and code structure, and outline implementation details for several specific architectures. We show selected accuracy and performance results from some of the largest high resolution cosmological simulations so far performed, including benchmarks evolving more than 3.6 trillion particles.« less

Simulating the universe(s): from cosmic bubble collisions to cosmological observables with numerical relativity

NASA Astrophysics Data System (ADS)

Wainwright, Carroll L.; Johnson, Matthew C.; Peiris, Hiranya V.; Aguirre, Anthony; Lehner, Luis; Liebling, Steven L.

2014-03-01

The theory of eternal inflation in an inflaton potential with multiple vacua predicts that our universe is one of many bubble universes nucleating and growing inside an ever-expanding false vacuum. The collision of our bubble with another could provide an important observational signature to test this scenario. We develop and implement an algorithm for accurately computing the cosmological observables arising from bubble collisions directly from the Lagrangian of a single scalar field. We first simulate the collision spacetime by solving Einstein's equations, starting from nucleation and ending at reheating. Taking advantage of the collision's hyperbolic symmetry, the simulations are performed with a 1+1-dimensional fully relativistic code that uses adaptive mesh refinement. We then calculate the comoving curvature perturbation in an open Friedmann-Robertson-Walker universe, which is used to determine the temperature anisotropies of the cosmic microwave background radiation. For a fiducial Lagrangian, the anisotropies are well described by a power law in the cosine of the angular distance from the center of the collision signature. For a given form of the Lagrangian, the resulting observational predictions are inherently statistical due to stochastic elements of the bubble nucleation process. Further uncertainties arise due to our imperfect knowledge about inflationary and pre-recombination physics. We characterize observational predictions by computing the probability distributions over four phenomenological parameters which capture these intrinsic and model uncertainties. This represents the first fully-relativistic set of predictions from an ensemble of scalar field models giving rise to eternal inflation, yielding significant differences from previous non-relativistic approximations. Thus, our results provide a basis for a rigorous confrontation of these theories with cosmological data.

Stochastic Order Redshift Technique (SORT): a simple, efficient and robust method to improve cosmological redshift measurements

NASA Astrophysics Data System (ADS)

Tejos, Nicolas; Rodríguez-Puebla, Aldo; Primack, Joel R.

2018-01-01

We present a simple, efficient and robust approach to improve cosmological redshift measurements. The method is based on the presence of a reference sample for which a precise redshift number distribution (dN/dz) can be obtained for different pencil-beam-like sub-volumes within the original survey. For each sub-volume we then impose that: (i) the redshift number distribution of the uncertain redshift measurements matches the reference dN/dz corrected by their selection functions and (ii) the rank order in redshift of the original ensemble of uncertain measurements is preserved. The latter step is motivated by the fact that random variables drawn from Gaussian probability density functions (PDFs) of different means and arbitrarily large standard deviations satisfy stochastic ordering. We then repeat this simple algorithm for multiple arbitrary pencil-beam-like overlapping sub-volumes; in this manner, each uncertain measurement has multiple (non-independent) 'recovered' redshifts which can be used to estimate a new redshift PDF. We refer to this method as the Stochastic Order Redshift Technique (SORT). We have used a state-of-the-art N-body simulation to test the performance of SORT under simple assumptions and found that it can improve the quality of cosmological redshifts in a robust and efficient manner. Particularly, SORT redshifts (zsort) are able to recover the distinctive features of the so-called 'cosmic web' and can provide unbiased measurement of the two-point correlation function on scales ≳4 h-1Mpc. Given its simplicity, we envision that a method like SORT can be incorporated into more sophisticated algorithms aimed to exploit the full potential of large extragalactic photometric surveys.

Chemical enrichment in cosmological, smoothed particle hydrodynamics simulations

NASA Astrophysics Data System (ADS)

Wiersma, Robert P. C.; Schaye, Joop; Theuns, Tom; Dalla Vecchia, Claudio; Tornatore, Luca

2009-10-01

We present an implementation of stellar evolution and chemical feedback for smoothed particle hydrodynamics simulations. We consider the timed release of individual elements by both massive (Type II supernovae and stellar winds) and intermediate-mass stars (Type Ia supernovae and asymptotic giant branch stars). We illustrate the results of our method using a suite of cosmological simulations that include new prescriptions for radiative cooling, star formation and galactic winds. Radiative cooling is implemented element-by-element, in the presence of an ionizing radiation background, and we track all 11 elements that contribute significantly to the radiative cooling. While all simulations presented here use a single set of physical parameters, we take specific care to investigate the robustness of the predictions of chemodynamical simulations with respect to the ingredients, the methods and the numerical convergence. A comparison of nucleosynthetic yields taken from the literature indicates that relative abundance ratios may only be reliable at the factor of 2 level, even for a fixed initial mass function. Abundances relative to iron are even more uncertain because the rate of Type Ia supernovae is not well known. We contrast two reasonable definitions of the metallicity of a resolution element and find that while they agree for high metallicities, there are large differences at low metallicities. We argue that the discrepancy is indicative of the lack of metal mixing caused by the fact that metals are stuck to particles. We argue that since this is a (numerical) sampling problem, solving it by using a poorly constrained physical process such as diffusion could have undesired consequences. We demonstrate that the two metallicity definitions result in redshift z = 0 stellar masses that can differ by up to a factor of 2, because of the sensitivity of the cooling rates to the elemental abundances. Finally, we use several 5123 particle simulations to investigate the

Computational complexity of the landscape II-Cosmological considerations

NASA Astrophysics Data System (ADS)

Denef, Frederik; Douglas, Michael R.; Greene, Brian; Zukowski, Claire

2018-05-01

We propose a new approach for multiverse analysis based on computational complexity, which leads to a new family of "computational" measure factors. By defining a cosmology as a space-time containing a vacuum with specified properties (for example small cosmological constant) together with rules for how time evolution will produce the vacuum, we can associate global time in a multiverse with clock time on a supercomputer which simulates it. We argue for a principle of "limited computational complexity" governing early universe dynamics as simulated by this supercomputer, which translates to a global measure for regulating the infinities of eternal inflation. The rules for time evolution can be thought of as a search algorithm, whose details should be constrained by a stronger principle of "minimal computational complexity". Unlike previously studied global measures, ours avoids standard equilibrium considerations and the well-known problems of Boltzmann Brains and the youngness paradox. We also give various definitions of the computational complexity of a cosmology, and argue that there are only a few natural complexity classes.

Supermassive Black Hole Binaries in High Performance Massively Parallel Direct N-body Simulations on Large GPU Clusters

NASA Astrophysics Data System (ADS)

Spurzem, R.; Berczik, P.; Zhong, S.; Nitadori, K.; Hamada, T.; Berentzen, I.; Veles, A.

2012-07-01

Astrophysical Computer Simulations of Dense Star Clusters in Galactic Nuclei with Supermassive Black Holes are presented using new cost-efficient supercomputers in China accelerated by graphical processing cards (GPU). We use large high-accuracy direct N-body simulations with Hermite scheme and block-time steps, parallelised across a large number of nodes on the large scale and across many GPU thread processors on each node on the small scale. A sustained performance of more than 350 Tflop/s for a science run on using simultaneously 1600 Fermi C2050 GPUs is reached; a detailed performance model is presented and studies for the largest GPU clusters in China with up to Petaflop/s performance and 7000 Fermi GPU cards. In our case study we look at two supermassive black holes with equal and unequal masses embedded in a dense stellar cluster in a galactic nucleus. The hardening processes due to interactions between black holes and stars, effects of rotation in the stellar system and relativistic forces between the black holes are simultaneously taken into account. The simulation stops at the complete relativistic merger of the black holes.

Towards Forming a Primordial Protostar in a Cosmological AMR Simulation

NASA Astrophysics Data System (ADS)

Turk, Matthew J.; Abel, Tom; O'Shea, Brian W.

2008-03-01

Modeling the formation of the first stars in the universe is a well-posed problem and ideally suited for computational investigation.We have conducted high-resolution numerical studies of the formation of primordial stars. Beginning with primordial initial conditions appropriate for a ΛCDM model, we used the Eulerian adaptive mesh refinement code (Enzo) to achieve unprecedented numerical resolution, resolving cosmological scales as well as sub-stellar scales simultaneously. Building on the work of Abel, Bryan and Norman (2002), we followed the evolution of the first collapsing cloud until molecular hydrogen is optically thick to cooling radiation. In addition, the calculations account for the process of collision-induced emission (CIE) and add approximations to the optical depth in both molecular hydrogen roto-vibrational cooling and CIE. Also considered are the effects of chemical heating/cooling from the formation/destruction of molecular hydrogen. We present the results of these simulations, showing the formation of a 10 Jupiter-mass protostellar core bounded by a strongly aspherical accretion shock. Accretion rates are found to be as high as one solar mass per year.

Cosmology and the weak interaction

NASA Technical Reports Server (NTRS)

Schramm, David N.

1989-01-01

The weak interaction plays a critical role in modern Big Bang cosmology. Two of its most publicized comological connections are emphasized: big bang nucleosynthesis and dark matter. The first of these is connected to the cosmological prediction of neutrine flavors, N(sub nu) is approximately 3 which in now being confirmed. The second is interrelated to the whole problem of galacty and structure formation in the universe. The role of the weak interaction both for dark matter candidates and for the problem of generating seeds to form structure is demonstrated.

Generalized teleparallel cosmology and initial singularity crossing

DOE Office of Scientific and Technical Information (OSTI.GOV)

Awad, Adel; Nashed, Gamal, E-mail: Adel.Awad@bue.edu.eg, E-mail: gglnashed@sci.asu.edu.eg

We present a class of cosmological solutions for a generalized teleparallel gravity with f ( T )= T +α̃ (− T ) {sup n} , where α̃ is some parameter and n is an integer or half-integer. Choosing α̃ ∼ G {sup n} {sup −1}, where G is the gravitational constant, and working with an equation of state p = w ρ, one obtains a cosmological solution with multiple branches. The dynamics of the solution describes standard cosmology at late times, but the higher-torsion correction changes the nature of the initial singularity from big bang to a sudden singularity. Themore » milder behavior of the sudden singularity enables us to extend timelike or lightlike curves, through joining two disconnected branches of solution at the singularity, leaving the singularity traversable. We show that this extension is consistent with the field equations through checking the known junction conditions for generalized teleparallel gravity. This suggests that these solutions describe a contracting phase a prior to the expanding phase of the universe.« less

Approximate Bayesian computation for forward modeling in cosmology

DOE Office of Scientific and Technical Information (OSTI.GOV)

Akeret, Joël; Refregier, Alexandre; Amara, Adam

Bayesian inference is often used in cosmology and astrophysics to derive constraints on model parameters from observations. This approach relies on the ability to compute the likelihood of the data given a choice of model parameters. In many practical situations, the likelihood function may however be unavailable or intractable due to non-gaussian errors, non-linear measurements processes, or complex data formats such as catalogs and maps. In these cases, the simulation of mock data sets can often be made through forward modeling. We discuss how Approximate Bayesian Computation (ABC) can be used in these cases to derive an approximation to themore » posterior constraints using simulated data sets. This technique relies on the sampling of the parameter set, a distance metric to quantify the difference between the observation and the simulations and summary statistics to compress the information in the data. We first review the principles of ABC and discuss its implementation using a Population Monte-Carlo (PMC) algorithm and the Mahalanobis distance metric. We test the performance of the implementation using a Gaussian toy model. We then apply the ABC technique to the practical case of the calibration of image simulations for wide field cosmological surveys. We find that the ABC analysis is able to provide reliable parameter constraints for this problem and is therefore a promising technique for other applications in cosmology and astrophysics. Our implementation of the ABC PMC method is made available via a public code release.« less

Smoot Group Cosmology

Science.gov Websites

Image: NASA WMAP George F. Smoot and John Mather share the 2006 Nobel prize "for their the Universe About Cosmology Planck Satellite Launched Cosmology Videos Professor George Smoot's group science goals regarding cosmology. George Smoot named Director of Korean Cosmology Institute The GRB

[Early Greek medicine and Plato's cosmology].

PubMed

Rhee, Kee-Bag

2004-06-01

The purpose of this paper is to show the influence of Early greek medicine on Plato's Cosmology. Alcmaeon holds that health depends on proportion (equality; isonomia) or proportioned mixture of opposing factors. This notion dominated nearly all greek medicine, and also influenced Plato's cosmology greatly. Generally each greek doctors believed that man consisted of opposing factors, though these are designated differently. Alcmaeon takes powers - hot and dry, cold and hot, vitter, sweet and the rest as those factors. On the other hand, Philistion of Locri adopts the four element theory of Empedocles. He conceives that human body as a mixture of the four elements, and health consists in proportion of these opposing four element, basically as Alcmaeon. This notion is accepted by Plato. Only Plato differs from Philistion in that he doesn't consider the four elements as the ultimate factors. In Timaeus Plato explain that the Demiourgos constructed the four elements through introducing 'proportion' into the primitive materials (the oppositives) by means of shapes and numbers. And Plato thinks that the cosmic body and soul was constructed basically in the same way as the four elements. This is true of the human body and soul. Also Plato explicates diseases from standpoint of proportion or symmetry. Moreover according to Philebus, the good states (i.e. 'health', 'music', 'season' etc) in the cosmos arises out of the right mixture of the limit and the unlimited. In the other word this mixture is proportioned mixture of the oppositives by aid of ratios. In short Plato believes that both the cosmos itself and the good states s proportioned mixture of the oppositives. Thus Plato' cosmology is fundamentally based upon Alcmaeon's or Philistion's concept of Health.

Statistical Issues in Galaxy Cluster Cosmology

NASA Technical Reports Server (NTRS)

Mantz, Adam

2013-01-01

The number and growth of massive galaxy clusters are sensitive probes of cosmological structure formation. Surveys at various wavelengths can detect clusters to high redshift, but the fact that cluster mass is not directly observable complicates matters, requiring us to simultaneously constrain scaling relations of observable signals with mass. The problem can be cast as one of regression, in which the data set is truncated, the (cosmology-dependent) underlying population must be modeled, and strong, complex correlations between measurements often exist. Simulations of cosmological structure formation provide a robust prediction for the number of clusters in the Universe as a function of mass and redshift (the mass function), but they cannot reliably predict the observables used to detect clusters in sky surveys (e.g. X-ray luminosity). Consequently, observers must constrain observable-mass scaling relations using additional data, and use the scaling relation model in conjunction with the mass function to predict the number of clusters as a function of redshift and luminosity.

Cosmological tests of the Hoyle-Narlikar conformal gravity

NASA Technical Reports Server (NTRS)

Canuto, V. M.; Narlikar, J. V.

1980-01-01

For the first time the Hoyle-Narlikar theory with creation of matter and a variable gravitational constant G, is subjected to the following cosmological tests: (1) the magnitude versus z relation, (2) the N(m) versus m relation for quasars, (3) the metric angular diameters versus z relation, (4) the isophotal angles versus z relation, (5) the log N-log S radio source count, and finally (6) the 3 K radiation. It is shown that the theory passes all these tests just as well as the standard cosmology, with the additional advantage that the geometry of the universe is uniquely determined, with a curvature parameter equal to zero. It is also interesting to note that the variability of G affects the log N-log S curve in a way similar to the density evolution introduced in standard cosmologies. The agreement with the data is therefore achieved without recourse to an ad hoc density evolution.

Structure formation by a fifth force: N-body versus linear simulations

NASA Astrophysics Data System (ADS)

Li, Baojiu; Zhao, Hongsheng

2009-08-01

We lay out the frameworks to numerically study the structure formation in both linear and nonlinear regimes in general dark-matter-coupled scalar field models, and give an explicit example where the scalar field serves as a dynamical dark energy. Adopting parameters of the scalar field which yield a realistic cosmic microwave background (CMB) spectrum, we generate the initial conditions for our N-body simulations, which follow the spatial distributions of the dark matter and the scalar field by solving their equations of motion using the multilevel adaptive grid technique. We show that the spatial configuration of the scalar field tracks well the voids and clusters of dark matter. Indeed, the propagation of scalar degree of freedom effectively acts as a fifth force on dark matter particles, whose range and magnitude are determined by the two model parameters (μ,γ), local dark matter density as well as the background value for the scalar field. The model behaves like the ΛCDM paradigm on scales relevant to the CMB spectrum, which are well beyond the probe of the local fifth force and thus not significantly affected by the matter-scalar coupling. On scales comparable or shorter than the range of the local fifth force, the fifth force is perfectly parallel to gravity and their strengths have a fixed ratio 2γ2 determined by the matter-scalar coupling, provided that the chameleon effect is weak; if on the other hand there is a strong chameleon effect (i.e., the scalar field almost resides at its effective potential minimum everywhere in the space), the fifth force indeed has suppressed effects in high density regions and shows no obvious correlation with gravity, which means that the dark-matter-scalar-field coupling is not simply equivalent to a rescaling of the gravitational constant or the mass of the dark matter particles. We show these spatial distributions and (lack of) correlations at typical redshifts (z=0,1,5.5) in our multigrid million-particle simulations

Post-Newtonian celestial dynamics in cosmology: Field equations

NASA Astrophysics Data System (ADS)

Kopeikin, Sergei M.; Petrov, Alexander N.

2013-02-01

Post-Newtonian celestial dynamics is a relativistic theory of motion of massive bodies and test particles under the influence of relatively weak gravitational forces. The standard approach for development of this theory relies upon the key concept of the isolated astronomical system supplemented by the assumption that the background spacetime is flat. The standard post-Newtonian theory of motion was instrumental in the explanation of the existing experimental data on binary pulsars, satellite, and lunar laser ranging, and in building precise ephemerides of planets in the Solar System. Recent studies of the formation of large-scale structures in our Universe indicate that the standard post-Newtonian mechanics fails to describe more subtle dynamical effects in motion of the bodies comprising the astronomical systems of larger size—galaxies and clusters of galaxies—where the Riemann curvature of the expanding Friedmann-Lemaître-Robertson-Walker universe interacts with the local gravitational field of the astronomical system and, as such, cannot be ignored. The present paper outlines theoretical principles of the post-Newtonian mechanics in the expanding Universe. It is based upon the gauge-invariant theory of the Lagrangian perturbations of cosmological manifold caused by an isolated astronomical N-body system (the Solar System, a binary star, a galaxy, and a cluster of galaxies). We postulate that the geometric properties of the background manifold are described by a homogeneous and isotropic Friedmann-Lemaître-Robertson-Walker metric governed by two primary components—the dark matter and the dark energy. The dark matter is treated as an ideal fluid with the Lagrangian taken in the form of pressure along with the scalar Clebsch potential as a dynamic variable. The dark energy is associated with a single scalar field with a potential which is hold unspecified as long as the theory permits. Both the Lagrangians of the dark matter and the scalar field are

The bar-halo interaction - I. From fundamental dynamics to revised N-body requirements

NASA Astrophysics Data System (ADS)

Weinberg, Martin D.; Katz, Neal

2007-02-01

A galaxy remains near equilibrium for most of its history. Only through resonances can non-axisymmetric features, such as spiral arms and bars, exert torques over large scales and change the overall structure of the galaxy. In this paper, we describe the resonant interaction mechanism in detail, derive explicit criteria for the particle number required to simulate these dynamical processes accurately using N-body simulations, and illustrate them with numerical experiments. To do this, we perform a direct numerical solution of perturbation theory, in short, by solving for each orbit in an ensemble and make detailed comparisons with N-body simulations. The criteria include: sufficient particle coverage in phase space near the resonance and enough particles to minimize gravitational potential fluctuations that will change the dynamics of the resonant encounter. These criteria are general in concept and can be applied to any dynamical interaction. We use the bar-halo interaction as our primary example owing to its technical simplicity and astronomical ubiquity. Some of our more surprising findings are as follows. First, the inner Lindblad like resonance, responsible for coupling the bar to the central halo cusp, requires more than equal-mass particles within the virial radius or inside the bar radius for a Milky Way like bar in a Navarro, Frenk & White profile. Secondly, orbits that linger near the resonance receive more angular momentum than orbits that move through the resonance quickly. Small-scale fluctuations present in state-of-the-art particle-particle simulations can knock orbits out of resonance, preventing them from lingering and, thereby, decrease the torque per orbit. This can be offset by the larger number of orbits affected by the resonance due to the diffusion. However, noise from orbiting substructure remains at least an order of magnitude too small to be of consequence. Applied to N-body simulations, the required particle numbers are sufficiently high

Virtual Petaflop Simulation: Parallel Potential Solvers and New Integrators for Gravitational Systems

NASA Technical Reports Server (NTRS)

Lake, George; Quinn, Thomas; Richardson, Derek C.; Stadel, Joachim

1999-01-01

"The orbit of any one planet depends on the combined motion of all the planets, not to mention the actions of all these on each other. To consider simultaneously all these causes of motion and to define these motions by exact laws allowing of convenient calculation exceeds, unless I am mistaken, the forces of the entire human intellect" -Isaac Newton 1687. Epochal surveys are throwing down the gauntlet for cosmological simulation. We describe three keys to meeting the challenge of N-body simulation: adaptive potential solvers, adaptive integrators and volume renormalization. With these techniques and a dedicated Teraflop facility, simulation can stay even with observation of the Universe. We also describe some problems in the formation and stability of planetary systems. Here, the challenge is to perform accurate integrations that retain Hamiltonian properties for 10(exp 13) timesteps.

Disentangling interacting dark energy cosmologies with the three-point correlation function

NASA Astrophysics Data System (ADS)

Moresco, Michele; Marulli, Federico; Baldi, Marco; Moscardini, Lauro; Cimatti, Andrea

2014-10-01

We investigate the possibility of constraining coupled dark energy (cDE) cosmologies using the three-point correlation function (3PCF). Making use of the CODECS N-body simulations, we study the statistical properties of cold dark matter (CDM) haloes for a variety of models, including a fiducial ΛCDM scenario and five models in which dark energy (DE) and CDM mutually interact. We measure both the halo 3PCF, ζ(θ), and the reduced 3PCF, Q(θ), at different scales (2 < r [h-1 Mpc ] < 40) and redshifts (0 ≤ z ≤ 2). In all cDE models considered in this work, Q(θ) appears flat at small scales (for all redshifts) and at low redshifts (for all scales), while it builds up the characteristic V-shape anisotropy at increasing redshifts and scales. With respect to the ΛCDM predictions, cDE models show lower (higher) values of the halo 3PCF for perpendicular (elongated) configurations. The effect is also scale-dependent, with differences between ΛCDM and cDE models that increase at large scales. We made use of these measurements to estimate the halo bias, that results in fair agreement with the one computed from the two-point correlation function (2PCF). The main advantage of using both the 2PCF and 3PCF is to break the bias-σ8 degeneracy. Moreover, we find that our bias estimates are approximately independent of the assumed strength of DE coupling. This study demonstrates the power of a higher order clustering analysis in discriminating between alternative cosmological scenarios, for both present and forthcoming galaxy surveys, such as e.g. Baryon Oscillation Spectroscopic Survey and Euclid.

Development and Evaluation of an Order-N Formulation for Multi-Flexible Body Space Systems

NASA Technical Reports Server (NTRS)

Ghosh, Tushar K.; Quiocho, Leslie J.

2013-01-01

This paper presents development of a generic recursive Order-N algorithm for systems with rigid and flexible bodies, in tree or closed-loop topology, with N being the number of bodies of the system. Simulation results are presented for several test cases to verify and evaluate the performance of the code compared to an existing efficient dense mass matrix-based code. The comparison brought out situations where Order-N or mass matrix-based algorithms could be useful.

Evaluating galactic habitability using high-resolution cosmological simulations of galaxy formation

NASA Astrophysics Data System (ADS)

Forgan, Duncan; Dayal, Pratika; Cockell, Charles; Libeskind, Noam

2017-01-01

We present the first model that couples high-resolution simulations of the formation of local group galaxies with calculations of the galactic habitable zone (GHZ), a region of space which has sufficient metallicity to form terrestrial planets without being subject to hazardous radiation. These simulations allow us to make substantial progress in mapping out the asymmetric three-dimensional GHZ and its time evolution for the Milky Way (MW) and Triangulum (M33) galaxies, as opposed to works that generally assume an azimuthally symmetric GHZ. Applying typical habitability metrics to MW and M33, we find that while a large number of habitable planets exist as close as a few kiloparsecs from the galactic centre, the probability of individual planetary systems being habitable rises as one approaches the edge of the stellar disc. Tidal streams and satellite galaxies also appear to be fertile grounds for habitable planet formation. In short, we find that both galaxies arrive at similar GHZs by different evolutionary paths, as measured by the first and third quartiles of surviving biospheres. For the MW, this interquartile range begins as a narrow band at large radii, expanding to encompass much of the Galaxy at intermediate times before settling at a range of 2-13 kpc. In the case of M33, the opposite behaviour occurs - the initial and final interquartile ranges are quite similar, showing gradual evolution. This suggests that Galaxy assembly history strongly influences the time evolution of the GHZ, which will affect the relative time lag between biospheres in different galactic locations. We end by noting the caveats involved in such studies and demonstrate that high-resolution cosmological simulations will play a vital role in understanding habitability on galactic scales, provided that these simulations accurately resolve chemical evolution.

Comparing cosmological hydrodynamic simulations with observations of high- redshift galaxy formation

NASA Astrophysics Data System (ADS)

Finlator, Kristian Markwart

We use cosmological hydrodynamic simulations to study the impact of outflows and radiative feedback on high-redshift galaxies. For outflows, we consider simulations that assume (i) no winds, (ii) a "constant-wind" model in which the mass-loading factor and outflow speed are constant, and (iii) "momentum-driven" winds in which both parameters vary smoothly with mass. In order to treat radiative feedback, we develop a moment-based radiative transfer technique that operates in both post-processing and coupled radiative hydrodynamic modes. We first ask how outflows impact the broadband spectral energy distributions (SEDs) of six observed reionization-epoch galaxies. Simulations reproduce five regardless of the outflow prescription, while the sixth suggests an unusually bursty star formation history. We conclude that (i) simulations broadly account for available constraints on reionization-epoch galaxies, (ii) individual SEDs do not constrain outflows, and (iii) SED comparisons efficiently isolate objects that challenge simulations. We next study how outflows impact the galaxy mass metallicity relation (MZR). Momentum-driven outflows uniquely reproduce observations at z = 2. In this scenario, galaxies obey two equilibria: (i) The rate at which a galaxy processes gas into stars and outflows tracks its inflow rate; and (ii) The gas enrichment rate owing to star formation balances the dilution rate owing to inflows. Combining these conditions indicates that the MZR is dominated by the (instantaneous) variation of outflows with mass, with more-massive galaxies driving less gas into outflows per unit stellar mass formed. Turning to radiative feedback, we use post-processing simulations to study the topology of reionization. Reionization begins in overdensities and then "leaks" directly into voids, with filaments reionizing last owing to their high density and low emissivity. This result conflicts with previous findings that voids ionize last. We argue that it owes to the

Astrophysical cosmology

NASA Astrophysics Data System (ADS)

Bardeen, J. M.

The last several years have seen a tremendous ferment of activity in astrophysical cosmology. Much of the theoretical impetus has come from particle physics theories of the early universe and candidates for dark matter, but what promise to be even more significant are improved direct observations of high z galaxies and intergalactic matter, deeper and more comprehensive redshift surveys, and the increasing power of computer simulations of the dynamical evolution of large scale structure. Upper limits on the anisotropy of the microwave background radiation are gradually getting tighter and constraining more severely theoretical scenarios for the evolution of the universe.

A generic multi-flex-body dynamics, controls simulation tool for space station

NASA Technical Reports Server (NTRS)

London, Ken W.; Lee, John F.; Singh, Ramen P.; Schubele, Buddy

1991-01-01

An order (n) multiflex body Space Station simulation tool is introduced. The flex multibody modeling is generic enough to model all phases of Space Station from build up through to Assembly Complete configuration and beyond. Multibody subsystems such as the Mobile Servicing System (MSS) undergoing a prescribed translation and rotation are also allowed. The software includes aerodynamic, gravity gradient, and magnetic field models. User defined controllers can be discrete or continuous. Extensive preprocessing of 'body by body' NASTRAN flex data is built in. A significant aspect, too, is the integrated controls design capability which includes model reduction and analytic linearization.

Flexible Multi-Body Spacecraft Simulator: Design, Construction, and Experiments

DTIC Science & Technology

2017-12-01

BODY SPACECRAFT SIMULATOR: DESIGN , CONSTRUCTION, AND EXPERIMENTS by Adam L. Atwood December 2017 Thesis Advisor: Mark Karpenko Second...TYPE AND DATES COVERED Master’s thesis 4. TITLE AND SUBTITLE FLEXIBLE MULTI-BODY SPACECRAFT SIMULATOR: DESIGN , CONSTRUCTION, AND EXPERIMENTS 5...spacecraft simulator for use in testing optimal control-based slew and maneuver designs . The simulator is modified from an earlier prototype, which

Measurement of the cross section for the 4He(α, n)7Be reaction as a possible solution to the cosmological lithium problem

NASA Astrophysics Data System (ADS)

Kawabata, T.; Furuno, T.; Ichikawa, M.; Iwasa, N.; Kanada-En'yo, Y.; Koshikawa, A.; Kubono, S.; Miyawaki, E.; Morimoto, T.; Murata, M.; Nanamura, T.; Nishimura, S.; Shikata, Y.; Takahashi, Y.; Takeda, T.; Tsumura, M.; Watanabe, K.

2017-06-01

The cross section for the 4He(α,n)7Be reaction was measured at low energies between Eα = 38.50 and 39.64 MeV motivated by the cosmological lithium problem. On the basis of the detailed balance principle, the cross section for the 7Be(n,α)4He reaction was obtained at Ec.m. = 0.20-0.81 MeV close to the Big Bang nucleosynthesis (BBN) energy window for the first time. The obtained cross sections are significantly smaller than the theoretical estimation widely used in the BBN calculations. The present results suggest the 7Be(n,α)4He reaction rate is not large enough to solve the cosmological lithium problem.

Reconciling Dwarf Galaxies with ΛCDM Cosmology: Simulating a Realistic Population of Satellites around a Milky Way-mass Galaxy

NASA Astrophysics Data System (ADS)

Wetzel, Andrew R.; Hopkins, Philip F.; Kim, Ji-hoon; Faucher-Giguère, Claude-André; Kereš, Dušan; Quataert, Eliot

2016-08-01

Low-mass “dwarf” galaxies represent the most significant challenges to the cold dark matter (CDM) model of cosmological structure formation. Because these faint galaxies are (best) observed within the Local Group (LG) of the Milky Way (MW) and Andromeda (M31), understanding their formation in such an environment is critical. We present first results from the Latte Project: the Milky Way on Feedback in Realistic Environments (FIRE). This simulation models the formation of an MW-mass galaxy to z=0 within ΛCDM cosmology, including dark matter, gas, and stars at unprecedented resolution: baryon particle mass of 7070 {M}⊙ with gas kernel/softening that adapts down to 1 {pc} (with a median of 25{--}60 {pc} at z=0). Latte was simulated using the GIZMO code with a mesh-free method for accurate hydrodynamics and the FIRE-2 model for star formation and explicit feedback within a multi-phase interstellar medium. For the first time, Latte self-consistently resolves the spatial scales corresponding to half-light radii of dwarf galaxies that form around an MW-mass host down to {M}{star}≳ {10}5 {M}⊙ . Latte’s population of dwarf galaxies agrees with the LG across a broad range of properties: (1) distributions of stellar masses and stellar velocity dispersions (dynamical masses), including their joint relation; (2) the mass-metallicity relation; and (3) diverse range of star formation histories, including their mass dependence. Thus, Latte produces a realistic population of dwarf galaxies at {M}{star}≳ {10}5 {M}⊙ that does not suffer from the “missing satellites” or “too big to fail” problems of small-scale structure formation. We conclude that baryonic physics can reconcile observed dwarf galaxies with standard ΛCDM cosmology.

An analytical method to simulate the H I 21-cm visibility signal for intensity mapping experiments

NASA Astrophysics Data System (ADS)

Sarkar, Anjan Kumar; Bharadwaj, Somnath; Marthi, Visweshwar Ram

2018-01-01

Simulations play a vital role in testing and validating H I 21-cm power spectrum estimation techniques. Conventional methods use techniques like N-body simulations to simulate the sky signal which is then passed through a model of the instrument. This makes it necessary to simulate the H I distribution in a large cosmological volume, and incorporate both the light-cone effect and the telescope's chromatic response. The computational requirements may be particularly large if one wishes to simulate many realizations of the signal. In this paper, we present an analytical method to simulate the H I visibility signal. This is particularly efficient if one wishes to simulate a large number of realizations of the signal. Our method is based on theoretical predictions of the visibility correlation which incorporate both the light-cone effect and the telescope's chromatic response. We have demonstrated this method by applying it to simulate the H I visibility signal for the upcoming Ooty Wide Field Array Phase I.

Cosmological singularity theorems and splitting theorems for N-Bakry-Émery spacetimes

NASA Astrophysics Data System (ADS)

Woolgar, Eric; Wylie, William

2016-02-01

We study Lorentzian manifolds with a weight function such that the N-Bakry-Émery tensor is bounded below. Such spacetimes arise in the physics of scalar-tensor gravitation theories, including Brans-Dicke theory, theories with Kaluza-Klein dimensional reduction, and low-energy approximations to string theory. In the "pure Bakry-Émery" N = ∞ case with f uniformly bounded above and initial data suitably bounded, cosmological-type singularity theorems are known, as are splitting theorems which determine the geometry of timelike geodesically complete spacetimes for which the bound on the initial data is borderline violated. We extend these results in a number of ways. We are able to extend the singularity theorems to finite N-values N ∈ (n, ∞) and N ∈ (-∞, 1]. In the N ∈ (n, ∞) case, no bound on f is required, while for N ∈ (-∞, 1] and N = ∞, we are able to replace the boundedness of f by a weaker condition on the integral of f along future-inextendible timelike geodesics. The splitting theorems extend similarly, but when N = 1, the splitting is only that of a warped product for all cases considered. A similar limited loss of rigidity has been observed in a prior work on the N-Bakry-Émery curvature in Riemannian signature when N = 1 and appears to be a general feature.

On the cosmological gravitational waves and cosmological distances

NASA Astrophysics Data System (ADS)

Belinski, V. A.; Vereshchagin, G. V.

2018-03-01

We show that solitonic cosmological gravitational waves propagated through the Friedmann universe and generated by the inhomogeneities of the gravitational field near the Big Bang can be responsible for increase of cosmological distances.

Seeking sterile neutrinos in Finslerian cosmology

NASA Astrophysics Data System (ADS)

Wang, Deng; Meng, Xin-He

2017-11-01

For the first time, to search for sterile neutrinos in the framework of Finler geometry, we constrain four cosmological models using the most stringent constraint we can provide so far. We find that the Finslerian massless sterile neutrino model can, respectively, give a better cosmological fit to data and alleviate the current H_0 tension more effectively than the other three models. For the Finslerian massless sterile neutrino model, we obtain the constraint N_eff=3.237^{+0.092}_{-0.185}, which is consistent with Δ N_eff > 0 at the 1.03σ confidence level (CL). This gives a very weak hint of massless sterile neutrinos and may imply the non-existence of massless sterile neutrinos in the Finslerian cosmological setting. For the Finslerian massive sterile neutrino model, we obtain the constraints N_eff=3.143^{+0.064}_{-0.066}, which favors Δ N_eff > 0 at the 1.47σ CL, and m_{ν , sterile}^eff < 0.121 eV at the 2σ CL which is much tighter than the Planck results. This very tight restriction appears to indicate the massive sterile neutrinos are also non-existent in the Finslerian scenarios. Consequently, one may conclude that the sterile neutrinos are possibly non-existent in the Finslerian universe. Our results are compatible with the recent results of the neutrino oscillation experiments implemented by the Daya Bay and MINOS collaborations and the cosmic ray one carried out by the IceCube collaboration.

Origin of chemically distinct discs in the Auriga cosmological simulations

NASA Astrophysics Data System (ADS)

Grand, Robert J. J.; Bustamante, Sebastián; Gómez, Facundo A.; Kawata, Daisuke; Marinacci, Federico; Pakmor, Rüdiger; Rix, Hans-Walter; Simpson, Christine M.; Sparre, Martin; Springel, Volker

2018-03-01

The stellar disc of the Milky Way shows complex spatial and abundance structure that is central to understanding the key physical mechanisms responsible for shaping our Galaxy. In this study, we use six very high resolution cosmological zoom-in simulations of Milky Way-sized haloes to study the prevalence and formation of chemically distinct disc components. We find that our simulations develop a clearly bimodal distribution in the [α/Fe]-[Fe/H] plane. We find two main pathways to creating this dichotomy, which operate in different regions of the galaxies: (a) an early (z > 1) and intense high-[α/Fe] star formation phase in the inner region (R ≲ 5 kpc) induced by gas-rich mergers, followed by more quiescent low-[α/Fe] star formation; and (b) an early phase of high-[α/Fe] star formation in the outer disc followed by a shrinking of the gas disc owing to a temporarily lowered gas accretion rate, after which disc growth resumes. In process (b), a double-peaked star formation history around the time and radius of disc shrinking accentuates the dichotomy. If the early star formation phase is prolonged (rather than short and intense), chemical evolution proceeds as per process (a) in the inner region, but the dichotomy is less clear. In the outer region, the dichotomy is only evident if the first intense phase of star formation covers a large enough radial range before disc shrinking occurs; otherwise, the outer disc consists of only low-[α/Fe] sequence stars. We discuss the implication that both processes occurred in the Milky Way.

Relativistic numerical cosmology with silent universes

NASA Astrophysics Data System (ADS)

Bolejko, Krzysztof

2018-01-01

Relativistic numerical cosmology is most often based either on the exact solutions of the Einstein equations, or perturbation theory, or weak-field limit, or the BSSN formalism. The silent universe provides an alternative approach to investigate relativistic evolution of cosmological systems. The silent universe is based on the solution of the Einstein equations in 1  +  3 comoving coordinates with additional constraints imposed. These constraints include: the gravitational field is sourced by dust and cosmological constant only, both rotation and magnetic part of the Weyl tensor vanish, and the shear is diagnosable. This paper describes the code simsilun (free software distributed under the terms of the reposi General Public License), which implements the equations of the silent universe. The paper also discusses applications of the silent universe and it uses the Millennium simulation to set up the initial conditions for the code simsilun. The simulation obtained this way consists of 16 777 216 worldlines, which are evolved from z  =  80 to z  =  0. Initially, the mean evolution (averaged over the whole domain) follows the evolution of the background ΛCDM model. However, once the evolution of cosmic structures becomes nonlinear, the spatial curvature evolves from ΩK =0 to ΩK ≈ 0.1 at the present day. The emergence of the spatial curvature is associated with ΩM and Ω_Λ being smaller by approximately 0.05 compared to the ΛCDM.

Cosmological Origins of Water

NASA Astrophysics Data System (ADS)

Gagliano, Alexander; Taylor, Morgan; Black, William; Smidt, Joseph; Wiggins, Brandon K.

2018-01-01

Recent models indicate that the sun's protoplanetary disk provided insufficient pathways for water formation, as evidenced by [D/H]H2O measurements in asteroids and Earth's oceans. It is therefore likely that the early universe contained sites conducive to water chemistry. This research tracks the timeline and abundance rates of water using cosmological simulations in Enzo. A 64 Mpc cube of space is evolved from z = 200 to z = 2. Simulations are then centered on a massive halo, and a 26-species reaction network is applied using operator split to track water formation rates. Density projection plots with metallicity contours predict regions of water formation, which are then compared to simulated abundances at both galactic and extragalactic scales. Observational signatures of formation sites are further discussed, and allow for additional validation of the simulations used.

On the stability of soliton and hairy black hole solutions of 𝔰𝔲(N) Einstein-Yang-Mills theory with a negative cosmological constant

NASA Astrophysics Data System (ADS)

Baxter, J. Erik; Winstanley, Elizabeth

2016-02-01

We investigate the stability of spherically symmetric, purely magnetic, soliton and black hole solutions of four-dimensional 𝔰𝔲(N) Einstein-Yang-Mills theory with a negative cosmological constant Λ. These solutions are described by N - 1 magnetic gauge field functions ωj. We consider linear, spherically symmetric, perturbations of these solutions. The perturbations decouple into two sectors, known as the sphaleronic and gravitational sectors. For any N, there are no instabilities in the sphaleronic sector if all the magnetic gauge field functions ωj have no zeros and satisfy a set of N - 1 inequalities. In the gravitational sector, we prove that there are solutions which have no instabilities in a neighbourhood of stable embedded 𝔰𝔲(2) solutions, provided the magnitude of the cosmological constant |" separators=" Λ | is sufficiently large.

Comparing Results of SPH/N-body Impact Simulations Using Both Solid and Rubble-pile Target Asteroids

NASA Astrophysics Data System (ADS)

Durda, Daniel D.; Bottke, W. F.; Enke, B. L.; Nesvorný, D.; Asphaug, E.; Richardson, D. C.

2006-09-01

We have been investigating the properties of satellites and the morphology of size-frequency distributions (SFDs) resulting from a suite of 160 SPH/N-body simulations of impacts into 100-km diameter parent asteroids (Durda et al. 2004, Icarus 170, 243-257; Durda et al. 2006, Icarus, in press). These simulations have produced many valuable insights into the outcomes of cratering and disruptive impacts but were limited to monolithic basalt targets. As a natural consequence of collisional evolution, however, many asteroids have undergone a series of battering impacts that likely have left their interiors substantially fractured, if not completely rubblized. In light of this, we have re-mapped the matrix of simulations using rubble-pile target objects. We constructed the rubble-pile targets by filling the interior of the 100-km diameter spherical shell (the target envelope) with randomly sized solid spheres in mutual contact. We then assigned full damage (which reduces tensile and shear stresses to zero) to SPH particles in the contacts between the components; the remaining volume is void space. The internal spherical components have a power-law distribution of sizes simulating fragments of a pre-shattered parent object. First-look analysis of the rubble-pile results indicate some general similarities to the simulations with the monolithic targets (e.g., similar trends in the number of small, gravitationally bound satellite systems as a function of impact conditions) and some significant differences (e.g., size of largest remnants and smaller debris affecting size frequency distributions of resulting families). We will report details of a more thorough analysis and the implications for collisional models of the main asteroid belt. This work is supported by the National Science Foundation, grant number AST0407045.

Cosmology on a Mesh

NASA Astrophysics Data System (ADS)

Gill, Stuart P. D.; Knebe, Alexander; Gibson, Brad K.; Flynn, Chris; Ibata, Rodrigo A.; Lewis, Geraint F.

2003-04-01

An adaptive multi grid approach to simulating the formation of structure from collisionless dark matter is described. MLAPM (Multi-Level Adaptive Particle Mesh) is one of the most efficient serial codes available on the cosmological "market" today. As part of Swinburne University's role in the development of the Square Kilometer Array, we are implementing hydrodynamics, feedback, and radiative transfer within the MLAPM adaptive mesh, in order to simulate baryonic processes relevant to the interstellar and intergalactic media at high redshift. We will outline our progress to date in applying the existing MLAPM to a study of the decay of satellite galaxies within massive host potentials.

First star formation in ultralight particle dark matter cosmology

NASA Astrophysics Data System (ADS)

Hirano, Shingo; Sullivan, James M.; Bromm, Volker

2018-01-01

The formation of the first stars in the high-redshift Universe is a sensitive probe of the small-scale, particle physics nature of dark matter (DM). We carry out cosmological simulations of primordial star formation in ultralight, axion-like particle DM cosmology, with masses of 10-22 and 10-21 eV, with de Broglie wavelengths approaching galactic scales (˜ kpc). The onset of star formation is delayed, and shifted to more massive host structures. For the lightest DM particle mass explored here, first stars form at z ˜ 7 in structures with ˜109 M⊙, compared to the standard minihalo environment within the Λ cold dark matter (ΛCDM) cosmology, where z ˜ 20-30 and ˜105-106 M⊙. Despite this greatly altered DM host environment, the thermodynamic behaviour of the metal-free gas as it collapses into the DM potential well asymptotically approaches a very similar evolutionary track. Thus, the fragmentation properties are predicted to remain the same as in ΛCDM cosmology, implying a similar mass scale for the first stars. These results predict intense starbursts in the axion cosmologies, which may be amenable to observations with the James Webb Space Telescope.

GRAPE-4: A special-purpose computer for gravitational N-body problems

DOE Office of Scientific and Technical Information (OSTI.GOV)

Makino, Junichiro; Taiji, Makoto; Ebisuzaki, Toshikazu

1995-12-01

We describe GRAPE-4, a special-purpose computer for gravitational N-body simulations. In gravitational N-body simulations, almost all computing time is spent for the calculation of interaction between particles. GRAPE-4 is a specialized hardware to calculate the interaction between particles. It is used with a general-purpose host computer that performs all calculations other than the force calculation. With this architecture, it is relatively easy to realize a massively parallel system. In 1991, we developed the GRAPE-3 system with the peak speed equivalent to 14.4 Gflops. It consists of 48 custom pipelined processors. In 1992 we started the development of GRAPE-4. The GRAPE-4more » system will consist of 1920 custom pipeline chips. Each chip has the speed of 600 Mflops, when operated on 30 MHz clock. A prototype system with two custom LSIs has been completed July 1994, and the full system is now under manufacturing.« less

zeldovich-PLT: Zel'dovich approximation initial conditions generator

NASA Astrophysics Data System (ADS)

Eisenstein, Daniel; Garrison, Lehman

2016-05-01

zeldovich-PLT generates Zel'dovich approximation (ZA) initial conditions (i.e. first-order Lagrangian perturbation theory) for cosmological N-body simulations, optionally applying particle linear theory (PLT) corrections.

A finite-element simulation of galvanic coupling intra-body communication based on the whole human body.

PubMed

Song, Yong; Zhang, Kai; Hao, Qun; Hu, Lanxin; Wang, Jingwen; Shang, Fuzhou

2012-10-09

Simulation based on the finite-element (FE) method plays an important role in the investigation of intra-body communication (IBC). In this paper, a finite-element model of the whole body model used for the IBC simulation is proposed and verified, while the FE simulation of the galvanic coupling IBC with different signal transmission paths has been achieved. Firstly, a novel finite-element method for modeling the whole human body is proposed, and a FE model of the whole human body used for IBC simulation was developed. Secondly, the simulations of the galvanic coupling IBC with the different signal transmission paths were implemented. Finally, the feasibility of the proposed method was verified by using in vivo measurements within the frequency range of 10 kHz-5 MHz, whereby some important conclusions were deduced. Our results indicate that the proposed method will offer significant advantages in the investigation of the galvanic coupling intra-body communication.

A Finite-Element Simulation of Galvanic Coupling Intra-Body Communication Based on the Whole Human Body

PubMed Central

Song, Yong; Zhang, Kai; Hao, Qun; Hu, Lanxin; Wang, Jingwen; Shang, Fuzhou

2012-01-01

Simulation based on the finite-element (FE) method plays an important role in the investigation of intra-body communication (IBC). In this paper, a finite-element model of the whole body model used for the IBC simulation is proposed and verified, while the FE simulation of the galvanic coupling IBC with different signal transmission paths has been achieved. Firstly, a novel finite-element method for modeling the whole human body is proposed, and a FE model of the whole human body used for IBC simulation was developed. Secondly, the simulations of the galvanic coupling IBC with the different signal transmission paths were implemented. Finally, the feasibility of the proposed method was verified by using in vivo measurements within the frequency range of 10 kHz–5 MHz, whereby some important conclusions were deduced. Our results indicate that the proposed method will offer significant advantages in the investigation of the galvanic coupling intra-body communication. PMID:23202010

Correlation between centre offsets and gas velocity dispersion of galaxy clusters in cosmological simulations

NASA Astrophysics Data System (ADS)

Li, Ming-Hua; Zhu, Weishan; Zhao, Dong

2018-05-01

The gas is the dominant component of baryonic matter in most galaxy groups and clusters. The spatial offsets of gas centre from the halo centre could be an indicator of the dynamical state of cluster. Knowledge of such offsets is important for estimate the uncertainties when using clusters as cosmological probes. In this paper, we study the centre offsets roff between the gas and that of all the matter within halo systems in ΛCDM cosmological hydrodynamic simulations. We focus on two kinds of centre offsets: one is the three-dimensional PB offsets between the gravitational potential minimum of the entire halo and the barycentre of the ICM, and the other is the two-dimensional PX offsets between the potential minimum of the halo and the iterative centroid of the projected synthetic X-ray emission of the halo. Haloes at higher redshifts tend to have larger values of rescaled offsets roff/r200 and larger gas velocity dispersion σ v^gas/σ _{200}. For both types of offsets, we find that the correlation between the rescaled centre offsets roff/r200 and the rescaled 3D gas velocity dispersion, σ _v^gas/σ _{200} can be approximately described by a quadratic function as r_{off}/r_{200} ∝ (σ v^gas/σ _{200} - k_2)2. A Bayesian analysis with MCMC method is employed to estimate the model parameters. Dependence of the correlation relation on redshifts and the gas mass fraction are also investigated.

Singularity free N-body simulations called 'Dynamic Universe Model' don't require dark matter

NASA Astrophysics Data System (ADS)

Naga Parameswara Gupta, Satyavarapu

For finding trajectories of Pioneer satellite (Anomaly), New Horizons satellite going to Pluto, the Calculations of Dynamic Universe model can be successfully applied. No dark matter is assumed within solar system radius. The effect on the masses around SUN shows as though there is extra gravitation pull toward SUN. It solves the Dynamics of Extra-solar planets like Planet X, satellite like Pioneer and NH for 3-Position, 3-velocity 3-accelaration for their masses, considering the complex situation of Multiple planets, Stars, Galaxy parts and Galaxy centre and other Galaxies Using simple Newtonian Physics. It already solved problems Missing mass in Galaxies observed by galaxy circular velocity curves successfully. Singularity free Newtonian N-body simulations Historically, King Oscar II of Sweden an-nounced a prize to a solution of N-body problem with advice given by Güsta Mittag-Leffler in 1887. He announced `Given a system of arbitrarily many mass points that attract each according to Newton's law, under the assumption that no two points ever collide, try to find a representation of the coordinates of each point as a series in a variable that is some known function of time and for all of whose values the series converges uniformly.'[This is taken from Wikipedia]. The announced dead line that time was1st June 1888. And after that dead line, on 21st January 1889, Great mathematician Poincaré claimed that prize. Later he himself sent a telegram to journal Acta Mathematica to stop printing the special issue after finding the error in his solution. Yet for such a man of science reputation is important than money. [ Ref Book `Celestial mechanics: the waltz of the planets' By Alessandra Celletti, Ettore Perozzi, page 27]. He realized that he has been wrong in his general stability result! But till now nobody could solve that problem or claimed that prize. Later all solutions resulted in singularities and collisions of masses, given by many people

Higgs cosmology

NASA Astrophysics Data System (ADS)

Rajantie, Arttu

2018-01-01

The discovery of the Higgs boson in 2012 and other results from the Large Hadron Collider have confirmed the standard model of particle physics as the correct theory of elementary particles and their interactions up to energies of several TeV. Remarkably, the theory may even remain valid all the way to the Planck scale of quantum gravity, and therefore it provides a solid theoretical basis for describing the early Universe. Furthermore, the Higgs field itself has unique properties that may have allowed it to play a central role in the evolution of the Universe, from inflation to cosmological phase transitions and the origin of both baryonic and dark matter, and possibly to determine its ultimate fate through the electroweak vacuum instability. These connections between particle physics and cosmology have given rise to a new and growing field of Higgs cosmology, which promises to shed new light on some of the most puzzling questions about the Universe as new data from particle physics experiments and cosmological observations become available. This article is part of the Theo Murphy meeting issue `Higgs cosmology'.

The Atacama Cosmology Telescope: Cosmological Parameters from the 2008 Power Spectrum

NASA Technical Reports Server (NTRS)

Dunkley, J.; Hlozek, R.; Sievers, J.; Acquaviva, V.; Ade, P. A. R.; Aguirre, P.; Amiri, M.; Appel, J. W.; Barrientos, L. F.; Battistelli, E. S.;

Hydrodynamic Simulations and Tomographic Reconstructions of the Intergalactic Medium

NASA Astrophysics Data System (ADS)

Stark, Casey William

discuss differences in their convergence behavior, their overall agreement, and the implications for cosmological constraints. In the second part of my thesis, I present a tomographic reconstruction method that allows us to make 3D maps of the IGM with Mpc resolution. In order to make reconstructions of large surveys computationally feasible, I developed a new Wiener Filter application with an algorithm specialized to our problem, which significantly reduces the space and time complexity compared to previous implementations. I explore two scientific applications of the maps: finding protoclusters by searching the maps for large, contiguous regions of low flux and finding cosmic voids by searching the maps for regions of high flux. Using a large N-body simulation, I identify and characterize both protoclusters and voids at z = 2.5, in the middle of the redshift range being mapped by ongoing surveys. I provide simple methods for identifying protocluster and void candidates in the tomographic flux maps, and then test them on mock surveys and reconstructions. I present forecasts for sample purity and completeness and other scientific applications of these large, high-redshift objects.

The Atacama Cosmology Telescope: Cosmological Parameters from Three Seasons of Data

NASA Technical Reports Server (NTRS)

Seivers, Jonathan L.; Hlozek, Renee A.; Nolta, Michael R.; Acquaviva, Viviana; Addison, Graeme E.; Ade, Peter A. R.; Aguirre, Paula; Amiri, Mandana; Appel, John W.; Barrientos, L. Felipe;

Cosmological singularity theorems and splitting theorems for N-Bakry-Émery spacetimes

DOE Office of Scientific and Technical Information (OSTI.GOV)

Woolgar, Eric, E-mail: ewoolgar@ualberta.ca; Wylie, William, E-mail: wwylie@syr.edu

We study Lorentzian manifolds with a weight function such that the N-Bakry-Émery tensor is bounded below. Such spacetimes arise in the physics of scalar-tensor gravitation theories, including Brans-Dicke theory, theories with Kaluza-Klein dimensional reduction, and low-energy approximations to string theory. In the “pure Bakry-Émery” N = ∞ case with f uniformly bounded above and initial data suitably bounded, cosmological-type singularity theorems are known, as are splitting theorems which determine the geometry of timelike geodesically complete spacetimes for which the bound on the initial data is borderline violated. We extend these results in a number of ways. We are able tomore » extend the singularity theorems to finite N-values N ∈ (n, ∞) and N ∈ (−∞, 1]. In the N ∈ (n, ∞) case, no bound on f is required, while for N ∈ (−∞, 1] and N = ∞, we are able to replace the boundedness of f by a weaker condition on the integral of f along future-inextendible timelike geodesics. The splitting theorems extend similarly, but when N = 1, the splitting is only that of a warped product for all cases considered. A similar limited loss of rigidity has been observed in a prior work on the N-Bakry-Émery curvature in Riemannian signature when N = 1 and appears to be a general feature.« less

Optimal Shape Design of Mail-Slot Nacelle on N3-X Hybrid Wing-Body Configuration

NASA Technical Reports Server (NTRS)

Kim, Hyoungjin; Liou, Meng-Sing

2013-01-01

System studies show that a N3-X hybrid wing-body aircraft with a turboelectric distributed propulsion system using a mail-slot inlet/nozzle nacelle can meet the environmental and performance goals for N+3 generation transports (three generations beyond the current air transport technology level) set by NASA's Subsonic Fixed Wing Project. In this study, a Navier-Stokes flow simulation of N3-X on hybrid unstructured meshes was conducted, including the mail-slot propulsor. The geometry of the mail-slot propulsor was generated by a CAD (Computer-Aided Design)-free shape parameterization. A novel body force model generation approach was suggested for a more realistic and efficient simulation of the flow turning, pressure rise and loss effects of the fan blades and the inlet-fan interactions. Flow simulation results of the N3-X demonstrates the validity of the present approach. An optimal Shape design of the mail-slot nacelle surface was conducted to reduce strength of shock waves and flow separations on the cowl surface.

The Dirac-Milne cosmology

NASA Astrophysics Data System (ADS)

Benoit-Lévy, Aurélien; Chardin, Gabriel

2014-05-01

We study an unconventional cosmology, in which we investigate the consequences that antigravity would pose to cosmology. We present the main characteristics of the Dirac-Milne Universe, a cosmological model where antimatter has a negative active gravitational mass. In this non-standard Universe, separate domains of matter and antimatter coexist at our epoch without annihilation, separated by a gravitationally induced depletion zone. We show that this cosmology does not require a priori the Dark Matter and Dark Energy components of the standard model of cosmology. Additionally, inflation becomes an unnecessary ingredient. Investigating this model, we show that the classical cosmological tests such as primordial nucleosynthesis, Type Ia supernovæ and Cosmic Microwave Background are surprisingly concordant.

BOOK REVIEW: Observational Cosmology Observational Cosmology

NASA Astrophysics Data System (ADS)

Howell, Dale Andrew

2013-04-01

Observational Cosmology by Stephen Serjeant fills a niche that was underserved in the textbook market: an up-to-date, thorough cosmology textbook focused on observations, aimed at advanced undergraduates. Not everything about the book is perfect - some subjects get short shrift, in some cases jargon dominates, and there are too few exercises. Still, on the whole, the book is a welcome addition. For decades, the classic textbooks of cosmology have focused on theory. But for every Sunyaev-Zel'dovich effect there is a Butcher-Oemler effect; there are as many cosmological phenomena established by observations, and only explained later by theory, as there were predicted by theory and confirmed by observations. In fact, in the last decade, there has been an explosion of new cosmological findings driven by observations. Some are so new that you won't find them mentioned in books just a few years old. So it is not just refreshing to see a book that reflects the new realities of cosmology, it is vital, if students are to truly stay up on a field that has widened in scope considerably. Observational Cosmology is filled with full-color images, and graphs from the latest experiments. How exciting it is that we live in an era where satellites and large experiments have gathered so much data to reveal astounding details about the origin of the universe and its evolution. To have all the latest data gathered together and explained in one book will be a revelation to students. In fact, at times it was to me. I've picked up modern cosmological knowledge through a patchwork of reading papers, going to colloquia, and serving on grant and telescope allocation panels. To go back and see them explained from square one, and summarized succinctly, filled in quite a few gaps in my own knowledge and corrected a few misconceptions I'd acquired along the way. To make room for all these graphs and observational details, a few things had to be left out. For one, there are few derivations

Supernova Cosmology Inference with Probabilistic Photometric Redshifts (SCIPPR)

NASA Astrophysics Data System (ADS)

Peters, Christina; Malz, Alex; Hlozek, Renée

2018-01-01

The Bayesian Estimation Applied to Multiple Species (BEAMS) framework employs probabilistic supernova type classifications to do photometric SN cosmology. This work extends BEAMS to replace high-confidence spectroscopic redshifts with photometric redshift probability density functions, a capability that will be essential in the era the Large Synoptic Survey Telescope and other next-generation photometric surveys where it will not be possible to perform spectroscopic follow up on every SN. We present the Supernova Cosmology Inference with Probabilistic Photometric Redshifts (SCIPPR) Bayesian hierarchical model for constraining the cosmological parameters from photometric lightcurves and host galaxy photometry, which includes selection effects and is extensible to uncertainty in the redshift-dependent supernova type proportions. We create a pair of realistic mock catalogs of joint posteriors over supernova type, redshift, and distance modulus informed by photometric supernova lightcurves and over redshift from simulated host galaxy photometry. We perform inference under our model to obtain a joint posterior probability distribution over the cosmological parameters and compare our results with other methods, namely: a spectroscopic subset, a subset of high probability photometrically classified supernovae, and reducing the photometric redshift probability to a single measurement and error bar.

Realistic Simulation for Body Area and Body-To-Body Networks.

PubMed

Alam, Muhammad Mahtab; Ben Hamida, Elyes; Ben Arbia, Dhafer; Maman, Mickael; Mani, Francesco; Denis, Benoit; D'Errico, Raffaele

2016-04-20

In this paper, we present an accurate and realistic simulation for body area networks (BAN) and body-to-body networks (BBN) using deterministic and semi-deterministic approaches. First, in the semi-deterministic approach, a real-time measurement campaign is performed, which is further characterized through statistical analysis. It is able to generate link-correlated and time-varying realistic traces (i.e., with consistent mobility patterns) for on-body and body-to-body shadowing and fading, including body orientations and rotations, by means of stochastic channel models. The full deterministic approach is particularly targeted to enhance IEEE 802.15.6 proposed channel models by introducing space and time variations (i.e., dynamic distances) through biomechanical modeling. In addition, it helps to accurately model the radio link by identifying the link types and corresponding path loss factors for line of sight (LOS) and non-line of sight (NLOS). This approach is particularly important for links that vary over time due to mobility. It is also important to add that the communication and protocol stack, including the physical (PHY), medium access control (MAC) and networking models, is developed for BAN and BBN, and the IEEE 802.15.6 compliance standard is provided as a benchmark for future research works of the community. Finally, the two approaches are compared in terms of the successful packet delivery ratio, packet delay and energy efficiency. The results show that the semi-deterministic approach is the best option; however, for the diversity of the mobility patterns and scenarios applicable, biomechanical modeling and the deterministic approach are better choices.

Realistic Simulation for Body Area and Body-To-Body Networks

PubMed Central

Alam, Muhammad Mahtab; Ben Hamida, Elyes; Ben Arbia, Dhafer; Maman, Mickael; Mani, Francesco; Denis, Benoit; D’Errico, Raffaele

2016-01-01

In this paper, we present an accurate and realistic simulation for body area networks (BAN) and body-to-body networks (BBN) using deterministic and semi-deterministic approaches. First, in the semi-deterministic approach, a real-time measurement campaign is performed, which is further characterized through statistical analysis. It is able to generate link-correlated and time-varying realistic traces (i.e., with consistent mobility patterns) for on-body and body-to-body shadowing and fading, including body orientations and rotations, by means of stochastic channel models. The full deterministic approach is particularly targeted to enhance IEEE 802.15.6 proposed channel models by introducing space and time variations (i.e., dynamic distances) through biomechanical modeling. In addition, it helps to accurately model the radio link by identifying the link types and corresponding path loss factors for line of sight (LOS) and non-line of sight (NLOS). This approach is particularly important for links that vary over time due to mobility. It is also important to add that the communication and protocol stack, including the physical (PHY), medium access control (MAC) and networking models, is developed for BAN and BBN, and the IEEE 802.15.6 compliance standard is provided as a benchmark for future research works of the community. Finally, the two approaches are compared in terms of the successful packet delivery ratio, packet delay and energy efficiency. The results show that the semi-deterministic approach is the best option; however, for the diversity of the mobility patterns and scenarios applicable, biomechanical modeling and the deterministic approach are better choices. PMID:27104537

Emergent cosmology revisited

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bag, Satadru; Sahni, Varun; Shtanov, Yuri

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

Higgs cosmology.

PubMed

Rajantie, Arttu

2018-03-06

The discovery of the Higgs boson in 2012 and other results from the Large Hadron Collider have confirmed the standard model of particle physics as the correct theory of elementary particles and their interactions up to energies of several TeV. Remarkably, the theory may even remain valid all the way to the Planck scale of quantum gravity, and therefore it provides a solid theoretical basis for describing the early Universe. Furthermore, the Higgs field itself has unique properties that may have allowed it to play a central role in the evolution of the Universe, from inflation to cosmological phase transitions and the origin of both baryonic and dark matter, and possibly to determine its ultimate fate through the electroweak vacuum instability. These connections between particle physics and cosmology have given rise to a new and growing field of Higgs cosmology, which promises to shed new light on some of the most puzzling questions about the Universe as new data from particle physics experiments and cosmological observations become available.This article is part of the Theo Murphy meeting issue 'Higgs cosmology'. © 2018 The Author(s).

Cosmology from a gauge induced gravity

NASA Astrophysics Data System (ADS)

Falciano, F. T.; Sadovski, G.; Sobreiro, R. F.; Tomaz, A. A.

2017-09-01

The main goal of the present work is to analyze the cosmological scenario of the induced gravity theory developed in previous works. Such a theory consists on a Yang-Mills theory in a four-dimensional Euclidian spacetime with { SO}(m,n) such that m+n=5 and m\\in {0,1,2} as its gauge group. This theory undergoes a dynamical gauge symmetry breaking via an Inönü-Wigner contraction in its infrared sector. As a consequence, the { SO}(m,n) algebra is deformed into a Lorentz algebra with the emergency of the local Lorentz symmetries and the gauge fields being identified with a vierbein and a spin connection. As a result, gravity is described as an effective Einstein-Cartan-like theory with ultraviolet correction terms and a propagating torsion field. We show that the cosmological model associated with this effective theory has three different regimes. In particular, the high curvature regime presents a de Sitter phase which tends towards a Λ CDM model. We argue that { SO}(m,n) induced gravities are promising effective theories to describe the early phase of the universe.

Cosmological evolution as squeezing: a toy model for group field cosmology

NASA Astrophysics Data System (ADS)

Adjei, Eugene; Gielen, Steffen; Wieland, Wolfgang

2018-05-01

We present a simple model of quantum cosmology based on the group field theory (GFT) approach to quantum gravity. The model is formulated on a subspace of the GFT Fock space for the quanta of geometry, with a fixed volume per quantum. In this Hilbert space, cosmological expansion corresponds to the generation of new quanta. Our main insight is that the evolution of a flat Friedmann–Lemaître–Robertson–Walker universe with a massless scalar field can be described on this Hilbert space as squeezing, familiar from quantum optics. As in GFT cosmology, we find that the three-volume satisfies an effective Friedmann equation similar to the one of loop quantum cosmology, connecting the classical contracting and expanding solutions by a quantum bounce. The only free parameter in the model is identified with Newton’s constant. We also comment on the possible topological interpretation of our squeezed states. This paper can serve as an introduction into the main ideas of GFT cosmology without requiring the full GFT formalism; our results can also motivate new developments in GFT and its cosmological application.

An experimental apparatus to simulate body-powered prosthetic usage: Development and preliminary evaluation.

PubMed

Gao, Fan; Rodriguez, Johanan; Kapp, Susan

2016-06-01

Harness fitting in the body-powered prosthesis remains more art than science due to a lack of consistent and quantitative evaluation. The aim of this study was to develop a mechanical, human-body-shaped apparatus to simulate body-powered upper limb prosthetic usage and evaluate its capability of quantitative examination of harness configuration. The apparatus was built upon a torso of a wooden mannequin and integrated major mechanical joints to simulate terminal device operation. Sensors were used to register cable tension, cable excursion, and grip force simultaneously. The apparatus allowed the scapula to move up to 127 mm laterally and the load cell can measure the cable tension up to 445 N. Our preliminary evaluation highlighted the needs and importance of investigating harness configurations in a systematic and controllable manner. The apparatus allows objective, systematic, and quantitative evaluation of effects of realistic harness configurations and will provide insightful and working knowledge on harness fitting in upper limb amputees using body-powered prosthesis. © The International Society for Prosthetics and Orthotics 2015.

Massive optimal data compression and density estimation for scalable, likelihood-free inference in cosmology

NASA Astrophysics Data System (ADS)

Alsing, Justin; Wandelt, Benjamin; Feeney, Stephen

2018-07-01

Many statistical models in cosmology can be simulated forwards but have intractable likelihood functions. Likelihood-free inference methods allow us to perform Bayesian inference from these models using only forward simulations, free from any likelihood assumptions or approximations. Likelihood-free inference generically involves simulating mock data and comparing to the observed data; this comparison in data space suffers from the curse of dimensionality and requires compression of the data to a small number of summary statistics to be tractable. In this paper, we use massive asymptotically optimal data compression to reduce the dimensionality of the data space to just one number per parameter, providing a natural and optimal framework for summary statistic choice for likelihood-free inference. Secondly, we present the first cosmological application of Density Estimation Likelihood-Free Inference (DELFI), which learns a parametrized model for joint distribution of data and parameters, yielding both the parameter posterior and the model evidence. This approach is conceptually simple, requires less tuning than traditional Approximate Bayesian Computation approaches to likelihood-free inference and can give high-fidelity posteriors from orders of magnitude fewer forward simulations. As an additional bonus, it enables parameter inference and Bayesian model comparison simultaneously. We demonstrate DELFI with massive data compression on an analysis of the joint light-curve analysis supernova data, as a simple validation case study. We show that high-fidelity posterior inference is possible for full-scale cosmological data analyses with as few as ˜104 simulations, with substantial scope for further improvement, demonstrating the scalability of likelihood-free inference to large and complex cosmological data sets.

Cosmological test using strong gravitational lensing systems

NASA Astrophysics Data System (ADS)

Yuan, C. C.; Wang, F. Y.

2015-09-01

As one of the probes of universe, strong gravitational lensing systems allow us to compare different cosmological models and constrain vital cosmological parameters. This purpose can be reached from the dynamic and geometry properties of strong gravitational lensing systems, for instance, time-delay Δτ of images, the velocity dispersion σ of the lensing galaxies and the combination of these two effects, Δτ/σ2. In this paper, in order to carry out one-on-one comparisons between ΛCDM universe and Rh = ct universe, we use a sample containing 36 strong lensing systems with the measurement of velocity dispersion from the Sloan Lens Advanced Camera for Surveys (SLACS) and Lens Structure and Dynamic survey (LSD) survey. Concerning the time-delay effect, 12 two-image lensing systems with Δτ are also used. In addition, Monte Carlo simulations are used to compare the efficiency of the three methods as mentioned above. From simulations, we estimate the number of lenses required to rule out one model at the 99.7 per cent confidence level. Comparing with constraints from Δτ and the velocity dispersion σ, we find that using Δτ/σ2 can improve the discrimination between cosmological models. Despite the independence tests of these methods reveal a correlation between Δτ/σ2 and σ, Δτ/σ2 could be considered as an improved method of σ if more data samples are available.

N-S/DSMC hybrid simulation of hypersonic flow over blunt body including wakes

NASA Astrophysics Data System (ADS)

Li, Zhonghua; Li, Zhihui; Li, Haiyan; Yang, Yanguang; Jiang, Xinyu

2014-12-01

A hybrid N-S/DSMC method is presented and applied to solve the three-dimensional hypersonic transitional flows by employing the MPC (modular Particle-Continuum) technique based on the N-S and the DSMC method. A sub-relax technique is adopted to deal with information transfer between the N-S and the DSMC. The hypersonic flows over a 70-deg spherically blunted cone under different Kn numbers are simulated using the CFD, DSMC and hybrid N-S/DSMC method. The present computations are found in good agreement with DSMC and experimental results. The present method provides an efficient way to predict the hypersonic aerodynamics in near-continuum transitional flow regime.

A whirling plane of satellite galaxies around Centaurus A challenges cold dark matter cosmology

NASA Astrophysics Data System (ADS)

Müller, Oliver; Pawlowski, Marcel S.; Jerjen, Helmut; Lelli, Federico

2018-02-01

The Milky Way and Andromeda galaxies are each surrounded by a thin plane of satellite dwarf galaxies that may be corotating. Cosmological simulations predict that most satellite galaxy systems are close to isotropic with random motions, so those two well-studied systems are often interpreted as rare statistical outliers. We test this assumption using the kinematics of satellite galaxies around the Centaurus A galaxy. Our statistical analysis reveals evidence for corotation in a narrow plane: Of the 16 Centaurus A satellites with kinematic data, 14 follow a coherent velocity pattern aligned with the long axis of their spatial distribution. In standard cosmological simulations, <0.5% of Centaurus A–like systems show such behavior. Corotating satellite systems may be common in the universe, challenging small-scale structure formation in the prevailing cosmological paradigm.

HI Fluctuations at Large Redshifts: III - Simulating the Signal Expected at GMRT

NASA Astrophysics Data System (ADS)

Bharadwaj, Somnath; Srikant, P. S.

2004-03-01

We simulate the distribution of neutral hydrogen (HI) at the redshifts z D 1:3 and 3:4 using a cosmological N-body simulation along with a prescription for assigning HI masses to the particles. The HI is distributed in clouds whose properties are consistent with those of the damped Lyman- absorption systems (DLAs) seen in quasar spectra. The clustering properties of these clouds are identical to those of the dark matter. We use this to simulate the redshifted HI emission expected at 610 MHzand 325 MHz, two of the observing bands at theGMRT. These are used to predict the correlations expected between the complex visibilities measured at different baselines and frequencies in radio-interferometric observations with the GMRT. The visibility correlations directly probe the power spectrum of HI fluctuations at the epoch when the HI emission originated, and this holds the possibility of using HI observations to study large-scale structures at high z.

The Atacama Cosmology Telescope: Dynamical Masses for 44 SZ-Selected Galaxy Clusters over 755 Square Degrees

NASA Technical Reports Server (NTRS)

Sifon, Cristobal; Battaglia, Nick; Hasselfield, Matthew; Menanteau, Felipe; Barrientos, L. Felipe; Bond, J. Richard; Crichton, Devin; Devlin, Mark J.; Dunner, Rolando; Hilton, Matt;

Dark Energy Survey Year 1 results: the impact of galaxy neighbours on weak lensing cosmology with IM3SHAPE

DOE PAGES

Samuroff, S.

2017-12-26

We use a suite of simulated images based on Year 1 of the Dark Energy Survey to explore the impact of galaxy neighbours on shape measurement and shear cosmology. The hoopoe image simulations include realistic blending, galaxy positions, and spatial variations in depth and PSF properties. Using the im3shape maximum-likelihood shape measurement code, we identify four mechanisms by which neighbours can have a non-negligible influence on shear estimation. These effects, if ignored, would contribute a net multiplicative bias ofmore » $$m \\sim 0.03 - 0.09$$ in the DES Y1 im3shape catalogue, though the precise impact will be dependent on both the measurement code and the selection cuts applied. This can be reduced to percentage level or less by removing objects with close neighbours, at a cost to the effective number density of galaxies $$n_\\mathrm{eff}$$ of 30%. We use the cosmological inference pipeline of DES Y1 to explore the cosmological implications of neighbour bias and show that omitting blending from the calibration simulation for DES Y1 would bias the inferred clustering amplitude $$S_8\\equiv \\sigma_8 (\\Omega _\\mathrm{m} /0.3)^{0.5}$$ by $$2 \\sigma$$ towards low values. Lastly, we use the hoopoe simulations to test the effect of neighbour-induced spatial correlations in the multiplicative bias. We find the impact on the recovered $$S_8$$ of ignoring such correlations to be subdominant to statistical error at the current level of precision.« less

Dark Energy Survey Year 1 results: the impact of galaxy neighbours on weak lensing cosmology with IM3SHAPE

DOE Office of Scientific and Technical Information (OSTI.GOV)

Samuroff, S.

We use a suite of simulated images based on Year 1 of the Dark Energy Survey to explore the impact of galaxy neighbours on shape measurement and shear cosmology. The hoopoe image simulations include realistic blending, galaxy positions, and spatial variations in depth and PSF properties. Using the im3shape maximum-likelihood shape measurement code, we identify four mechanisms by which neighbours can have a non-negligible influence on shear estimation. These effects, if ignored, would contribute a net multiplicative bias ofmore » $$m \\sim 0.03 - 0.09$$ in the DES Y1 im3shape catalogue, though the precise impact will be dependent on both the measurement code and the selection cuts applied. This can be reduced to percentage level or less by removing objects with close neighbours, at a cost to the effective number density of galaxies $$n_\\mathrm{eff}$$ of 30%. We use the cosmological inference pipeline of DES Y1 to explore the cosmological implications of neighbour bias and show that omitting blending from the calibration simulation for DES Y1 would bias the inferred clustering amplitude $$S_8\\equiv \\sigma_8 (\\Omega _\\mathrm{m} /0.3)^{0.5}$$ by $$2 \\sigma$$ towards low values. Lastly, we use the hoopoe simulations to test the effect of neighbour-induced spatial correlations in the multiplicative bias. We find the impact on the recovered $$S_8$$ of ignoring such correlations to be subdominant to statistical error at the current level of precision.« less

PoMiN: A Post-Minkowskian N-Body Solver

NASA Astrophysics Data System (ADS)

Feng, Justin; Baumann, Mark; Hall, Bryton; Doss, Joel; Spencer, Lucas; Matzner, Richard

2018-05-01

PoMiN is a lightweight N-body code based on the Post-Minkowskian N-body Hamiltonian of Ledvinka, Schafer, and Bicak, which includes General Relativistic effects up to first order in Newton's constant G, and all orders in the speed of light c. PoMiN is a single file written in C and uses a fourth-order Runge-Kutta integration scheme. PoMiN has also been written to handle an arbitrary number of particles (both massive and massless) with a computational complexity that scales as O(N^2).

Computer simulation of plasma and N-body problems

NASA Technical Reports Server (NTRS)

Harries, W. L.; Miller, J. B.

1975-01-01

The following FORTRAN language computer codes are presented: (1) efficient two- and three-dimensional central force potential solvers; (2) a three-dimensional simulator of an isolated galaxy which incorporates the potential solver; (3) a two-dimensional particle-in-cell simulator of the Jeans instability in an infinite self-gravitating compressible gas; and (4) a two-dimensional particle-in-cell simulator of a rotating self-gravitating compressible gaseous system of which rectangular coordinate and superior polar coordinate versions were written.

Cluster cosmology with next-generation surveys.

NASA Astrophysics Data System (ADS)

Ascaso, B.

2017-03-01

The advent of next-generation surveys will provide a large number of cluster detections that will serve the basis for constraining cos mological parameters using cluster counts. The main two observational ingredients needed are the cluster selection function and the calibration of the mass-observable relation. In this talk, we present the methodology designed to obtain robust predictions of both ingredients based on realistic cosmological simulations mimicking the following next-generation surveys: J-PAS, LSST and Euclid. We display recent results on the selection functions for these mentioned surveys together with others coming from other next-generation surveys such as eROSITA, ACTpol and SPTpol. We notice that the optical and IR surveys will reach the lowest masses between 0.3simulations, obtaining similar scatter to other observational results limited to higher redshifts. Finally, we describe the technique that we are developing to perform a Fisher Matrix analysis to provide cosmological constraints for the considered next-generation surveys and introduce very preliminary results.

Stability analysis in tachyonic potential chameleon cosmology

DOE Office of Scientific and Technical Information (OSTI.GOV)

Farajollahi, H.; Salehi, A.; Tayebi, F.

2011-05-01

We study general properties of attractors for tachyonic potential chameleon scalar-field model which possess cosmological scaling solutions. An analytic formulation is given to obtain fixed points with a discussion on their stability. The model predicts a dynamical equation of state parameter with phantom crossing behavior for an accelerating universe. We constrain the parameters of the model by best fitting with the recent data-sets from supernovae and simulated data points for redshift drift experiment generated by Monte Carlo simulations.

A METHOD TO EXTRACT THE REDSHIFT DISTORTION {beta} PARAMETER IN CONFIGURATION SPACE FROM MINIMAL COSMOLOGICAL ASSUMPTIONS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tocchini-Valentini, Domenico; Barnard, Michael; Bennett, Charles L.

2012-10-01

We present a method to extract the redshift-space distortion {beta} parameter in configuration space with a minimal set of cosmological assumptions. We show that a novel combination of the observed monopole and quadrupole correlation functions can remove efficiently the impact of mild nonlinearities and redshift errors. The method offers a series of convenient properties: it does not depend on the theoretical linear correlation function, the mean galaxy density is irrelevant, only convolutions are used, and there is no explicit dependence on linear bias. Analyses based on dark matter N-body simulations and Fisher matrix demonstrate that errors of a few percentmore » on {beta} are possible with a full-sky, 1 (h {sup -1} Gpc){sup 3} survey centered at a redshift of unity and with negligible shot noise. We also find a baryonic feature in the normalized quadrupole in configuration space that should complicate the extraction of the growth parameter from the linear theory asymptote, but that does not have a major impact on our method.« less

BMS in cosmology

NASA Astrophysics Data System (ADS)

Kehagias, A.; Riotto, A.

2016-05-01

Symmetries play an interesting role in cosmology. They are useful in characterizing the cosmological perturbations generated during inflation and lead to consistency relations involving the soft limit of the statistical correlators of large-scale structure dark matter and galaxies overdensities. On the other hand, in observational cosmology the carriers of the information about these large-scale statistical distributions are light rays traveling on null geodesics. Motivated by this simple consideration, we study the structure of null infinity and the associated BMS symmetry in a cosmological setting. For decelerating Friedmann-Robertson-Walker backgrounds, for which future null infinity exists, we find that the BMS transformations which leaves the asymptotic metric invariant to leading order. Contrary to the asymptotic flat case, the BMS transformations in cosmology generate Goldstone modes corresponding to scalar, vector and tensor degrees of freedom which may exist at null infinity and perturb the asymptotic data. Therefore, BMS transformations generate physically inequivalent vacua as they populate the universe at null infinity with these physical degrees of freedom. We also discuss the gravitational memory effect when cosmological expansion is taken into account. In this case, there are extra contribution to the gravitational memory due to the tail of the retarded Green functions which are supported not only on the light-cone, but also in its interior. The gravitational memory effect can be understood also from an asymptotic point of view as a transition among cosmological BMS-related vacua.

BMS in cosmology

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kehagias, A.; Riotto, A.; Center for Astroparticle Physics

Symmetries play an interesting role in cosmology. They are useful in characterizing the cosmological perturbations generated during inflation and lead to consistency relations involving the soft limit of the statistical correlators of large-scale structure dark matter and galaxies overdensities. On the other hand, in observational cosmology the carriers of the information about these large-scale statistical distributions are light rays traveling on null geodesics. Motivated by this simple consideration, we study the structure of null infinity and the associated BMS symmetry in a cosmological setting. For decelerating Friedmann-Robertson-Walker backgrounds, for which future null infinity exists, we find that the BMS transformationsmore » which leaves the asymptotic metric invariant to leading order. Contrary to the asymptotic flat case, the BMS transformations in cosmology generate Goldstone modes corresponding to scalar, vector and tensor degrees of freedom which may exist at null infinity and perturb the asymptotic data. Therefore, BMS transformations generate physically inequivalent vacua as they populate the universe at null infinity with these physical degrees of freedom. We also discuss the gravitational memory effect when cosmological expansion is taken into account. In this case, there are extra contribution to the gravitational memory due to the tail of the retarded Green functions which are supported not only on the light-cone, but also in its interior. The gravitational memory effect can be understood also from an asymptotic point of view as a transition among cosmological BMS-related vacua.« less

Philosophical Roots of Cosmology

NASA Astrophysics Data System (ADS)

Ivanovic, M.

2008-10-01

We shall consider the philosophical roots of cosmology in the earlier Greek philosophy. Our goal is to answer the question: Are earlier Greek theories of pure philosophical-mythological character, as often philosophers cited it, or they have scientific character. On the bases of methodological criteria, we shall contend that the latter is the case. In order to answer the question about contemporary situation of the relation philosophy-cosmology, we shall consider the next question: Is contemporary cosmology completely independent of philosophical conjectures? The answer demands consideration of methodological character about scientific status of contemporary cosmology. We also consider some aspects of the relation contemporary philosophy-cosmology.

Did I Say Cosmology? On Modern Cosmologies and Ancient World-views

NASA Astrophysics Data System (ADS)

Iwaniszewski, S.

2009-08-01

The modern cosmology that emerged from observational astronomy in 16th century Europe meant a radical break-away from earlier conceptions of the world. While all ancient and nonwestern worldviews usually describe a multidimensional reality in which diverse environmental, economic, sociopolitical and ideological factors intersect, modern cosmologies espouse the vision of a radically different universe which is completely dehumanized, ethically indifferent and universally valid. Despite these differences cosmology and worldview tend to be used interchangeably to depict ancient and nonwestern worldviews.Any correspondences which can be found between different parts of ancient and/or nonwestern worldviews and modern cosmologies tend to transfer modern conceptions to the premodern world. Ignoring ancient cultural contexts, we risk imposing modern cosmological concepts on past worldview categories. While we have to describe ancient astronomies in our own terms, our ultimate goal is to understand them on their own terms.

Second-order Cosmological Perturbations Engendered by Point-like Masses

DOE Office of Scientific and Technical Information (OSTI.GOV)

Brilenkov, Ruslan; Eingorn, Maxim, E-mail: ruslan.brilenkov@gmail.com, E-mail: maxim.eingorn@gmail.com

2017-08-20

In the ΛCDM framework, presenting nonrelativistic matter inhomogeneities as discrete massive particles, we develop the second‐order cosmological perturbation theory. Our approach relies on the weak gravitational field limit. The derived equations for the second‐order scalar, vector, and tensor metric corrections are suitable at arbitrary distances, including regions with nonlinear contrasts of the matter density. We thoroughly verify fulfillment of all Einstein equations, as well as self‐consistency of order assignments. In addition, we achieve logical positive results in the Minkowski background limit. Feasible investigations of the cosmological back-reaction manifestations by means of relativistic simulations are also outlined.

Sample variance in weak lensing: How many simulations are required?

DOE Office of Scientific and Technical Information (OSTI.GOV)

Petri, Andrea; May, Morgan; Haiman, Zoltan

Constraining cosmology using weak gravitational lensing consists of comparing a measured feature vector of dimension N b with its simulated counterpart. An accurate estimate of the N b × N b feature covariance matrix C is essential to obtain accurate parameter confidence intervals. When C is measured from a set of simulations, an important question is how large this set should be. To answer this question, we construct different ensembles of N r realizations of the shear field, using a common randomization procedure that recycles the outputs from a smaller number N s ≤ N r of independent ray-tracing N-bodymore » simulations. We study parameter confidence intervals as a function of (N s, N r) in the range 1 ≤ N s ≤ 200 and 1 ≤ N r ≲ 105. Previous work [S. Dodelson and M. D. Schneider, Phys. Rev. D 88, 063537 (2013)] has shown that Gaussian noise in the feature vectors (from which the covariance is estimated) lead, at quadratic order, to an O(1/N r) degradation of the parameter confidence intervals. Using a variety of lensing features measured in our simulations, including shear-shear power spectra and peak counts, we show that cubic and quartic covariance fluctuations lead to additional O(1/N 2 r) error degradation that is not negligible when N r is only a factor of few larger than N b. We study the large N r limit, and find that a single, 240 Mpc/h sized 512 3-particle N-body simulation (N s = 1) can be repeatedly recycled to produce as many as N r = few × 10 4 shear maps whose power spectra and high-significance peak counts can be treated as statistically independent. Lastly, a small number of simulations (N s = 1 or 2) is sufficient to forecast parameter confidence intervals at percent accuracy.« less

Sample variance in weak lensing: How many simulations are required?

DOE PAGES

Petri, Andrea; May, Morgan; Haiman, Zoltan

2016-03-24

Constraining cosmology using weak gravitational lensing consists of comparing a measured feature vector of dimension N b with its simulated counterpart. An accurate estimate of the N b × N b feature covariance matrix C is essential to obtain accurate parameter confidence intervals. When C is measured from a set of simulations, an important question is how large this set should be. To answer this question, we construct different ensembles of N r realizations of the shear field, using a common randomization procedure that recycles the outputs from a smaller number N s ≤ N r of independent ray-tracing N-bodymore » simulations. We study parameter confidence intervals as a function of (N s, N r) in the range 1 ≤ N s ≤ 200 and 1 ≤ N r ≲ 105. Previous work [S. Dodelson and M. D. Schneider, Phys. Rev. D 88, 063537 (2013)] has shown that Gaussian noise in the feature vectors (from which the covariance is estimated) lead, at quadratic order, to an O(1/N r) degradation of the parameter confidence intervals. Using a variety of lensing features measured in our simulations, including shear-shear power spectra and peak counts, we show that cubic and quartic covariance fluctuations lead to additional O(1/N 2 r) error degradation that is not negligible when N r is only a factor of few larger than N b. We study the large N r limit, and find that a single, 240 Mpc/h sized 512 3-particle N-body simulation (N s = 1) can be repeatedly recycled to produce as many as N r = few × 10 4 shear maps whose power spectra and high-significance peak counts can be treated as statistically independent. Lastly, a small number of simulations (N s = 1 or 2) is sufficient to forecast parameter confidence intervals at percent accuracy.« less

Cosmological perturbations in teleparallel Loop Quantum Cosmology

DOE Office of Scientific and Technical Information (OSTI.GOV)

Haro, Jaime, E-mail: jaime.haro@upc.edu

2013-11-01

Cosmological perturbations in Loop Quantum Cosmology (LQC) are usually studied incorporating either holonomy corrections, where the Ashtekar connection is replaced by a suitable sinus function in order to have a well-defined quantum analogue, or inverse-volume corrections coming from the eigenvalues of the inverse-volume operator. In this paper we will develop an alternative approach to calculate cosmological perturbations in LQC based on the fact that, holonomy corrected LQC in the flat Friedmann-Lemaître-Robertson-Walker (FLRW) geometry could be also obtained as a particular case of teleparallel F(T) gravity (teleparallel LQC). The main idea of our approach is to mix the simple bounce providedmore » by holonomy corrections in LQC with the non-singular perturbation equations given by F(T) gravity, in order to obtain a matter bounce scenario as a viable alternative to slow-roll inflation. In our study, we have obtained an scale invariant power spectrum of cosmological perturbations. However, the ratio of tensor to scalar perturbations is of order 1, which does not agree with the current observations. For this reason, we suggest a model where a transition from the matter domination to a quasi de Sitter phase is produced in order to enhance the scalar power spectrum.« less

One dark matter mystery: halos in the cosmic web

NASA Astrophysics Data System (ADS)

Gaite, Jose

2015-01-01

The current cold dark matter cosmological model explains the large scale cosmic web structure but is challenged by the observation of a relatively smooth distribution of matter in galactic clusters. We consider various aspects of modeling the dark matter around galaxies as distributed in smooth halos and, especially, the smoothness of the dark matter halos seen in N-body cosmological simulations. We conclude that the problems of the cold dark matter cosmology on small scales are more serious than normally admitted.

Computer simulation of multigrid body dynamics and control

NASA Technical Reports Server (NTRS)

Swaminadham, M.; Moon, Young I.; Venkayya, V. B.

1990-01-01

The objective is to set up and analyze benchmark problems on multibody dynamics and to verify the predictions of two multibody computer simulation codes. TREETOPS and DISCOS have been used to run three example problems - one degree-of-freedom spring mass dashpot system, an inverted pendulum system, and a triple pendulum. To study the dynamics and control interaction, an inverted planar pendulum with an external body force and a torsional control spring was modeled as a hinge connected two-rigid body system. TREETOPS and DISCOS affected the time history simulation of this problem. System state space variables and their time derivatives from two simulation codes were compared.

A hydrodynamic approach to cosmology - Methodology

NASA Technical Reports Server (NTRS)

Cen, Renyue

1992-01-01

The present study describes an accurate and efficient hydrodynamic code for evolving self-gravitating cosmological systems. The hydrodynamic code is a flux-based mesh code originally designed for engineering hydrodynamical applications. A variety of checks were performed which indicate that the resolution of the code is a few cells, providing accuracy for integral energy quantities in the present simulations of 1-3 percent over the whole runs. Six species (H I, H II, He I, He II, He III) are tracked separately, and relevant ionization and recombination processes, as well as line and continuum heating and cooling, are computed. The background radiation field is simultaneously determined in the range 1 eV to 100 keV, allowing for absorption, emission, and cosmological effects. It is shown how the inevitable numerical inaccuracies can be estimated and to some extent overcome.

Testing the conditional mass function of dark matter haloes against numerical N-body simulations

NASA Astrophysics Data System (ADS)

Tramonte, D.; Rubiño-Martín, J. A.; Betancort-Rijo, J.; Dalla Vecchia, C.

2017-05-01

We compare the predicted conditional mass function (CMF) of dark matter haloes from two theoretical prescriptions against numerical N-body simulations, both in overdense and underdense regions and at different Eulerian scales ranging from 5 to 30 h-1 Mpc. In particular, we consider in detail a locally implemented rescaling of the unconditional mass function (UMF) already discussed in the literature, and also a generalization of the standard rescaling method described in the extended Press-Schechter formalism. First, we test the consistency of these two rescalings by verifying the normalization of the CMF at different scales, and showing that none of the proposed cases provides a normalized CMF. In order to satisfy the normalization condition, we include a modification in the rescaling procedure. After this modification, the resulting CMF generally provides a better description of numerical results. We finally present an analytical fit to the ratio between the CMF and the UMF (also known as the matter-to-halo bias function) in underdense regions, which could be of special interest to speed up the computation of the halo abundance when studying void statistics. In this case, the CMF prescription based on the locally implemented rescaling provides a slightly better description of the numerical results when compared to the standard rescaling.

An order (n) algorithm for the dynamics simulation of robotic systems

NASA Technical Reports Server (NTRS)

Chun, H. M.; Turner, J. D.; Frisch, Harold P.

1989-01-01

The formulation of an Order (n) algorithm for DISCOS (Dynamics Interaction Simulation of Controls and Structures), which is an industry-standard software package for simulation and analysis of flexible multibody systems is presented. For systems involving many bodies, the new Order (n) version of DISCOS is much faster than the current version. Results of the experimental validation of the dynamics software are also presented. The experiment is carried out on a seven-joint robot arm at NASA's Goddard Space Flight Center. The algorithm used in the current version of DISCOS requires the inverse of a matrix whose dimension is equal to the number of constraints in the system. Generally, the number of constraints in a system is roughly proportional to the number of bodies in the system, and matrix inversion requires O(p exp 3) operations, where p is the dimension of the matrix. The current version of DISCOS is therefore considered an Order (n exp 3) algorithm. In contrast, the Order (n) algorithm requires inversion of matrices which are small, and the number of matrices to be inverted increases only linearly with the number of bodies. The newly-developed Order (n) DISCOS is currently capable of handling chain and tree topologies as well as multiple closed loops. Continuing development will extend the capability of the software to deal with typical robotics applications such as put-and-place, multi-arm hand-off and surface sliding.

Cosmology and particle physics

NASA Technical Reports Server (NTRS)

Turner, Michael S.

1988-01-01

The interplay between cosmology and elementary particle physics is discussed. The standard cosmology is reviewed, concentrating on primordial nucleosynthesis and discussing how the standard cosmology has been used to place constraints on the properties of various particles. Baryogenesis is discussed, showing how a scenario in which the B-, C-, and CP-violating interactions in GUTs provide a dynamical explanation for the predominance of matter over antimatter and for the present baryon-to-photon ratio. It is shown how the very early dynamical evolution of a very weakly coupled scalar field which is initially displaced from the minimum of its potential may explain a handful of very fundamental cosmological facts which are not explained by the standard cosmology.

Robust three-body water simulation model

NASA Astrophysics Data System (ADS)

Tainter, C. J.; Pieniazek, P. A.; Lin, Y.-S.; Skinner, J. L.

2011-05-01

The most common potentials used in classical simulations of liquid water assume a pairwise additive form. Although these models have been very successful in reproducing many properties of liquid water at ambient conditions, none is able to describe accurately water throughout its complicated phase diagram. The primary reason for this is the neglect of many-body interactions. To this end, a simulation model with explicit three-body interactions was introduced recently [R. Kumar and J. L. Skinner, J. Phys. Chem. B 112, 8311 (2008), 10.1021/jp8009468]. This model was parameterized to fit the experimental O-O radial distribution function and diffusion constant. Herein we reparameterize the model, fitting to a wider range of experimental properties (diffusion constant, rotational correlation time, density for the liquid, liquid/vapor surface tension, melting point, and the ice Ih density). The robustness of the model is then verified by comparing simulation to experiment for a number of other quantities (enthalpy of vaporization, dielectric constant, Debye relaxation time, temperature of maximum density, and the temperature-dependent second and third virial coefficients), with good agreement.

Smoot Cosmology Group

Science.gov Websites

Links We bring the universe to you! University of California Berkeley Cosmology Group Lawrence Computational Cosmology Center Institute for Nuclear & Particle Astrophysics Supernova Acceleration Probe

Systematic Biases in Weak Lensing Cosmology with the Dark Energy Survey

DOE Office of Scientific and Technical Information (OSTI.GOV)

Samuroff, Simon

This thesis sets out a practical guide to applying shear measurements as a cosmological tool. We first present one of two science-ready galaxy shape catalogues from Year 1 of the Dark Energy Survey (DES Y1), which covers 1500 square degrees in four bandsmore » $griz$, with a median redshift of $0.59$. We describe the shape measurement process implemented by the DES Y1 imshape catalogue, which contains 21.9 million high-quality $r$-band bulge/disc fits. In Chapter 3 a new suite of image simulations, referred to as Hoopoe, are presented. The Hoopoe dataset is tailored to DES Y1 and includes realistic blending, spatial masks and variation in the point spread function. We derive shear corrections, which we show are robust to changes in calibration method, galaxy binning and variance within the simulated dataset. Sources of systematic uncertainty in the simulation-based shear calibration are discussed, leading to a final estimate of the $$1\\sigma$$ uncertainties in the residual multiplica tive bias after calibration of 0.025. Chapter 4 describes an extension of the analysis on the Hoopoe simulations into a detailed investigation of the impact of galaxy neighbours on shape measurement and shear cosmology. Four mechanisms by which neighbours can have a non-negligible influence on shear measurement are identified. These effects, if ignored, would contribute a net multiplicative bias of $$m \\sim 0.03 - 0.09$$ in DES Y1, though the precise impact will depend on both the measurement code and the selection cuts applied. We use the cosmological inference pipeline of DES Y1 to explore the cosmological implications of neighbour bias and show that omitting blending from the calibration simulation for DES Y1 would bias the inferred clustering amplitude $$S_8 \\equiv \\sigma_8 (\\omegam /0.3)^{0.5}$$ by $$1.5 \\sigma$$ towards low values. Finally, we use the Hoopoe simulations to test the effect of neighbour-induced spatial correlations in the multiplicative bias. We find the

An accurate halo model for fitting non-linear cosmological power spectra and baryonic feedback models

NASA Astrophysics Data System (ADS)

Mead, A. J.; Peacock, J. A.; Heymans, C.; Joudaki, S.; Heavens, A. F.

2015-12-01

We present an optimized variant of the halo model, designed to produce accurate matter power spectra well into the non-linear regime for a wide range of cosmological models. To do this, we introduce physically motivated free parameters into the halo-model formalism and fit these to data from high-resolution N-body simulations. For a variety of Λ cold dark matter (ΛCDM) and wCDM models, the halo-model power is accurate to ≃ 5 per cent for k ≤ 10h Mpc-1 and z ≤ 2. An advantage of our new halo model is that it can be adapted to account for the effects of baryonic feedback on the power spectrum. We demonstrate this by fitting the halo model to power spectra from the OWLS (OverWhelmingly Large Simulations) hydrodynamical simulation suite via parameters that govern halo internal structure. We are able to fit all feedback models investigated at the 5 per cent level using only two free parameters, and we place limits on the range of these halo parameters for feedback models investigated by the OWLS simulations. Accurate predictions to high k are vital for weak-lensing surveys, and these halo parameters could be considered nuisance parameters to marginalize over in future analyses to mitigate uncertainty regarding the details of feedback. Finally, we investigate how lensing observables predicted by our model compare to those from simulations and from HALOFIT for a range of k-cuts and feedback models and quantify the angular scales at which these effects become important. Code to calculate power spectra from the model presented in this paper can be found at https://github.com/alexander-mead/hmcode.

Cosmological constraints from the convergence 1-point probability distribution

NASA Astrophysics Data System (ADS)

Patton, Kenneth; Blazek, Jonathan; Honscheid, Klaus; Huff, Eric; Melchior, Peter; Ross, Ashley J.; Suchyta, Eric

2017-11-01

We examine the cosmological information available from the 1-point probability density function (PDF) of the weak-lensing convergence field, utilizing fast L-PICOLA simulations and a Fisher analysis. We find competitive constraints in the Ωm-σ8 plane from the convergence PDF with 188 arcmin2 pixels compared to the cosmic shear power spectrum with an equivalent number of modes (ℓ < 886). The convergence PDF also partially breaks the degeneracy cosmic shear exhibits in that parameter space. A joint analysis of the convergence PDF and shear 2-point function also reduces the impact of shape measurement systematics, to which the PDF is less susceptible, and improves the total figure of merit by a factor of 2-3, depending on the level of systematics. Finally, we present a correction factor necessary for calculating the unbiased Fisher information from finite differences using a limited number of cosmological simulations.

Dynamical analysis on f(R, G) cosmology

NASA Astrophysics Data System (ADS)

Santos da Costa, S.; Roig, F. V.; Alcaniz, J. S.; Capozziello, S.; De Laurentis, M.; Benetti, M.

2018-04-01

We use a dynamical system approach to study the cosmological viability of f(R, G) gravity theories. The method consists of formulating the evolution equations as an autonomous system of ordinary differential equations, using suitable variables. The formalism is applied to a class of models in which f(R, G)\\propto RnG1-n and its solutions and corresponding stability are analysed in detail. New accelerating solutions that can be attractors in the phase space are found. We also find that this class of models does not exhibit a matter-dominated epoch, a solution which is inconsistent with current cosmological observations.

Field Theoretical Methods in Cosmology

NASA Astrophysics Data System (ADS)

Singh, Anupam

1995-01-01

To optimally utilize all the exciting cosmological data coming in we need to sharpen also the theoretical tools available to cosmologists. One such indispensible tool to understand hot big bang cosmology is finite temperature field theory. We review and summarise the efforts made by us to use finite temperature field theory to address issues of current interest to cosmologists. An introduction to both the real time and the imaginary time formalisms is provided. The imaginary time formalism is illustrated by applying it to understand the interesting possibility of late Time Phase Transitions. Recent observations of the space distribution of quasars indicate a very notable peak in space density at a redshift of 2 to 3. It is pointed out that this may be the result of a phase transition which has a critical temperature of roughly a few meV (in the cosmological units, h = c = k = 1), which is natural in the context of massive neutrinos. In fact, the neutrino masses required for quasar production and those required to solve the solar neutrino problem by the MSW mechanism are consistent with each other. As a bonus, the cosmological constant implied by this model may also help resolve the discrepancy between the recently measured value of the Hubble Constant and the age of the universe. We illustrate the real time formalism by studying one of the most important time-dependent and non-equilibrium phenomena associated with phase transitions. The non-equilibrium dynamics of the first stage of the reheating process, that is dissipation via particle production is studied in scalar field theories. We show that a complete understanding of the mechanism of dissipation via particle production requires a non-perturbative resummation. We then study a Hartree approximation and clearly exhibit dissipative effects related to particle production. The effect of dissipation by Goldstone bosons is studied non-perturbatively in the large N limit in an O(N) theory. We also place our work in

THE INNER STRUCTURE OF DWARF-SIZED HALOS IN WARM AND COLD DARK MATTER COSMOLOGIES

DOE Office of Scientific and Technical Information (OSTI.GOV)

González-Samaniego, A.; Avila-Reese, V.; Colín, P.

2016-03-10

By means of N-body + hydrodynamic zoom-in simulations we study the evolution of the inner dark matter and stellar mass distributions of central dwarf galaxies formed in halos of virial masses M{sub v} = (2–3) × 10{sup 10} h{sup −1} M{sub ⊙} at z = 0, both in a warm dark matter (WDM) and cold dark matter (CDM) cosmology. The half-mode mass in the WDM power spectrum of our simulations is M{sub f} = 2 × 10{sup 10} h{sup −1} M{sub ⊙}. In the dark matter (DM) only simulations halo density profiles are well described by the Navarro–Frenk–White parametric fit in both cosmologies, though the WDM halos have concentrations lower bymore » factors of 1.5–2.0 than their CDM counterparts. In the hydrodynamic simulations, the effects of baryons significantly flatten the inner density, velocity dispersion, and pseudo phase space density profiles of the WDM halos but not of the CDM ones. The density slope, measured at ≈0.02R{sub v}, α{sub 0.02}, becomes shallow in periods of 2–5 Gyr in the WDM runs. We explore whether this flattening process correlates with the global star formation (SF), M{sub s}/M{sub v} ratio, gas outflow, and internal specific angular momentum histories. We do not find any clear trends, but when α{sub 0.02} is shallower than −0.5, M{sub s}/M{sub v} is always between 0.25% and 1%. We conclude that the main reason for the formation of the shallow core is the presence of strong gas mass fluctuations inside the inner halo, which are a consequence of the feedback driven by a very bursty and sustained SF history in shallow gravitational potentials. Our WDM halos, which assemble late and are less concentrated than the CDM ones, obey these conditions. There are also (rare) CDM systems with extended mass assembly histories that obey these conditions and form shallow cores. The dynamical heating and expansion processes behind the DM core flattening apply also to the stars in such a way that the stellar age and metallicity gradients

An Advanced N -body Model for Interacting Multiple Stellar Systems

DOE Office of Scientific and Technical Information (OSTI.GOV)

Brož, Miroslav

We construct an advanced model for interacting multiple stellar systems in which we compute all trajectories with a numerical N -body integrator, namely the Bulirsch–Stoer from the SWIFT package. We can then derive various observables: astrometric positions, radial velocities, minima timings (TTVs), eclipse durations, interferometric visibilities, closure phases, synthetic spectra, spectral energy distribution, and even complete light curves. We use a modified version of the Wilson–Devinney code for the latter, in which the instantaneous true phase and inclination of the eclipsing binary are governed by the N -body integration. If all of these types of observations are at one’s disposal,more » a joint χ {sup 2} metric and an optimization algorithm (a simplex or simulated annealing) allow one to search for a global minimum and construct very robust models of stellar systems. At the same time, our N -body model is free from artifacts that may arise if mutual gravitational interactions among all components are not self-consistently accounted for. Finally, we present a number of examples showing dynamical effects that can be studied with our code and we discuss how systematic errors may affect the results (and how to prevent this from happening).« less

^{7}Be(n,α)^{4}He Reaction and the Cosmological Lithium Problem: Measurement of the Cross Section in a Wide Energy Range at n_TOF at CERN.

PubMed

Barbagallo, M; Musumarra, A; Cosentino, L; Maugeri, E; Heinitz, S; Mengoni, A; Dressler, R; Schumann, D; Käppeler, F; Colonna, N; Finocchiaro, P; Ayranov, M; Damone, L; Kivel, N; Aberle, O; Altstadt, S; Andrzejewski, J; Audouin, L; Bacak, M; Balibrea-Correa, J; Barros, S; Bécares, V; Bečvář, F; Beinrucker, C; Berthoumieux, E; Billowes, J; Bosnar, D; Brugger, M; Caamaño, M; Calviani, M; Calviño, F; Cano-Ott, D; Cardella, R; Casanovas, A; Castelluccio, D M; Cerutti, F; Chen, Y H; Chiaveri, E; Cortés, G; Cortés-Giraldo, M A; Cristallo, S; Diakaki, M; Domingo-Pardo, C; Dupont, E; Duran, I; Fernandez-Dominguez, B; Ferrari, A; Ferreira, P; Furman, W; Ganesan, S; García-Rios, A; Gawlik, A; Glodariu, T; Göbel, K; Gonçalves, I F; González-Romero, E; Griesmayer, E; Guerrero, C; Gunsing, F; Harada, H; Heftrich, T; Heyse, J; Jenkins, D G; Jericha, E; Katabuchi, T; Kavrigin, P; Kimura, A; Kokkoris, M; Krtička, M; Leal-Cidoncha, E; Lerendegui, J; Lederer, C; Leeb, H; Lo Meo, S; Lonsdale, S J; Losito, R; Macina, D; Marganiec, J; Martínez, T; Massimi, C; Mastinu, P; Mastromarco, M; Mazzone, A; Mendoza, E; Milazzo, P M; Mingrone, F; Mirea, M; Montesano, S; Nolte, R; Oprea, A; Pappalardo, A; Patronis, N; Pavlik, A; Perkowski, J; Piscopo, M; Plompen, A; Porras, I; Praena, J; Quesada, J; Rajeev, K; Rauscher, T; Reifarth, R; Riego-Perez, A; Rout, P; Rubbia, C; Ryan, J; Sabate-Gilarte, M; Saxena, A; Schillebeeckx, P; Schmidt, S; Sedyshev, P; Smith, A G; Stamatopoulos, A; Tagliente, G; Tain, J L; Tarifeño-Saldivia, A; Tassan-Got, L; Tsinganis, A; Valenta, S; Vannini, G; Variale, V; Vaz, P; Ventura, A; Vlachoudis, V; Vlastou, R; Vollaire, J; Wallner, A; Warren, S; Weigand, M; Weiß, C; Wolf, C; Woods, P J; Wright, T; Žugec, P

2016-10-07

The energy-dependent cross section of the ^{7}Be(n,α)^{4}He reaction, of interest for the so-called cosmological lithium problem in big bang nucleosynthesis, has been measured for the first time from 10 meV to 10 keV neutron energy. The challenges posed by the short half-life of ^{7}Be and by the low reaction cross section have been overcome at n_TOF thanks to an unprecedented combination of the extremely high luminosity and good resolution of the neutron beam in the new experimental area (EAR2) of the n_TOF facility at CERN, the availability of a sufficient amount of chemically pure ^{7}Be, and a specifically designed experimental setup. Coincidences between the two alpha particles have been recorded in two Si-^{7}Be-Si arrays placed directly in the neutron beam. The present results are consistent, at thermal neutron energy, with the only previous measurement performed in the 1960s at a nuclear reactor. The energy dependence reported here clearly indicates the inadequacy of the cross section estimates currently used in BBN calculations. Although new measurements at higher neutron energy may still be needed, the n_TOF results hint at a minor role of this reaction in BBN, leaving the long-standing cosmological lithium problem unsolved.

The Atacama Cosmology Telescope: A Measurement of the Thermal Sunyaev-Zel'dovich Effect Using the Skewness of the CMB Temperature Distribution

NASA Technical Reports Server (NTRS)

Wilson, Michael J.; Sherwin, Blake D.; Hill, J. Collin; Addison, Graeme; Battaglia, Nick; Bond, J. Richard; Das, Sudeep; Devlin, Mark J.; Dunkley, Joanna; Duenner, Rolando;

Holographic signatures of cosmological singularities.

PubMed

Engelhardt, Netta; Hertog, Thomas; Horowitz, Gary T

2014-09-19

To gain insight into the quantum nature of cosmological singularities, we study anisotropic Kasner solutions in gauge-gravity duality. The dual description of the bulk evolution towards the singularity involves N=4 super Yang-Mills theory on the expanding branch of deformed de Sitter space and is well defined. We compute two-point correlators of Yang-Mills operators of large dimensions using spacelike geodesics anchored on the boundary. The correlators show a strong signature of the singularity around horizon scales and decay at large boundary separation at different rates in different directions. More generally, the boundary evolution exhibits a process of particle creation similar to that in inflation. This leads us to conjecture that information on the quantum nature of cosmological singularities is encoded in long-wavelength features of the boundary wave function.

New central configurations of the (n + 1) -body problem

NASA Astrophysics Data System (ADS)

Fernandes, Antonio Carlos; Garcia, Braulio Augusto; Llibre, Jaume; Mello, Luis Fernando

2018-01-01

In this article we study central configurations of the (n + 1) -body problem. For the planar (n + 1) -body problem we study central configurations performed by n ≥ 2 bodies with equal masses at the vertices of a regular n-gon and one body with null mass. We also study spatial central configurations considering n bodies with equal masses at the vertices of a regular polyhedron and one body with null mass.

Supernova Cosmology Project

Science.gov Websites

, 2014 The Supernova Cosmology Project and High-Z Team share the 2015 Breakthrough Prize in Fundamental Perlmutter, leader of the international Supernova Cosmology Project, and principal investigator of the

Dark-ages Reionization and Galaxy Formation Simulation - XIV. Gas accretion, cooling, and star formation in dwarf galaxies at high redshift

NASA Astrophysics Data System (ADS)

Qin, Yuxiang; Duffy, Alan R.; Mutch, Simon J.; Poole, Gregory B.; Geil, Paul M.; Mesinger, Andrei; Wyithe, J. Stuart B.

2018-06-01

We study dwarf galaxy formation at high redshift (z ≥ 5) using a suite of high-resolution, cosmological hydrodynamic simulations and a semi-analytic model (SAM). We focus on gas accretion, cooling, and star formation in this work by isolating the relevant process from reionization and supernova feedback, which will be further discussed in a companion paper. We apply the SAM to halo merger trees constructed from a collisionless N-body simulation sharing identical initial conditions to the hydrodynamic suite, and calibrate the free parameters against the stellar mass function predicted by the hydrodynamic simulations at z = 5. By making comparisons of the star formation history and gas components calculated by the two modelling techniques, we find that semi-analytic prescriptions that are commonly adopted in the literature of low-redshift galaxy formation do not accurately represent dwarf galaxy properties in the hydrodynamic simulation at earlier times. We propose three modifications to SAMs that will provide more accurate high-redshift simulations. These include (1) the halo mass and baryon fraction which are overestimated by collisionless N-body simulations; (2) the star formation efficiency which follows a different cosmic evolutionary path from the hydrodynamic simulation; and (3) the cooling rate which is not well defined for dwarf galaxies at high redshift. Accurate semi-analytic modelling of dwarf galaxy formation informed by detailed hydrodynamical modelling will facilitate reliable semi-analytic predictions over the large volumes needed for the study of reionization.

7Be(n,α) and 7Be(n,p) cross-section measurement for the cosmological lithium problem at the n_TOF facility at CERN

NASA Astrophysics Data System (ADS)

Barbagallo, M.; Colonna, N.; Aberle, O.; Andrzejewski, J.; Audouin, L.; Bécares, V.; Bacak, M.; Balibrea, J.; Barros, S.; Bečvář, F.; Beinrucker, C.; Berthoumieux, E.; Billowes, J.; Bosnar, D.; Brugger, M.; Caamaño, M.; Calviño, F.; Calviani, M.; Cano-Ott, D.; Cardella, R.; Casanovas, A.; Castelluccio, D. M.; Cerutti, F.; Chen, Y. H.; Chiaveri, E.; Cortés, G.; Cortés-Giraldo, M. A.; Cosentino, L.; Damone, L. A.; Diakaki, M.; Domingo-Pardo, C.; Dressler, R.; Dupont, E.; Durán, I.; Fernández-Domínguez, B.; Ferrari, A.; Ferreira, P.; Finocchiaro, P.; Furman, V.; Göbel, K.; García, A. R.; Gawlik, A.; Glodariu, T.; Gonçalves, I. F.; González, E.; Goverdovski, A.; Griesmayer, E.; Guerrero, C.; Gunsing, F.; Harada, H.; Heftrich, T.; Heinitz, S.; Heyse, J.; Jenkins, D. G.; Jericha, E.; Käppeler, F.; Kadi, Y.; Katabuchi, T.; Kavrigin, P.; Ketlerov, V.; Khryachkov, V.; Kimura, A.; Kivel, N.; Kokkoris, M.; Krtička, M.; Leal-Cidoncha, E.; Lederer, C.; Leeb, H.; Lerendegui-Marco, J.; Meo, S. Lo; Lonsdale, S. J.; Losito, R.; Macina, D.; Marganiec, J.; Martínez, T.; Massimi, C.; Mastinu, P.; Mastromarco, M.; Matteucci, F.; Maugeri, E. A.; Mendoza, E.; Mengoni, A.; Milazzo, P. M.; Mingrone, F.; Mirea, M.; Montesano, S.; Musumarra, A.; Nolte, R.; Oprea, A.; Patronis, N.; Pavlik, A.; Perkowski, J.; Porras, J. I.; Praena, J.; Quesada, J. M.; Rajeev, K.; Rauscher, T.; Reifarth, R.; Riego-Perez, A.; Rout, P. C.; Rubbia, C.; Ryan, J. A.; Sabaté-Gilarte, M.; Saxena, A.; Schillebeeckx, P.; Schmidt, S.; Schumann, D.; Sedyshev, P.; Smith, A. G.; Stamatopoulos, A.; Tagliente, G.; Tain, J. L.; Tarifeño-Saldivia, A.; Tassan-Got, L.; Tsinganis, A.; Valenta, S.; Vannini, G.; Variale, V.; Vaz, P.; Ventura, A.; Vlachoudis, V.; Vlastou, R.; Wallner, A.; Warren, S.; Weigand, M.; Weiss, C.; Wolf, C.; Woods, P. J.; Wright, T.; Žugec, P.

2017-09-01

The Cosmological Lithium Problem refers to the large discrepancy between the abundance of primordial 7Li predicted by the standard theory of Big Bang Nucleosynthesis and the value inferred from the so-called "Spite plateau" in halo stars. A possible explanation for this longstanding puzzle in Nuclear Astrophysics is related to the incorrect estimation of the destruction rate of 7Be, which is responsible for the production of 95% of primordial Lithium. While charged-particle induced reactions have mostly been ruled out, data on the 7Be(n,α) and 7Be(n,p) reactions are scarce or completely missing, so that a large uncertainty still affects the abundance of 7Li predicted by the standard theory of Big Bang Nucleosynthesis. Both reactions have been measured at the n_TOF facility at CERN, providing for the first time data in a wide neutron energy range.

The Atacama Cosmology Telescope: cosmological parameters from three seasons of data

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sievers, Jonathan L.; Appel, John William; Hlozek, Renée A.

2013-10-01

We present constraints on cosmological and astrophysical parameters from high-resolution microwave background maps at 148 GHz and 218 GHz made by the Atacama Cosmology Telescope (ACT) in three seasons of observations from 2008 to 2010. A model of primary cosmological and secondary foreground parameters is fit to the map power spectra and lensing deflection power spectrum, including contributions from both the thermal Sunyaev-Zeldovich (tSZ) effect and the kinematic Sunyaev-Zeldovich (kSZ) effect, Poisson and correlated anisotropy from unresolved infrared sources, radio sources, and the correlation between the tSZ effect and infrared sources. The power ℓ{sup 2}C{sub ℓ}/2π of the thermal SZmore » power spectrum at 148 GHz is measured to be 3.4±1.4  μK{sup 2} at ℓ = 3000, while the corresponding amplitude of the kinematic SZ power spectrum has a 95% confidence level upper limit of 8.6  μK{sup 2}. Combining ACT power spectra with the WMAP 7-year temperature and polarization power spectra, we find excellent consistency with the LCDM model. We constrain the number of effective relativistic degrees of freedom in the early universe to be N{sub eff} = 2.79±0.56, in agreement with the canonical value of N{sub eff} = 3.046 for three massless neutrinos. We constrain the sum of the neutrino masses to be Σm{sub ν} < 0.39 eV at 95% confidence when combining ACT and WMAP 7-year data with BAO and Hubble constant measurements. We constrain the amount of primordial helium to be Y{sub p} = 0.225±0.034, and measure no variation in the fine structure constant α since recombination, with α/α{sub 0} = 1.004±0.005. We also find no evidence for any running of the scalar spectral index, dn{sub s}/dln k = −0.004±0.012.« less

PoMiN: A Post-Minkowskian N-body Solver

NASA Astrophysics Data System (ADS)

Feng, Justin; Baumann, Mark; Hall, Bryton; Doss, Joel; Spencer, Lucas; Matzner, Richard

2018-06-01

In this paper, we introduce PoMiN, a lightweight N-body code based on the post-Minkowskian N-body Hamiltonian of Ledvinka et al., which includes general relativistic effects up to first order in Newton’s constant G, and all orders in the speed of light c. PoMiN is written in C and uses a fourth-order Runge–Kutta integration scheme. PoMiN has also been written to handle an arbitrary number of particles (both massive and massless), with a computational complexity that scales as O(N 2). We describe the methods we used to simplify and organize the Hamiltonian, and the tests we performed (convergence, conservation, and analytical comparison tests) to validate the code.

Constraining Cosmological Models with Different Observations

NASA Astrophysics Data System (ADS)

Wei, J. J.

2016-07-01

cosmological probes, perhaps even out to redshifts much greater (z≫2) than those accessible using SNe Ia. However, the currently available sample of SNe Ia is still quite small. Our simulations have shown that if SLSNe Ic can be commonly detected in the future, they have the potential of greatly refining the measurement of cosmological parameters, particularly the parameter w_{de} of the dark energy equation of state. In Chapter 3, we focus on GRB cosmology. We firstly use GRBs as standard candles in constructing the Hubble diagram at redshifts beyond the current reach of SNe Ia observations. Then we measure high-z star formation rate (SFR) using GRBs. We confirm that the latest Swift sample of GRBs reveals an increasing evolution in the GRB rate relative to SFR at high redshifts. The observed discrepancy between the GRB rate and the SFR may be eliminated by assuming a cosmic evolution in metallicity. Assuming that the SFR and GRB rate are related via an evolving metallicity, we find that the GRB data constrain the slope of the high-z SFR to be -2.41_{-2.09}^{+1.87}. In addition, first stars can only form in structures that are suitably dense, which can be parameterized by the minimum dark matter halo mass M_{min}. M_{min} must play an important role in star formation. We can constrain M_{min}<10^{12.5} M_{⊙} at 68% confidence level from the GRB data. In Chapter 4, we assemble a catalog of 69 strong gravitational lensing systems, and carefully introduce how to constrain cosmological parameters using these important data. We find that both ΛCDM and the R_{h}=ct Universe account for the lens observations quite well, though the precision of these measurements does not appear to be good enough to favor one model over the other. In Chapters 5 and 6, we use measurements of the galaxy-cluster angular diameter distances and 32 age measurements of passively evolving galaxies to test and compare the standard model (ΛCDM) and the R_{h}=ct Universe, respectively. We show that both

Building cosmological frozen stars

NASA Astrophysics Data System (ADS)

Kastor, David; Traschen, Jennie

2017-02-01

Janis-Newman-Winicour (JNW) solutions generalize Schwarzschild to include a massless scalar field. While they share the familiar infinite redshift feature of Schwarzschild, they suffer from the presence of naked singularities. Cosmological versions of JNW spacetimes were discovered some years ago, in the most general case, by Fonarev. Fonarev solutions are also plagued by naked singularities, but have the virtue, unlike e.g. Schwarzschild-deSitter, of being dynamical. Given that exact dynamical cosmological black hole solutions are scarce, Fonarev solutions merit further study. We show how Fonarev solutions can be obtained via generalized dimensional reduction from simpler static vacuum solutions. These results may lead towards constructions of actual dynamical cosmological black holes. In particular, we note that cosmological versions of extremal charged dilaton black holes are known. JNW spacetimes represent a different limiting case of the family of charged dilaton black holes, which have been important in the context of string theory, and better understanding their cosmological versions of JNW spacetimes thus provides a second data point towards finding cosmological versions of the entire family.

Cosmological Models and Stability

NASA Astrophysics Data System (ADS)

Andersson, Lars

Principles in the form of heuristic guidelines or generally accepted dogma play an important role in the development of physical theories. In particular, philosophical considerations and principles figure prominently in the work of Albert Einstein. As mentioned in the talk by Jiří Bičák at this conference, Einstein formulated the equivalence principle, an essential step on the road to general relativity, during his time in Prague 1911-1912. In this talk, I would like to discuss some aspects of cosmological models. As cosmology is an area of physics where "principles" such as the "cosmological principle" or the "Copernican principle" play a prominent role in motivating the class of models which form part of the current standard model, I will start by comparing the role of the equivalence principle to that of the principles used in cosmology. I will then briefly describe the standard model of cosmology to give a perspective on some mathematical problems and conjectures on cosmological models, which are discussed in the later part of this paper.

The Cosmological Lithium Problem and the Measurement of the 7Be(n, α) Reaction at n_TOF-CERN

NASA Astrophysics Data System (ADS)

Musumarra, Agatino; Barbagallo, Massimo

A possible explanation of the so-called "Cosmological Lithium Problem", an important unsolved problem in Nuclear Astrophysics, involves large systematic uncertainties in the cross-sections of reactions leading to the destruction of 7Be during the Big-Bang Nucleosynthesis (BBN). Among these reactions, the 7Be(n, α) is the most uncertain. So far, only a single measurement with thermal neutrons has been performed. Therefore, BBN calculations had to rely on rather uncertain theoretical extrapolations. The short half-life of 7Be (53.29 d) and the low cross section have prevented, up to now, to obtain experimental data at keV neutron energies typical for BBN studies. We have measured for the first time at n_TOF the 7Be(n, α) reaction in a wide neutron energy range, from thermal up to 10 keV. This measurement has been performed, at the new beam line (EAR2) of the Neutron-Time-Of-Flight facility n_TOF at CERN. The two α-particles, emitted back-to-back in the reaction, have been detected by mean of sandwiches of silicon detectors and, by exploiting the coincidence technique, we were able to suppress the large γ and n-induced background. The 7Be isotope production and purification has been performed by PSI-Zurich Switzerland.

Cosmological Particle Data Compression in Practice

NASA Astrophysics Data System (ADS)

Zeyen, M.; Ahrens, J.; Hagen, H.; Heitmann, K.; Habib, S.

2017-12-01

In cosmological simulations trillions of particles are handled and several terabytes of unstructured particle data are generated in each time step. Transferring this data directly from memory to disk in an uncompressed way results in a massive load on I/O and storage systems. Hence, one goal of domain scientists is to compress the data before storing it to disk while minimizing the loss of information. To prevent reading back uncompressed data from disk, this can be done in an in-situ process. Since the simulation continuously generates data, the available time for the compression of one time step is limited. Therefore, the evaluation of compression techniques has shifted from only focusing on compression rates to include run-times and scalability.In recent years several compression techniques for cosmological data have become available. These techniques can be either lossy or lossless, depending on the technique. For both cases, this study aims to evaluate and compare the state of the art compression techniques for unstructured particle data. This study focuses on the techniques available in the Blosc framework with its multi-threading support, the XZ Utils toolkit with the LZMA algorithm that achieves high compression rates, and the widespread FPZIP and ZFP methods for lossy compressions.For the investigated compression techniques, quantitative performance indicators such as compression rates, run-time/throughput, and reconstruction errors are measured. Based on these factors, this study offers a comprehensive analysis of the individual techniques and discusses their applicability for in-situ compression. In addition, domain specific measures are evaluated on the reconstructed data sets, and the relative error rates and statistical properties are analyzed and compared. Based on this study future challenges and directions in the compression of unstructured cosmological particle data were identified.

Investigating cluster astrophysics and cosmology with cross-correlation of the thermal Sunyaev-Zel'dovich effect and weak lensing

NASA Astrophysics Data System (ADS)

Osato, Ken; Flender, Samuel; Nagai, Daisuke; Shirasaki, Masato; Yoshida, Naoki

2018-03-01

Recent detections of the cross-correlation of the thermal Sunyaev-Zel'dovich (tSZ) effect and weak gravitational lensing (WL) enable unique studies of cluster astrophysics and cosmology. In this work, we present constraints on the amplitude of the non-thermal pressure fraction in galaxy clusters, α0, and the amplitude of the matter power spectrum, σ8, using measurements of the tSZ power spectrum from Planck, and the tSZ-WL cross-correlation from Planck and the Red Cluster Sequence Lensing Survey. We fit the data to a semi-analytic model with the covariance matrix using N-body simulations. We find that the tSZ power spectrum alone prefers σ8 ˜ 0.85 and a large fraction of non-thermal pressure (α0 ˜ 0.2-0.3). The tSZ-WL cross-correlation on the other hand prefers a significantly lower σ8 ˜ 0.6 and low α0 ˜ 0.05. We show that this tension can be mitigated by allowing for a steep slope in the stellar mass-halo mass relation, which would cause a reduction of the gas in low-mass haloes. In such a model, the combined data prefer σ8 ˜ 0.7 and α0 ˜ 0.2, consistent with predictions from hydrodynamical simulations.

Cosmological constraints from the convergence 1-point probability distribution

DOE Office of Scientific and Technical Information (OSTI.GOV)

Patton, Kenneth; Blazek, Jonathan; Honscheid, Klaus

2017-06-29

Here, we examine the cosmological information available from the 1-point probability density function (PDF) of the weak-lensing convergence field, utilizing fast l-picola simulations and a Fisher analysis. We find competitive constraints in the Ωm–σ8 plane from the convergence PDF with 188 arcmin 2 pixels compared to the cosmic shear power spectrum with an equivalent number of modes (ℓ < 886). The convergence PDF also partially breaks the degeneracy cosmic shear exhibits in that parameter space. A joint analysis of the convergence PDF and shear 2-point function also reduces the impact of shape measurement systematics, to which the PDF is lessmore » susceptible, and improves the total figure of merit by a factor of 2–3, depending on the level of systematics. Finally, we present a correction factor necessary for calculating the unbiased Fisher information from finite differences using a limited number of cosmological simulations.« less

Evaluation of an Interactive Undergraduate Cosmology Curriculum

NASA Astrophysics Data System (ADS)

White, Aaron; Coble, Kimberly A.; Martin, Dominique; Hayes, Patrycia; Targett, Tom; Cominsky, Lynn R.

2018-06-01

The Big Ideas in Cosmology is an immersive set of web-based learning modules that integrates text, figures, and visualizations with short and long interactive tasks as well as labs that allow students to manipulate and analyze real cosmological data. This enables the transformation of general education astronomy and cosmology classes from primarily lecture and book-based courses to a format that builds important STEM skills, while engaging those outside the field with modern discoveries and a more realistic sense of practices and tools used by professional astronomers. Over two semesters, we field-tested the curriculum in general education cosmology classes at a state university in California [N ~ 80]. We administered pre- and post-instruction multiple-choice and open-ended content surveys as well as the CLASS, to gauge the effectiveness of the course and modules. Questions addressed included the structure, composition, and evolution of the universe, including students’ reasoning and “how we know.”Module development and evaluation was supported by NASA ROSES E/PO Grant #NNXl0AC89G, the Illinois Space Grant Consortium, the Fermi E/PO program, Sonoma State University’s Space Science Education and Public Outreach Group, and San Francisco State University. The modules are published by Great River Learning/Kendall-Hunt.

Problem of time in slightly inhomogeneous cosmology

NASA Astrophysics Data System (ADS)

Anderson, Edward

2016-07-01

The problem of time (PoT) is a multi-faceted conceptual incompatibility between various areas of Theoretical Physics. While usually stated as between GR and QM, in fact 8/9ths of it is already present at the classical level. Thus we adopt a ‘top-down’ classical and then quantum approach. I consider a local resolution to the PoT that is Machian, which was previously realized for relational triangle and minisuperspace models. This resolution has three levels: classical, semiclassical and combined semiclassical-histories-records. This article’s specific model is a slightly inhomogeneous cosmology considered for now at the classical level. This is motivated by how the inhomogeneous fluctuations that underlie structure formation—galaxies and CMB hotspots—might have been seeded by quantum cosmological fluctuations, as magnified by some inflationary mechanism. In particular, I consider the perturbations about {{{S}}}3 case of this involving up to second order, which has a number of parallels with the Halliwell-Hawking model but has a number of conceptual differences and useful upgrades. The article’s main features are that the elimination part of the model’s thin sandwich is straightforward, but the modewise split of the constraints fail to be first-class constraints. Thus the elimination part only arises as an intermediate geometry between superspace and Riem. The reduced geometries have surprising singularities influenced by the matter content of the Universe, though the N-body problem anticipates these with its collinear singularities. I also give a ‘basis set’ of Kuchař beables for this model arena.

Running of the spectrum of cosmological perturbations in string gas cosmology

NASA Astrophysics Data System (ADS)

Brandenberger, Robert; Franzmann, Guilherme; Liang, Qiuyue

2017-12-01

We compute the running of the spectrum of cosmological perturbations in string gas cosmology, making use of a smooth parametrization of the transition between the early Hagedorn phase and the later radiation phase. We find that the running has the same sign as in simple models of single scalar field inflation. Its magnitude is proportional to (1 -ns) (ns being the slope index of the spectrum), and it is thus parametrically larger than for inflationary cosmology, where it is proportional to (1 -ns)2 .

Moving-mesh cosmology: characteristics of galaxies and haloes

NASA Astrophysics Data System (ADS)

Kereš, Dušan; Vogelsberger, Mark; Sijacki, Debora; Springel, Volker; Hernquist, Lars

2012-09-01

We discuss cosmological hydrodynamic simulations of galaxy formation performed with the new moving-mesh code AREPO, which promises higher accuracy compared with the traditional smoothed particle hydrodynamics (SPH) technique that has been widely employed for this problem. In this exploratory study, we deliberately limit the complexity of the physical processes followed by the code for ease of comparison with previous calculations, and include only cooling of gas with a primordial composition, heating by a spatially uniform ultraviolet background, and a simple subresolution model for regulating star formation in the dense interstellar medium. We use an identical set of physics in corresponding simulations carried out with the well-tested SPH code GADGET, adopting also the same high-resolution gravity solver. We are thus able to compare both simulation sets on an object-by-object basis, allowing us to cleanly isolate the impact of different hydrodynamical methods on galaxy and halo properties. In accompanying papers, Vogelsberger et al. and Sijacki et al., we focus on an analysis of the global baryonic statistics predicted by the simulation codes, and complementary idealized simulations that highlight the differences between the hydrodynamical schemes. Here we investigate their influence on the baryonic properties of simulated galaxies and their surrounding haloes. We find that AREPO leads to significantly higher star formation rates for galaxies in massive haloes and to more extended gaseous discs in galaxies, which also feature a thinner and smoother morphology than their GADGET counterparts. Consequently, galaxies formed in AREPO have larger sizes and higher specific angular momentum than their SPH correspondents. Interestingly, the more efficient cooling flows in AREPO yield higher densities and lower entropies in halo centres compared to GADGET, whereas the opposite trend is found in halo outskirts. The cooling differences leading to higher star formation rates

Higgs cosmology

PubMed Central

2018-01-01

The discovery of the Higgs boson in 2012 and other results from the Large Hadron Collider have confirmed the standard model of particle physics as the correct theory of elementary particles and their interactions up to energies of several TeV. Remarkably, the theory may even remain valid all the way to the Planck scale of quantum gravity, and therefore it provides a solid theoretical basis for describing the early Universe. Furthermore, the Higgs field itself has unique properties that may have allowed it to play a central role in the evolution of the Universe, from inflation to cosmological phase transitions and the origin of both baryonic and dark matter, and possibly to determine its ultimate fate through the electroweak vacuum instability. These connections between particle physics and cosmology have given rise to a new and growing field of Higgs cosmology, which promises to shed new light on some of the most puzzling questions about the Universe as new data from particle physics experiments and cosmological observations become available. This article is part of the Theo Murphy meeting issue ‘Higgs cosmology’. PMID:29358352

Some new concepts in the n-body and 3-body problems

NASA Astrophysics Data System (ADS)

Kyrala, A.

1982-06-01

A new approach to the n-body problem in terms of an rms particle velocity and a harmonic mean particle separation has been constructed by using averaging procedures formulated in terms of a single parameter. A systematic classification of escape and collision processes by means of specific polynomials, which can be used somewhat like generating functions, is introduced. For n-body problems with non-null total angular momentum, an rms angular momentum is defined which together with a harmonic mean centroidal moment of inertia characterizes the rotational kinetic energy. Finally, a graphical construction for the equipotentials of the three-body problem is given and attention is drawn to the use of the apex, defined as the point of least average separation, in this problem. It is supposed that the n-bodies interact with one another via the Newtonian potential in an inertial system.

Power law cosmology model comparison with CMB scale information

NASA Astrophysics Data System (ADS)

Tutusaus, Isaac; Lamine, Brahim; Blanchard, Alain; Dupays, Arnaud; Zolnierowski, Yves; Cohen-Tanugi, Johann; Ealet, Anne; Escoffier, Stéphanie; Le Fèvre, Olivier; Ilić, Stéphane; Pisani, Alice; Plaszczynski, Stéphane; Sakr, Ziad; Salvatelli, Valentina; Schücker, Thomas; Tilquin, André; Virey, Jean-Marc

2016-11-01

Despite the ability of the cosmological concordance model (Λ CDM ) to describe the cosmological observations exceedingly well, power law expansion of the Universe scale radius, R (t )∝tn, has been proposed as an alternative framework. We examine here these models, analyzing their ability to fit cosmological data using robust model comparison criteria. Type Ia supernovae (SNIa), baryonic acoustic oscillations (BAO) and acoustic scale information from the cosmic microwave background (CMB) have been used. We find that SNIa data either alone or combined with BAO can be well reproduced by both Λ CDM and power law expansion models with n ˜1.5 , while the constant expansion rate model (n =1 ) is clearly disfavored. Allowing for some redshift evolution in the SNIa luminosity essentially removes any clear preference for a specific model. The CMB data are well known to provide the most stringent constraints on standard cosmological models, in particular, through the position of the first peak of the temperature angular power spectrum, corresponding to the sound horizon at recombination, a scale physically related to the BAO scale. Models with n ≥1 lead to a divergence of the sound horizon and do not naturally provide the relevant scales for the BAO and the CMB. We retain an empirical footing to overcome this issue: we let the data choose the preferred values for these scales, while we recompute the ionization history in power law models, to obtain the distance to the CMB. In doing so, we find that the scale coming from the BAO data is not consistent with the observed position of the first peak of the CMB temperature angular power spectrum for any power law cosmology. Therefore, we conclude that when the three standard probes (SNIa, BAO, and CMB) are combined, the Λ CDM model is very strongly favored over any of these alternative models, which are then essentially ruled out.

Cosmology with decaying cosmological constant—exact solutions and model testing

DOE Office of Scientific and Technical Information (OSTI.GOV)

Szydłowski, Marek; Stachowski, Aleksander, E-mail: marek.szydlowski@uj.edu.pl, E-mail: aleksander.stachowski@uj.edu.pl

We study dynamics of Λ(t) cosmological models which are a natural generalization of the standard cosmological model (the ΛCDM model). We consider a class of models: the ones with a prescribed form of Λ(t)=Λ{sub bare}+α{sup 2}/t{sup 2}. This type of a Λ(t) parametrization is motivated by different cosmological approaches. We interpret the model with running Lambda (Λ(t)) as a special model of an interacting cosmology with the interaction term −dΛ(t)/dt in which energy transfer is between dark matter and dark energy sectors. For the Λ(t) cosmology with a prescribed form of Λ(t) we have found the exact solution in themore » form of Bessel functions. Our model shows that fractional density of dark energy Ω{sub e} is constant and close to zero during the early evolution of the universe. We have also constrained the model parameters for this class of models using the astronomical data such as SNIa data, BAO, CMB, measurements of H(z) and the Alcock-Paczyński test. In this context we formulate a simple criterion of variability of Λ with respect to t in terms of variability of the jerk or sign of estimator (1−Ω{sub m},0−Ω{sub Λ,0}). The case study of our model enable us to find an upper limit α{sup 2} < 0.012 (2σ C.L.) describing the variation from the cosmological constant while the LCDM model seems to be consistent with various data.« less

The growth of central and satellite galaxies in cosmological smoothed particle hydrodynamics simulations

NASA Astrophysics Data System (ADS)

Simha, Vimal; Weinberg, David H.; Davé, Romeel; Gnedin, Oleg Y.; Katz, Neal; Kereš, Dušan

2009-10-01

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to 7 × 109Msolar. Most friends-of-friends haloes in the simulation have a distinct central galaxy, typically 2-5 times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive haloes form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5-1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Typically, satellites in our simulation are 0.1-0.2 mag bluer than in models that assume no gas accretion on to satellites after a halo merger. Since z = 1, 27 per cent of central galaxies (above 3 × 1010Msolar) and 22 per cent of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain `central' objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.

The nature of massive transition galaxies in CANDELS, GAMA and cosmological simulations

NASA Astrophysics Data System (ADS)

Pandya, Viraj; Brennan, Ryan; Somerville, Rachel S.; Choi, Ena; Barro, Guillermo; Wuyts, Stijn; Taylor, Edward N.; Behroozi, Peter; Kirkpatrick, Allison; Faber, Sandra M.; Primack, Joel; Koo, David C.; McIntosh, Daniel H.; Kocevski, Dale; Bell, Eric F.; Dekel, Avishai; Fang, Jerome J.; Ferguson, Henry C.; Grogin, Norman; Koekemoer, Anton M.; Lu, Yu; Mantha, Kameswara; Mobasher, Bahram; Newman, Jeffrey; Pacifici, Camilla; Papovich, Casey; van der Wel, Arjen; Yesuf, Hassen M.

2017-12-01

We explore observational and theoretical constraints on how galaxies might transition between the 'star-forming main sequence' (SFMS) and varying 'degrees of quiescence' out to z = 3. Our analysis is focused on galaxies with stellar mass M* > 1010 M⊙, and is enabled by GAMA and CANDELS observations, a semi-analytic model (SAM) of galaxy formation, and a cosmological hydrodynamical 'zoom in' simulation with momentum-driven AGN feedback. In both the observations and the SAM, transition galaxies tend to have intermediate Sérsic indices, half-light radii, and surface stellar mass densities compared to star-forming and quiescent galaxies out to z = 3. We place an observational upper limit on the average population transition time-scale as a function of redshift, finding that the average high-redshift galaxy is on a 'fast track' for quenching whereas the average low-redshift galaxy is on a 'slow track' for quenching. We qualitatively identify four physical origin scenarios for transition galaxies in the SAM: oscillations on the SFMS, slow quenching, fast quenching, and rejuvenation. Quenching time-scales in both the SAM and the hydrodynamical simulation are not fast enough to reproduce the quiescent population that we observe at z ∼ 3. In the SAM, we do not find a clear-cut morphological dependence of quenching time-scales, but we do predict that the mean stellar ages, cold gas fractions, SMBH (supermassive black hole) masses and halo masses of transition galaxies tend to be intermediate relative to those of star-forming and quiescent galaxies at z < 3.

Axion cosmology

NASA Astrophysics Data System (ADS)

Marsh, David J. E.

2016-07-01

Axions comprise a broad class of particles that can play a major role in explaining the unknown aspects of cosmology. They are also well-motivated within high energy physics, appearing in theories related to CP-violation in the standard model, supersymmetric theories, and theories with extra-dimensions, including string theory, and so axion cosmology offers us a unique view onto these theories. I review the motivation and models for axions in particle physics and string theory. I then present a comprehensive and pedagogical view on the cosmology and astrophysics of axion-like particles, starting from inflation and progressing via BBN, the CMB, reionization and structure formation, up to the present-day Universe. Topics covered include: axion dark matter (DM); direct and indirect detection of axions, reviewing existing and future experiments; axions as dark radiation; axions and the cosmological constant problem; decays of heavy axions; axions and stellar astrophysics; black hole superradiance; axions and astrophysical magnetic fields; axion inflation, and axion DM as an indirect probe of inflation. A major focus is on the population of ultralight axions created via vacuum realignment, and its role as a DM candidate with distinctive phenomenology. Cosmological observations place robust constraints on the axion mass and relic density in this scenario, and I review where such constraints come from. I next cover aspects of galaxy formation with axion DM, and ways this can be used to further search for evidence of axions. An absolute lower bound on DM particle mass is established. It is ma > 10-24eV from linear observables, extending to ma ≳ 10-22eV from non-linear observables, and has the potential to reach ma ≳ 10-18eV in the future. These bounds are weaker if the axion is not all of the DM, giving rise to limits on the relic density at low mass. This leads to the exciting possibility that the effects of axion DM on structure formation could one day be detected

Tachyon with an inverse power-law potential in a braneworld cosmology

NASA Astrophysics Data System (ADS)

Bilić, Neven; Domazet, Silvije; Djordjevic, Goran S.

2017-08-01

We study a tachyon cosmological model based on the dynamics of a 3-brane in the bulk of the second Randall-Sundrum model extended to more general warp functions. A well known prototype of such a generalization is the bulk with a selfinteracting scalar field. As a consequence of a generalized bulk geometry the cosmology on the observer brane is modified by the scale dependent four-dimensional gravitational constant. In particular, we study a power law warp factor which generates an inverse power-law potential V\\propto \\varphi-n of the tachyon field φ. We find a critical power n cr that divides two subclasses with distinct asymptotic behaviors: a dust universe for n>n_cr and a quasi de Sitter universe for 0.

Sociology of Modern Cosmology

NASA Astrophysics Data System (ADS)

López-Corredoira, M.

2009-08-01

Certain results of observational cosmology cast critical doubt on the foundations of standard cosmology but leave most cosmologists untroubled. Alternative cosmological models that differ from the Big Bang have been published and defended by heterodox scientists; however, most cosmologists do not heed these. This may be because standard theory is correct and all other ideas and criticisms are incorrect, but it is also to a great extent due to sociological phenomena such as the ``snowball effect'' or ``groupthink''. We might wonder whether cosmology, the study of the Universe as a whole, is a science like other branches of physics or just a dominant ideology.

A whirling plane of satellite galaxies around Centaurus A challenges cold dark matter cosmology.

PubMed

Müller, Oliver; Pawlowski, Marcel S; Jerjen, Helmut; Lelli, Federico

2018-02-02

The Milky Way and Andromeda galaxies are each surrounded by a thin plane of satellite dwarf galaxies that may be corotating. Cosmological simulations predict that most satellite galaxy systems are close to isotropic with random motions, so those two well-studied systems are often interpreted as rare statistical outliers. We test this assumption using the kinematics of satellite galaxies around the Centaurus A galaxy. Our statistical analysis reveals evidence for corotation in a narrow plane: Of the 16 Centaurus A satellites with kinematic data, 14 follow a coherent velocity pattern aligned with the long axis of their spatial distribution. In standard cosmological simulations, <0.5% of Centaurus A-like systems show such behavior. Corotating satellite systems may be common in the universe, challenging small-scale structure formation in the prevailing cosmological paradigm. Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.

Linear response to long wavelength fluctuations using curvature simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Baldauf, Tobias; Zaldarriaga, Matias; Seljak, Uroš

2016-09-01

We study the local response to long wavelength fluctuations in cosmological N -body simulations, focusing on the matter and halo power spectra, halo abundance and non-linear transformations of the density field. The long wavelength mode is implemented using an effective curved cosmology and a mapping of time and distances. The method provides an alternative, more direct, way to measure the isotropic halo biases. Limiting ourselves to the linear case, we find generally good agreement between the biases obtained from the curvature method and the traditional power spectrum method at the level of a few percent. We also study the responsemore » of halo counts to changes in the variance of the field and find that the slope of the relation between the responses to density and variance differs from the naïve derivation assuming a universal mass function by approximately 8–20%. This has implications for measurements of the amplitude of local non-Gaussianity using scale dependent bias. We also analyze the halo power spectrum and halo-dark matter cross-spectrum response to long wavelength fluctuations and derive second order halo bias from it, as well as the super-sample variance contribution to the galaxy power spectrum covariance matrix.« less

UNITY: Confronting Supernova Cosmology's Statistical and Systematic Uncertainties in a Unified Bayesian Framework

NASA Astrophysics Data System (ADS)

Rubin, D.; Aldering, G.; Barbary, K.; Boone, K.; Chappell, G.; Currie, M.; Deustua, S.; Fagrelius, P.; Fruchter, A.; Hayden, B.; Lidman, C.; Nordin, J.; Perlmutter, S.; Saunders, C.; Sofiatti, C.; Supernova Cosmology Project, The

2015-11-01

While recent supernova (SN) cosmology research has benefited from improved measurements, current analysis approaches are not statistically optimal and will prove insufficient for future surveys. This paper discusses the limitations of current SN cosmological analyses in treating outliers, selection effects, shape- and color-standardization relations, unexplained dispersion, and heterogeneous observations. We present a new Bayesian framework, called UNITY (Unified Nonlinear Inference for Type-Ia cosmologY), that incorporates significant improvements in our ability to confront these effects. We apply the framework to real SN observations and demonstrate smaller statistical and systematic uncertainties. We verify earlier results that SNe Ia require nonlinear shape and color standardizations, but we now include these nonlinear relations in a statistically well-justified way. This analysis was primarily performed blinded, in that the basic framework was first validated on simulated data before transitioning to real data. We also discuss possible extensions of the method.

Cosmological imprints of frozen-in light sterile neutrinos

DOE Office of Scientific and Technical Information (OSTI.GOV)

Roland, Samuel B.; Shakya, Bibhushan, E-mail: rolandsa@umich.edu, E-mail: bshakya@umich.edu

We investigate observable cosmological aspects of sterile neutrino dark matter produced via the freeze-in mechanism. The study is performed in a framework that admits many cosmologically interesting variations: high temperature production via annihilation processes from higher dimensional operators or low temperature production from decays of a scalar, with the decaying scalar in or out of equilibrium with the thermal bath, in supersymmetric or non-supersymmetric setups, thus allowing us to both extract generic properties and highlight features unique to particular variations. We find that while such sterile neutrinos are generally compatible with all cosmological constraints, interesting scenarios can arise where darkmore » matter is cold, warm, or hot, has nontrivial momentum distributions, or provides contributions to the effective number of relativistic degrees of freedom N {sub eff} during Big Bang nucleosynthesis large enough to be probed by future measurements.« less

Simulations of the WFIRST Supernova Survey and Forecasts of Cosmological Constraints

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hounsell, R.; Scolnic, D.; Foley, R. J.

The Wide Field InfraRed Survey Telescope (WFIRST) was the highest rankedlarge space-based mission of the 2010 New Worlds, New Horizons decadal survey.It is now a NASA mission in formulation with a planned launch in themid-2020's. A primary mission objective is to precisely constrain the nature ofdark energy through multiple probes, including Type Ia supernovae (SNe Ia).Here, we present the first realistic simulations of the WFIRST SN survey basedon current hardware specifications and using open-source tools. We simulate SNlight curves and spectra as viewed by the WFIRST wide-field channel (WFC)imager and integral field channel (IFC) spectrometer, respectively. We examine11 survey strategiesmore » with different time allocations between the WFC and IFC,two of which are based upon the strategy described by the WFIRST ScienceDefinition Team, which measures SN distances exclusively from IFC data. Wepropagate statistical and, crucially, systematic uncertainties to predict thedark energy task force figure of merit (DETF FoM) for each strategy. Theincrease in FoM values with SN search area is limited by the overhead times foreach exposure. For IFC-focused strategies the largest individual systematicuncertainty is the wavelength-dependent calibration uncertainty, whereas forWFC-focused strategies, it is the intrinsic scatter uncertainty. We find thatthe best IFC-focused and WFC-exclusive strategies have comparable FoM values.Even without improvements to other cosmological probes, the WFIRST SN surveyhas the potential to increase the FoM by more than an order of magnitude fromthe current values. Although the survey strategies presented here have not beenfully optimized, these initial investigations are an important step in thedevelopment of the final hardware design and implementation of the WFIRSTmission.« less

The SLUGGS survey: a comparison of total-mass profiles of early-type galaxies from observations and cosmological simulations, to ˜4 effective radii

NASA Astrophysics Data System (ADS)

Bellstedt, Sabine; Forbes, Duncan A.; Romanowsky, Aaron J.; Remus, Rhea-Silvia; Stevens, Adam R. H.; Brodie, Jean P.; Poci, Adriano; McDermid, Richard; Alabi, Adebusola; Chevalier, Leonie; Adams, Caitlin; Ferré-Mateu, Anna; Wasserman, Asher; Pandya, Viraj

2018-06-01

We apply the Jeans Anisotropic Multi-Gaussian Expansion dynamical modelling method to SAGES Legacy Unifying Globulars and GalaxieS (SLUGGS) survey data of early-type galaxies in the stellar mass range 1010 < M*/M⊙ < 1011.6 that cover a large radial range of 0.1-4.0 effective radii. We combine SLUGGS and ATLAS3D data sets to model the total-mass profiles of a sample of 21 fast-rotator galaxies, utilizing a hyperparameter method to combine the two independent data sets. The total-mass density profile slope values derived for these galaxies are consistent with those measured in the inner regions of galaxies by other studies. Furthermore, the total-mass density slopes (γtot) appear to be universal over this broad stellar mass range, with an average value of γtot = -2.24 ± 0.05 , i.e. slightly steeper than isothermal. We compare our results to model galaxies from the Magneticum and EAGLE cosmological hydrodynamic simulations, in order to probe the mechanisms that are responsible for varying total-mass density profile slopes. The simulated-galaxy slopes are shallower than the observed values by ˜0.3-0.5, indicating that the physical processes shaping the mass distributions of galaxies in cosmological simulations are still incomplete. For galaxies with M* > 1010.7 M⊙ in the Magneticum simulations, we identify a significant anticorrelation between total-mass density profile slopes and the fraction of stellar mass formed ex situ (i.e. accreted), whereas this anticorrelation is weaker for lower stellar masses, implying that the measured total-mass density slopes for low-mass galaxies are less likely to be determined by merger activity.

Integrable scalar cosmologies. II. Can they fit into Gauged Extended Supergravity or be encoded in N=1 superpotentials?

NASA Astrophysics Data System (ADS)

Fré, P.; Sorin, A. S.; Trigiante, M.

2014-04-01

The question whether the integrable one-field cosmologies classified in a previous paper by Fré, Sagnotti and Sorin can be embedded as consistent one-field truncations into Extended Gauged Supergravity or in N=1 supergravity gauged by a superpotential without the use of D-terms is addressed in this paper. The answer is that such an embedding is very difficult and rare but not impossible. Indeed, we were able to find two examples of integrable models embedded in supergravity in this way. Both examples are fitted into N=1 supergravity by means of a very specific and interesting choice of the superpotential W(z). The question whether there are examples of such an embedding in Extended Gauged Supergravity remains open. In the present paper, relying on the embedding tensor formalism we classified all gaugings of the N=2 STU model, confirming, in the absence on hypermultiplets, the uniqueness of the stable de Sitter vacuum found several years ago by Fré, Trigiante and Van Proeyen and excluding the embedding of any integrable cosmological model. A detailed analysis of the space of exact solutions of the first supergravity-embedded integrable cosmological model revealed several new features worth an in-depth consideration. When the scalar potential has an extremum at a negative value, the Universe necessarily collapses into a Big Crunch notwithstanding its spatial flatness. The causal structure of these Universes is quite different from that of the closed, positive curved, Universe: indeed, in this case the particle and event horizons do not coincide and develop complicated patterns. The cosmological consequences of this unexpected mechanism deserve careful consideration. The Cartan fieldshi associated with the Cartan generators of the Lie algebra G, whose number equals the rank r of G/H. For instance, in models associated with toroidal or orbifold compactifications, fields of this type are generically interpreted as radii of the underlying multi-tori. The axion fieldsb

Direct measurement of the 7Be(n, α)4 He reaction cross sections for the cosmological Li problem

NASA Astrophysics Data System (ADS)

Kawabata, Takahiro; Fujikawa, Yuki; Furuno, Tatsuya; Goto, Tatsuya; Hashimoto, Toshikazu; Ichikawa, Masaya; Itoh, Makoto; Iwasa, Naohito; Kanada-En'yo, Yoshiko; Koshikawa, Ami; Kubono, Shigeru; Miyawaki, Eisuke; Mizuno, Masatoshi; Mizutani, Keigo; Morimoto, Takahiro; Murata, Motoki; Nanamura, Takuya; Nishimura, Shunji; Nanamura, Takuya; Okamoto, Shintaro; Sakaguchi, Yuichi; Sakata, Itsushi; Sakaue, Akane; Sawada, Ryo; Shikata, Yuki; Takahashi, Yu; Takechi, Daiki; Takeda, Tomoya; Takimoto, Chisato; Tsumura, Miho; Watanabe, Ken; Yoshida, Sota

2017-11-01

The cross sections of the 7Be(n, α)4He reaction for p-wave neutrons were experimentally determined at Ec.m. = 0.20-0.81 MeV close to the Big Bang nucleosynthesis (BBN) energy window for the first time on the basis of the detailed balance principle by measuring the time-reverse reaction. The obtained cross sections are much larger than the cross sections for s-wave neutrons inferred from the recent measurement at the n_TOF facility in CERN, but significantly smaller than the theoretical estimation widely used in the BBN calculations. The present results suggest the 7Be(n, α)4 He reaction rate is not large enough to solve the cosmological lithium problem

On the Minimal Accuracy Required for Simulating Self-gravitating Systems by Means of Direct N-body Methods

NASA Astrophysics Data System (ADS)

Portegies Zwart, Simon; Boekholt, Tjarda

2014-04-01

The conservation of energy, linear momentum, and angular momentum are important drivers of our physical understanding of the evolution of the universe. These quantities are also conserved in Newton's laws of motion under gravity. Numerical integration of the associated equations of motion is extremely challenging, in particular due to the steady growth of numerical errors (by round-off and discrete time-stepping and the exponential divergence between two nearby solutions. As a result, numerical solutions to the general N-body problem are intrinsically questionable. Using brute force integrations to arbitrary numerical precision we demonstrate empirically that ensembles of different realizations of resonant three-body interactions produce statistically indistinguishable results. Although individual solutions using common integration methods are notoriously unreliable, we conjecture that an ensemble of approximate three-body solutions accurately represents an ensemble of true solutions, so long as the energy during integration is conserved to better than 1/10. We therefore provide an independent confirmation that previous work on self-gravitating systems can actually be trusted, irrespective of the intrinsically chaotic nature of the N-body problem.

Cosmology solved? Maybe

NASA Astrophysics Data System (ADS)

Turner, Michael S.

1999-03-01

For two decades the hot big-bang model as been referred to as the standard cosmology - and for good reason. For just as long cosmologists have known that there are fundamental questions that are not answered by the standard cosmology and point to a grander theory. The best candidate for that grander theory is inflation + cold dark matter. It holds that the Universe is flat, that slowly moving elementary particles left over from the earliest moments provide the cosmic infrastructure, and that the primeval density inhomogeneities that seed all the structure arose from quantum fluctuations. There is now prima facie evidence that supports two basic tenets of this paradigm. An avalanche of high-quality cosmological observations will soon make this case stronger or will break it. Key questions remain to be answered; foremost among them are: identification and detection of the cold dark matter particles and elucidation of the dark-energy component. These are exciting times in cosmology!

Linear perturbations in spherically symmetric dust cosmologies including a cosmological constant

NASA Astrophysics Data System (ADS)

Meyer, Sven; Bartelmann, Matthias

2017-12-01

We study the dynamical behaviour of gauge-invariant linear perturbations in spherically symmetric dust cosmologies including a cosmological constant. In contrast to spatially homogeneous FLRW models, the reduced degree of spatial symmetry causes a non-trivial dynamical coupling of gauge-invariant quantities already at first order perturbation theory and the strength and influence of this coupling on the spacetime evolution is investigated here. We present results on the underlying dynamical equations augmented by a cosmological constant and integrate them numerically. We also present a method to derive cosmologically relevant initial variables for this setup. Estimates of angular power spectra for each metric variable are computed and evaluated on the central observer's past null cone. By comparing the full evolution to the freely evolved initial profiles, the coupling strength will be determined for a best fit radially inhomogeneous patch obtained in previous works (see [1]). We find that coupling effects are not noticeable within the cosmic variance limit and can therefore safely be neglected for a relevant cosmological scenario. On the contrary, we find very strong coupling effects in a best fit spherical void model matching the distance redshift relation of SNe which is in accordance with previous findings using parametric void models.

N-Body Simulations of Planetary Accretion Around M Dwarf Stars

NASA Astrophysics Data System (ADS)

Ogihara, Masahiro; Ida, Shigeru

2009-07-01

We have investigated planetary accretion from planetesimals in terrestrial planet regions inside the ice line around M dwarf stars through N-body simulations including tidal interactions with disk gas. Because of low luminosity of M dwarfs, habitable zones (HZs) are located in inner regions (~0.1 AU). In the close-in HZ, type-I migration and the orbital decay induced by eccentricity damping are efficient according to the high disk gas density in the small orbital radii. Since the orbital decay is terminated around the disk inner edge and the disk edge is close to the HZ, the protoplanets accumulated near the disk edge affect formation of planets in the HZ. Ice lines are also in relatively inner regions at ~0.3 AU. Due to the small orbital radii, icy protoplanets accrete rapidly and undergo type-I migration before disk depletion. The rapid orbital decay, the proximity of the disk inner edge, and large amount of inflow of icy protoplanets are characteristic in planetary accretion in terrestrial planet regions around M dwarfs. In the case of full efficiency of type-I migration predicted by the linear theory, we found that protoplanets that migrate to the vicinity of the host star undergo close scatterings and collisions, and four to six planets eventually remain in mutual mean-motion resonances and their orbits have small eccentricities (lsim0.01) and they are stable both before and after disk gas decays. In the case of slow migration, the resonant capture is so efficient that densely packed ~40 small protoplanets remain in mutual mean-motion resonances. In this case, they start orbit crossing, after the disk gas decays and eccentricity damping due to tidal interaction with gas is no more effective. Through merging of the protoplanets, several planets in widely separated non-resonant orbits with relatively large eccentricities (~0.05) are formed. Thus, the final orbital configurations (separations, resonant or non-resonant, eccentricity, and distribution) of the

Cosmological abundance of the QCD axion coupled to hidden photons

NASA Astrophysics Data System (ADS)

Kitajima, Naoya; Sekiguchi, Toyokazu; Takahashi, Fuminobu

2018-06-01

We study the cosmological evolution of the QCD axion coupled to hidden photons. For a moderately strong coupling, the motion of the axion field leads to an explosive production of hidden photons by tachyonic instability. We use lattice simulations to evaluate the cosmological abundance of the QCD axion. In doing so, we incorporate the backreaction of the produced hidden photons on the axion dynamics, which becomes significant in the non-linear regime. We find that the axion abundance is suppressed by at most O (102) for the decay constant fa =1016GeV, compared to the case without the coupling. For a sufficiently large coupling, the motion of the QCD axion becomes strongly damped, and as a result, the axion abundance is enhanced. Our results show that the cosmological upper bound on the axion decay constant can be relaxed by a few hundred for a certain range of the coupling to hidden photons.

The Naked Truth: The Face and Body Sensitive N170 Response Is Enhanced for Nude Bodies

PubMed Central

Hietanen, Jari K.; Nummenmaa, Lauri

2011-01-01

Recent event-related potential studies have shown that the occipitotemporal N170 component - best known for its sensitivity to faces - is also sensitive to perception of human bodies. Considering that in the timescale of evolution clothing is a relatively new invention that hides the bodily features relevant for sexual selection and arousal, we investigated whether the early N170 brain response would be enhanced to nude over clothed bodies. In two experiments, we measured N170 responses to nude bodies, bodies wearing swimsuits, clothed bodies, faces, and control stimuli (cars). We found that the N170 amplitude was larger to opposite and same-sex nude vs. clothed bodies. Moreover, the N170 amplitude increased linearly as the amount of clothing decreased from full clothing via swimsuits to nude bodies. Strikingly, the N170 response to nude bodies was even greater than that to faces, and the N170 amplitude to bodies was independent of whether the face of the bodies was visible or not. All human stimuli evoked greater N170 responses than did the control stimulus. Autonomic measurements and self-evaluations showed that nude bodies were affectively more arousing compared to the other stimulus categories. We conclude that the early visual processing of human bodies is sensitive to the visibility of the sex-related features of human bodies and that the visual processing of other people's nude bodies is enhanced in the brain. This enhancement is likely to reflect affective arousal elicited by nude bodies. Such facilitated visual processing of other people's nude bodies is possibly beneficial in identifying potential mating partners and competitors, and for triggering sexual behavior. PMID:22110574

Cosmological ``Truths''

NASA Astrophysics Data System (ADS)

Bothun, Greg

2011-10-01

Ever since Aristotle placed us, with certainty, in the Center of the Cosmos, Cosmological models have more or less operated from a position of known truths for some time. As early as 1963, for instance, it was ``known'' that the Universe had to be 15-17 billion years old due to the suspected ages of globular clusters. For many years, attempts to determine the expansion age of the Universe (the inverse of the Hubble constant) were done against this preconceived and biased notion. Not surprisingly when more precise observations indicated a Hubble expansion age of 11-13 billion years, stellar models suddenly changed to produce a new age for globular cluster stars, consistent with 11-13 billion years. Then in 1980, to solve a variety of standard big bang problems, inflation was introduced in a fairly ad hoc manner. Inflation makes the simple prediction that the net curvature of spacetime is zero (i.e. spacetime is flat). The consequence of introducing inflation is now the necessary existence of a dark matter dominated Universe since the known baryonic material could comprise no more than 1% of the necessary energy density to make spacetime flat. As a result of this new cosmological ``truth'' a significant world wide effort was launched to detect the dark matter (which obviously also has particle physics implications). To date, no such cosmological component has been detected. Moreover, all available dynamical inferences of the mass density of the Universe showed in to be about 20% of that required for closure. This again was inconsistent with the truth that the real density of the Universe was the closure density (e.g. Omega = 1), that the observations were biased, and that 99% of the mass density had to be in the form of dark matter. That is, we know the universe is two component -- baryons and dark matter. Another prevailing cosmological truth during this time was that all the baryonic matter was known to be in galaxies that populated our galaxy catalogs. Subsequent

A charge optimized many-body potential for titanium nitride (TiN).

PubMed

Cheng, Y-T; Liang, T; Martinez, J A; Phillpot, S R; Sinnott, S B

2014-07-02

This work presents a new empirical, variable charge potential for TiN systems in the charge-optimized many-body potential framework. The potential parameters were determined by fitting them to experimental data for the enthalpy of formation, lattice parameters, and elastic constants of rocksalt structured TiN. The potential does a good job of describing the fundamental physical properties (defect formation and surface energies) of TiN relative to the predictions of first-principles calculations. This potential is used in classical molecular dynamics simulations to examine the interface of fcc-Ti(0 0 1)/TiN(0 0 1) and to characterize the adsorption of oxygen atoms and molecules on the TiN(0 0 1) surface. The results indicate that the potential is well suited to model TiN thin films and to explore the chemistry associated with their oxidation.

Simulations of the Formation and Evolution of X-ray Clusters

NASA Astrophysics Data System (ADS)

Bryan, G. L.; Klypin, A.; Norman, M. L.

1994-05-01

We describe results from a set of Omega = 1 Cold plus Hot Dark Matter (CHDM) and Cold Dark Matter (CDM) simulations. We examine the formation and evolution of X-ray clusters in a cosmological setting with sufficient numbers to perform statistical analysis. We find that CDM, normalized to COBE, seems to produce too many large clusters, both in terms of the luminosity (dn/dL) and temperature (dn/dT) functions. The CHDM simulation produces fewer clusters and the temperature distribution (our numerically most secure result) matches observations where they overlap. The computed cluster luminosity function drops below observations, but we are almost surely underestimating the X-ray luminosity. Because of the lower fluctuations in CHDM, there are only a small number of bright clusters in our simulation volume; however we can use the simulated clusters to fix the relation between temperature and velocity dispersion, allowing us to use collisionless N-body codes to probe larger length scales with correspondingly brighter clusters. The hydrodynamic simulations have been performed with a hybrid particle-mesh scheme for the dark matter and a high resolution grid-based piecewise parabolic method for the adiabatic gas dynamics. This combination has been implemented for massively parallel computers, allowing us to achive grids as large as 512(3) .

Is cosmology a historical science?

NASA Astrophysics Data System (ADS)

Grignon, Claude

2012-06-01

To explain the formation and the evolution of the Universe, cosmology settles universal laws. In this respect, cosmology belongs to the category of the nomothetic sciences, which write and think in mathematics. But cosmology is also akin to the historical sciences; like archaeology, geology or the biology of evolution, cosmology infers history from the vestiges of the past; moreover, it is not an experimental but an observational science. Due to this ambivalence, cosmology confronts divergent epistemological options. Nomothetic and historical sciences use indeed different, even opposite conceptions of such fundamental notions as time, causality and chance. Is it possible to make the history of the Universe intelligible without referring to the narrative conception of history congruent with the course of the historical world?.

On under-determination in cosmology

NASA Astrophysics Data System (ADS)

Butterfield, Jeremy

2014-05-01

I discuss how modern cosmology illustrates under-determination of theoretical hypotheses by data, in ways that are different from most philosophical discussions. I emphasise cosmology's concern with what data could in principle be collected by a single observer (Section 2); and I give a broadly sceptical discussion of cosmology's appeal to the cosmological principle as a way of breaking the under-determination (Section 3). I confine most of the discussion to the history of the observable universe from about one second after the Big Bang, as described by the mainstream cosmological model: in effect, what cosmologists in the early 1970s dubbed the 'standard model', as elaborated since then. But in the closing Section 4, I broach some questions about times earlier than one second.

Cosmological Entropy Bounds

NASA Astrophysics Data System (ADS)

Brustein, R.

I review some basic facts about entropy bounds in general and about cosmological entropy bounds. Then I review the causal entropy bound, the conditions for its validity and its application to the study of cosmological singularities. This article is based on joint work with Gabriele Veneziano and subsequent related research.

Real- and redshift-space halo clustering in f(R) cosmologies

NASA Astrophysics Data System (ADS)

Arnalte-Mur, Pablo; Hellwing, Wojciech A.; Norberg, Peder

2017-05-01

We present two-point correlation function statistics of the mass and the haloes in the chameleon f(R) modified gravity scenario using a series of large-volume N-body simulations. Three distinct variations of f(R) are considered (F4, F5 and F6) and compared to a fiducial Λ cold dark matter (ΛCDM) model in the redshift range z ∈ [0, 1]. We find that the matter clustering is indistinguishable for all models except for F4, which shows a significantly steeper slope. The ratio of the redshift- to real-space correlation function at scales >20 h-1 Mpc agrees with the linear General Relativity (GR) Kaiser formula for the viable f(R) models considered. We consider three halo populations characterized by spatial abundances comparable to that of luminous red galaxies and galaxy clusters. The redshift-space halo correlation functions of F4 and F5 deviate significantly from ΛCDM at intermediate and high redshift, as the f(R) halo bias is smaller than or equal to that of the ΛCDM case. Finally, we introduce a new model-independent clustering statistic to distinguish f(R) from GR: the relative halo clustering ratio - R. The sampling required to adequately reduce the scatter in R will be available with the advent of the next-generation galaxy redshift surveys. This will foster a prospective avenue to obtain largely model-independent cosmological constraints on this class of modified gravity models.

An Introduction to Galaxies and Cosmology

NASA Astrophysics Data System (ADS)

Jones, Mark H.; Lambourne, Robert J. A.; Serjeant, Stephen

2015-01-01

Introduction; 1. The Milky Way - our galaxy; 2. Normal galaxies; 3. Active galaxies; 4. The spatial distribution of galaxies; 5. Introducing cosmology - the science of the Universe; 6. Big bang cosmology - the evolving Universe; 7. Observational cosmology - measuring the Universe; 8. Questioning cosmology - outstanding problems about the Universe; Answers and comments; Appendix; Glossary; Further reading; Acknowledgements; Figure references; Index.

cosmoabc: Likelihood-free inference for cosmology

NASA Astrophysics Data System (ADS)

Ishida, Emille E. O.; Vitenti, Sandro D. P.; Penna-Lima, Mariana; Trindade, Arlindo M.; Cisewski, Jessi; M.; de Souza, Rafael; Cameron, Ewan; Busti, Vinicius C.

2015-05-01

Approximate Bayesian Computation (ABC) enables parameter inference for complex physical systems in cases where the true likelihood function is unknown, unavailable, or computationally too expensive. It relies on the forward simulation of mock data and comparison between observed and synthetic catalogs. cosmoabc is a Python Approximate Bayesian Computation (ABC) sampler featuring a Population Monte Carlo variation of the original ABC algorithm, which uses an adaptive importance sampling scheme. The code can be coupled to an external simulator to allow incorporation of arbitrary distance and prior functions. When coupled with the numcosmo library, it has been used to estimate posterior probability distributions over cosmological parameters based on measurements of galaxy clusters number counts without computing the likelihood function.

Philosophical aspects of modern cosmology

NASA Astrophysics Data System (ADS)

Zinkernagel, Henrik

2014-05-01

Cosmology is the attempt to understand in scientific terms the structure and evolution of the universe as a whole. This ambition has been with us since the ancient Greeks, even if the developments in modern cosmology have provided a picture of the universe dramatically different from that of Pythagoras, Plato and Aristotle. The cosmological thinking of these figures, e.g. the belief in uniform circular motion of the heavens, was closely related to their philosophical ideas, and it shaped the field of cosmology at least up to the times of Copernicus and Kepler.

Multibody Simulation Software Testbed for Small-Body Exploration and Sampling

NASA Technical Reports Server (NTRS)

Acikmese, Behcet; Blackmore, James C.; Mandic, Milan

2011-01-01

G-TAG is a software tool for the multibody simulation of a spacecraft with a robotic arm and a sampling mechanism, which performs a touch-and-go (TAG) maneuver for sampling from the surface of a small celestial body. G-TAG utilizes G-DYN, a multi-body simulation engine described in the previous article, and interfaces to controllers, estimators, and environmental forces that affect the spacecraft. G-TAG can easily be adapted for the analysis of the mission stress cases to support the design of a TAG system, as well as for comprehensive Monte Carlo simulations to analyze and evaluate a particular TAG system design. Any future small-body mission will benefit from using G-TAG, which has already been extensively used in Comet Odyssey and Galahad Asteroid New Frontiers proposals.

Galactic Angular Momentum in Cosmological Zoom-in Simulations. I. Disk and Bulge Components and the Galaxy-Halo Connection

NASA Astrophysics Data System (ADS)

Sokołowska, Aleksandra; Capelo, Pedro R.; Fall, S. Michael; Mayer, Lucio; Shen, Sijing; Bonoli, Silvia

2017-02-01

We investigate the angular momentum evolution of four disk galaxies residing in Milky-Way-sized halos formed in cosmological zoom-in simulations with various sub-grid physics and merging histories. We decompose these galaxies, kinematically and photometrically, into their disk and bulge components. The simulated galaxies and their components lie on the observed sequences in the j *-M * diagram, relating the specific angular momentum and mass of the stellar component. We find that galaxies in low-density environments follow the relation {j}* \\propto {M}* α past major mergers, with α ˜ 0.6 in the case of strong feedback, when bulge-to-disk ratios are relatively constant, and α ˜ 1.4 in the other cases, when secular processes operate on shorter timescales. We compute the retention factors (I.e., the ratio of the specific angular momenta of stars and dark matter) for both disks and bulges and show that they vary relatively slowly after averaging over numerous but brief fluctuations. For disks, the retention factors are usually close to unity, while for bulges, they are a few times smaller. Our simulations therefore indicate that galaxies and their halos grow in a quasi-homologous way.

Systematic simulations of modified gravity: chameleon models

DOE Office of Scientific and Technical Information (OSTI.GOV)

Brax, Philippe; Davis, Anne-Christine; Li, Baojiu

2013-04-01

In this work we systematically study the linear and nonlinear structure formation in chameleon theories of modified gravity, using a generic parameterisation which describes a large class of models using only 4 parameters. For this we have modified the N-body simulation code ecosmog to perform a total of 65 simulations for different models and parameter values, including the default ΛCDM. These simulations enable us to explore a significant portion of the parameter space. We have studied the effects of modified gravity on the matter power spectrum and mass function, and found a rich and interesting phenomenology where the difference withmore » the ΛCDM paradigm cannot be reproduced by a linear analysis even on scales as large as k ∼ 0.05 hMpc{sup −1}, since the latter incorrectly assumes that the modification of gravity depends only on the background matter density. Our results show that the chameleon screening mechanism is significantly more efficient than other mechanisms such as the dilaton and symmetron, especially in high-density regions and at early times, and can serve as a guidance to determine the parts of the chameleon parameter space which are cosmologically interesting and thus merit further studies in the future.« less

Possibility of determining cosmological parameters from measurements of gravitational waves emitted by coalescing, compact binaries

NASA Astrophysics Data System (ADS)

Marković, Dragoljub

1993-11-01

We explore the feasibility of using LIGO and/or VIRGO gravitational-wave measurements of coalescing, neutron-star-neutron-star (NS-NS) binaries and black-hole-neutron-star (BH-NS) binaries at cosmological distances to determine the cosmological parameters of our Universe. From the observed gravitational waveforms one can infer, as direct observables, the luminosity distance D of the source and the binary's two ``redshifted masses,'' M'1≡M1(1+z) and M'2≡M2(1+z), where Mi are the actual masses and z≡Δλ/λ is the binary's cosmological redshift. Assuming that the NS mass spectrum is sharply peaked about 1.4Msolar, as binary pulsar and x-ray source observations suggest, the redshift can be estimated as z=M'NS/1.4Msolar-1. The actual distance-redshift relation D(z) for our Universe is strongly dependent on its cosmological parameters [the Hubble constant H0, or h0≡H0/100 km s-1Mpc-1, the mean mass density ρm, or density parameter Ω0≡(8π/3H20)ρm, and the cosmological constant, Λ, or λ0≡Λ/(3H20)], so by a statistical study of (necessarily noisy) measurements of D and z for a large number of binaries, one can deduce the cosmological parameters. The various noise sources that will plague such a cosmological study are discussed and estimated, and the accuracies of the inferred parameters are determined as functions of the detectors' noise characteristics, the number of binaries observed, and the neutron-star mass spectrum. The dominant source of error is the detectors' intrinsic noise, though stochastic gravitational lensing of the waves by intervening matter might significantly influence the inferred cosmological constant λ0, when the detectors reach ``advanced'' stages of development. The estimated errors of parameters inferred from BH-NS measurements can be described by the following rough analytic fits: Δh0/h0~=0.02(N/h0)(τR)-1/2 (for N/h0<~2), where N is the detector's noise level (strain/Hz) in units of the ``advanced LIGO'' noise level, R is the

Simulations of the pairwise kinematic Sunyaev-Zel'dovich signal

DOE PAGES

Flender, Samuel; Bleem, Lindsey; Finkel, Hal; ...

2016-05-26

The pairwise kinematic Sunyaev–Zel'dovich (kSZ) signal from galaxy clusters is a probe of their line of sight momenta, and thus a potentially valuable source of cosmological information. In addition to the momenta, the amplitude of the measured signal depends on the properties of the intracluster gas and observational limitations such as errors in determining cluster centers and redshifts. In this work, we simulate the pairwise kSZ signal of clusters atmore » $$z\\lt 1$$, using the output from a cosmological N-body simulation and including the properties of the intracluster gas via a model that can be varied in post-processing. We find that modifications to the gas profile due to star formation and feedback reduce the pairwise kSZ amplitude of clusters by $$\\sim 50\\%$$, relative to the naive "gas traces mass" assumption. We demonstrate that miscentering can reduce the overall amplitude of the pairwise kSZ signal by up to 10%, while redshift errors can lead to an almost complete suppression of the signal at small separations. We confirm that a high-significance detection is expected from the combination of data from current generation, high-resolution cosmic microwave background experiments, such as the South Pole Telescope, and cluster samples from optical photometric surveys, such as the Dark Energy Survey. As a result, we forecast that future experiments such as Advanced ACTPol in conjunction with data from the Dark Energy Spectroscopic Instrument will yield detection significances of at least $$20\\sigma $$, and up to $$57\\sigma $$ in an optimistic scenario.« less

Numerical Results for a Polytropic Cosmology Interpreted as a Dust Universe Producing Gravitational Waves

NASA Astrophysics Data System (ADS)

Klapp, J.; Cervantes-Cota, J.; Chauvet, P.

1990-11-01

RESUMEN. A nivel cosmol6gico pensamos que se ha estado prodticiendo radiaci6n gravitacional en cantidades considerables dentro de las galaxias. Si los eventos prodnctores de radiaci6n gravitatoria han venido ocurriendo desde Ia epoca de Ia formaci6n de las galaxias, cuando menos, sus efectos cosmol6gicos pueden ser tomados en cuenta con simplicidad y elegancia al representar la producci6n de radiaci6n y, por consiguiente, su interacci6n con materia ordinaria fenomenol6gicamente a trave's de una ecuaci6n de estado politr6pica, como lo hemos mostrado en otros trabajos. Presentamos en este articulo resultados nunericos de este modelo. ABSTRACT A common believe in cosmology is that gravitational radiation in considerable quantities is being produced within the galaxies. Ifgravitational radiation production has been running since the galaxy formation epoch, at least, its cosmological effects can be assesed with simplicity and elegance by representing the production of radiation and, therefore, its interaction with ordinary matter phenomenologically through a polytropic equation of state as shown already elsewhere. We present in this paper the numerical results of such a model. K words: COSMOLOGY - GRAVITATION

A THOUSAND SHADOWS OF ANDROMEDA: ROTATING PLANES OF SATELLITES IN THE MILLENNIUM-II COSMOLOGICAL SIMULATION

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ibata, Rodrigo A.; Martin, Nicolas F.; Ibata, Neil G.

In a recent contribution, Bahl and Baumgardt investigated the incidence of planar alignments of satellite galaxies in the Millennium-II simulation and concluded that vast, thin planes of dwarf galaxies, similar to that observed in the Andromeda galaxy (M31), occur frequently by chance in Λ-cold dark matter cosmology. However, their analysis did not capture the essential fact that the observed alignment is simultaneously radially extended, yet thin, and kinematically unusual. With the caveat that the Millennium-II simulation may not have sufficient mass resolution to identify confidently simulacra of low-luminosity dwarf galaxies, we re-examine that simulation for planar structures, using the samemore » method as employed by Ibata et al. on the real M31 satellites. We find that 0.04% of host galaxies display satellite alignments that are at least as extreme as the observations, when we consider their extent, thickness, and number of members rotating in the same sense. We further investigate the angular momentum properties of the co-planar satellites, and find that the median of the specific angular momentum derived from the line-of-sight velocities in the real M31 structure (1.3 × 10{sup 4} km s{sup –1} kpc) is very high compared to systems drawn from the simulations. This analysis confirms that it is highly unlikely that the observed structure around the Andromeda galaxy is due to a chance occurrence. Interestingly, the few extreme systems that are similar to M31 arise from the accretion of a massive sub-halo with its own spatially concentrated entourage of orphan satellites.« less

Graviton fluctuations erase the cosmological constant

NASA Astrophysics Data System (ADS)

Wetterich, C.

2017-10-01

Graviton fluctuations induce strong non-perturbative infrared renormalization effects for the cosmological constant. The functional renormalization flow drives a positive cosmological constant towards zero, solving the cosmological constant problem without the need to tune parameters. We propose a simple computation of the graviton contribution to the flow of the effective potential for scalar fields. Within variable gravity, with effective Planck mass proportional to the scalar field, we find that the potential increases asymptotically at most quadratically with the scalar field. The solutions of the derived cosmological equations lead to an asymptotically vanishing cosmological "constant" in the infinite future, providing for dynamical dark energy in the present cosmological epoch. Beyond a solution of the cosmological constant problem, our simplified computation also entails a sizeable positive graviton-induced anomalous dimension for the quartic Higgs coupling in the ultraviolet regime, substantiating the successful prediction of the Higgs boson mass within the asymptotic safety scenario for quantum gravity.

Was Newtonian cosmology really inconsistent?

NASA Astrophysics Data System (ADS)

Vickers, Peter

This paper follows up a debate as to the consistency of Newtonian cosmology. Whereas Malament [(1995). Is Newtonian cosmology really inconsistent? Philosophy of Science 62, 489-510] has shown that Newtonian cosmology is not inconsistent, to date there has been no analysis of Norton's claim [(1995). The force of Newtonian cosmology: Acceleration is relative. Philosophy of Science 62, 511-522.] that Newtonian cosmology was inconsistent prior to certain advances in the 1930s, and in particular prior to Seeliger's seminal paper of Seeliger [(1895). Über das Newton'sche Gravitationsgesetz. Astronomische Nachrichten 137 (3273), 129-136.] In this paper I agree that there are assumptions, Newtonian and cosmological in character, and relevant to the real history of science, which are inconsistent. But there are some important corrections to make to Norton's account. Here I display for the first time the inconsistencies-four in total-in all their detail. Although this extra detail shows there to be several different inconsistencies, it also goes some way towards explaining why they went unnoticed for 200 years.

Warps and waves in the stellar discs of the Auriga cosmological simulations

NASA Astrophysics Data System (ADS)

Gómez, Facundo A.; White, Simon D. M.; Grand, Robert J. J.; Marinacci, Federico; Springel, Volker; Pakmor, Rüdiger

2017-03-01

Recent studies have revealed an oscillating asymmetry in the vertical structure of the Milky Way's disc. Here, we analyse 16 high-resolution, fully cosmological simulations of the evolution of individual Milky Way-sized galaxies, carried out with the magnetohydrodynamic code AREPO. At redshift zero, about 70 per cent of our galactic discs show strong vertical patterns, with amplitudes that can exceed 2 kpc. Half of these are typical 'integral sign' warps. The rest are oscillations similar to those observed in the Milky Way. Such structures are thus expected to be common. The associated mean vertical motions can be as large as 30 km s-1. Cold disc gas typically follows the vertical patterns seen in the stars. These perturbations have a variety of causes: close encounters with satellites, distant fly-bys of massive objects, accretion of misaligned cold gas from halo infall or from mergers. Tidally induced vertical patterns can be identified in both young and old stellar populations, whereas those originating from cold gas accretion are seen mainly in the younger populations. Galaxies with regular or at most weakly perturbed discs are usually, but not always, free from recent interactions with massive companions, although we have one case where an equilibrium compact disc reforms after a merger.

Cosmology. A first course

NASA Astrophysics Data System (ADS)

Lachieze-Rey, Marc

This book delivers a quantitative account of the science of cosmology, designed for a non-specialist audience. The basic principles are outlined using simple maths and physics, while still providing rigorous models of the Universe. It offers an ideal introduction to the key ideas in cosmology, without going into technical details. The approach used is based on the fundamental ideas of general relativity such as the spacetime interval, comoving coordinates, and spacetime curvature. It provides an up-to-date and thoughtful discussion of the big bang, and the crucial questions of structure and galaxy formation. Questions of method and philosophical approaches in cosmology are also briefly discussed. Advanced undergraduates in either physics or mathematics would benefit greatly from use either as a course text or as a supplementary guide to cosmology courses.

Properties of Galaxies and Groups: Nature versus Nurture

NASA Astrophysics Data System (ADS)

Niemi, Sami-Matias

2011-09-01

Due to the inherently nonlinear nature of gravity cosmological N-body simulations have become an invaluable tool when the growth of structure is being studied and modelled closer to the present epoch. Large simulations with high dynamical range have made it possible to model the formation and growth of cosmic structure with unprecedented accuracy. Moreover, galaxies, the basic building blocks of the Universe, can also be modelled in cosmological context. However, despite all the simulations and successes in recent decades, there are still many unanswered questions in the field of galaxy formation and evolution. One of the longest standing issue being the significance of the formation place and thus initial conditions to a galaxy's evolution in respect to environment, often formulated simply as "nature versus nurture" like in human development and psychology. Unfortunately, our understanding of galaxy evolution in different environments is still limited, albeit, for example, the morphology-density relation has shown that the density of the galaxy's local environment can affect its properties. Consequently, the environment should play a role in galaxy evolution, however despite the efforts, the exact role of the galaxy's local environment to its evolution remains open. This thesis introduction discusses briefly the background cosmology, cosmological N-body simulations and semi-analytical models. The second part is reserved for groups of galaxies, whether they are gravitationally bound, and what this may imply for galaxy evolution. The third part of the thesis concentrates on describing results of a case study of isolated field elliptical galaxies. The final chapter discusses another case study of luminous infra-red galaxies.

ICE-COLA: towards fast and accurate synthetic galaxy catalogues optimizing a quasi-N-body method

NASA Astrophysics Data System (ADS)

Izard, Albert; Crocce, Martin; Fosalba, Pablo

2016-07-01

Next generation galaxy surveys demand the development of massive ensembles of galaxy mocks to model the observables and their covariances, what is computationally prohibitive using N-body simulations. COmoving Lagrangian Acceleration (COLA) is a novel method designed to make this feasible by following an approximate dynamics but with up to three orders of magnitude speed-ups when compared to an exact N-body. In this paper, we investigate the optimization of the code parameters in the compromise between computational cost and recovered accuracy in observables such as two-point clustering and halo abundance. We benchmark those observables with a state-of-the-art N-body run, the MICE Grand Challenge simulation. We find that using 40 time-steps linearly spaced since zI ˜ 20, and a force mesh resolution three times finer than that of the number of particles, yields a matter power spectrum within 1 per cent for k ≲ 1 h Mpc-1 and a halo mass function within 5 per cent of those in the N-body. In turn, the halo bias is accurate within 2 per cent for k ≲ 0.7 h Mpc-1 whereas, in redshift space, the halo monopole and quadrupole are within 4 per cent for k ≲ 0.4 h Mpc-1. These results hold for a broad range in redshift (0 < z < 1) and for all halo mass bins investigated (M > 1012.5 h-1 M⊙). To bring accuracy in clustering to one per cent level we study various methods that re-calibrate halo masses and/or velocities. We thus propose an optimized choice of COLA code parameters as a powerful tool to optimally exploit future galaxy surveys.

Astrometric cosmology .

NASA Astrophysics Data System (ADS)

Lattanzi, M. G.

The accurate measurement of the motions of stars in our Galaxy can provide access to the cosmological signatures in the disk and halo, while astrometric experiments from within our Solar System can uniquely probe possible deviations from General Relativity. This article will introduce to the fact that astrometry has the potential, thanks also to impressive technological advancements, to become a key player in the field of local cosmology. For example, accurate absolute kinematics at the scale of the Milky Way can, for the first time in situ, account for the predictions made by the cold dark matter model for the Galactic halo, and eventually map out the distribution of dark matter, or other formation mechanisms, required to explain the signatures recently identified in the old component of the thick disk. Final notes dwell on to what extent Gaia can fulfill the expectations of astrometric cosmology and on what must instead be left to future, specifically designed, astrometric experiments.

Inhomogeneous anisotropic cosmology

DOE PAGES

Kleban, Matthew; Senatore, Leonardo

2016-10-12

In homogeneous and isotropic Friedmann-Robertson-Walker cosmology, the topology of the universe determines its ultimate fate. If the Weak Energy Condition is satisfied, open and flat universes must expand forever, while closed cosmologies can recollapse to a Big Crunch. A similar statement holds for homogeneous but anisotropic (Bianchi) universes. Here in this paper, we prove that arbitrarily inhomogeneous and anisotropic cosmologies with "flat'' (including toroidal) and "open'' (including compact hyperbolic) spatial topology that are initially expanding must continue to expand forever at least in some region at a rate bounded from below by a positive number, despite the presence of arbitrarilymore » large density fluctuations and/or the formation of black holes. Because the set of 3-manifold topologies is countable, a single integer determines the ultimate fate of the universe, and, in a specific sense, most 3-manifolds are "flat" or "open". Our result has important implications for inflation: if there is a positive cosmological constant (or suitable inflationary potential) and initial conditions for the inflaton, cosmologies with "flat'' or "open" topology must expand forever in some region at least as fast as de Sitter space, and are therefore very likely to begin inflationary expansion eventually, regardless of the scale of the inflationary energy or the spectrum and amplitude of initial inhomogeneities and gravitational waves. Our result is also significant for numerical general relativity, which often makes use of periodic (toroidal) boundary conditions.« less

Inhomogeneous anisotropic cosmology

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kleban, Matthew; Senatore, Leonardo

In homogeneous and isotropic Friedmann-Robertson-Walker cosmology, the topology of the universe determines its ultimate fate. If the Weak Energy Condition is satisfied, open and flat universes must expand forever, while closed cosmologies can recollapse to a Big Crunch. A similar statement holds for homogeneous but anisotropic (Bianchi) universes. Here in this paper, we prove that arbitrarily inhomogeneous and anisotropic cosmologies with "flat'' (including toroidal) and "open'' (including compact hyperbolic) spatial topology that are initially expanding must continue to expand forever at least in some region at a rate bounded from below by a positive number, despite the presence of arbitrarilymore » large density fluctuations and/or the formation of black holes. Because the set of 3-manifold topologies is countable, a single integer determines the ultimate fate of the universe, and, in a specific sense, most 3-manifolds are "flat" or "open". Our result has important implications for inflation: if there is a positive cosmological constant (or suitable inflationary potential) and initial conditions for the inflaton, cosmologies with "flat'' or "open" topology must expand forever in some region at least as fast as de Sitter space, and are therefore very likely to begin inflationary expansion eventually, regardless of the scale of the inflationary energy or the spectrum and amplitude of initial inhomogeneities and gravitational waves. Our result is also significant for numerical general relativity, which often makes use of periodic (toroidal) boundary conditions.« less

Inhomogeneous anisotropic cosmology

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kleban, Matthew; Senatore, Leonardo; Kavli Institute for Particle Astrophysics and Cosmology, Stanford University and SLAC,2575 Sand Hill Road, M/S 29, Menlo Park, CA 94025

In homogeneous and isotropic Friedmann-Robertson-Walker cosmology, the topology of the universe determines its ultimate fate. If the Weak Energy Condition is satisfied, open and flat universes must expand forever, while closed cosmologies can recollapse to a Big Crunch. A similar statement holds for homogeneous but anisotropic (Bianchi) universes. Here, we prove that arbitrarily inhomogeneous and anisotropic cosmologies with “flat” (including toroidal) and “open” (including compact hyperbolic) spatial topology that are initially expanding must continue to expand forever at least in some region at a rate bounded from below by a positive number, despite the presence of arbitrarily large density fluctuationsmore » and/or the formation of black holes. Because the set of 3-manifold topologies is countable, a single integer determines the ultimate fate of the universe, and, in a specific sense, most 3-manifolds are “flat” or “open”. Our result has important implications for inflation: if there is a positive cosmological constant (or suitable inflationary potential) and initial conditions for the inflaton, cosmologies with “flat” or “open” topology must expand forever in some region at least as fast as de Sitter space, and are therefore very likely to begin inflationary expansion eventually, regardless of the scale of the inflationary energy or the spectrum and amplitude of initial inhomogeneities and gravitational waves. Our result is also significant for numerical general relativity, which often makes use of periodic (toroidal) boundary conditions.« less

Constraining the Intergalactic and Circumgalactic Media with Lyman-Alpha Absorption

NASA Astrophysics Data System (ADS)

Sorini, Daniele; Onorbe, Jose; Hennawi, Joseph F.; Lukic, Zarija

2018-01-01

Lyman-alpha (Ly-a) absorption features detected in quasar spectra in the redshift range 0cosmology. In the first part of my thesis, I overcome numerical challenges posed by cosmological hydrodynamic simulations by developing a novel semi-analytic technique to predict various statistics of Ly-a absorption in the IGM with large N-body cosmological simulations. The technique developed is more accurate than previous attempts in the literature, and can be applied to Gpc-scale N-body simulations, allowing an accurate investigation of the Ly-a absorption on unprecedentedly large scales. In the second part of my thesis, I compare predictions of state-of-the-art hydrodynamic cosmological simulations with observations of the mean Ly-a absorption around foreground quasars, damped Ly-a absorbers, and Lyman-break galaxies, at different transverse distances (~20kpc-20Mpc) from background quasars. Far from galaxies >2Mpc, the simulations asymptotically match the observations, because the ΛCDM model successfully describes the ambient IGM. This represents a critical advantage of studying the mean absorption profile. However, significant differences between the simulations, and between simulations and observations are present on scales 20kpc-2Mpc, illustrating the challenges of accurately modeling and resolving galaxy formation physics. It is noteworthy that these differences are observed as far out as ~2Mpc, indicating that the `sphere-of-influence' of galaxies could extend to approximately ~20 times the halo virial radius (~100kpc). Current observations are very precise on these scales and can thus strongly discriminate between different galaxy formation models. I demonstrate that the Ly-a absorption profile is primarily sensitive to the underlying temperature-density relationship of diffuse gas around galaxies, and argue that it

The Evolution of Dwarf Galaxy Satellites with Different Dark Matter Density Profiles in the ErisMod Simulations. I. The Early Infalls

NASA Astrophysics Data System (ADS)

Tomozeiu, Mihai; Mayer, Lucio; Quinn, Thomas

2016-02-01

We present the first simulations of tidal stirring of dwarf galaxies in the Local Group carried out in a fully cosmological context. We use the ErisDARK cosmological simulation of a Milky Way (MW)-sized galaxy to identify some of the most massive subhalos (Mvir > 108 M⊙) that fall into the main host before z = 2. Subhalos are replaced before infall with extremely high-resolution models of dwarf galaxies comprising a faint stellar disk embedded in a dark matter halo. The set of models contains cuspy halos as well as halos with “cored” profiles (with the cusp coefficient γ = 0.6) consistent with recent results of hydrodynamical simulations of dwarf galaxy formation. The simulations are then run to z = 0 with as many as 54 million particles and resolutions as small as ∼4 pc using the new parallel N-body code ChaNGa. The stellar components of all satellites are significantly affected by tidal stirring, losing stellar mass, and undergoing a morphological transformation toward a pressure supported spheroidal system. However, while some remnants with cuspy halos maintain significant rotational flattening and disk-like features, all the shallow halo models achieve vrot/σ⋆ < 0.5 and round shapes typical of dSph satellites of the MW and M31. Mass loss is also enhanced in the latter, and remnants can reach luminosities and velocity dispersions as low as those of ultra-faint dwarfs.

Comparison and Historical Evolution of Ancient Greek Cosmological Ideas and Mathematical Models

NASA Astrophysics Data System (ADS)

Pinotsis, Antonios D.

2005-12-01

We present a comparative study of the cosmological ideas and mathematical models in ancient Greece. We show that the heliocentric system introduced by Aristarchus of Samos was the outcome of much intellectual activity. Many Greek philosophers, mathematicians and astronomers such as Anaximander, Philolaus, Hicetas, Ecphantus and Heraclides of Pontus contributed to this. Also, Ptolemy was influenced by the cosmological model of Heraclides of Pontus for the explanation of the apparent motions of Mercury and Venus. Apollonius, who wrote the definitive work on conic sections, introduced the theory of eccentric circles and implemented them together with epicycles instead of considering that the celestial bodies travel in elliptic orbits. This is due to the deeply rooted belief that the orbits of the celestial bodies were normal circular motions around the Earth, which was still. There was also a variety of important ideas which are relevant to modern science. We present the ideas of Plato that are consistent with modern relativity theories, as well as Aristarchus' estimations of the size of the Universe in comparison with the size of the planetary system. As a first approximation, Hipparchus' theory of eccentric circles was equivalent to the first two laws of Kepler. The significance of the principle of independence and superposition of motions in the formulation of ancient cosmological models is also clarified.

Simulating the cold dark matter-neutrino dipole with TianNu

DOE PAGES

Inman, Derek; Yu, Hao-Ran; Zhu, Hong-Ming; ...

2017-04-20

Measurements of neutrino mass in cosmological observations rely on two-point statistics that are hindered by significant degeneracies with the optical depth and galaxy bias. The relative velocity effect between cold dark matter and neutrinos induces a large scale dipole in the matter density field and may be able to provide orthogonal constraints to standard techniques. In this paper, we numerically investigate this dipole in the TianNu simulation, which contains cold dark matter and 50 meV neutrinos. We first compute the dipole using a new linear response technique where we treat the displacement caused by the relative velocity as a phasemore » in Fourier space and then integrate the matter power spectrum over redshift. Then, we compute the dipole numerically in real space using the simulation density and velocity fields. We find excellent agreement between the linear response and N-body methods. Finally, utilizing the dipole as an observational tool requires two tracers of the matter distribution that are differently biased with respect to the neutrino density.« less

Parameterized post-Newtonian cosmology

NASA Astrophysics Data System (ADS)

Sanghai, Viraj A. A.; Clifton, Timothy

2017-03-01

Einstein’s theory of gravity has been extensively tested on solar system scales, and for isolated astrophysical systems, using the perturbative framework known as the parameterized post-Newtonian (PPN) formalism. This framework is designed for use in the weak-field and slow-motion limit of gravity, and can be used to constrain a large class of metric theories of gravity with data collected from the aforementioned systems. Given the potential of future surveys to probe cosmological scales to high precision, it is a topic of much contemporary interest to construct a similar framework to link Einstein’s theory of gravity and its alternatives to observations on cosmological scales. Our approach to this problem is to adapt and extend the existing PPN formalism for use in cosmology. We derive a set of equations that use the same parameters to consistently model both weak fields and cosmology. This allows us to parameterize a large class of modified theories of gravity and dark energy models on cosmological scales, using just four functions of time. These four functions can be directly linked to the background expansion of the universe, first-order cosmological perturbations, and the weak-field limit of the theory. They also reduce to the standard PPN parameters on solar system scales. We illustrate how dark energy models and scalar-tensor and vector-tensor theories of gravity fit into this framework, which we refer to as ‘parameterized post-Newtonian cosmology’ (PPNC).

Effects of 17β-estradiol on emissions of greenhouse gases in simulative natural water body.

PubMed

Ruan, Aidong; Zhao, Ying; Liu, Chenxiao; Zong, Fengjiao; Yu, Zhongbo

2015-05-01

Environmental estrogens are widely spread across the world and are increasingly thought of as serious contaminators. The present study looks at the influence of different concentrations of 17β-estradiol on greenhouse gas emissions (CO2 , CH4 , and N2 O) in simulated systems to explore the relationship between environmental estrogen-pollution and greenhouse gas emissions in natural water bodies. The present study finds that 17β-estradiol pollution in simulated systems has significant promoting effects on the emissions of CH4 and CO2 , although no significant effects on N2 O emissions. The present study indicates that 17β-estradiol has different effects on the different elements cycles; the mechanism of microbial ecology is under review. © 2015 SETAC.

Cosmological constraints on interacting light particles

DOE Office of Scientific and Technical Information (OSTI.GOV)

Brust, Christopher; Cui, Yanou; Sigurdson, Kris, E-mail: cbrust@perimeterinstitute.ca, E-mail: yanou.cui@ucr.edu, E-mail: krs@phas.ubc.ca

2017-08-01

Cosmological observations are becoming increasingly sensitive to the effects of light particles in the form of dark radiation (DR) at the time of recombination. The conventional observable of effective neutrino number, N {sub eff}, is insufficient for probing generic, interacting models of DR. In this work, we perform likelihood analyses which allow both free-streaming effective neutrinos (parametrized by N {sub eff}) and interacting effective neutrinos (parametrized by N {sub fld}). We motivate an alternative parametrization of DR in terms of N {sub tot} (total effective number of neutrinos) and f {sub fs} (the fraction of effective neutrinos which are free-streaming),more » which is less degenerate than using N {sub eff} and N {sub fld}. Using the Planck 2015 likelihoods in conjunction with measurements of baryon acoustic oscillations (BAO), we find constraints on the total amount of beyond the Standard Model effective neutrinos (both free-streaming and interacting) of Δ N {sub tot} < 0.39 at 2σ. In addition, we consider the possibility that this scenario alleviates the tensions between early-time and late-time cosmological observations, in particular the measurements of σ{sub 8} (the amplitude of matter power fluctuations at 8 h {sup −1} Mpc), finding a mild preference for interactions among light species. We further forecast the sensitivities of a variety of future experiments, including Advanced ACTPol (a representative CMB Stage-III experiment), CMB Stage-IV, and the Euclid satellite. This study is relevant for probing non-standard neutrino physics as well as a wide variety of new particle physics models beyond the Standard Model that involve dark radiation.« less

Differentiating G-inflation from string gas cosmology using the effective field theory approach

DOE Office of Scientific and Technical Information (OSTI.GOV)

He, Minxi; Liu, Junyu; Lu, Shiyun

A characteristic signature of String Gas Cosmology is primordial power spectra for scalar and tensor modes which are almost scale-invariant but with a red tilt for scalar modes but a blue tilt for tensor modes. This feature, however, can also be realized in the so-called G-inflation model, in which Horndeski operators are introduced which leads to a blue tensor tilt by softly breaking the Null Energy Condition. In this article we search for potential observational differences between these two cosmologies by performing detailed perturbation analyses based on the Effective Field Theory approach. Our results show that, although both two modelsmore » produce blue tilted tensor perturbations, they behave differently in three aspects. Firstly, String Gas Cosmology predicts a specific consistency relation between the index of the scalar modes n {sub s} and that of tensor ones n {sub t} , which is hard to be reproduced by G-inflation. Secondly, String Gas Cosmology typically predicts non-Gaussianities which are highly suppressed on observable scales, while G-inflation gives rise to observationally large non-Gaussianities because the kinetic terms in the action become important during inflation. However, after finely tuning the model parameters of G-inflation it is possible to obtain a blue tensor spectrum and negligible non-Gaussianities with a degeneracy between the two models. This degeneracy can be broken by a third observable, namely the scale dependence of the nonlinearity parameter, which vanishes for G-inflation but has a blue tilt in the case of String Gas Cosmology. Therefore, we conclude that String Gas Cosmology is in principle observationally distinguishable from the single field inflationary cosmology, even allowing for modifications such as G-inflation.« less

Planck 2015 results: XIII. Cosmological parameters

DOE PAGES

Ade, P. A. R.; Aghanim, N.; Arnaud, M.; ...

2016-09-20

Here, this paper presents cosmological results based on full-mission Planck observations of temperature and polarization anisotropies of the cosmic microwave background (CMB) radiation. Our results are in very good agreement with the 2013 analysis of the Planck nominal-mission temperature data, but with increased precision. The temperature and polarization power spectra are consistent with the standard spatially-flat 6-parameter ΛCDM cosmology with a power-law spectrum of adiabatic scalar perturbations (denoted “base ΛCDM” in this paper). From the Planck temperature data combined with Planck lensing, for this cosmology we find a Hubble constant, H 0 = (67.8 ± 0.9) km s -1Mpc -1, a matter density parameter Ω m = 0.308 ± 0.012, and a tilted scalar spectral index with n s = 0.968 ± 0.006, consistent with the 2013 analysis. Note that in this abstract we quote 68% confidence limits on measured parameters and 95% upper limits on other parameters. We present the first results of polarization measurements with the Low Frequency Instrument at large angular scales. Combined with the Planck temperature and lensing data, these measurements give a reionization optical depth of τ = 0.066 ± 0.016, corresponding to a reionization redshift of z re= 8.8more » $$+1.7\\atop{-1.4}$$. These results are consistent with those from WMAP polarization measurements cleaned for dust emission using 353-GHz polarization maps from the High Frequency Instrument. We find no evidence for any departure from base ΛCDM in the neutrino sector of the theory; for example, combining Planck observations with other astrophysical data we find N eff = 3.15 ± 0.23 for the effective number of relativistic degrees of freedom, consistent with the value N eff = 3.046 of the Standard Model of particle physics. The sum of neutrino masses is constrained to Σ m ν < 0.23 eV. The spatial curvature of our Universe is found to be very close to zero, with | Ω K | < 0.005. Adding a tensor component as a

Planck 2015 results: XIII. Cosmological parameters

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ade, P. A. R.; Aghanim, N.; Arnaud, M.

Here, this paper presents cosmological results based on full-mission Planck observations of temperature and polarization anisotropies of the cosmic microwave background (CMB) radiation. Our results are in very good agreement with the 2013 analysis of the Planck nominal-mission temperature data, but with increased precision. The temperature and polarization power spectra are consistent with the standard spatially-flat 6-parameter ΛCDM cosmology with a power-law spectrum of adiabatic scalar perturbations (denoted “base ΛCDM” in this paper). From the Planck temperature data combined with Planck lensing, for this cosmology we find a Hubble constant, H 0 = (67.8 ± 0.9) km s -1Mpc -1, a matter density parameter Ω m = 0.308 ± 0.012, and a tilted scalar spectral index with n s = 0.968 ± 0.006, consistent with the 2013 analysis. Note that in this abstract we quote 68% confidence limits on measured parameters and 95% upper limits on other parameters. We present the first results of polarization measurements with the Low Frequency Instrument at large angular scales. Combined with the Planck temperature and lensing data, these measurements give a reionization optical depth of τ = 0.066 ± 0.016, corresponding to a reionization redshift of z re= 8.8more » $$+1.7\\atop{-1.4}$$. These results are consistent with those from WMAP polarization measurements cleaned for dust emission using 353-GHz polarization maps from the High Frequency Instrument. We find no evidence for any departure from base ΛCDM in the neutrino sector of the theory; for example, combining Planck observations with other astrophysical data we find N eff = 3.15 ± 0.23 for the effective number of relativistic degrees of freedom, consistent with the value N eff = 3.046 of the Standard Model of particle physics. The sum of neutrino masses is constrained to Σ m ν < 0.23 eV. The spatial curvature of our Universe is found to be very close to zero, with | Ω K | < 0.005. Adding a tensor component as a

Anisotropic Weyl symmetry and cosmology

NASA Astrophysics Data System (ADS)

Moon, Taeyoon; Oh, Phillial; Sohn, Jongsu

2010-11-01

We construct an anisotropic Weyl invariant theory in the ADM formalism and discuss its cosmological consequences. It extends the original anisotropic Weyl invariance of Hořava-Lifshitz gravity using an extra scalar field. The action is invariant under the anisotropic transformations of the space and time metric components with an arbitrary value of the critical exponent z. One of the interesting features is that the cosmological constant term maintains the anisotropic symmetry for z = -3. We also include the cosmological fluid and show that it can preserve the anisotropic Weyl invariance if the equation of state satisfies P = zρ/3. Then, we study cosmology of the Einstein-Hilbert-anisotropic Weyl (EHaW) action including the cosmological fluid, both with or without anisotropic Weyl invariance. The correlation of the critical exponent z and the equation of state parameter bar omega provides a new perspective of the cosmology. It is also shown that the EHaW action admits a late time accelerating universe for an arbitrary value of z when the anisotropic conformal invariance is broken, and the anisotropic conformal scalar field is interpreted as a possible source of dark energy.