Search NASA⌕ Search

SEARCH · Search NASA

Results for “density”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Enhanced particle transport events approaching the density limit of the J-TEXT tokamak

Enhanced particle transport events are discovered and analyzed as the density limit of the J-TEXT tokamak is approached. Edge shear layer collapse is observed and the ratio of Reynolds power to turbulence production decreases. Simultaneously, the divergence of turbulence internal energy flux (i.e. turbulence spreading) increases, indicating that shear layer collapse triggers an outward spreading event. Studies of correlations show that the enhanced particle transport events are quasi-coherent, and manifested primarily in density fluctuations which exhibit positive skewness. Furthermore, electron adiabaticity emerges as the critical parameter which signals transport event onset. For α < 0.35 as density approaches the Greenwald density, both turbulence spreading and density fluctuations rise rapidly. Taken together, these results elucidate the connections between edge shear layer, density fluctuations, particle transport events, turbulence spreading and plasma edge cooling as the density limit is approached.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Observation and analysis of lower-hybrid-current-drive density limit in EAST

Lower hybrid current drive (LHCD) experiments with line-averaged density up to ~5.1 × 10 19 m –3 were performed in EAST L-mode plasmas. When the line-averaged density rises above a critical value, the hard x-ray (HXR) emission falls to the noise level, indicating that the LHCD density limit is encountered. The experimental results show that the LHCD density limit can be increased with higher wave source frequency (ƒ 0 ) and higher magnetic field (B t ). Although a higher LHCD density limit is obtained by a higher magnetic field for both 2.45 GHz and 4.6 GHz waves, the results show a stronger dependence on the magnetic field for the 4.6 GHz case. Analysis suggests that, for normal operation with a relatively low magnetic field (1.6 T $\leqslant$ B t $\leqslant$ 2.5 T) on EAST, the dominant mechanisms responsible for the LHCD density limit are different between the 2.45 GHz and 4.6 GHz waves. The wave accessibility plays a more significant role during 4.6 GHz LHCD experiments, while parasitic losses due to parametric decay instability (PDIs) dominate the accessibility issue in the 2.45 GHz case. Collisional loss in the scrape-off layer (SOL) may explain the 4.6 GHz result when combined with the accessibility limit at high density and low temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The role of isotope mass on neutral fueling and density pedestal structure in the DIII-D tokamak

Abstract Experimental measurements on DIII-D of hydrogen neutral penetration lengths ( λ n 0 ) on the high field side (HFS) are longer by a factor of 2 than for deuterium consistent with the thermal velocity ratio for neutrals at the same temperature ( v t h H / v t h D = 2 ) . This ratio is constant for both low and high pedestal electron density. At low pedestal density ( n e ∼ 4 × 10 19 m − 3 ) , the neutral penetration length is greater than the density pedestal width for both isotopes, and the additional 40% increase of neutral penetration in hydrogen widens the pedestal by the same amount. As the density pedestal height increases ( n e ∼ 6 × 10 19 m − 3 ) , the neutral penetration lengths drop below the density pedestal widths for both isotopes, and the increased penetration of hydrogen has no increased effect on the pedestal width compared to deuterium. Extrapolating to future reactor-relevant high electron density pedestals, the isotope-mass change in neutral fueling on the HFS from the deepest neutral penetration of hydrogen, to the shortest neutral penetration of tritium will be negligible (0.2–0.4 cm) in comparison to estimates of the density pedestal width (6–8.5 cm).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Density limit in peeling-limited pedestals at and above the Greenwald value in DIII-D high poloidal beta plasmas

