Search NASA⌕ Search

SEARCH · Search NASA

Results for “density model”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Plasma Properties in the Earth's Magnetosheath Near the Subsolar Magnetopause: Implications for Geocoronal Density Estimates

Combined in situ ion measurements and remote sensing of energetic neutral atoms are used to determine the geocoronal Hydrogen density at large (∼10 R E ) distances from the Earth. This method for determining the geocoronal density requires global magnetospheric modeling. Observations in the Earth's subsolar magnetosheath from the Magnetospheric Multiscale mission are used to determine the accuracy of using global models to predict the geocoronal density. On average, gas dynamic and magnetohydrodynamic (MHD) models and observations are in reasonable agreement, with differences <25%. In addition, the MHD model subsolar magnetopause is about 0.5 R E sunward of the observed location. However, variations around averages are large (up to a factor of 2), indicating that global models introduce relatively large uncertainties in geocoronal density estimates. Finally, the critical ion flux in the Interstellar Boundary Explorer IBEX‐Hi energy range is often minimally affected by fluctuations of a factor of 2 in the density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Atomistic Insights into Nucleation and Formation of Hexagonal Boron Nitride on Nickel from First-Principles-Based Reactive Molecular Dynamics Simulations

Atomistic-scale insights into the growth of a continuous, atomically thin hexagonal boron nitride (hBN) lattice from elemental boron and nitrogen on Ni substrates were obtained from multiscale modeling combining density functional theory (DFT) and reactive molecular dynamics. The quantum mechanical calculations focused on the adsorption and reaction energetics for the hBN building-block species, i.e., atomic B, N, BxNy (x, y = 1, 2), on Ni(111) and Ni(211), and the diffusion pathways of elemental B and N on these slab model surfaces and in the sublayer. B can diffuse competitively on both the surface and in the sublayer, while N diffuses strictly on the substrate surface. The DFT data were then used to generate a classical description of the Ni–B and Ni–N pair interactions within the formulation of the reactive force field, ReaxFF. Using the potential developed from this work, the elementary nucleation and growth process of an hBN monolayer structure from elemental B and N is shown at the atomistic scale. The nucleation initiates from the growth of linear BN chains, which evolve into branched and then hexagonal lattices. Subsequent DFT calculations confirmed the structure evolution energetically and validate the self-consistency of this multiscale modeling framework. On the basis of this framework, the fundamental aspects regarding crystal quality and the role of temperature and substrates used during hBN growth can also be understood.

Liu, Song↗
DOE OSTI2019

Bridging the Homogeneous-Heterogeneous Divide: Modeling Spin for Reactivity in Single Atom Catalysis

Single atom catalysts (SACs) are emergent catalytic materials that have the promise of merging the scalability of heterogeneous catalysts with the high activity and atom economy of homogeneous catalysts. Computational, first-principles modeling can provide essential insight into SAC mechanism and active site configuration, where the sub-nm-scale environment can challenge even the highest-resolution experimental spectroscopic techniques. Nevertheless, the very properties that make SACs attractive in catalysis, such as localized d electrons of the isolated transition metal center, make them challenging to study with conventional computational modeling using density functional theory (DFT). For example, Fe/N-doped graphitic SACs have exhibited spin-state dependent reactivity that remains poorly understood. However, spin-state ordering in DFT is very sensitive to the nature of the functional approximation chosen. In this work, we develop accurate benchmarks from correlated wavefunction theory (WFT) for relevant octahedral complexes. We use those benchmarks to evaluate optimal DFT functional choice for predicting spin state ordering in small octahedral complexes as well as models of pyridinic and pyrrolic nitrogen environments expected in larger SACs. Using these guidelines, we determine Fe/N-doped graphene SAC model properties and reactivity as well as their sensitivities to DFT functional choice. Finally, we conclude with broad recommendations for computational modeling of open-shell transition metal single-atom catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Enhanced Pair-Density-Wave Vertices in a Bilayer Hubbard Model at Half Filling

Motivated by the pair-density-wave (PDW) state found in the one-dimensional Kondo-Heisenberg chain, we report on a determinant quantum Monte Carlo study of pair fields for a two-dimensional half-filled Hubbard layer coupled to an itinerant, noninteracting layer with one electron per site. In a specific range of interlayer hopping, the pairing vertex associated with PDW order becomes more attractive than that for uniform d-wave pairing, although both remain subdominant to the leading antiferromagnetic correlations at half filling. In conclusion, our result sheds light on where one potentially may find a PDW state in such a model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Phonon anharmonicity in silicon from 100 to 1500 K

