Search NASA⌕ Search

SEARCH · Search NASA

Results for “evolution”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Evolution of highly multimodal Rayleigh–Taylor instabilities

Rayleigh–Taylor (RT) instabilities are important fluid instabilities that arise in inertial confinement fusion (ICF) capsule implosions, and many other contexts. Multi-mode coupling is observed in experiments and plays a substantial role in material mix from RT instabilities. In this work, we study the evolution of highly multimodal perturbations (power law distribution) that approximate those found at manufactured material interfaces. We use simulations of over 2000 different perturbations in the LANL code xRAGE to identify distinct phases in the processes of bubble growth and bubble merger which can be visualized in a 2D phase portrait with clear regimes of mode growth and decay. Our results show that the dynamic evolution of the instability strongly depends on the mode of the perturbations and mode interactions. The merger process accelerates bubble growth. A non-Markovian region and a transition of the instability from: (1) initial exponential growth to (2) linear growth and to (3) quadratic growth and asymptotic behavior, are clearly captured in the phase space. We have developed a quantitative model of bubble growth that reproduces the dynamic behavior of ensembles of perturbations. Implications for ICF capsules designed for robustness against instabilities are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Analytical simulation of meander morphology from equilibrium to long-term evolution: Impacts of channel geometry and vegetation-induced coarsening

This study introduces an innovative approach to modeling meandering river morphology, integrating and investigating the effects of geometric characteristics and vegetation-induced channel coarsening. The developed comprehensive framework combines several advanced techniques: Genetic Programming for refining the scour factor of transverse bed slope, a Leaf Area Index (LAI)-enhanced analytical model for quantifying vegetative flow resistance, and an upstream-weighted moving average method for efficient approximation of the convolution integral in meander migration calculations. The model is validated against both an idealized Kinoshita meander and a natural bend of the Tumen River (China) in equilibrium, demonstrating its robustness across diverse scales and conditions. The model's ability to simulate the long-term evolution, including cutoff events, provides valuable insight for river management strategies. The current findings demonstrate that channel geometry, particularly width-to-depth ratio, plays a dominant role in meander evolution, with wider channels prone to more complex and rapid morphological changes. Vegetation effects are most pronounced in channels with moderate width-to-depth ratios, where they can significantly influence migration rates and bed topography. A combination of channel widening and deepening, coupled with strategic vegetation management, can effectively enhance navigability while maintaining channel stability in the studied Tumen River reach. Sensitivity analyses highlight the complex interplay between hydraulic conditions, sediment characteristics, and vegetation in shaping river morphology. This research advances understanding of the multifaceted nature of meandering river systems and offers practical tools for informed decision-making in river engineering and environmental management, particularly in the context of climate change and increasing anthropogenic pressures on fluvial ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantification of grain boundary effects on the geometrically necessary dislocation density evolution and strain hardening of polycrystalline Mg 4Al using in situ tensile testing in scanning electron microscope and HR-EBSD

In situ tensile testing in a scanning electron microscope (SEM) in conjunction with high-resolution electron backscatter diffraction (HR-EBSD) under load was used to characterize the evolution of geometrically necessary dislocation (GND) densities at individual grain boundaries as a function of applied strain in a polycrystalline Mg—4Al alloy. The increase in GND density was investigated at plastic strains of 0 %, 0.6 %, 2.2 %, 3.3 % from the area including 76 grains and correlated with (i) geometric compatibility between slip systems across grain boundaries, and (ii) plastic incompatibility. We develop expressions for the grain boundary GND density evolution as a function of plastic strain and plastic incompatibility, from which uniaxial tensile stress-strain response of polycrystalline Mg—4Al are computed and compared with experimental measurement. The findings in this study contribute to understanding the mechanisms governing the strain hardening response of single-phase polycrystalline alloys and more reliable prediction of mechanical behaviors in diverse microstructures.

Geometrically necessary dislocations↗
DOE OSTI2025

Twin nucleation and growth in hexagonal close-packed metals: The role of slip-mediated plasticity on twin embryo formation and evolution