High pressure, peeling limited pedestals with pedestal normalized beta β N,ped >2 and pedestal top density n ped at or above the Greenwald density nG have been achieved in DIII-D high poloidal beta plasmas, with high global normalized beta β N >3 and energy confinement H 98 ~1.2-1.7. Higher β N allows higher pedestal density above the Greenwald value and higher pedestal pressure, even with a low injected torque. MHD modeling confirms that the experimental profiles lie near the peeling-mode unstable boundary with high normalized pressure gradient βMHD and high edge current density. Experimental analysis and stability calculations indicate that the high poloidal beta with strong Shafranov shift, high βMHD and weak/negative magnetic shear improves the pedestal stability by decoupling the peeling and ballooning modes and stabilizing the ballooning modes, thus facilitating access to the second stable region of peeling-ballooning mode. The access to the second stable peeling-ballooning stability region opens Super-H-like channels without extremely strong shaping or strong torque injection. The high-pressure peeling pedestal allows the pedestal density to go beyond the Greenwald limit with strong ExB shear maintained: pedestal pressure increases with pedestal density even when n ped >n G , until reaching the ideal MHD instability boundary, where giant ELMs occur. The giant ELMs are dominated by a strong n=1 component and cause a large reduction of the edge pressure, but a negligible change of the core pressure, consistent with kink/peeling-mode induced instability. The pedestal recovers from the collapse and typically sustains a high baseline density, around the Greenwald limit, during the whole discharge duration. Experiments also found that internal transport barriers and n ped ~0.9nG, peeling limited pedestals could be simultaneously achieved in high β N plasmas, while an internal feedback between ITB strength and pedestal performance is found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

COMPOSITIONS AND ORIGINS OF OUTER PLANET SYSTEMS: INSIGHTS FROM THE ROCHE CRITICAL DENSITY

We consider the Roche critical density ({rho}{sub Roche}), the minimum density of an orbiting object that, at a given distance from its planet, is able to hold itself together by self-gravity. It is directly related to the more familiar ''Roche limit,'' the distance from a planet at which a strengthless orbiting object of given density is pulled apart by tides. The presence of a substantial ring requires that transient clumps have an internal density less than {rho}{sub Roche}. Conversely, in the presence of abundant material for accretion, an orbiting object with density greater than {rho}{sub Roche} will grow. Comparing the {rho}{sub Roche} values at which the Saturn and Uranus systems transition rapidly from disruption-dominated (rings) to accretion-dominated (moons), we infer that the material composing Uranus' rings is likely more rocky, as well as less porous, than that composing Saturn's rings. From the high values of {rho}{sub Roche} at the innermost ring moons of Jupiter and Neptune, we infer that those moons may be composed of denser material than expected, or more likely that they are interlopers that formed farther from their planets and have since migrated inward, now being held together by internal material strength. Finally, the ''Portia group'' of eight closely packed Uranian moons has an overall surface density similar to that of Saturn's A ring. Thus, it can be seen as an accretion-dominated ring system, of similar character to the standard ring systems except that its material has a characteristic density greater than the local {rho}{sub Roche}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Generalization of the Kohn-Sham system that can represent arbitrary one-electron density matrices

Density functional theory is currently the most widely applied method in electronic structure theory. The Kohn-Sham method, based on a fictitious system of noninteracting particles, is the workhorse of the theory. The particular form of the Kohn-Sham wave function admits only idempotent one-electron density matrices whereas wave functions of correlated electrons in post-Hartree-Fock methods invariably have fractional occupation numbers. Here we show that by generalizing the orbital concept and introducing a suitable dot product as well as a probability density, a noninteracting system can be chosen that can represent the one-electron density matrix of any system, even one with fractional occupation numbers. This fictitious system ensures that the exact electron density is accessible within density functional theory. It can also serve as the basis for reduced density matrix functional theory. Moreover, to aid the analysis of the results the orbitals may be assigned energies from a mean-field Hamiltonian. This produces energy levels that are akin to Hartree-Fock orbital energies such that conventional analyses based on Koopmans' theorem are available. Lastly, this system is convenient in formalisms that depend on creation and annihilation operators as they are trivially applied to single-determinant wave functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ab initio electronic density in solids by many-body plane-wave auxiliary-field quantum Monte Carlo calculations