Inelastic neutron scattering was performed on silicon powder to measure the phonon density of states (DOS) from 100 to 1500 K. The mean fractional energy shifts with temperature of the modes were $\langle$Δε i /ε i ΔT$\rangle$=₋0.07, giving a mean isobaric Grüneisen parameter of +6.95±0.67, which is significantly different from the isothermal parameter of +0.98. These large effects are beyond the predictions from quasiharmonic models using density functional theory or experimental data, demonstrating large effects from phonon anharmonicity. At 1500 K the anharmonicity contributes 0.15k B /atom to the vibrational entropy, compared to 0.03k B /atom from quasiharmonicity. Lastly, excellent agreement was found between the entropy from phonon DOS measurements and the reference NIST-JANAF thermodynamic entropy from calorimetric measurements.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Protostellar formation in rotating interstellar clouds. IV Nonisothermal collapse

Radiative transfer in the Eddington approximation is included in a multidimensional, self-gravitational, hydrodynamical computer code. Details of the numerical solution and thermodynamic relations are given. Comparison calculations with previous spherically symmetrical models of protostellar collapse are used to validate the basic approach and the artifices which allow the explicit hydrodynamics code to follow the accretion of gas onto a quasi-equilibrium core. A series of axisymmetric models is used to investigate the importance of rotation in collapsing clouds, as the initial amount of angular momentum is lowered, with an emphasis on the possible formation of rings. Rings readily form even in the nonisothermal regime except for very low initial angular momenta; even these clouds may experience ring formation prior to reaching stellar densities. The models imply that other effects (such as gravitational torques or turbulent viscosity) may be necesary to avoid binary formation and thus result in a presolar nebula consistent with the assumptions of either Safronov or Cameron.

Boss, A. P.↗
DOE OSTI2015

Density-dependent electron transport and precise modeling of GaN high electron mobility transistors

We report on the direct measurement of two-dimensional sheet charge density dependence of electron transport in AlGaN/GaN high electron mobility transistors (HEMTs). Pulsed IV measurements established increasing electron velocities with decreasing sheet charge densities, resulting in saturation velocity of 1.9 × 10{sup 7 }cm/s at a low sheet charge density of 7.8 × 10{sup 11 }cm{sup −2}. An optical phonon emission-based electron velocity model for GaN is also presented. It accommodates stimulated longitudinal optical (LO) phonon emission which clamps the electron velocity with strong electron-phonon interaction and long LO phonon lifetime in GaN. A comparison with the measured density-dependent saturation velocity shows that it captures the dependence rather well. Finally, the experimental result is applied in TCAD-based device simulator to predict DC and small signal characteristics of a reported GaN HEMT. Good agreement between the simulated and reported experimental results validated the measurement presented in this report and established accurate modeling of GaN HEMTs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Hidden molecules in planetary nebulae: New detections of HCN and HCO{sup +} from a multi-object survey

Searches for HCN and HCO{sup +} have been conducted toward 17 planetary nebulae (PNs) in the age range 800 to 13,000 years using the facilities of the Arizona Radio Observatory (ARO). For both molecules, observations of the J=1→0 transition near 88–89 GHz were carried out with the ARO 12 m, including measurements with the new ALMA prototype antenna, while the J=3→2 lines near 265–267 GHz were sought with the ARO Sub-Millimeter Telescope (SMT). HCN and HCO{sup +} were newly detected in 13 of the 17 target sources in at least one transition. Nine PNs were common to both molecules: Hb5, K3-17, K3-58, M1-7, M4-14, M3-28, M3-55, NGC 2440, and K4-47, while HCO{sup +} was also identified in K3-83 and M2-9, and HCN in K3-45 and NGC 6772. From radiative transfer modeling, column densities for HCN and HCO{sup +} in these sources were determined to be N{sub tot}(HCN) ∼ 0.2–27 × 10{sup 13} cm{sup −2} and N{sub tot}(HCO{sup +}) ∼ 0.3–8.7 × 10{sup 13} cm{sup −2}. Gas densities, assumed to be in clumped regions, were established to be n(H{sub 2}) ∼ 0.1–5.2 × 10{sup 6} cm{sup −3}. Fractional abundances, relative to H{sub 2}, for both molecules were found to be f(HCN) ∼ 0.1–9.1 × 10{sup −7} and f(HCO{sup +}) ∼ 0.04–7.4 × 10{sup −7}. The abundances of both species were found to remain relatively constant with nebular age over a 10,000 year time span, in contrast to predictions of chemical models. The HCN/HCO{sup +} ratio varied from 17 to <0.2, and roughly correlates with the C/O ratio. Polyatomic molecules appear to be common constituents of PNs.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2002

