Search NASA⌕ Search

SEARCH · Search NASA

Results for “zero energy state”

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 397 records · Page 22

Variability in Wind Energy Generation across the Contiguous United States

ERA5 provides high-resolution, high-quality hourly wind speeds at 100 m and is a unique resource for quantifying temporal variability in likely wind-derived power production across the United States. Gross capacity factors (CF) in seven independent system operators (ISOs) are estimated using the location and rated power of each wind turbine, a simplified power curve, and ERA5 output from 1979 to 2018. Excluding the California ISO, the marginal probability of a calm (zero power production) is less than 0.1 in any ERA5 grid cell. When a calm occurs, the mean co-occurrence across wind-turbine-containing grid cells ranges from 0.38 to 0.39 for ISOs in the Midwest and central plains [Midcontinent (or Midwest) ISO (MISO), Southwest Power Pool (SPP), and the Electric Reliability Council of Texas (ERCOT) region], increasing to 0.54–0.58 for ISOs in the eastern United States [Pennsylvania–New Jersey–Maryland interconnection (PJM), New York ISO (NYISO), and New England ISO (NEISO)]. Periods with low gross CF have a median duration of ≤6 h, except in California, and are most likely during summer. Additionally, gross CF exhibit highest variance at periods of 1 day in ERCOT and SPP; on synoptic scales in MISO, NEISO, and NYISO; and on interannual time scales in PJM. This implies differences in optimal strategies for ensuring resilience of supply. Theoretical scenarios show adding wind energy capacity near existing wind farms is advantageous even in areas with high existing installed capacity (IC), while expanding into areas with lower IC is more beneficial to reducing ramps and the probability of gross CF falling below 20%. These results emphasize the benefits of large balancing areas and aggregation in reducing wind power variability and the likelihood of wind droughts.

17 WIND ENERGY↗

Approximate Quantum Codes From Long Wormholes

We discuss families of approximate quantum error correcting codes which arise as the nearly-degenerate ground states of certain quantum many-body Hamiltonians composed of non-commuting terms. For exact codes, the conditions for error correction can be formulated in terms of the vanishing of a two-sided mutual information in a low-temperature thermofield double state. We consider a notion of distance for approximate codes obtained by demanding that this mutual information instead be small, and we evaluate this mutual information for the SYK model and for a family of low-rank SYK models. After an extrapolation to nearly zero temperature, we find that both kinds of models produce fermionic codes with constant rate as the number, N , of fermions goes to infinity. For SYK, the distance scales as N 1 / 2 , and for low-rank SYK, the distance can be arbitrarily close to linear scaling, e.g. N .99 , while maintaining a constant rate. We also consider an analog of the no low-energy trivial states property which we dub the no low-energy adiabatically accessible states property and show that these models do have low-energy states that can be prepared adiabatically in a time that does not scale with system size N . We discuss a holographic model of these codes in which the large code distance is a consequence of the emergence of a long wormhole geometry in a simple model of quantum gravity.

Physics↗

Electric Vehicle Infrastructure Projection Tool (EVI-Pro)

This presentation includes results from the second California installment of EVI-Pro. Pursuant to California Assembly Bill 2127, evolving market and technology conditions warrant updating this statewide infrastructure assessment at least every two years. EVI-Pro has been updated to consider California's 2030 goal to have 5 million zero emission vehicles (ZEVs) on the road by 2030. CEC and NREL with the support of UC Davis and other state agencies, have set out to refine EVI-Pro to reflect increasing PEV market share, evolving vehicle and charging technology, and observed charging behavior.

ADVANCED PROPULSION SYSTEMS,ENERGY STORAGE↗

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗

Demystifying group-4 polyolefin hydrogenolysis catalysis: Gaseous propane hydrogenolysis mechanism over the same catalysts

