Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reproducibility”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Dust Radiative Transfer Modeling of the Infrared Ring around the Magnetar SGR 1900+14

A peculiar infrared ring-like structure was discovered by Spitzer around the strongly magnetized neutron star SGR 1900+14. This infrared (IR) structure was suggested to be due to a dust-free cavity, produced by the Soft Gamma-ray Repeaters (SGRs) Giant Flare occurring in 1998, and kept illuminated by surrounding stars. Using a 3D dust radiative transfer code, we aimed to reproduce the emission morphology and the integrated emission flux of this structure assuming different spatial distributions and densities for the dust, and different positions for the illuminating stars. We found that a dust-free ellipsoidal cavity can reproduce the shape, flux, and spectrum of the ring-like IR emission, provided that the illuminating stars are inside the cavity and that the interstellar medium has high gas density ( n {sub H} ∼ 1000 cm{sup −3}). We further constrain the emitting region to have a sharp inner boundary and to be significantly extended in the radial direction, possibly even just a cavity in a smooth molecular cloud. We discuss possible scenarios for the formation of the dustless cavity and the particular geometry that allows it to be IR-bright.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Constraints on the Optical Depth of Galaxy Groups and Clusters

Future data from galaxy redshift surveys, combined with high-resolutions maps of the cosmic microwave background, will enable measurements of the pairwise kinematic Sunyaev–Zel’dovich (kSZ) signal with unprecedented statistical significance. This signal probes the matter-velocity correlation function, scaled by the average optical depth (τ) of the galaxy groups and clusters in the sample, and is thus of fundamental importance for cosmology. However, in order to translate pairwise kSZ measurements into cosmological constraints, external constraints on τ are necessary. In this work, we present a new model for the intracluster medium, which takes into account star formation, feedback, non-thermal pressure, and gas cooling. Our semi-analytic model is computationally efficient and can reproduce results of recent hydrodynamical simulations of galaxy cluster formation. We calibrate the free parameters in the model using recent X-ray measurements of gas density profiles of clusters, and gas masses of groups and clusters. Our observationally calibrated model predicts the average τ{sub 500} (i.e., the integrated τ within a disk of size R {sub 500}) to better than 6% modeling uncertainty (at 95% confidence level). If the remaining uncertainties associated with other astrophysical uncertainties and X-ray selection effects can be better understood, our model for the optical depth should break the degeneracy between optical depth and cluster velocity in the analysis of future pairwise kSZ measurements and improve cosmological constraints with the combination of upcoming galaxy and CMB surveys, including the nature of dark energy, modified gravity, and neutrino mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Light curve and SED modeling of the gamma-ray binary 1FGL J1018.6–5856: Constraints on the orbital geometry and relativistic flow

We present broadband spectral energy distributions and light curves of the gamma-ray binary 1FGL J1018.6–5856 measured in the X-ray and the gamma-ray bands. We find that the orbital modulation in the low-energy gamma-ray band is similar to that in the X-ray band, suggesting a common spectral component. However, above a GeV the orbital light curve changes significantly. We suggest that the GeV band contains significant flux from a pulsar magnetosphere, while the X-ray to TeV light curves are dominated by synchrotron and Compton emission from an intrabinary shock (IBS). We find that a simple one-zone model is inadequate to explain the IBS emission, but that beamed Synchrotron-self Compton radiation from adiabatically accelerated plasma in the shocked pulsar wind can reproduce the complex multiband light curves, including the variable X-ray spike coincident with the gamma-ray maximum. Furthermore, the model requires an inclination of ~50° and an orbital eccentricity of ~0.35, consistent with the limited constraints from existing optical observations. This picture motivates searches for pulsations from the energetic young pulsar powering the wind shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Deflection of the Cartwheel CME: ForeCAT Results