We present accurate many-body results of the electronic densities in several solid materials, including Si, NaCl, and Cu. These results are obtained using the ab initio auxiliary-field quantum Monte Carlo (AFQMC) method working in a plane-wave basis with norm-conserving, multiple-projector pseudopotentials. AFQMC has been shown to be an excellent many-body total energy method. Computation of observables and correlation functions other than the ground-state energy requires back-propagation, whose adaption and implementation in the plane-wave basis AFQMC framework are discussed in the present paper. This development allows us to compute correlation functions, electronic densities, and interatomic forces, paving the way for geometry optimizations and calculations of thermodynamic properties in solids. Finite supercell size effects are considerably more subtle in the many-body framework than in independent-electron calculations. We analyze the convergence of the electronic density, and obtain best estimates for the thermodynamic limit. The densities from several typical density functionals are benchmarked against our near-exact results. The electronic densities we obtained can also be used to help construct improved density functionals.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Imposing correct jellium response is key to predict the density response by orbital-free DFT

Orbital-free density functional theory constitutes a computationally highly effective tool for modeling electronic structures of systems ranging from room-temperature materials to warm dense matter. Its accuracy critically depends on the employed kinetic energy (KE) density functional, which has to be supplied as an external input. In this work we consider several nonlocal and Laplacian-level KE functionals and use an external harmonic perturbation to compute the static density response at T = 0 K in the linear and beyond-linear response regimes. We test for the satisfaction of exact conditions in the limit of uniform densities and for how approximate KE functionals reproduce the density response of realistic materials (e.g., Al and Si) against the Kohn-Sham DFT reference, which employs the exact KE. The results illustrate that several functionals violate exact conditions in the uniform electron gas (UEG) limit. We find a strong correlation between the accuracy of the KE functionals in the UEG limit and in the strongly inhomogeneous case. This empirically demonstrates the importance of imposing the limit of UEG response for uniform densities and validates the use of the Lindhard function in the formulation of kernels for nonlocal functionals. This conclusion is substantiated by additional calculations for bulk aluminum (Al) with a face-centered cubic (fcc) lattice and silicon (Si) with an fcc lattice, body-centered cubic (bcc) lattice, and semiconducting crystal diamond state. The analysis of fcc Al, and fcc as well as bcc Si data follows closely the conclusions drawn for the UEG, allowing us to extend our conclusions to realistic systems that are subject to density inhomogeneities induced by ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Skyrme pseudopotentials at next-to-next-to-leading order: Construction of local densities and first symmetry-breaking calculations

There is an ongoing quest to improve on the spectroscopic quality of nuclear energy density functionals (EDFs) of the Skyrme type through extensions of its traditional form. One direction for such activities is the inclusion of terms of higher order in gradients in the EDF. We report on exploratory symmetry-breaking calculations performed for an extension of the Skyrme EDF that includes central terms with four gradients at next-to-next-to-leading order (N2LO) and for which the high-quality parametrization SN2LO1 has been constructed recently. Up to now, the investigation of such functionals with higher-order terms was limited to infinite matter and spherically symmetric configurations of singly and doubly magic nuclei. We address here nuclei and phenomena that require us to consider axial and nonaxial deformation, both for reflection-symmetric and also reflection-asymmetric shapes, as well as the breaking of time-reversal invariance. Achieving these calculations demanded a number of formal developments. These all resulted from the formulation of the N2LO EDF requiring the introduction of new local densities with additional gradients that are not present in the EDF at NLO. Their choice is not unique, but can differ in the way the gradients are coupled. While designing a numerical implementation of N2LO EDFs in Cartesian three-dimensional coordinate-space representation, we have developed a novel definition and a new unifying notation for normal and pair densities that contain gradients at arbitrary order. Besides having mnemonic advantages, the new notation allows for the easy identification of redundancies and reducibilities in a given set of local densities, and the new definition makes it straightforward to construct densities that automatically adopt the symmetries of the many-body state they are constructed from. The resulting scheme resolves several issues with some of the choices that have been made for local densities in the past, in particular when breaking time-reversal symmetry. Guided by general practical considerations, we propose an alternative form of the N2LO contribution to the Skyrme EDF that is built from a different set of densities. It has exactly the same physics content, but is much more efficient to handle in formal discussions and, compared to the original formulation, leads to a substantial reduction of computational cost and memory requirements in deformed codes. As representative examples for the performance of SN2LO1, we have chosen the ground states of even-even Kr and Nd isotopes, the fission barrier of 240 Pu as well as the superdeformed rotational band of 194 Hg. Overall, for the nuclei and phenomena studied here, the SN2LO1 parametrization does not yet present a systematic improvement over standard NLO parametrizations. This finding calls for improved fit protocols that better discriminate between NLO and N2LO terms and better exploit the unique features of the additional degrees of freedom offered by the latter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Density- and wavefunction-normalized Cartesian spherical harmonics for l ≤ 20