The Velocity Distribution of Isolated Radio Pulsars

We infer the velocity distribution of radio pulsars based on large-scale 0.4 GHz pulsar surveys. We do so by modelling evolution of the locations, velocities, spins, and radio luminosities of pulsars; calculating pulsed flux according to a beaming model and random orientation angles of spin and beam; applying selection effects of pulsar surveys; and comparing model distributions of measurable pulsar properties with survey data using a likelihood function. The surveys analyzed have well-defined characteristics and cover approx. 95% of the sky. We maximize the likelihood in a 6-dimensional space of observables P, dot-P, DM, absolute value of b, mu, F (period, period derivative, dispersion measure, Galactic latitude, proper motion, and flux density). The models we test are described by 12 parameters that characterize a population's birth rate, luminosity, shutoff of radio emission, birth locations, and birth velocities. We infer that the radio beam luminosity (i) is comparable to the energy flux of relativistic particles in models for spin-driven magnetospheres, signifying that radio emission losses reach nearly 100% for the oldest pulsars; and (ii) scales approximately as E(exp 1/2) which, in magnetosphere models, is proportional to the voltage drop available for acceleration of particles. We find that a two-component velocity distribution with characteristic velocities of 90 km/ s and 500 km/ s is greatly preferred to any one-component distribution; this preference is largely immune to variations in other population parameters, such as the luminosity or distance scale, or the assumed spin-down law. We explore some consequences of the preferred birth velocity distribution: (1) roughly 50% of pulsars in the solar neighborhood will escape the Galaxy, while approx. 15% have velocities greater than 1000 km/ s (2) observational bias against high velocity pulsars is relatively unimportant for surveys that reach high Galactic absolute value of z distances, but is severe for spatially bounded surveys; (3) an important low-velocity population exists that increases the fraction of neutron stars retained by globular clusters and is consistent with the number of old objects that accrete from the interstellar medium; (4) under standard assumptions for supernova remnant expansion and pulsar spin-down, approx. 10% of pulsars younger than 20 kyr will appear to lie outside of their host remnants. Finally, we comment on the ramifications of our birth velocity distribution for binary survival and the population of inspiraling binary neutron stars relevant to some GRB models and potential sources for LIGO.

Arzoumanian, Z.↗
DOE OSTI2014

Class I methanol (CH{sub 3}OH) maser conditions near supernova remnants

We present results from calculations of the physical conditions necessary for the occurrence of 36.169 (4{sub –1}-3{sub 0} E), 44.070 (7{sub 0}-6{sub 1} A {sup +}), 84.521 (5{sub –1}-4{sub 0} E), and 95.169 (8{sub 0}-7{sub 1} A {sup +}) GHz methanol (CH{sub 3}OH) maser emission lines near supernova remnants (SNRs), using the MOLPOP-CEP program. The calculations show that given a sufficient methanol abundance, methanol maser emission arises over a wide range of densities and temperatures, with optimal conditions at n ∼ 10{sup 4}-10{sup 6} cm{sup –3} and T > 60 K. The 36 GHz and 44 GHz transitions display more significant maser optical depths compared to the 84 GHz and 95 GHz transitions over the majority of physical conditions. It is also shown that line ratios are an important and applicable probe of the gas conditions. The line ratio changes are largely a result of the E-type transitions becoming quenched faster at increasing densities. The modeling results are discussed using recent observations of CH{sub 3}OH and hydroxyl (OH) masers near the SNRs G1.4–0.1, W28, and Sgr A East.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Unusual Exchange Couplings and Intermediate Temperature Weyl State in Co 3 Sn 2 S 2