Twinning is a key deformation mechanism in hexagonal close-packed (hcp) metals, which are typified by a lack of easily activated slip systems that can accommodate a general state of loading. Pragmatically, the nucleation and evolution of twin domains occurs concomitantly with slip, such that the eventual twin network is conditioned by both external and internal stresses resulting, among others, from the evolution of dislocations. However our understanding of the interplay between dislocations, and twin nucleation and stability remains limited. This work focuses on elucidating the influence of dislocation-mediated plasticity on the formation, growth, and stability of twin embryos. First, a new mesoscale spectral crystal plasticity-twinning framework is developed and used to quantify the change in the free energy landscape following the nucleation {$10\bar{1}2$} of twins in Mg. Representative twin morphologies are modeled under conditions of limited and profuse slip activity, which are emulative of small- and bulk-scale samples, respectively. Then, the driving traction profiles around twin embryos are investigated via a sharp interface approach to obtain insights into how concurrent slip-mediated plasticity can influence the growth/stabilization of nanometric twins. The driving traction profiles are further utilized to determine the stability of twin embryos post loading. The initial dislocation density in the samples, and within the different domains (i.e., twin vs. parent), is seen to have a significant effect on the twin nucleation stress. Namely, the activation of high levels of concomitant plasticity in the parent grain is seen to significantly drive the formation of nanometric twin nuclei at stresses as low as 250 MPa. Further, the sharp interface analysis reveals that profuse plasticity in the parent grain simultaneously alters the forward and back stresses, such that the magnitude/polarity of the driving tractions become increasingly favorable for nanometric twin growth when slip is active. Finally, prior plasticity in the parent grain is seen to result in favorable driving tractions for nanometric twin growth, even at applied stresses as low as ~ 100 MPa. In conclusion, these results are in stark contrast to a case without any dislocations, wherein applied stresses as high as ~ 500 MPa are necessary to grow the twin domains.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effects of Zn and Ca on the deformation texture evolution of Mg–Zn–Ca alloys at elevated temperatures

This work utilized crystal plasticity simulations and experimental thermo-mechanical processing to examine how Zn and Ca content affects texture evolution in Mg–Zn–Ca alloys during simulated hot rolling. Four different compositions were studied: unalloyed Mg and three Mg alloys of ZX0p50 (Mg- 0.5 wt% Zn- 0.1 wt% Ca), ZX30 (Mg-3 wt% Zn- 0.1 wt% Ca), and ZX31 (Mg-3 wt% Zn- 0.3 wt% Ca). Multi-pass Gleeble experiments simulated hot rolling using plane strain conditions. Unalloyed Mg developed a strong basal texture, while the addition of Zn and Ca significantly weakened the texture in ZX30 and ZX31; this reduction was not seen in ZX0p50. Notably, ZX31 exhibited a split basal texture aligned with the rolling direction. Crystal plasticity simulations revealed that the weaker basal texture in ZX30 and ZX31 resulted from reduced activation of basal and twinning modes compared to unalloyed Mg and Zx0p50. Increased Zn and Ca content raised the ratios of basal to pyramidal critical resolved sheer stress (CRSS) and twin to pyramidal CRSS, enhancing pyramidal slip activity and weakening the basal texture. Extension twinning played a critical role in texture development during the first deformation pass. In unalloyed Mg and ZX0p50, twinning led to a strong basal texture, while ZX30 and ZX31 had weak basal textures that increased only slightly with further deformation. The split basal poles simulated in ZX31 were consistent with the experimental findings, highlighting the interplay between alloy composition and texture evolution.

Crystal plasticity finite element↗
DOE OSTI2025

Residual strain and microstructure evolution in 316L stainless steel fabricated by hybrid additive/subtractive manufacturing