The widely used pseudoatom formalism in experimental X-ray charge-density studies makes use of real spherical harmonics when describing the angular component of aspherical deformations of the atomic electron density in molecules and crystals. The analytical form of the density-normalized Cartesian spherical harmonic functions for up to l ≤ 7 and the corresponding normalization coefficients were reported previously by Paturle & Coppens. It was shown that the analytical form for normalization coefficients is available primarily for l ≤ 4. Only in very special cases it is possible to derive an analytical representation of the normalization coefficients for 4 < l ≤ 7. In most cases for l > 4 the density normalization coefficients were calculated numerically to within seven significant figures. In this study we review the literature on the density-normalized spherical harmonics, clarify the existing notations, use the Paturle–Coppens method in the Wolfram Mathematica software to derive the Cartesian spherical harmonics for l ≤ 20 and determine the density normalization coefficients to 35 significant figures, and computer-generate a Fortran90 code. The article primarily targets researchers who work in the field of experimental X-ray electron density, but may be of some use to all who are interested in Cartesian spherical harmonics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

eFESTA: Ensemble Feature Exploration with Surface Density Estimates

We propose surface density estimate (SDE) to model the spatial distribution of surface features-isosurfaces, ridge surfaces, and streamsurfaces-in 3D ensemble simulation data. The inputs of SDE computation are surface features represented as polygon meshes, and no field datasets are required (e.g., scalar fields or vector fields). The SDE is defined as the kernel density estimate of the infinite set of points on the input surfaces and is approximated by accumulating the surface densities of triangular patches. We also propose an algorithm to guide the selection of a proper kernel bandwidth for SDE computation. An ensemble Feature Exploration method based on Surface densiTy EstimAtes (eFESTA) is then proposed to extract and visualize the major trends of ensemble surface features. For an ensemble of surface features, each surface is first transformed into a density field based on its contribution to the SDE, and the resulting density fields are organized into a hierarchical representation based on the pairwise distances between them. The hierarchical representation is then used to guide visual exploration of the density fields as well as the underlying surface features. We demonstrate the application of our method using isosurface in ensemble scalar fields, Lagrangian coherent structures in uncertain unsteady flows, and streamsurfaces in ensemble fluid flows.

Computational modeling↗
DOE OSTI2019

Effects of group size and group density on trade–offs in resource selection by a group–territorial central–place foraging woodpecker