Understanding magnetism and its possible correlations to topological properties has emerged to the forefront as a difficult topic in studying magnetic Weyl semimetals. Co 3 Sn 2 S 2 is a newly discovered magnetic Weyl semimetal with a kagome lattice of cobalt ions and has triggered intense interest for rich fantastic phenomena. In this study, we report the magnetic exchange couplings of Co 3 Sn 2 S 2 using inelastic neutron scattering and two density functional theory (DFT) based methods: constrained magnetism and multiple-scattering Green’s function methods. Co 3 Sn 2 S 2 exhibits highly anisotropic magnon dispersions and linewidths below T C , and paramagnetic excitations above T C . The spin-wave spectra in the ferromagnetic ground state is well described by the dominant third-neighbor “across-hexagon” J d model. Our density functional theory calculations reveal that both the symmetry-allowed 120° antiferromagnetic orders support Weyl points in the intermediate temperature region, with distinct numbers and the locations of Weyl points. Our study highlights the important role Co 3 Sn 2 S 2 can play in advancing our understanding of kagome physics and exploring the interplay between magnetism and band topology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2018

Alternative Explanations for Extreme Supersolar Iron Abundances Inferred from the Energy Spectrum of Cygnus X-1

Here we study a 1-200 keV energy spectrum of the black hole binary Cygnus X-1 taken with NuSTAR and Suzaku. This is the first report of a NuSTAR observation of Cyg X-1 in the intermediate state, and the observation was taken during the part of the binary orbit where absorption due to the companion's stellar wind is minimal. The spectrum includes a multi-temperature thermal disk component, a cutoff power-law component, and relativistic and nonrelativistic reflection components. Our initial fits with publicly available constant density reflection models (relxill and reflionx) lead to extremely high iron abundances (greater than 9.96 and 10.6 minus 0.9 (sub) to plus 1.6 (sup) times solar, respectively). Although supersolar iron abundances have been reported previously for Cyg X-1, our measurements are much higher and such variability is almost certainly unphysical. Using a new version of reflionx that we modified to make the electron density a free parameter, we obtain better fits to the spectrum even with solar iron abundances. We report on how the higher density (n (sub e) equal to 3.98 minus 0.25 (sub) to plus 0.12 (sup) times 10 (sup 20) per cubic centimeter) impacts other parameters such as the inner radius and inclination of the disk.

Tomsick, John A.↗
DOE OSTI2014

Calculated and Measured Dose Distributions in an Anthropomorphic RANDO Phantom - Paper 105

Patient-specific quality assurance (QA) is a major component of the Intensity Modulated Radiation Therapy (IMRT) clinical process. The objective of this often time-intensive measurement process is to verify that the prescribed dose is accurately delivered to the target volume while sparing surrounding healthy tissue. The IMRT dose distributions can be measured using a wide variety of techniques and measurement devices including film (radiographic and radiochromic), ion chamber arrays, diode arrays, and Electronic Portal Imaging Devices (EPIDs). Of all the available modalities, film is one of the most versatile and most frequently utilized measurement systems. In rotation delivery techniques (such as volume modulated arc therapy and helical tomotherapy), composite dose distributions are obtained by irradiating a test phantom containing dosimeter(s). When film is utilized as the dosimeter, it can be placed in coronal, sagittal, or transverse planes inside the phantom. Because the film is irradiated using the rotational delivery sequence used to treat the patient, the dose given to specific anatomical structures (such as the parotids, spinal cord, bladder, rectum etc...) can be evaluated. Thomas et al (Med Phys 32: 3793-3800) recently reported the results of IMRT verifications for ten patients planned and treated with helical tomotherapy. This study reported that nearly 92% of the film pixels within the analyzed region of interest were within a gamma criterion of 2 mm and 2%. Like most tomotherapy users, Thomas et al used a test phantom made of homogeneous water equivalent material. Heterogeneous phantoms that simulate the different tissue densities in the human body do exist, but are infrequently used for patient-specific IMRT QA. The potential advantage of using heterogeneous phantoms is that the measured dose distributions will more adequately reflect the clinical dose distributions in the patients. The use of homogenous phantoms during IMRT quality assurance can potentially mask errors that arise from limitations in the treatment planning system's modeling of density corrections. The purpose of this study was to evaluate the dosimetric accuracy of a helical tomotherapy system using an anthropomorphic RANDO phantom (The Phantom Laboratory, Salem, NY) for patient-specific QA. The results of this study will allow a benchmark to be established for dosimetric accuracy in heterogeneous materials. Performing patient-specific QA is an integral part of the IMRT patient treatment process. Test phantoms are typically homogeneous, and do test the treatment planning systems ability to account for heterogeneity corrections. In this study, an anthropomorphic RANDO phantom was used to evaluate the doses calculated by a helical tomotherapy system. The dosimetric agreement between calculated and measured was worse than published results in homogenous phantoms, especially for target volumes located in the lung. Further investigation will be needed to determine if these errors are due to film response at heterogeneous tissue interfaces or from actual dose calculation errors. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