A kinetic/mechanistic investigation of gaseous propane hydrogenolysis over the single-site heterogeneous polyolefin depolymerization catalysts AlS/ZrNp 2 and AlS/HfNp 2 (AlS = sulfated alumina, Np = neopentyl), is use to probe intrinsic catalyst properties without the complexities introduced by time- and viscosity-dependent polymer medium effects. In a polymer-free automated plug-flow catalytic reactor, propane hydrogenolysis turnover frequencies approach 3,000 h −1 at 150 °C. Both catalysts exhibit approximately linear relationships between rate and [H 2 ] at substoichiometric [H 2 ] with rate law orders of 0.66 ± 0.09 and 0.48 ± 0.07 for Hf and Zr, respectively; at higher [H 2 ], the rates approach zero-order in [H 2 ]. Reaction orders in [C 3 H 8 ] and [catalyst] are essentially zero-order under all conditions, with the former implying rapid, irreversible alkane binding/activation. This rate law, activation parameter, and DFT energy span analysis support a scenario in which [H 2 ] is pivotal in one of two plausible and competing rate-determining transition states—bimolecular metal-alkyl bond hydrogenolysis vs. unimolecular β-alkyl elimination. The Zr and Hf catalyst activation parameters, ΔH ‡ = 16.8 ± 0.2 kcal mol −1 and 18.2 ± 0.6 kcal mol −1 , respectively, track the relative turnover frequencies, while ΔS ‡ = −19.1 ± 0.8 and −16.7 ± 1.4 cal mol −1 K −1 , respectively, imply highly organized transition states. These catalysts maintain activity up to 200 °C, while time-on-stream data indicate multiday activities with an extrapolated turnover number ~92,000 at 150 °C for the Zr catalyst. This methodology is attractive for depolymerization catalyst discovery and process optimization.

03 NATURAL GAS↗

Assessing Nutrient and Carbon Responses to Agricultural Conservation Practices in Two Midwest Watersheds

The United States is undertaking efforts to transition to a low-carbon economy in response to the heightened impacts of climate change, which are largely attributed to decades of carbon-intensive development. A concerted effort by the energy sector to achieve net-zero carbon emissions is aimed at strategically reducing the carbon footprint of the energy supply chain, particularly in the production of feedstocks for biofuels. The production of biofuels relies heavily on land use and management practices in biomass cultivation. Among other factors, soil organic carbon (SOC) plays a crucial role in the biofuel carbon cycle, impacting land productivity, greenhouse gas emissions, and water quality. The process of carbon drawdown during plant growth and storage helps to mitigate the release of carbon dioxide (CO 2 ) into the atmosphere. Another key parameter is the release of nitrous oxide (N 2 O) — a potent greenhouse gas with a global warming potential 273 times that of CO 2 — from soil. While prioritizing low-carbon production, sustainable bioenergy also necessitates improved water quality and ecosystem services. Agricultural conservation practices can contribute to enhanced soil health and reduce nutrient and sediment loss into streams. However, studies that explore the broader impact of these practices on soil carbon storage and N 2 O emissions at a watershed scale are limited. Specifically, our understanding of how low-carbon feedstock production and management affect nutrient cycle dynamics is incomplete. This study assesses watershed responses to land management practices in two agriculturally dominant watersheds: the Raccoon River watershed and the Southfork of Iowa River watershed. The study focuses on carbon and nutrient dynamics, characterizing spatial and temporal variations in SOC, N 2 O, and nutrient loadings to examine the relationship between nutrient and carbon responses and agricultural conservation practices. The study employs the newly developed Soil Water Analysis Tool for Carbon (SWAT-C) model, calibrated using 20 years’ of climate and water monitoring data, to simulate and evaluate two agricultural conservation practices: no-till and crop residue harvest with cover crop planting. We compared the calibrated SWAT-C model with a historical baseline model, which allowed us to analyze various aspects of the watershed, including stream flow, suspended sediments, nitrogen, phosphorus, organic carbon, SOC at different depths, and N 2 O emissions from topsoil. Our SWAT-C modeling results indicate that, compared with results from the historical baseline model, no-till practices are positively correlated with SOC accumulation, reduced N 2 O emissions, and decreased soil loss. We also observed no change or slightly increased nitrogen and phosphorus loss to water bodies in the watersheds compared with the baseline. Conversely, crop residue harvest with cover crop planting improved water quality by reducing nutrient and soil losses but also increased N 2 O emissions. This SWAT-C modeling study establishes the groundwork for further smaller-scale and/or sub-basin-level assessments of nitrogen, phosphorus, and carbon cycles. It provides valuable science-based insights to inform policy decisions — particularly in the context of transitioning to biofuel feedstock production in these watersheds — by addressing concerns about nutrient exports to local water bodies and, ultimately, the Mississippi River.

