Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COSMOLOGICAL MODELS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,747 records · Page 97

Cosmic Dawn (CoDa): the first radiation-hydrodynamics simulation of reionization and galaxy formation in the Local Universe

Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for cold dark matter. Reionization modelling requires simulating volumes large enough [~ (100 Mpc) 3 ] to sample reionization ‘patchiness’, while resolving millions of galaxy sources above ~10 8 M ⊙ combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the Local Universe today. Cosmic Dawn (‘CoDa’) is the first such fully coupled, radiation-hydrodynamics simulation of reionization of the Local Universe. Our new hybrid CPU-GPU code, RAMSES-CUDATON, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc) 3 with 4096 3 particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 × 109 M ⊙ , decreasing the abundance of faint galaxies around M AB1600 = [-10, -12]. For most of reionization, star formation was dominated by haloes between 10 10 –10 11 M ⊙ , so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. Intergalactic filaments display sheathed structures, with hot envelopes surrounding cooler cores, but do not self-shield, unlike regions denser than 100 $\langle$ρ$\rangle$.

Astronomy & Astrophysics↗

Milky Way Satellite Census. III. Constraints on Dark Matter Properties from Observations of Milky Way Satellite Galaxies

We perform a comprehensive study of Milky Way (MW) satellite galaxies to constrain the fundamental properties of dark matter (DM). This analysis fully incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and marginalizes over uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW disk. Our results are consistent with the cold, collisionless DM paradigm and yield the strongest cosmological constraints to date on particle models of warm, interacting, and fuzzy dark matter. At $95\%$ confidence, we report limits on (i) the mass of thermal relic warm DM, $m_{\rm WDM} > 6.5\ \mathrm{keV}$ (free-streaming length, $\lambda_{\rm{fs}} \lesssim 10\,h^{-1}\ \mathrm{kpc}$), (ii) the velocity-independent DM-proton scattering cross section, $\sigma_{0} < 8.8\times 10^{-29}\ \mathrm{cm}^{2}$ for a $100\ \mathrm{MeV}$ DM particle mass (DM-proton coupling, $c_p \lesssim (0.3\ \mathrm{GeV})^{-2}$), and (iii) the mass of fuzzy DM, $m_{\phi}> 2.9 \times 10^{-21}\ \mathrm{eV}$ (de Broglie wavelength, $\lambda_{\rm{dB}} \lesssim 0.5\ \mathrm{kpc}$). These constraints are complementary to other observational and laboratory constraints on DM properties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Toward Transparent Data Management in Multi-layer Storage Hierarchy for HPC Systems

Upcoming exascale high performance computing (HPC) systems are expected to comprise multi-tier storage hierarchy, and thus will necessitate innovative storage and I/O mechanisms. Traditional disk and block-based interfaces and file systems face severe challenges in utilizing capabilities of storage hierarchies due to the lack of hierarchy support and semantic interfaces. Object-based and semantically-rich data abstractions for scientific data management on large scale systems offer a sustainable solution to these challenges. Such data abstractions can also simplify users involvement in data movement. Here, we take the first steps of realizing such an object abstraction and explore storage mechanisms for these objects to enhance I/O performance, especially for scientific applications. We explore how an object-based interface can facilitate next generation scalable computing systems by presenting the mapping of data I/O from two real world HPC scientific use cases: a plasma physics simulation code (VPIC) and a cosmology simulation code (HACC). Our storage model stores data objects in different physical organizations to support data movement across layers of memory/storage hierarchy. Our implementation sclaes well to 16K parallel processes, and compared to the state of the art, such as MPI-IO and HDF5, our object-based data abstractions and data placement strategy in multi-level storage hierarchy achieves up to 7 X I/O performance improvement for scientific data.

97 MATHEMATICS AND COMPUTING↗

SMART U(1)$$_X$$: standard model with axion, right handed neutrinos, two Higgs doublets and U(1)$$_X$$ gauge symmetry

Abstract To address five fundamental shortcomings of the Standard Model (SM) of particle physics and cosmology, we propose a phenomenologically viable framework based on a $$U(1)_X \times U(1)_{PQ}$$ U ( 1 ) X × U ( 1 ) PQ extension of the SM, that we call “SMART U(1) $$_X$$ X ”. The $$U(1)_X$$ U ( 1 ) X gauge symmetry is a well-known generalization of the $$U(1)_{B-L}$$ U ( 1 ) B - L symmetry and $$U(1)_{PQ}$$ U ( 1 ) PQ is the global Peccei–Quinn (PQ) symmetry. Three right handed neutrinos are added to cancel $$U(1)_X$$ U ( 1 ) X related anomalies, and they play a crucial role in understanding the observed neutrino oscillations and explaining the observed baryon asymmetry in the universe via leptogenesis. Implementation of PQ symmetry helps resolve the strong CP problem and also provides axion as a compelling dark matter (DM) candidate. The $$U(1)_X$$ U ( 1 ) X gauge symmetry enables us to implement the inflection-point inflation scenario with $$H_{inf} \lesssim 2 \times 10^{7}$$ H inf ≲ 2 × 10 7 GeV, where $$H_{inf}$$ H inf is the value of Hubble parameter during inflation. This is crucial to overcome a potential axion domain wall problem as well as the axion isocurvature problem. The SMART U(1) $$_X$$ X framework can be successfully implemented in the presence of SU (5) grand unification, as we briefly show.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