The extreme thermal gradients in additive manufacturing (AM) result in residual strains. Subsequent machining operations can relieve the AM strains but can introduce additional strains. AM and machining on the same platform, also called interleaving, during hybrid manufacturing has been shown to improve productivity. However, the impacts of these operations on the resulting residual strains are not well understood. In this investigation, neutron diffraction was used to measure residual strains in hybrid manufactured 316L stainless steel where blown-powder laser directed energy deposition (DED), and machining operations were interleaved in 1-inch intervals. We observe that the dominant texture developed during AM, resulting from different spot sizes and deposition parameters, significantly affects residual strain evolution. A dominant <011> texture along the build direction does not show any strain localization during interleaved deposition and machining whereas a dominant <001> texture shows strong tensile strain localization consistent with the pauses for machining at the end of AM step. Using coolant during machining can further alter the residual strain evolution. Our work shows that a combination of AM parameters and machining strategy can be used to influence the nature and magnitude of residual strains in components fabricated using convergent manufacturing technologies.

Hybrid manufacturing↗
DOE OSTI2025

Influence of implantation temperature and He implantation-induced defects on morphological evolution of co-deposited Cu-Mo nanocomposites

Here, we investigate the effect of high-temperature helium (He) implantation on microstructural evolution in physical-vapor-co-deposited nanocomposite thin films of copper (Cu) and molybdenum (Mo). The microstructure morphologies of He-implanted and He-free domains are characterized using transmission electron microscopy and statistical analysis. High implantation temperatures (500°C and 750°C) lead to coarsening of Cu and Mo domains and their eventual reorientation. The microstructure evolution in He-implanted and He-free domains is comparable, indicating that implantation-induced defects do not accelerate the coarsening of the nanocomposite as compared to annealing alone. This observation contrasts with previously reported effects of implantation-induced defects on single-phase nanocrystalline metals, which include enhancement of grain growth by increasing self-diffusivity or its inhibition by pinning of grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Porosity, swelling, and composition evolution in high-burnup monolithic U-Mo fuel

The microstructural progression of very high-burnup (>8 × 10 21 fissions/cm 3 ) monolithic uranium-molybdenum (U-10wt.%Mo) was analyzed, providing crucial insights into the behavior of post-recrystallized nuclear fuel, where scant data exists. Three focused ion beam cuboids sourced from a fuel plate with varying local burnups of 8.86 × 10 21 , 9.05 × 10 21 , and 9.36 × 10 21 fissions/cm 3 were characterized. The porosity and composition of the samples were evaluated to characterize the evolution of the microstructure as a function of fission density and different locations on the fuel plate, while simultaneously isolating plate-specific parameters such as Zr diffusion barrier thickness, hot-isostatic-pressing conditions, enrichment, and reactor conditions. The porosity was segmented, and the three-dimensional distribution of the porosity was extrapolated from the two-dimensional segmentation. The composition was assessed and quantified using energy-dispersive X-ray spectroscopy areal mapping. The porosity fraction increased as a function of the burnup from 27.77±0.51, 35.12±1.54, and 37.71±0.44 % for 8.86 × 10 21 , 9.05 × 10 21 , and 9.36 × 10 21 fissions/cm 3 , respectively. When compared to literature, the porosity volume fraction plateaus at burnups greater than 6 × 10 21 fissions/cm 3 , while the pore size grows linearly as a function of fission density. The number of large pores increased in number density as a function of burnup, while the smallest pores (<0.3 µm) increased up to 9.05 × 10 21 fissions/cm 3 , followed by a decrease at 9.36 × 10 21 fissions/cm 3 . The delamination and cracking in the fuel plate propagated through an interconnected porosity sublayer identified ∼5 µm from the diffusion barrier. The local swelling of the specimens was within or near the prediction bounds of the Robinson-Williams model for local swelling. The fission products, strontium, barium, cerium, and cesium, precipitated into the pores, while neodymium accumulated adjacent to the pores. Furthermore, these findings have direct implications for the development of fuel performance codes and the accurate documentation of the microstructure evolution in high burnup U-Mo, thus enhancing the safety and efficiency of nuclear fuel usage.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Helium co-implantation effects on dislocation loop evolution during in situ ion irradiations of a novel additive-manufactured Fe-9Cr steel alloy