54 ENVIRONMENTAL SCIENCES↗

Mitigating urban climate-energy feedback with citywide building-integrated photovoltaics implementation

Urban heat islands (UHIs) intensify cooling demands, whose waste heat further exacerbates UHI, creating a self-reinforcing feedback loop. Building-integrated photovoltaics (BIPV) can break this cycle by generating electricity and altering urban climate-energy interactions, yet its net impact remains unresolved. Using a cross-scale modeling framework for 118,521 buildings in Hong Kong, we demonstrate that BIPV glazing reduces cooling demand through improved insulation and lower solar heat gain, while inducing diurnal thermal asymmetry of daytime warming (up to +1.9°C) and nighttime cooling (up to −2.5°C). This microclimate regulation indirectly contributes to a 0.6% reduction in citywide building energy use. Cumulative savings from electricity generation and reduced energy use reach 4.7% (2,060.2 GWh) and 10.1% (4,401.3 GWh) under low- and high-coverage deployment, respectively, with building-level savings of −20.1% ± 11.7% (mean ± SD) at high coverage. We establish BIPV as a configurable climate-energy regulator that enables net-zero, heat-resilient planning and quantifies deployment trade-offs in cooling-dominated cities worldwide.

Building-integrated photovoltaics (BIPV)↗

Zero-Emission Transit Bus Needs Assessment

The transition to zero-emissions vehicles (ZEVs) in public transit has gained traction due to significant federal investments from the Bipartisan Infrastructure Law (BIL) and the Inflation Reduction Act (IRA). This needs assessment, commissioned by the Joint Office of Energy and Transportation and conducted by researchers at the Idaho National Laboratory, explores the current state of electrification in transit agencies, identifying barriers to implementation, potential funding sources, and operational considerations necessary for a successful transition. The assessment involved qualitative interviews with representatives from 19 transit service providers across diverse geographic regions. Key findings highlight the challenges related to bus facilities and operations, which require careful planning for charging infrastructure and maintenance capabilities to accommodate battery electric buses (BEBs) and hydrogen fuel cell buses (HFCBs). Agencies reported operational hurdles due to the shorter range of BEBs compared to diesel buses, necessitating revised scheduling and routing strategies. Despite these challenges, many agencies expressed optimism about their capacity to adapt. Funding availability emerged as a critical factor influencing the transition to ZEVs. While agencies welcomed increased financial support, particularly from the Low or No Emission Grant Program (Lo-No), concerns about the sustainability of this funding and the ongoing operational costs were prevalent. The need for a comprehensive funding inventory was underscored to ensure transit agencies are aware of all available resources. Technological constraints were significant barriers to ZEV adoption. The limited range of BEBs was frequently cited as a concern, leading to operational challenges and reliability issues. Agencies reported difficulties in sourcing replacement parts, which exacerbated downtime and maintenance challenges. Workforce development and training were identified as pivotal for a successful transition. Many agencies rely heavily on manufacturers for technician training, highlighting the need for scalable training programs that equip staff with the necessary skills to maintain electric powertrains effectively. This assessment offers actionable recommendations for the Joint Office, including enhancing outreach to transit agencies, developing resources for effective utility partnerships, and facilitating comprehensive training programs. Establishing a zero-emission bus evaluation program to track performance metrics such as cost, range, and reliability could provide valuable insights for transit agencies. The needs assessment provides a detailed examination of the challenges and opportunities facing transit agencies in their transition to zero-emissions bus fleets. By addressing these issues through targeted support, stakeholders can collaboratively work towards a cleaner, more sustainable public transportation system that benefits all communities.

33 - ADVANCED PROPULSION SYSTEMS↗

Comparative techno-economic assessment of osmotically-assisted reverse osmosis and batch-operated vacuum-air-gap membrane distillation for high-salinity water desalination