Trade-offs in resource selection by central-place foragers are driven by the need to balance the benefits of selecting resources against the costs of travel from the central place. For group-territorial central-place foraging birds, trade-offs in resource selection are likely to be complicated by a competitive advantage for larger groups at high group density that may limit accessibility of high-quality distant resources to small groups. We used the group-territorial, central-place foraging Red-cockaded Woodpecker Leuconotopicus borealis (RCW) as a case study to test predictions that increases in group density lead to differences in foraging distances and resource selection for groups of different sizes. We used GPS tracking and LiDAR-derived habitat data to model effects of group size on foraging distances and selection for high-quality pines (≥ 35.6 cm diameter at breast height (dbh)) and lower quality pines (25.4–35.6 cm dbh) by RCW groups across low ( n = 14), moderate ( n = 10) and high group density ( n = 10) conditions. At low and moderate group density, all RCW groups selected distant high-quality pines in addition to those near the central place because competition for resources was low. In contrast, at high group density, larger groups travelled further to select high-quality pines, whereas smaller groups selected high-quality pines only when they were close to the central place and, conversely, were more likely to select lower quality pines at greater distances from the central place. Selection for high-quality pines only when close to the cavity tree cluster at high group density is important to long-term fitness of small RCW groups because it allows them to maximize benefits from both territorial defence and selecting high-quality resources while minimizing costs of competition. Furthermore these relationships suggest that intraspecific competition at high group density entails substantive costs to smaller groups of territorial central-place foragers by limiting accessibility of distant high-quality foraging resources.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

MO-F-CAMPUS-J-04: Tissue Segmentation-Based MR Electron Density Mapping Method for MR-Only Radiation Treatment Planning of Brain

Purpose: Automatically derive electron density of tissues using MR images and generate a pseudo-CT for MR-only treatment planning of brain tumours. Methods: 20 stereotactic radiosurgery (SRS) patients’ T1-weighted MR images and CT images were retrospectively acquired. First, a semi-automated tissue segmentation algorithm was developed to differentiate tissues with similar MR intensities and large differences in electron densities. The method started with approximately 12 slices of manually contoured spatial regions containing sinuses and airways, then air, bone, brain, cerebrospinal fluid (CSF) and eyes were automatically segmented using edge detection and anatomical information including location, shape, tissue uniformity and relative intensity distribution. Next, soft tissues - muscle and fat were segmented based on their relative intensity histogram. Finally, intensities of voxels in each segmented tissue were mapped into their electron density range to generate pseudo-CT by linearly fitting their relative intensity histograms. Co-registered CT was used as a ground truth. The bone segmentations of pseudo-CT were compared with those of co-registered CT obtained by using a 300HU threshold. The average distances between voxels on external edges of the skull of pseudo-CT and CT in three axial, coronal and sagittal slices with the largest width of skull were calculated. The mean absolute electron density (in Hounsfield unit) difference of voxels in each segmented tissues was calculated. Results: The average of distances between voxels on external skull from pseudo-CT and CT were 0.6±1.1mm (mean±1SD). The mean absolute electron density differences for bone, brain, CSF, muscle and fat are 78±114 HU, and 21±8 HU, 14±29 HU, 57±37 HU, and 31±63 HU, respectively. Conclusion: The semi-automated MR electron density mapping technique was developed using T1-weighted MR images. The generated pseudo-CT is comparable to that of CT in terms of anatomical position of tissues and similarity of electron density assignment. This method can allow MR-only treatment planning.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Multifactorial genetic control and magnesium levels govern the production of a Streptomyces antibiotic with unusual cell density dependence

Streptomyces bacteria are renowned both for their antibiotic production capabilities and for their cryptic metabolic potential. Their metabolic repertoire is subject to stringent genetic control, with many of the associated biosynthetic gene clusters being repressed by the conserved nucleoid-associated protein Lsr2. In an effort to stimulate new antibiotic production in wild Streptomyces isolates, we leveraged the activity of an Lsr2 knockdown construct and successfully enhanced antibiotic production in the wild Streptomyces isolate WAC07094. We determined that this new activity stemmed from increased levels of the angucycline-like family member saquayamycin. Saquayamycin has both antibiotic and anti-cancer activities, and intriguingly, beyond Lsr2-mediated repression, we found saquayamycin production was also suppressed at high density on solid or in liquid growth media; its levels were greatest in low-density cultures. This density-dependent control was exerted at the level of the cluster-situated regulatory gene sqnR and was mediated in part through the activity of the PhoRP two-component regulatory system, where deleting phoRP led to both constitutive antibiotic production and sqnR expression. This suggests that PhoP functions to repress the expression of sqnR at high cell density. We further discovered that magnesium supplementation could alleviate this density dependence, although its action was independent of PhoP. Finally, we revealed that the nitrogen-responsive regulators GlnR and AfsQ1 could relieve the repression exerted by Lsr2 and PhoP. Intriguingly, we found that this low density-dependent production of saquayamycin was not unique to WAC07094; saquayamycin production by another wild isolate also exhibited low-density activation, suggesting that this spatial control may serve an important ecological function in their native environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