QCD constraints on the equation of state for compact stars

In recent years, there have been several successful attempts to constrain the equation of state of neutron star matter using input from low-energy nuclear physics and observational data. We demonstrate that significant further restrictions can be placed by additionally requiring the pressure to approach that of deconfined quark matter at high densities. Remarkably, the new constraints turn out to be highly insensitive to the amount - or even presence - of quark matter inside the stars. In this framework, we also present a simple effective equation of state for cold quark matter that consistently incorporates the effects of interactions and furthermore includes a built-in estimate of the inherent systematic uncertainties. This goes beyond the MIT bag model description in a crucial way, yet leads to an equation of state that is equally straightforward to use.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Toroidal modeling of plasma flow damping and density pump-out by RMP during ELM mitigation in HL-2A

Reduction of both the plasma density and toroidal flow speed, due to application of the predominantly n = 1 (n is the toroidal mode number) resonant magnetic perturbation (RMP) for controlling the edge localized mode in the HL-2A tokamak, is numerically investigated utilizing the quasi-linear initial-value code MARS-Q (Liu et al 2013 Phys. Plasmas 20 042503). Simulation results reveal that the neoclassical toroidal viscosity (NTV) due to three dimensional fields plays the key role in modifying the plasma momentum and particle transport in the HL-2A discharge. By comparing the modeling results with the measured density pump-out in the experiment, the electron NTV particle flux model, in combination with the free-boundary condition for the axisymmetric change of the density at the plasma edge, is found to yield the best agreement in terms of both the pump-out level and the overall time scale. Further sensitivity studies show that the simulated density pump-out level is reasonably robust against variations in the model assumptions, including the particle diffusion model and the non-ambipolar versus ambipolar NTV particle flux. The latter however affects the time scale for reaching the steady state solution. Finally, it is found that the plasma edge-peeling response, the NTV torque, as well as the plasma momentum and particle transport, all are sensitive to the toroidal phase difference between the upper and lower rows of the RMP coil currents in HL-2A, with the 30 ∘ coil phasing producing the minimal side effects on the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Model-independent measurement of the charge density distribution along an Fe atom probe needle using off-axis electron holography without mean inner potential effects

The one-dimensional charge density distribution along an electrically biased Fe atom probe needle is measured using a model-independent approach based on off-axis electron holography in the transmission electron microscope. Both the mean inner potential and the magnetic contribution to the phase shift are subtracted by taking differences between electron-optical phase images recorded with different voltages applied to the needle. The measured one-dimensional charge density distribution along the needle is compared with a similar result obtained using model-based fitting of the phase shift surrounding the needle. On the assumption of cylindrical symmetry, it is then used to infer the three-dimensional electric field and electrostatic potential around the needle with ∼10 nm spatial resolution, without needing to consider either the influence of the perturbed reference wave or the extension of the projected potential outside the field of view of the electron hologram. The present study illustrates how a model-independent approach can be used to measure local variations in charge density in a material using electron holography in the presence of additional contributions to the phase, such as those arising from changes in mean inner potential and specimen thickness.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Shear viscosity from perturbative quantum chromodynamics to the hadron resonance gas at finite baryon, strangeness, and electric charge densities

Through model-to-data comparisons from heavy-ion collisions, it has been shown that the quark gluon plasma has an extremely small shear viscosity at vanishing densities. At large baryon densities, significantly less is known about the nature of the shear viscosity from quantum chromodynamics (QCD). Within heavy-ion collisions, there are three conserved charges: baryon number (B), strangeness (S), and electric charge (Q). Here we calculate the shear viscosity in two limits using perturbative QCD (pQCD) and an excluded-volume hadron resonance gas at finite BSQ densities. We then develop a framework that interpolates between these two limits such that shear viscosity is possible to calculate across a wide range of finite BSQ densities. We find that the pQCD and hadron resonance gas calculations have different BSQ density dependencies such that a rather nontrivial shear viscosity appears at finite densities.

Phenomenology↗