New developments in pressure- and thermally driven membrane desalination technologies offer the potential to cost-effectively treat high-salinity waters, especially when powered by low-cost solar electricity and thermal energy. This paper presents a comparative techno-economic assessment of the state-of-the-art most promising pressure- and thermally driven membrane technologies for high recovery desalination, namely, osmotically-assisted reverse osmosis (OARO) and batch-operated vacuum-air-gap membrane distillation (batch V-AGMD), to produce potable water while concentrating brine within the range of 140–290 g/L TDS for minimum-liquid-discharge (MLD) and zero-liquid-discharge (ZLD) applications. It is shown that both OARO and batch V-AGMD can treat feedwater and brines with TDS in the range of 30–125 g/L with corresponding fresh water recovery rates of 85–25%. When low cost solar electricity and thermal energy are used, the resulting levelized cost of water (LCOW) from OARO is in the range of 0.70–6.28 $/m 3 , and that from batch-V-AGMD is in the range of 1.74–2.77 $/m 3 . OARO is more cost-effective than batch V-AGMD when feedwater salinity is below 70 g/L and recovery below 75%, whereas batch V-AGMD is more cost-effective at higher recovery rates and salinity levels. Finally, the sensitivity of this comparison on energy prices and module costs is discussed.

42 ENGINEERING↗

Carbon management technology pathways for reaching a U.S. Economy-Wide net-Zero emissions goal

The Carbon Management Study Group of the 37 th Energy Modeling Forum (EMF 37) designed seven scenarios to explore the role of three potentially key technology suites – point source carbon dioxide capture and storage (PSCCS), direct air capture of carbon dioxide (DACCS), and hydrogen systems (H 2 ) – in shaping the broader technology pathways to reaching net-zero carbon dioxide (CO 2 ) emissions in United States by 2050. Each scenario was run by up to 13 models participating in the EMF 37 study. Results show that carbon dioxide removal technologies were consistently a major part of successful pathways to net-zero U.S. CO 2 emissions in 2050. Achieving this net-zero CO 2 goal without any form of carbon dioxide capture and storage was found to be impossible for most models; some models also found it impossible to reach net-zero without DACCS. The marginal cost of achieving net-zero CO 2 emissions in 2050 was between two and 10 times higher without PSCCS and/or DACCS available. The carbon price at which DACCS was deployed as a backstop technology depended upon the assumed cost at which DACCS was available at scale. Carbon prices were between $\$$250 and $\$$500 per ton CO 2 when DACCS deployed as a backstop. The average CO 2 capture rate across all models in 2050 in the central net-zero scenario was 1.3 GtCO 2 /year, which implies a substantial upscaling of capacity to move and store CO 2 . Finally, hydrogen sensitivity scenarios showed that H 2 typically constituted a relatively small share of the overall U.S. energy system; however, H 2 deployed in applications that are considered hard to decarbonize, facilitating transition towards net-zero emissions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Experimental study of the 34 m Cl beam production at intermediate energies

Here the isomeric content of a 34 Cl beam produced in the intermediate-energy projectile fragmentation of a 150 MeV/u 36 Ar beam on a 3 mm-thick Be target was studied. β-delayed γ-ray spectroscopy was used to measure the population of 34 Cl fragments in the ground vs. isomeric states at zero degrees relative to the incoming primary beam for four different momentum settings of the fragment separator near the predicted central velocity of these fragments, as well as, at two non-zero-degree settings for one momentum setting. Of the settings explored, which excluded rigidities within 0.5% of the value predicted to maximize total 34 Cl yield due to unreacted primary beam, the maximum rate for the production of Cl was found at a rigidity setting 0.75% below the predicted peak 34 Cl yield. The maximum population of the isomeric state relative to the ground state was observed at a rigidity 1.25% below the predicted maximum 34 Cl yield. Studies such as this are important in generating the understanding needed for producing isomer-enriched rare-isotope beams.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Reactivity-initiated accidents in two pressurized water reactor high burnup core designs