We analyze the Cartwheel coronal mass ejection's (CME; 2008 April 9) trajectory in the low corona with the ForeCAT model. This complex event presented a significant rotation in the low corona and a reversal in its original latitude direction. We successfully reproduce the observed CME’s trajectory (latitude and longitude deflection) and speed. Through a χ{sup 2} test, we are able to constrain the CME’s mass to (2.3−3.0) × 10{sup 14} g and the CME’s initial shape. We are able to constrain the expansion of the CME as well: the angular width linearly increases until 2.1 R{sub ⊙}, and is constant afterward. In order to match the observed latitude, we include a non-radial initial speed of −42 km s{sup −1}. Despite allowing the CME to rotate in the model, the magnetic forces of the solar background are not able to reproduce the observed rotation. We suggest that the complex reversal in latitude and the significant rotation of the Cartwheel CME can be justified with an asymmetrical reconnection event that ejected the CME non-radially and also initiated its rotation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Nebular Continuum and Line Emission in Stellar Population Synthesis Models

Accounting for nebular emission when modeling galaxy spectral energy distributions (SEDs) is important, as both line and continuum emissions can contribute significantly to the total observed flux. In this work, we present a new nebular emission model integrated within the Flexible Stellar Population Synthesis code that computes the line and continuum emission for complex stellar populations using the photoionization code Cloudy. The self-consistent coupling of the nebular emission to the matched ionizing spectrum produces emission line intensities that correctly scale with the stellar population as a function of age and metallicity. This more complete model of galaxy SEDs will improve estimates of global gas properties derived with diagnostic diagrams, star formation rates based on H α , and physical properties derived from broadband photometry. Our models agree well with results from other photoionization models and are able to reproduce observed emission from H ii regions and star-forming galaxies. Our models show improved agreement with the observed H ii regions in the Ne iii/O ii plane and show satisfactory agreement with He ii emission from z = 2 galaxies, when including rotating stellar models. Models including post-asymptotic giant branch stars are able to reproduce line ratios consistent with low-ionization emission regions. The models are integrated into current versions of FSPS and include self-consistent nebular emission predictions for MIST and Padova+Geneva evolutionary tracks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Advection of Potential Temperature in the Atmosphere of Irradiated Exoplanets: A Robust Mechanism to Explain Radius Inflation

The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a two-dimensional steady-state atmospheric circulation model, the validity of which is assessed by comparison to three-dimensional calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occurring in the Earth’s atmosphere or oceans. In the deep atmosphere, the vanishing heating flux forces the atmospheric structure to converge to a hotter adiabat than the one obtained with 1D calculations, implying a larger radius for the planet. Not only do the calculations reproduce the observed radius of HD 209458b, but also reproduce the observed correlation between radius inflation and irradiation for transiting planets. Vertical advection of potential temperature induced by non-uniform atmospheric heating thus provides a robust mechanism to explain the inflated radii of irradiated hot Jupiters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Y-type Brown Dwarfs: Estimates of Mass and Age from New Astrometry, Homogenized Photometry, and Near-infrared Spectroscopy

The survey of the mid-infrared sky by the Wide-field Infrared Survey Explorer ( WISE ) led to the discovery of extremely cold, low-mass brown dwarfs, classified as Y dwarfs, which extend the T class to lower temperatures. Twenty-four Y dwarfs are known at the time of writing. Here we present improved parallaxes for four of these, determined using Spitzer images. We give new photometry for four late-type T and three Y dwarfs and new spectra of three Y dwarfs, obtained at Gemini Observatory. We also present previously unpublished photometry taken from HST , ESO, Spitzer , and WISE archives of 11 late-type T and 9 Y dwarfs. The near-infrared data are put onto the same photometric system, forming a homogeneous data set for the coolest brown dwarfs. We compare recent models to our photometric and spectroscopic data set. We confirm that nonequilibrium atmospheric chemistry is important for these objects. Nonequilibrium cloud-free models reproduce well the near-infrared spectra and mid-infrared photometry for the warmer Y dwarfs with 425 ≤ T {sub eff} (K) ≤ 450. A small amount of cloud cover may improve the model fits in the near-infrared for the Y dwarfs with 325 ≤ T {sub eff} (K) ≤ 375. Neither cloudy nor cloud-free models reproduce the near-infrared photometry for the T {sub eff} = 250 K Y dwarf W0855. We use the mid-infrared region, where most of the flux originates, to constrain our models of W0855. We find that W0855 likely has a mass of 1.5–8 Jupiter masses and an age of 0.3–6 Gyr. The Y dwarfs with measured parallaxes are within 20 pc of the Sun and have tangential velocities typical of the thin disk. The metallicities and ages we derive for the sample are generally solar-like. We estimate that the known Y dwarfs are 3 to 20 Jupiter-mass objects with ages of 0.6–8.5 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Intrinsic Far-infrared Continua of Type-1 Quasars