This study examines the impact of helium injection rate on dislocation loop evolution during ion irradiation of a novel, additively manufactured (AM) 9Cr reduced-activation ferritic/martensitic (RAFM) steel. Through in situ transmission electron microscopy ion irradiations and post-irradiation characterization, the effects of simultaneous 1 MeV Kr²⁺ and 12 keV He⁺ ion irradiation at 400°C to 5 dpa with varying helium levels were examined. The findings reveal that helium co-injection significantly increases dislocation loop density, primarily through the formation of small defect features exhibiting contrast and projected morphologies consistent with a⟨100⟩ loops (∼1 nm), while suppressing larger a/2⟨111⟩ loops. The findings underscore the complex interplay among continuously increasing helium content, temperature, and alloy composition in the evolution of dislocation loops under irradiation in RAFM steels. The study highlights the need for further investigation into helium-mediated mechanisms, such as proposed interactions between self-interstitial atoms and helium atoms, to enhance the performance of RAFM steels in advanced nuclear reactors.

ANA2↗
DOE OSTI2025

Amorphous electrocatalysts for oxygen and hydrogen evolution reactions: Advances in hydrogen production

The electrochemical splitting of water into oxygen and hydrogen is fundamental for renewable energy storage and conversion. The development of cost-effective and highly efficient electrocatalysts remains essential for industrial-scale implementation of this technology. Recent advances have highlighted the superior activity, stability and structural adaptability of amorphous electrocatalysts compared to their crystalline counterparts. This review critically examines synthesis strategies, characterisation techniques, and the electrochemical performance of amorphous materials for both oxygen evolution (OER) and hydrogen evolution (HER) reactions. Key factors influencing catalytic efficiency, including electronic structure and surface chemistry, are discussed in detail and contextualised with established literature. The review also highlights the critical role of enthalpic contributions in governing reaction energetics and catalyst performance, which aids in understanding and optimising electrocatalytic efficiency. Notably, ongoing research continues to reveal that amorphous catalysts consistently deliver improved performance in water-splitting applications, highlighting their growing relevance in electrocatalysis. The rationale for employing amorphous catalysts in water splitting is articulated, emphasising their unique advantages. By integrating recent findings and outlining future research directions, this review underscores the pivotal role of amorphous materials in advancing sustainable hydrogen production and identifies promising avenues for catalyst innovation.

Amorphous catalysts↗
DOE OSTI2025

Design of a separate effects MiniFuel irradiation experiment investigating microstructure evolution in high burnup UO 2

The microstructural evolution of UO 2 fuel pellets during commercial operation in light water reactors (LWRs) is known to vary significantly across the pellet radius due to spatial variations in local temperature and burnup. The primary obstacle to extending LWR refueling cycles to 24-month intervals is the susceptibility of certain high burnup fuel microstructures to fuel fragmentation, relocation, and dispersal (FFRD) during a loss of coolant accident (LOCA). Although FFRD of the high burnup structure in the rim region of a pellet is well studied, the fine fragmentation that has been observed in a second region, near the midradius of the pellet (termed the “dark zone”) following mock LOCA testing of high burnup commercial fuel rods is less understood. This paper describes the design, analysis, and execution of a separate effects MiniFuel irradiation experiment that aims to identify the specific temperature and burnup regimes under which FFRD-susceptible dark zone microstructures form. The small disc specimens (3 mm diameter by ∼0.3 mm thick) enable more precise control of the relatively uniform temperature and burnup conditions. A total of 42 specimens were fabricated with typical LWR fuel densities (∼96%–98% of theoretical density) and grain sizes (∼12 μm) and are being irradiated over a range of temperatures (600°C–1000°C) and discharge burnups (50–72 MWd/kg-U) that bound the midradius region of high burnup LWR fuel. Fuel specimens with identical 235 U enrichments were inserted in two irradiation locations in the High Flux Isotope Reactor and are currently undergoing irradiation to further evaluate the impact of rate effects (fission rate, time at temperature) on the microstructural evolution. The fuel fabrication and the thermal and neutronic simulations used for designing the experiment are detailed in this paper. A secondary objective of the experiment is to observe fission gas release (FGR) under the various irradiation conditions, and this work provides first-order predictions of FGR from all fuel specimens. The insights gained from these experiments will inform future high burnup core designs that could minimize the formation of susceptible microstructures and ultimately enable 24-month refueling cycles while minimizing the fraction of the fuel susceptible to FFRD.