Here, this paper presents a safety analysis of two proposed core loadings for 24-month Pressurized Water Reactor (PWR) fuel cycles. This analysis focuses on reactivity-initiated accidents (RIAs) and evaluates core safety performance impacts of rod-averaged burnup limits up to 75 GWD/MTU and less than 7 % enriched UO 2 . The capabilities of Polaris, PARCS, and RELAP5-3D are leveraged to evaluate the core neutronic and thermal–hydraulic behavior for normal-operation, uncontrolled control rod withdrawal (CRW) transients, and control rod ejection (CRE) accidents. The two core designs are compared to identify features of realistic high burnup/extended enrichment core design approaches which have significant safety impact, identify experimental data needs for high-fidelity predictive modeling, and provide recommendations for future high burnup core designs. The first core design evaluated in this study was developed by Southern Nuclear Company and used an ZrB 2 Integral Fuel Burnable Absorber (IFBA) and B 4 C Wet-Annular Burnable Absorber (WABA)-based burnable poison strategy. The second core design assessed in this work used a Gd 2 O 3 -doped UO 2 burnable poison, similar to that used in boiling water reactors or French PWRs. Results indicate that fuel thermal limits are maintained for limiting CRW and hot full power (HFP) CRE transients. Cladding failure is predicted for the highest energy deposition rods in each core during limiting hot zero power (HZP) CRE accidents (where maximum radially averaged enthalpy exceeds 120 cal/g), though licensing may be permissible with a limited number of failed rods. While concerns exist regarding high critical boron concentration during steady state for the IFBA core and large plenum pressures for the gadolinia core design, the analysis demonstrates adequate safety performance during limiting RIA accident scenarios for two representative high burnup core designs. Design changes limiting plenum pressures and implementation of accident tolerant fuel (ATF) cladding features which minimize hydriding and susceptibility to pellet-cladding mechanical interaction (PCMI) are recommended for future high burnup fuel concepts. To support the technical basis for burnup limit increases, high-fidelity fuel performance models are needed to address physical effects not considered in this analysis, and high burnup irradiated fuel tests are required to extend applicability of the fuel failure limits and validate existing and future models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Anomalous transport due to Weyl fermions in the chiral antiferromagnets Mn 3 X, X = Sn, Ge

The recent discoveries of strikingly large zero-field Hall and Nernst effects in antiferromagnets Mn 3 X(X = Sn, Ge) have brought the study of magnetic topological states to the forefront of condensed matter research and technological innovation. These effects are considered fingerprints of Weyl nodes residing near the Fermi energy, promoting Mn 3 X(X = Sn, Ge) as a fascinating platform to explore the elusive magnetic Weyl fermions. In this review, we provide recent updates on the insights drawn from experimental and theoretical studies of Mn 3 X(X = Sn, Ge) by combining previous reports with our new, comprehensive set of transport measurements of high-quality Mn 3 Sn and Mn 3 Ge single crystals. In particular, we report magnetotransport signatures specific to chiral anomalies in Mn 3 Ge and planar Hall effect in Mn 3 Sn, which have not yet been found in earlier studies. The results summarized here indicate the essential role of magnetic Weyl fermions in producing the large transverse responses in the absence of magnetization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Extending a scaling equation of state to QCD

Whether quantum chromodynamics (QCD) exhibits a phase transition at finite temperature and density is an open question. It is important for hydrodynamic modeling of heavy ion collisions and neutron-star mergers. Lattice QCD simulations have definitively shown that the transition from hadrons to quarks and gluons is a crossover when the baryon chemical potential is zero or small. Here, we combine the parametric scaling equation of state, usually associated with the three-dimensional Ising model, with a background equation of state based on a smooth crossover from hadrons to quarks and gluons. Comparison to experimental data from the Beam Energy Scan II at the Relativistic Heavy Ion Collider or in heavy ion experiments at other accelerators may allow the critical exponents and amplitudes in the scaling equation of state to be determined for QCD if a critical point exists.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Finite-range pairing in nuclear density functional theory

Pairing correlations are ubiquitous in low-energy states of atomic nuclei. To incorporate them within nuclear density functional theory, often used for global computations of nuclear properties, pairing functionals that generate nucleonic pair densities and pairing fields are introduced. Many pairing functionals currently used can be traced back to zero-range nucleon-nucleon interactions. Unfortunately, such functionals are plagued by deficiencies that become apparent in large model spaces that contain unbound single-particle (continuum) states. In particular, the underlying computational schemes diverge as the single-particle space increases, and the results depend on how marginally occupied states are incorporated. These problems become more pronounced for pairing functionals that contain gradient-density dependence, such as in the Fayans functional. To remedy this, finite-range pairing functionals are introduced. In this study, this is done by folding the pair density with Gaussians. Here, we show that a folding radius of about 1 fm offers the best compromise between quality and stability, and substantially reduces the pathological behavior in different numerical applications.