The range of currently proposed active galactic nucleus (AGN) far-infrared templates results in uncertainties in retrieving host galaxy information from infrared observations and also undermines constraints on the outer part of the AGN torus. We discuss how to test and reconcile these templates. Physically, the fraction of the intrinsic AGN IR-processed luminosity compared with that from the central engine should be consistent with the dust-covering factor. In addition, besides reproducing the composite spectral energy distributions (SEDs) of quasars, a correct AGN IR template combined with an accurate library of star-forming galaxy templates should be able to reproduce the IR properties of the host galaxies, such as the luminosity-dependent SED shapes and aromatic feature strengths. We develop tests based on these expected behaviors and find that the shape of the AGN intrinsic far-IR emission drops off rapidly starting at ∼20 μ m and can be matched by an Elvis et al.-like template with a minor modification. Despite the variations in the near- to mid-IR bands, AGNs in quasars and Seyfert galaxies have remarkably similar intrinsic far-IR SEDs at λ ∼ 20–100 μ m, suggesting a similar emission character of the outermost region of the circumnuclear torus. The variations of the intrinsic AGN IR SEDs among the type-1 quasar population can be explained by the changing relative strengths of four major dust components with similar characteristic temperatures, and there is evidence for compact AGN-heated dusty structures at sub-kiloparsec scales in the far-IR.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Thermonuclear Bursts with Short Recurrence Times from Neutron Stars Explained by Opacity-driven Convection

Thermonuclear flashes of hydrogen and helium accreted onto neutron stars produce the frequently observed Type I X-ray bursts. It is the current paradigm that almost all material burns in a burst, after which it takes hours to accumulate fresh fuel for the next burst. In rare cases, however, bursts are observed with recurrence times as short as minutes. We present the first one-dimensional multi-zone simulations that reproduce this phenomenon. Bursts that ignite in a relatively hot neutron star envelope leave a substantial fraction of the fuel unburned at shallow depths. In the wake of the burst, convective mixing events driven by opacity bring this fuel down to the ignition depth on the observed timescale of minutes. There, unburned hydrogen mixes with the metal-rich ashes, igniting to produce a subsequent burst. We find burst pairs and triplets, similar to the observed instances. Our simulations reproduce the observed fraction of bursts with short waiting times of ∼30%, and demonstrate that short recurrence time bursts are typically less bright and of shorter duration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

PTF11kx: A Type Ia Supernova with Hydrogen Emission Persisting after 3.5 Years