AMReX: Block-structured adaptive mesh refinement for multiphysics applications

Block-structured adaptive mesh refinement (AMR) provides the basis for the temporal and spatial discretization strategy for a number of Exascale Computing Project applications in the areas of accelerator design, additive manufacturing, astrophysics, combustion, cosmology, multiphase flow, and wind plant modeling. AMReX is a software framework that provides a unified infrastructure with the functionality needed for these and other AMR applications to be able to effectively and efficiently utilize machines from laptops to exascale architectures. AMR reduces the computational cost and memory footprint compared to a uniform mesh while preserving accurate descriptions of different physical processes in complex multiphysics algorithms. AMReX supports algorithms that solve systems of partial differential equations in simple or complex geometries and those that use particles and/or particle–mesh operations to represent component physical processes. In this article, we will discuss the core elements of the AMReX framework such as data containers and iterators as well as several specialized operations to meet the needs of the application projects. In addition, we will highlight the strategy that the AMReX team is pursuing to achieve highly performant code across a range of accelerator-based architectures for a variety of different applications.

Zhang, Weiqun↗

Isothermal Bondi Accretion in Jaffe and Hernquist Galaxies with a Central Black Hole: Fully Analytical Solutions

One of the most active fields of research of modern-day astrophysics is that of massive black hole formation and coevolution with the host galaxy. In these investigations, ranging from cosmological simulations, to semi-analytical modeling, to observational studies, the Bondi solution for accretion on a central point-mass is widely adopted. In this work we generalize the classical Bondi accretion theory to take into account the effects of the gravitational potential of the host galaxy, and of radiation pressure in the optically thin limit. Then, we present the fully analytical solution, in terms of the Lambert–Euler W -function, for isothermal accretion in Jaffe and Hernquist galaxies with a central black hole. The flow structure is found to be sensitive to the shape of the mass profile of the host galaxy. These results and the formulae that are provided, most importantly, the one for the critical accretion parameter, allow for a direct evaluation of all flow properties, and are then useful for the abovementioned studies. As an application, we examine the departure from the true mass accretion rate of estimates obtained using the gas properties at various distances from the black hole, under the hypothesis of classical Bondi accretion. An overestimate is obtained from regions close to the black hole, and an underestimate outside a few Bondi radii; the exact position of the transition between the two kinds of departure depends on the galaxy model.

79 ASTRONOMY AND ASTROPHYSICS↗

The formation and evolution of domain walls

Domain walls are sheet-like defects produced when the low energy vacuum has isolated degenerate minima. The researchers' computer code follows the evolution of a scalar field, whose dynamics are determined by its Lagrangian density. The topology of the scalar field determines the evolution of the domain walls. This approach treats both wall dynamics and reconnection. The researchers investigated not only potentials that produce single domain walls, but also potentials that produce a network of walls and strings. These networks arise in axion models where the U(1) Peccei-Quinn symmetry is broken into Z sub N discrete symmetries. If N equals 1, the walls are bounded by strings and the network quickly disappears. For N greater than 1, the network of walls and strings behaved qualitatively just as the wall network shown in the figures given here. This both confirms the researchers' pessimistic view that domain walls cannot play an important role in the formation of large scale structure and implies that axion models with multiple minimum can be cosmologically disastrous.

Press, William H.↗

Cold dark matter and degree-scale cosmic microwave background anisotropy statistics after COBE

We conduct a Monte Carlo simulation of the cosmic microwave background (CMB) anisotropy in the UCSB South Pole 1991 degree-scale experiment. We examine cold dark matter cosmology with large-scale structure seeded by the Harrison-Zel'dovich hierarchy of Gaussian-distributed primordial inhomogeneities normalized to the COBE-DMR measurement of large-angle CMB anisotropy. We find it statistically implausible (in the sense of low cumulative probability F lower than 5 percent, of not measuring a cosmological delta-T/T signal) that the degree-scale cosmological CMB anisotropy predicted in such models could have escaped a detection at the level of sensitivity achieved in the South Pole 1991 experiment.

Gorski, Krzysztof M.↗

Galaxy formation in Lambda greater than 0 Friedmann models: Consequences for the number counts versus redshift test

We study the effect of the cosmological constant Lambda on galaxy formation using a simple spherical top-hat overdensity model. We consider models with Omega(sub 0) = 0.2, lambda(sub 0) = 0, and Omega(sub 0) = 0.2, lambda(sub 0) = 0.8 (where Omega(sub 0) is the density parameter, and lambda(sub 0) identically equal Lambda/3 H(sub 0 exp 2) where H(sub 0) is the Hubble constant). We adjust the initial power spectrum amplitude so that both models reproduce the same large-scale structures. The galaxy formation era in the lambda(sub 0) = 0 model occurs early (z approximately 6) and is very short, whereas in the lambda(sub 0) = 0.8 model the galaxy formation era starts later (z approximately 4), and last much longer, possibly all the way to the present. Consequently, galaxies at low redshift (z less than 1) are significantly more evolved in the lambda(sub 0) = 0 model than in the lambda(sub 0) = 0.8 model. This result implies that previous attempts to determine Lambda using the number counts versus redshift test are probably unreliable.