On the density of states of disordered epitaxial graphene

The study is concerned with two types of disordered epitaxial graphene: (i) graphene with randomly located carbon vacancies and (ii) structurally amorphous graphene. The former type is considered in the coherent potential approximation, and for the latter type, a model of the density of states is proposed. The effects of two types of substrates, specifically, metal and semiconductor substrates are taken into account. The specific features of the density of states of epitaxial graphene at the Dirac point and the edges of the continuous spectrum are analyzed. It is shown that vacancies in epitaxial graphene formed on the metal substrate bring about logarithmic nulling of the density of states of graphene at the Dirac point and the edges of the continuous spectrum. If the Dirac point corresponds to the middle of the band gap of the semiconductor substrate, the linear trend of the density of states to zero in the vicinity of the Dirac point in defect-free graphene transforms into a logarithmic decrease in the presence of vacancies. In both cases, the graphene-substrate interaction is assumed to be weak (quasi-free graphene). In the study of amorphous epitaxial graphene, a simple model of free amorphous graphene is proposed as the initial model, in which account is taken of the nonzero density of states at the Dirac point, and then the interaction of the graphene sheet with the substrate is taken into consideration. It is shown that, near the Dirac point, the quadratic behavior of the density of states of free amorphous graphene transforms into a linear dependence for amorphous epitaxial graphene. In the study, the density of states of free graphene corresponds to the low-energy approximation of the electron spectrum.

36 MATERIALS SCIENCE↗
DOE OSTI2013

High Energy Density Physics and applications with a state-of-the art compact x-pinch

Recent advances in technology has made possible to create matter with extremely high energy density (energy densities and pressure exceeding 10 11 J/m 3 and 1 Mbar respectively). The field is new and complex. The basic question for high energy density physics (HEDP) is how does matter behave under extreme conditions of temperature, pressure, density and electromagnetic radiation? The conditions for studying HEDP are normally produced using high intensity short pulse laser, x-rays, particle beams and pulsed power z-pinches. Most of these installations occupy a large laboratory floor space and require a team consisting of a large number of scientists and engineers. This limits the number of experiments that can be performed to explore and understand the complex physics. A novel way of studying HEDP is with a compact x-pinch in university scale laboratory. The x-pinch is a configuration in which a pulsed current is passed through two or more wires placed between the electrodes making the shape of the letter ‘X’. Extreme conditions of magnetic field (> 200 MGauss for less than 1 ns), temperature (1 keV) and density (~ 10 22 cm -3 ) are produced at the cross-point, where two wires make contact. Further, supersonic jets are produced on either side of the cross-point. The physics of the formation of the plasma at the cross-point is complex. It is not clear what role radiation plays in the formation of high energy density plasma (>> 10 11 J/m 3 ) at the cross-point. Nor it is understood how the supersonic jets are formed. Present numerical codes do not contain complex physics that can take into account some of these aspects. Indeed, a comprehensive experimental study could answer some of the questions, which are relevant to wide-ranging fields such as inertial confinement fusion, astrophysical plasmas, high intensity laser plasma interactions and radiation physics. The main aim of the proposal was to increase the fundamental understanding of high energy density physics and particularly address the key issues associated with x-pinches, which include radiation transport, energetic particle transport, supersonic jet formation, using state-of-the-art compact pulsed power drivers. All the primary objectives of the proposed work were met. These objectives include: • Understanding of the fundamental physics of hot and dense plasma formation, implosion to less than 1 µm size due to the radiation enhanced collapse and energetic electron heating, • Study of the jet formation mechanism, which is of interest due to the astrophysical jets and deposition of energy by energetic electrons in jets, • Characterization of an x-pinch as a point x-ray source for the phase contrast radiography of beryllium cryogenic targets for the National Ignition Facility (NIF) experiments. The work carried out included a strong educational component involving both undergraduate and graduate students. Several undergraduate students from University of California San Diego participated in this project. A post-doctoral fellow, Dr. Simon Bott and two graduate students, David Haas and Erik Shipton contributed to every aspect of this project. The success of the project can be judged from the fact that fifteen peer-reviewed papers were published in high quality journals. In addition several presentations were made to a number of scientific meetings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Connecting the dots: Tracking galaxy evolution using constant cumulative number density at 3 ≤ z ≤ 7