The optical transient PTF11kx exhibited both the characteristic spectral features of Type Ia supernovae (SNe Ia) and the signature of ejecta interacting with circumstellar material (CSM) containing hydrogen, indicating the presence of a nondegenerate companion. In this paper, we present an optical spectrum at 1342 days after peak from Keck Observatory, in which the broad component of Hα emission persists with a similar profile as in early-time observations. We also present Spitzer IRAC detections obtained 1237 and 1818 days after peak, and an upper limit from Hubble Space Telescope ultraviolet imaging at 2133 days. We interpret our late-time observations in the context of published results—and reinterpret the early-time observations—in order to constrain the CSM's physical parameters and to compare to theoretical predictions for recurrent-nova systems. We find that the CSM's radial extent may be several times the distance between the star and the CSM's inner edge, and that the CSM column density may be two orders of magnitude lower than previous estimates. We show that the Hα luminosity decline is similar to other SNe with CSM interaction and demonstrate how our infrared photometry is evidence for newly formed, collisionally heated dust. We create a model for PTF11kx's late-time CSM interaction and find that X-ray reprocessing by photoionization and recombination cannot reproduce the observed Hα luminosity, suggesting that the X-rays are thermalized and that Hα radiates from collisional excitation. In conclusion, we discuss the implications of our results regarding the progenitor scenario and the geometric properties of the CSM for the PTF11kx system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Discovery of Extended Main-sequence Turnoffs in Four Young Massive Clusters in the Magellanic Clouds

An increasing number of young massive clusters (YMCs) in the Magellanic Clouds have been found to exhibit bimodal or extended main sequences (MSs) in their color–magnitude diagrams (CMDs). These features are usually interpreted in terms of a coeval stellar population with different stellar rotational rates, where the blue and red MS stars are populated by non- (or slowly) and rapidly rotating stellar populations, respectively. However, some studies have shown that an age spread of several million years is required to reproduce the observed wide turnoff regions in some YMCs. Here we present the ultraviolet–visual CMDs of four Large and Small Magellanic Cloud YMCs, NGC 330, NGC 1805, NGC 1818, and NGC 2164, based on high-precision Hubble Space Telescope photometry. We show that they all exhibit extended main-sequence turnoffs (MSTOs). The importance of age spreads and stellar rotation in reproducing the observations is investigated. The observed extended MSTOs cannot be explained by stellar rotation alone. Adopting an age spread of 35–50 Myr can alleviate this difficulty. We conclude that stars in these clusters are characterized by ranges in both their ages and rotation properties, but the origin of the age spread in these clusters remains unknown.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Trident: A Universal Tool for Generating Synthetic Absorption Spectra from Astrophysical Simulations

Hydrodynamical simulations are increasingly able to accurately model physical systems on stellar, galactic, and cosmological scales; however, the utility of these simulations is often limited by our ability to directly compare them with the data sets produced by observers: spectra, photometry, etc. To address this problem, we have created trident, a Python-based open-source tool for post-processing hydrodynamical simulations to produce synthetic absorption spectra and related data. trident can (i) create absorption-line spectra for any trajectory through a simulated data set mimicking both background quasar and down-the-barrel configurations; (ii) reproduce the spectral characteristics of common instruments like the Cosmic Origins Spectrograph; (iii) operate across the ultraviolet, optical, and infrared using customizable absorption-line lists; (iv) trace simulated physical structures directly to spectral features; (v) approximate the presence of ion species absent from the simulation outputs; (vi) generate column density maps for any ion; and (vii) provide support for all major astrophysical hydrodynamical codes. trident was originally developed to aid in the interpretation of observations of the circumgalactic medium and intergalactic medium, but it remains a general tool applicable in other contexts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Physical Origin of Long Gas Depletion Times in Galaxies

We present a model that explains why galaxies form stars on a timescale significantly longer than the timescales of processes governing the evolution of interstellar gas. We show that gas evolves from a non-star-forming to a star-forming state on a relatively short timescale, and thus the rate of this evolution does not limit the star formation rate (SFR). Instead, the SFR is limited because only a small fraction of star-forming gas is converted into stars before star-forming regions are dispersed by feedback and dynamical processes. Thus, gas cycles into and out of a star-forming state multiple times, which results in a long timescale on which galaxies convert gas into stars. Our model does not rely on the assumption of equilibrium and can be used to interpret trends of depletion times with the properties of observed galaxies and the parameters of star formation and feedback recipes in simulations. In particular, the model explains how feedback self-regulates the SFR in simulations and makes it insensitive to the local star formation efficiency. We illustrate our model using the results of an isolated L {sub *}-sized galaxy simulation that reproduces the observed Kennicutt–Schmidt relation for both molecular and atomic gas. Interestingly, the relation for molecular gas is almost linear on kiloparsec scales, although a nonlinear relation is adopted in simulation cells. We discuss how a linear relation emerges from non-self-similar scaling of the gas density PDF with the average gas surface density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