Nuclear density functional theory↗

The Energy in Modular (EMOD) Buildings Method: A Guide to Energy-Efficient Design for Industrialized Construction of Modular Buildings

Industrialized construction has immense potential to address the growing need globally to build and upgrade the building stock to be affordable, energy-efficient, and resilient. It can also help achieve the United States' goal of a 50% reduction in U.S. greenhouse gas (GHG) emissions by 2030. Despite this potential, and the ever-increasing push for electrification and decarbonization of households in the United States, industrialized construction has not yet been leveraged specifically to help address these challenges and accelerate the pathway to meet these goals. The National Renewable Energy Laboratory (NREL) aims to claim this missed opportunity by focusing on delivering affordable, grid-efficient net-zero energy (NZE) modular buildings for underserved communities to ensure an equitable transition to the future of clean energy, accelerate decarbonization of the built environment, and support the development of a high-productivity construction and energy efficiency workforce. The Energy in Modular (EMOD) method is our approach to designing, producing, and delivering affordable, net-zero energy, low-carbon, and healthier buildings at scale. The following energy efficiency strategies are part of the scope of this guide: envelope thermal control, envelope infiltration control, mechanical, electrical, and plumbing systems, smart controls, and solar plus storage. We draw synergies between design for manufacturing and assembly, process optimization, retrofit technologies, and digitization. Our goal is to influence the improvement and production of buildings to increase performance, enhance energy efficiency, and reduce GHG emissions. This guide documents the research and development efforts initiated by a set of design objectives to "modularize" a set of energy efficiency and low-carbon strategies into a housing unit while preserving and enhancing energy efficiency benefits and decarbonization pathways. This guide is intended to serve as a framework for housing developers, housing agencies, architects, energy experts, and process engineers or factory operator personnel who are critical to today's modular builder teams. This guide focuses on specific energy efficiency strategies, decarbonization pathways, and associated processes as part of NREL's research efforts. Stakeholders may substitute other means, methods, and technologies for the ones evaluated in this study.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Charged-Particle Bound States in Periodic Boxes

We consider the binding energy of a two-body system with a repulsive Coulomb interaction in a finite periodic volume. We define the finite-volume Coulomb potential as the usual Coulomb potential, except that the distance is defined the shortest separation between the two bodies in the periodic volume. We investigate this problem in one and three-dimensional periodic boxes and derive the asymptotic behavior of the volume dependence for bound states with zero angular momentum in terms of Whittaker functions. We benchmark our results against numerical calculations and show how the method can be used to extract asymptotic normalization coefficients for charged-particle bound states. Furthermore, the results we derive here have immediate applications for calculations of atomic nuclei in finite periodic volumes for the case where the leading finite-volume correction is associated with two charged clusters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Anharmonic phonons, thermal expansion, and nuclear quantum effects in Zn

The anharmonic behavior of phonons and thermal expansion of hexagonal zinc were studied from 15 to 690 K by inelastic neutron scattering (INS) and ab initio simulations. Phonon spectra were measured for Q-points covering the full Brillouin zone, giving the phonon density of states (DOS) and dispersions along high-symmetry directions. The dispersions were sharp at 15 K, but diffuse intensity was observed at energies above them. The dispersions broadened with temperature T and the diffuse intensity became stronger. This diffuse intensity appeared in all INS measurements and simulations, except for classical molecular dynamics at 15 K. The temperature-dependent effective potential (TDEP) method, which included the nuclear quantum effect from zero-point vibrational dynamics, was used to calculate the free energy and thermal expansion. For T < 100 K nuclear quantum effects were important for obtaining the correct negative thermal expansion, and path integral molecular dynamics (PIMD) was particularly effective for obtaining the negative thermal expansion in the basal plane. A Heisenberg-Langevin model for interacting phonons coupled to a thermal bath was able to reproduce the shape and intensity of the diffuse spectral features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