Martel, Hugo↗

Petascale Computing: Impact on Future NASA Missions

This slide presentation reviews NASA's use of a new super computer, called Columbia, capable of operating at 62 Tera Flops. This computer is the 4th fastest computer in the world. This computer will serve all mission directorates. The applications that it would serve are: aerospace analysis and design, propulsion subsystem analysis, climate modeling, hurricane prediction and astrophysics and cosmology.

Brooks, Walter↗

Using Black Hole Mergers to Explore Structure Formation

Observations of gravitational waves from massive black hole mergers will open a new window into the era of structure formation in the early universe. Past efforts have concentrated on calculating merger rates using different physical assumptions, resulting in merger rate estimates that span a wide range (0.1 - 1 0A4 mergers/year). We develop a semi-analytical, phenomenological model of massive black hole mergers that includes plausible combinations of several physical parameters, which we then turn around to determine how well observations with the Laser Interferometer Space Antenna (LISA) will be able to enhance our understanding of the universe during the critical z approx. 5 - 30 epoch. Our approach involves generating synthetic LISA observable data (total BH masses, BH mass ratios, redshifts, merger rates), which are then analyzed using a Markov Chain Monte Carlo method, thus finding constraints for the physical parameters of the mergers. We find that our method works well at estimating merger parameters and that the number of merger events is a key discriminant among models, therefore making our method robust against observational uncertainties. Our approach can also be extended to more physically-driven models and more general problems in cosmology.

Alicea-Munoz, E.↗

Using Black Hole Mergers to Explore Structure Formation

Observations of gravitational waves from massive black hole mergers will open a new window into the era of structure formation in the early universe. Past efforts have concentrated on calculating merger rates using different physical assumptions, resulting in merger rate estimates that span a wide range (0.1 - 10(exp 4) mergers/year). We develop a semi-analytical, phenomenological model of massive black hole mergers that includes plausible combinations of several physical parameters, which we then turn around to determine how well observations with the Laser Interferometer Space Antenna (LISA) will be able to enhance our understanding of the universe during the critical z approximately equal to 5-30 epoch. Our approach involves generating synthetic LISA observable data (total BH masses, BH mass ratios, redshifts, merger rates), which are then analyzed using a Markov Chain Monte Carlo method, thus finding constraints for the physical parameters of the mergers. We find that our method works well at estimating merger parameters and that the number of merger events is a key discriminant among models, therefore making our method robust against observational uncertainties. Our approach can also be extended to more physically-driven models and more general problems in cosmology. This work is supported in part by the Cooperative Education Program at NASA/GSFC.

Alicea-Munoz, E.↗

OBSERVABLE DEVIATIONS FROM HOMOGENEITY IN AN INHOMOGENEOUS UNIVERSE

How does inhomogeneity affect our interpretation of cosmological observations? It has long been wondered to what extent the observable properties of an inhomogeneous universe differ from those of a corresponding Friedmann–Lemaître–Robertson–Walker (FLRW) model, and how the inhomogeneities affect that correspondence. Here, we use numerical relativity to study the behavior of light beams traversing an inhomogeneous universe, and construct the resulting Hubble diagrams. The universe that emerges exhibits an average FLRW behavior, but inhomogeneous structures contribute to deviations in observables across the observer’s sky. We also investigate the relationship between angular diameter distance and the angular extent of a source, finding deviations that grow with source redshift. These departures from FLRW are important path-dependent effects, with implications for using real observables in an inhomogeneous universe such as our own.

79 ASTRONOMY AND ASTROPHYSICS↗

Cosmological constraints on Lorentz violating dark energy

The role of Lorentz invariance as a fundamental symmetry of nature has been lately reconsidered in different approaches to quantum gravity. It is thus natural to study whether other puzzles of physics may be solved within these proposals. This may be the case for the cosmological constant problem. Indeed, it has been shown that breaking Lorentz invariance provides Lagrangians that can drive the current acceleration of the universe without experiencing large corrections from ultraviolet physics. In this work, we focus on the simplest model of this type, called ΘCDM, and study its cosmological implications in detail. At the background level, this model cannot be distinguished from ΛCDM. The differences appear at the level of perturbations. We show that in ΘCDM, the spectrum of CMB anisotropies and matter fluctuations may be affected by a rescaling of the gravitational constant in the Poisson equation, by the presence of extra contributions to the anisotropic stress, and finally by the existence of extra clustering degrees of freedom. To explore these modifications accurately, we modify the Boltzmann code class. We then use the parameter inference code Monte Python to confront ΘCDM with data from WMAP-7, SPT and WiggleZ. We obtain strong bounds on the parameters accounting for deviations from ΛCDM. In particular, we find that the discrepancy between the gravitational constants appearing in the Poisson and Friedmann equations is constrained at the level of 1.8%.

79 ASTRONOMY AND ASTROPHYSICS↗