High-energy Gamma Rays from the Milky Way: Three-dimensional Spatial Models for the Cosmic-Ray and Radiation Field Densities in the Interstellar Medium

High-energy γ -rays of interstellar origin are produced by the interaction of cosmic-ray (CR) particles with the diffuse gas and radiation fields in the Galaxy. The main features of this emission are well understood and are reproduced by existing CR propagation models employing 2D galactocentric cylindrically symmetrical geometry. However, the high-quality data from instruments like the Fermi Large Area Telescope reveal significant deviations from the model predictions on few to tens of degrees scales, indicating the need to include the details of the Galactic spiral structure and thus requiring 3D spatial modeling. In this paper, the high-energy interstellar emissions from the Galaxy are calculated using the new release of the GALPROP code employing 3D spatial models for the CR source and interstellar radiation field (ISRF) densities. Three models for the spatial distribution of CR sources are used that are differentiated by their relative proportion of input luminosity attributed to the smooth disk or spiral arms. Two ISRF models are developed based on stellar and dust spatial density distributions taken from the literature that reproduce local near- to far-infrared observations. The interstellar emission models that include arms and bulges for the CR source and ISRF densities provide plausible physical interpretations for features found in the residual maps from high-energy γ -ray data analysis. The 3D models for CR and ISRF densities provide a more realistic basis that can be used for the interpretation of the nonthermal interstellar emissions from the Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Importance of Preventive Feedback: Inference from Observations of the Stellar Masses and Metallicities of Milky Way Dwarf Galaxies

Dwarf galaxies are known to have remarkably low star formation efficiency due to strong feedback. Adopting the dwarf galaxies of the Milky Way (MW) as a laboratory, we explore a flexible semi-analytic galaxy formation model to understand how the feedback processes shape the satellite galaxies of the MW. Using Markov Chain Monte Carlo, we exhaustively search a large parameter space of the model and rigorously show that the general wisdom of strong outflows as the primary feedback mechanism cannot simultaneously explain the stellar mass function and the mass–metallicity relation of the MW satellites. An extended model that assumes that a fraction of baryons is prevented from collapsing into low-mass halos in the first place can be accurately constrained to simultaneously reproduce those observations. The inference suggests that two different physical mechanisms are needed to explain the two different data sets. In particular, moderate outflows with weak halo mass dependence are needed to explain the mass–metallicity relation, and prevention of baryons falling into shallow gravitational potentials of low-mass halos (e.g., “pre-heating”) is needed to explain the low stellar mass fraction for a given subhalo mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Sensitivity to Thermonuclear Reaction Rates in Modeling the Abundance Anomalies of NGC 2419

Abundance anomalies in globular clusters provide strong evidence for multiple stellar populations within each cluster. These populations are usually interpreted as distinct generations, with the currently observed second-generation stars having formed in part from the ejecta of massive, first-generation “polluter” stars, giving rise to the anomalous abundance patterns. The precise nature of the polluters and their enrichment mechanism are still unclear. Even so, the chemical abundances measured in second-generation stars within the globular cluster NGC 2419 provide insight into this puzzling process. Previous work used Monte Carlo nuclear reaction network calculations to constrain the temperature–density conditions that could reproduce the observed abundances, thereby placing robust limits on the origins of the polluter material. The effect of individual reaction rates on these conditions has not been studied, however. Thus, we perform an exhaustive sensitivity study on the nuclear physics input to determine which reactions have the greatest impact on these predictions. We find that the {sup 30}Si(p,γ){sup 31}P, {sup 37}Ar(p,γ){sup 38}K, {sup 38}Ar(p,γ){sup 39}K, and {sup 39}K(p,γ){sup 40}Ca reactions are all critical in determining the temperature–density conditions, and ultimately, the origins of the polluter material. We conclude with recommendations for future experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Chemical Abundance Structure of the Inner Milky Way: A Signature of “Upside-down” Disk Formation