Using the cosmological smoothed particle hydrodynamical code GADGET-3, we make a realistic assessment of the technique of using constant cumulative number density as a tracer of galaxy evolution at high redshift. We find that over a redshift range of 3 ≤ z ≤ 7 one can on average track the growth of the stellar mass of a population of galaxies selected from the same cumulative number density bin to within ∼0.20 dex. Over the stellar mass range that we probe (10{sup 10.4}⩽M{sub s}/M{sub ⊙}⩽10{sup 10.8} at z = 3 and 10{sup 8.5}⩽M{sub s}/M{sub ⊙}⩽10{sup 9.6} at z = 7), one can reduce this bias by selecting galaxies based on an evolving cumulative number density. We find that this cumulative number density evolution exhibits a trend toward higher values which can be quantified by simple linear formulations of −0.10Δz for descendants and 0.12Δz for progenitors. Utilizing such an evolving cumulative number density increases the accuracy of descendant/progenitor tracking by a factor of ∼2. This result is in excellent agreement, within 0.10 dex, with abundance matching results over the same redshift range. However, we find that our more physically realistic cosmological hydrodynamic simulations produce a much larger scatter in descendant/progenitor stellar masses than previous studies, particularly when tracking progenitors. This large scatter makes the application of either the constant cumulative number density or evolving cumulative number density technique limited to average stellar masses of populations only, as the diverse mass assembly histories caused by stochastic physical processes such as gas accretion or mergers lead to an even larger scatter in other physical properties such as metallicity and star formation rate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

NITROGEN CAN LIMIT OVERSTORY TREE GROWTH FOLLOWING EXTREME STAND DENSITY INCREASE IN A PONDEROSA PINE FOREST

Extreme stand density increases have occurred in ponderosa pine forests throughout the western U.S. since the early 20th Century, with adverse implications for growth, physiological functioning, and mortality risk. Identifying primary stressors on large, old overstory trees in dense forests can inform management decisions to promote resilience and survival. We tested the impact of stand density increase on overstory treering growth, and the relative influence of water and nitrogen, in an old-growth ponderosa pine forest subject to variable density increase. We measured annual tree-ring growth and carbon discrimination in trees before stand density increased, in a climatically-similar period post-density increase, and in recent transition to drought. We expected density-driven water stress to drive reduced tree-ring growth in overstory trees in dense stands. We found reduced growth and higher mortality in dense stands, but nitrogen rather than water constrained growth, as determined by carbon isotope discrimination in tree rings, leaf nitrogen concentration, and soil nitrogen supply. In dense stands, less available nitrogen limited photosynthetic rate, leading to reduced assimilation of intracellular 13C and higher discrimination with low tree-ring growth and a reduced relationship with climate. This unexpected result illustrates that a variety of limiting factors can influence forest dynamics, as density-driven nitrogen limitation interacts with water stress to influence tree growth and physiological functioning.

dendroecology, tree-ring width, limiting factors, ↗