FFRD↗
DOE OSTI2024

Evolution of the nuclear spin-orbit splitting explored via the 32 Si( d,p ) 33 Si reaction using SOLARIS

The spin-orbit splitting between neutron 1p orbitals at 33 Si has been deduced using the single-neutron-adding (d,p) reaction in inverse kinematics with a beam of 32 Si, a long-lived radioisotope. Reaction products were analyzed by the newly implemented SOLARIS spectrometer at the reaccelerated-beam facility at the National Superconducting Cyclotron Laboratory. The measurements show reasonable agreement with shell-model calculations that incorporate modern cross-shell interactions, but they contradict the prediction of proton density depletion based on relativistic mean-field theory. The evolution of the neutron 1p-shell orbitals is systematically studied using the present and existing data in the isotonic chains of N = 17, 19, and 21. In each case, a smooth decrease in the separation of the 1p 3/2 - 1p 1/2 orbitals is seen as the respective p-orbitals approach zero binding, suggesting that the finite nuclear potential strongly influences the evolution of nuclear structure in this region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Understanding parton evolution in matter from renormalization group analysis

We perform a renormalization group (RG) analysis of collinear hadron production in deep inelastic scattering on nuclei. We consider the limit where the parent parton energy E is large, while the medium opacity remains small. We identify the fixed order and leading enhanced medium contributions to the semi-inclusive cross sections and derive RG equations that resum multiple emissions near the endpoints of the splitting functions at first order in opacity. These evolution equations treat the same type of radiation enhancement in matter as the modified Dokshitzer-Gribov-Lipatov-Altarelli-Parisi approach, but differ in the way one regulates the collinear divergences. They provide a unique analytic insight into the problem of resummation and a faster and more efficient path to phenomenology. The new RG evolution framework is applied to study fragmentation in eA reactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Foundational insights into the mechanical and molecular evolution of porcine skin gelatin during gelation

Gelatin is a widely used material in biomedical fields, particularly in regenerative medicine and tissue engineering, due to its biocompatibility and versatile properties. While prior research has explored methods to enhance gelatin's mechanical strength and stability, fundamental studies on gelatin, specifically its curing process, mechanical stiffness, and chemical evolution during gelation, remain limited. This study uses ultrasonic testing and Fourier Transform Infrared Spectroscopy (FTIR) to examine gelatin's stiffness and molecular changes during gelation. Samples of 175 and 300 Porcine Skin Bloom Strength Gelatin at concentrations of 2% and 6% (w/v) were analyzed. Through transmission ultrasonic testing helped identify key transition points in gelation, with higher concentrations exhibiting delayed transitions. FTIR revealed that C-N bond formation peaks early while N-H bond deformation persists. A correlation emerged between sound speed and peak absorbance, suggesting that changes in molecular mobility may contribute to the observed sound speed behavior during periods of active bond formation. However, as gelation continues, fewer bonding components may be available, potentially decreasing molecular movement and contributing to the observed increase in sound speed. These findings provide insights into gelatin's mechanical and chemical evolution, offering a framework for improved control over its gelation kinetics. Swept-Frequency Acoustic Interferometry (SFAI) was performed at the end of the curing process to measure the sound speed, enabling the calculation of the bulk moduli of the gelatin samples. The combined use of ultrasonic and FTIR testing provides a non-destructive method for characterizing gelatin and other biomaterials. This approach advances understanding of gelatin curing behavior and supports the development of safer biomaterials with tailored mechanical properties for various applications such as tissue engineering and regenerative medicine.

Biomaterials↗
DOE OSTI2026

Thermal Gradient Effects on Redox Evolution and Volatility-Driven Fractionation in Ternary U/Ce/Cs Condensates