We present a model for the [α/Fe]--[Fe/H] distribution of stars in the inner Galaxy, 3 kpc<R<5 kpc, measured as a function of vertical distance |z| from the midplane by Hayden et al. (H15). Motivated by an “upside-down” scenario for thick disk formation, in which the thickness of the star-forming gas layer contracts as the stellar mass of the disk grows, we combine one-zone chemical evolution with a simple prescription in which the scale-height of the stellar distribution drops linearly from z{sub h}=0.8 kpc to z{sub h}=0.2 kpc over a timescale t {sub c}, remaining constant thereafter. We assume a linear-exponential star formation history, M-dot {sub ∗}(t)∝te{sup −t/t{sub s}{sub f}}. With a star formation efficiency timescale τ{sub ∗}=M{sub g}(t)/ M-dot {sub ∗}(t)=2 Gyr, an outflow mass-loading factor η= M-dot {sub out}(t)/ M-dot {sub ∗}(t)=1.5, t{sub sf}=3 Gyr, and t{sub c}=2.5 Gyr, the model reproduces the observed locus of inner disk stars in [α/Fe]--[Fe/H] and the metallicity distribution functions (MDFs) measured by H15 at |z|=0--0.5 kpc, 0.5--1 kpc, and 1--2 kpc. Substantial changes to model parameters lead to disagreement with the H15 data; for example, models with t{sub c}=1 Gyr or t{sub sf}=1 Gyr fail to match the observed MDF at high-|z| and low-|z|, respectively. The inferred scale-height evolution, with z {sub h}(t) dropping on a timescale t{sub c}∼t{sub sf} at large lookback times, favors upside-down formation over dynamical heating of an initially thin stellar population as the primary mechanism regulating disk thickness. The failure of our short-t {sub c} models suggests that any model in which thick disk formation is a discrete event will not reproduce the continuous dependence of the MDF on |z| found by H15. Our scenario for the evolution of the inner disk can be tested by future measurements of the |z|-distribution and the age–metallicity distribution at R=3--5 kpc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Particle Acceleration and Fractional Transport in Turbulent Reconnection

We consider a large-scale environment of turbulent reconnection that is fragmented into a number of randomly distributed unstable current sheets (UCSs), and we statistically analyze the acceleration of particles within this environment. We address two important cases of acceleration mechanisms when particles interact with the UCS: (a) electric field acceleration and (b) acceleration by reflection at contracting islands. Electrons and ions are accelerated very efficiently, attaining an energy distribution of power-law shape with an index 1–2, depending on the acceleration mechanism. The transport coefficients in energy space are estimated from test-particle simulation data, and we show that the classical Fokker–Planck (FP) equation fails to reproduce the simulation results when the transport coefficients are inserted into it and it is solved numerically. The cause for this failure is that the particles perform Levy flights in energy space, while the distributions of the energy increments exhibit power-law tails. We then use the fractional transport equation (FTE) derived by Isliker et al., whose parameters and the order of the fractional derivatives are inferred from the simulation data, and solving the FTE numerically, we show that the FTE successfully reproduces the kinetic energy distribution of the test particles. We discuss in detail the analysis of the simulation data and the criteria that allow one to judge the appropriateness of either an FTE or a classical FP equation as a transport model.

79 ASTRONOMY AND ASTROPHYSICS↗