Understanding how thermal history influences redox evolution and chemical fractionation is essential for characterizing high-temperature condensation in complex materials, including nuclear debris. Here, we tested the hypothesis that distinct thermal regimes in a plasma flow reactor influence redox pathways and elemental partitioning in ternary U/Ce/Cs systems. A configurable plasma flow reactor was modified with an external tube furnace to impose two distinct thermal gradients: continuous ambient cooling and a furnace-assisted thermal hold-up near 1400 K followed by rapid cooling. Transmission electron microscopy characterized phase identity, morphology, and nanoscale element distributions, while inductively coupled plasma-mass spectrometry quantified bulk elemental ratios. Across both thermal regimes, uranium and cerium condensed as UO 2 and CeO 2 as dominant refractory oxide products. Uranium partially oxidized to α-UO 3 during extended ambient cooling, while furnace-assisted hold-up preserved UO 2 and produced partial reduction of cerium to Ce 2 O 3 . Cesium remained volatile upstream and condensed later in the reactor, forming Cs 2 O and Cs-uranate phases with the highest incorporation after thermal hold-up. Bulk ICP-MS measurements supported these observations. U/Ce ratios remained comparatively stable and Cs displayed delayed and apparent transient enrichment that matched the nanoscale measurements. This integrated approach provides a quantitative method for linking thermal gradients to redox evolution and volatility-driven fractionation. These results show how the plasma flow reactor can identify where equilibrium descriptions remain adequate and where kinetic effects from residence time and temperature history must be considered when interpreting condensation behavior in multicomponent systems.

and nuclear chemistry↗
DOE OSTI2025

Synthesis and Characterization of [Ni(H 2 O)(7-P Ph 2 N ArSO3 ) 2 ](NaBF 4 ) for Light-Driven Quantum Dot-Catalyst Hydrogen Evolution

Light-driven hydrogen generation from water is a route to carbon-neutral fuels. However, the integrated light absorbers and catalysts must be compatible and functional under the same conditions. A sulfate-functionalized complex, [Ni(H 2 O)(7-P Ph 2 N ArSO3 ) 2 ](NaBF 4 ), where P Ph 2 N ArSO3 = 4-(3,6-diphenyl-1,3,6-azadiphosphepan-1-yl)benzenesulfonate), was synthesized and characterized. A solid-state structure from single-crystal X-ray crystallography is also reported. The electrocatalytic activity for hydrogen evolution of this complex in 7:3 H 2 O:CH 3 CN was characterized by cyclic voltammetry and controlled potential electrolysis, with a maximum observed rate of 83 s –1 . Density functional theory was used to characterize the favorable binding of the sulfate functionalities to cadmium sulfide quantum dots. Experimental assembly of these catalysts onto cadmium sulfide quantum dots was successful. In conclusion, both fast electron transfer and light-driven pH-dependent hydrogen evolution were observed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamic Evolution of Mass and Physical Properties of Atmospheric Organic Aerosol under Solar Irradiance

Organic aerosol (OA) particles constitute a substantial fraction of sub-micron particulate mass in the atmosphere and play a critical role in climate system. OA undergoes dynamic aging processes in the atmosphere, with photolytic aging induced by ultraviolet solar irradiance being an important yet poorly characterized mechanism. Knowledge gaps persist regarding the role of volatility transformations during photolytic aging on the OA mass decay kinetics and the evolution of climate-relevant properties, such as hygroscopicity, hindering the model evaluation of OA spatiotemporal distributions and atmospheric budgets. In this study, we conduct isothermal photolytic aging experiments on both laboratory-generated secondary organic aerosols and ambient-collected particles from urban Atlanta, utilizing a high-sensitivity Quartz Crystal Microbalance. Our results reveal that photolytic aging reduces 40–66% of the low-volatility OA mass with lifetimes ranging from 8 to 200 hours under solar irradiance, and 44–92% of the photolytic mass loss is through slow evaporation of semi- or intermediate-volatile products, kinetically limited by their volatility. We observe up to ±50% changes in OA hygroscopicity with the transformation of fresh OA to photo-recalcitrant low-volatility products, associated with changes in oxygen-to-carbon ratio and molecular weight. A kinetic model incorporating photolytic volatility transformation provides a cohesive explanation for the observed photolysis-induced changes in mass, volatility, and hygroscopicity. Our results can help constrain model representation of the dynamic evolutions of mass and climate-relevant properties during photolytic aging processes of the ambient OA, improving our understanding of OA atmospheric behavior and climate impact.

Bai, Bin↗
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