Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RANDOM PHASE APPROXIMATION”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Investigation of the isoscalar monopole response in the proton-rich nucleus 14 O

Deuteron inelastic scattering on 14 O was measured in inverse kinematics using an active-target time projection chamber and a magnetic spectrograph. The experimental technique enabled precise measurements of deuteron recoiling particles in coincidence with beam-like fragments detected in the spectrograph focal plane. The double differential cross section was reconstructed for scattering angles of 3–6 degrees and excitation energies up to 26 MeV. The monopole strength distribution was obtained from the data using a multipole decomposition analysis. The results were compared to quasiparticle random-phase approximation (QRPA) and generator coordinate method (GCM) calculations. The QRPA calculation accurately describes experimental data in the energy range of 13 to 26 MeV. GCM calculations assuming a 12 C (g . s .) + p + p cluster configuration were used to determine the 0 + strength in 14 O below 13 MeV. The monopole transition strength of these cluster states provides a good description of the experimental distribution in the 9–11 MeV region, while the 0$^{+}_{2}$ state accounts for only a small fraction of the experimental strength around 6 MeV.

Active target↗
DOE OSTI2022

Redox Defect Thermochemistry of FeAl 2 O 4 Hercynite in Water Splitting from First-Principles Methods

Solar thermochemical hydrogen (STCH) production is a promising route to produce fuels from sunlight via high-temperature water splitting. However, efficient and technologically viable implementations of this process only allow a narrow window of thermodynamic boundary conditions that can be used to cycle the system, thus limiting the design space for suitable metal oxide redox mediators. An oxygen defect redox mechanism can contribute a favorable reduction entropy to expand this window, and computational evaluation of materials with high oxygen defect entropies could play a pivotal role in guiding the discovery and design of suitable oxides. This study employs first-principles calculations to investigate the redox mediating defect mechanism of the STCH candidate material, hercynite (FeAl 2 O 4 ). We compare the results of total energy calculations from density functional theory (DFT) with beyond-DFT approaches, including hybrid functionals and the random phase approximation, which are among the most advanced methods currently feasible for supercell defect calculations. Using the predicted formation energies, we perform thermodynamic modeling of FeAl2O4 reduction and oxidation via free energy minimization that incorporates ideal gas, configurational, and vibrational entropy contributions evaluated within the quasi-harmonic approximation. Special attention is devoted to understanding interactions among co-existing defects, such as the association of pairs and complexes of O vacancies and cation antisite defects, and the effect of mutually compensating defect charges. Our results corroborate the notion that the details of defect interactions can be decisive for the viability of hydrogen production within the desirable STCH process window.

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

Quantum Embedding Theory for Strongly Correlated States in Materials

Quantum embedding theories are promising approaches to investigate strongly correlated electronic states of active regions of large-scale molecular or condensed systems. Notable examples are spin defects in semiconductors and insulators. We present a detailed derivation of a quantum embedding theory recently introduced, which is based on the definition of effective Hamiltonians. The effect of the environment on a chosen active space is accounted for through screened Coulomb interactions evaluated using density functional theory. Importantly, the random phase approximation is not required, and the evaluation of virtual electronic orbitals is circumvented with algorithms previously developed in the context of calculations based on many-body perturbation theory. In addition, we generalize the quantum embedding theory to active spaces composed of orbitals that are not eigenstates of Kohn–Sham Hamiltonians. Finally, we report results for spin defects in semiconductors.

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

SCF Framework, HF Stability, and RPA Correlation for Jordan–Wigner-Transformed Spin Hamiltonians on Arbitrary Coupling Topologies

Mapping spins to fermions via the Jordan–Wigner (JW) transformation can render mean-field (Hartree–Fock, HF) descriptions effective for strongly correlated spin systems. As established in recent work, the application of such approaches is not limited by the nonlocal structure of JW strings or by site ordering because string operators can be absorbed into Thouless rotations of a Slater determinant, and the variational optimization of a unitary Lie-algebraic similarity transformation removes any ordering dependence. Leveraging these ideas, we develop a self-consistent field (SCF) scheme that expresses the mean-field energy as a functional of the single-particle density matrix, providing an alternative to gradient-based optimization of Thouless parameters. We derive the analytical orbital Hessian to diagnose HF stability and compute the ground-state correlation energy through the random-phase approximation (RPA). Benchmark results for the XXZ and J 1 –J 2 model on one- and two-dimensional lattices demonstrate that RPA significantly improves mean-field accuracy.

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

Metastable Metallic Phase of a Bilayer Blue Phosphorene Induced by Interlayer Bonding and Intralayer Charge Redistributions

We have carried out diffusion Monte Carlo calculations for an A 1 B -1 -stacked bilayer blue phosphorene to find that it undergoes a semiconductor-metal transition as the interlayer distance decreases. While the most stable bilayer structure is a semiconducting one with two monolayers coupled through weak van der Waals interaction, the metallic bilayer at a shorter interlayer distance is found to be only metastable. This is in contrast to a recent theoretical prediction based on random phase approximation that the metallic phase would be the most stable bilayer configuration of blue phosphorene. Our analysis of charge density distributions reveals that the metastable metallic phase is induced by interlayer chemical bonding and intralayer charge redistributions. This study enriches our understanding of interlayer binding of a blue phosphorene and contributes to establishment of correct energetic order between its different phases, which will be essential in devising an experimental pathway for a metallic phosphorene.

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

High-Performance Semiempirical Excited-State Molecular Dynamics Powered by Graphics Processing Units

Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.

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

Scaling Theory of Neutral Sequence-Specific Polyampholytes

The primary sequence of charged monomers in a polyampholyte has a profound effect on its conformational properties. In this work, theory and simulations are used to predict how the sequence influences the single-chain behavior of globally neutral polyampholytes (PAs) under salt-free conditions. We consider PAs with Markov statistics of charges, where charge “blockiness” is defined by the correlation λ along the chain. These PAs cover a wide spectrum of primary sequences ranging from alternating (λ = –1) to ideally random (λ = 0) to diblock (λ → 1). In a theta solvent, sufficiently long PAs of any primary sequence form globules, but their internal structure, as well as the position and sharpness of the coil–globule transition, crucially depends on λ. Further, using scaling arguments, we demonstrate that the excess charge of the globule correlation volume (blob) increases with increasing λ, thereby inducing stronger Coulomb attractions between neighboring blobs of opposite charge and leading to higher density PA globules. Scaling predictions are supported by the results of the random phase approximation (RPA) and coarse-grained molecular dynamics simulations. A generalization to good and poor solvents is presented that enables the construction of the resulting scaling diagram of PA conformations in terms of solvent quality and charge blockiness, thereby providing a comprehensive molecular understanding of the role of the PA primary sequence in the chain structure.

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

Electrostatic Interactions Control the Nanostructure of Conjugated Polyelectrolyte–Polymeric Ionic Liquid Blends

Polyelectrolyte complexation offers unique opportunities to compatibilize polymers with very different backbone chemistries and to control the morphology of the resulting blend via electrostatic manipulation. In this study, we demonstrate the ability to formulate homogeneous complexes of a conjugated polyelectrolyte with a polymeric ionic liquid, utilizing the electrostatic attraction among their oppositely charged side chains. Variation of electrostatic parameters, such as counterion concentration or polymer charge fraction, tunes the morphology of these polymer complexes from homogeneously disordered blend to weakly structured microemulsion where the local ordering arises from backbone-immiscibility-induced microphase segregation. Our experimental observations are in qualitative agreement with both field-theoretic simulation and random-phase approximation calculations. Simulated morphology snapshots suggest and experimental evidence also indicates that the microphase-segregated complex likely takes on a cocontinuous microemulsion structure. Our findings show that ionic interactions are an effective pathway to compatibilize polymers at macroscopic length scales while achieving controlled nanostructures in these ionic blends. Such systems have great potential for engineering the nanostructure of polymers to tailor applications such as nanofiltration, catalysis, and energy storage, where local ordering can enhance the physical properties of an otherwise macroscopically homogeneous structure.

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

The Influence of Charge Correlation and Ion Solvation on the Phase Behavior of Single-Ion Conducting Polymer Blend Electrolytes Using SAXS/SANS

Single-ion conducting polymer blends (SICPBs) have demonstrated exceptional electrochemical performance as solid-state battery electrolytes; however, their nanoscale morphology and thermodynamic behavior remain unexplored. In this work, we investigate blends composed of deuterated poly(ethylene oxide) and poly[lithium sulfonyl(trifluoromethane sulfonyl)imide methacrylate], dPEO/P(LiMTFSI), and report the first experimental study of the nanostructures of charge-neutral polymer blends using small-angle neutron scattering (SANS) and small-angle X-ray scattering (SAXS). Despite the macroscopic miscibility indicated by a single glass-transition temperature, SANS and SAXS results reveal disordered, charge-correlated nanostructures that are strongly influenced by blend composition and temperature. At low concentrations of charge polymer, the scattering is dominated by concentration fluctuations, and the random phase approximation is applied to extract values of the Flory–Huggins interaction parameter, χ SC . At higher charged polymer content, concentration fluctuations are suppressed, and a correlation model is used to characterize the nanostructures of the charge correlations. We find that the structures of the charge correlations are highly dependent on blend composition─consistent with predictions from Sing’s self-consistent field theory-liquid state models. Understanding these features is essential for uncovering the ion transport mechanism that leads to improved electrochemical performance previously reported in SICPB systems.

25 ENERGY STORAGE↗
DOE OSTI2025

Data-Efficient Methods for Determining Flory–Huggins χ Parameters in Multicomponent Polymer Formulations

Polymer formulations are essential in diverse applications including personal care products, coatings, paints, adhesives, and plastic materials. Designing these formulations requires navigating large, complex design spaces, where phase and self-assembly behavior critically impact performance. The Flory–Huggins χ parameter, which quantifies segmental miscibility, is widely used to parametrize the excess free energy of mixing in formulation models. In this work, we introduce two data-efficient, top-down methods for estimating χ parameters using the Random Phase Approximation (RPA): (i) Boundary Nonlinear Regression (Boundary-NLR), which fits theoretical spinodal boundaries to experimental phase boundaries, and (ii) Surrogate Model Inverse Parameter Estimation (SMIPE), which uses a Gaussian Process Classifier to fit sparse phase maps via a surrogate model. Both methods allow rapid parametrization of polymer field-theoretic models without the need for additional experiments. We evaluate these approaches on data sets involving polymer–solvent–nonsolvent ternary mixtures and block copolymer–solvent systems, demonstrating their robustness to experimental noise and their relevance for real-world formulation design.

copolymers↗
DOE OSTI2026

Mechanisms of Metal Additive-Induced Ordering During SNIPS Membrane Formation

Isoporous membranes can be fabricated by combining self-assembly with nonsolvent induced phase separation (SNIPS) using an amphiphilic block copolymer like polystyrene-b-poly(4-vinylpyridine) (SV). Poly(4-vinylpyridine) (V) is known to complex with metal salts, which are hypothesized to stabilize solution ordering and preserve structure during casting. We explored how the molar ratio of metal additive to the poly(4-vinylpyridine) block affected the final membrane morphology via scanning electron microscopy (SEM). Dynamic light scattering (DLS), small-angle X-ray scattering (SAXS), and in situ grazing-incidence SAXS were used to track changes in solution ordering and chain conformation as a function of the molar ratio of the additive to the V block. Additives induced aggregation, promoted the formation of more compact conformations in solution, and facilitated micelle ordering onto lattices at optimal ratios. Furthermore, these experimental results were supported by random phase approximation calculations, which helped explain how the thermodynamic order–disorder transition shifts with additive binding strength. Stronger additive–polymer interactions reduced the block copolymer volume fraction required for ordering in solution, allowing ordered domains to form at lower polymer concentrations.

Additives↗
DOE OSTI2021

Direct Structural Evidence for Interfacial Gradients in Asymmetric Polymer Nanocomposite Blends

Understanding the complex structure of polymer blends filled with nanoparticles (NPs) is key to design their macroscopic properties. In this work, the spatial distribution of hydrogenated (H) and deuterated (D) polymer chains asymmetric in mass is studied by small-angle neutron scattering. Depending on the chain mass, a qualitatively new large-scale organization of poly(vinyl acetate) chains beyond the random-phase approximation is evidenced in nanocomposites with attractive polymer–silica interactions. The silica is found to systematically induce bulk segregation. Only with long H-chains, a strong scattering signature is observed in the q range of the NP size: it is the sign of interfacial isotopic enrichment, that is, of contrasted polymer shells close to the NP surface. A quantitative model describing both the bulk segregation and the interfacial gradient (over ca. 10–20 nm depending on the NP size) is developed, showing that both are of comparable strength. In all cases, NP surfaces trap the polymer blend in a non-equilibrium state, with preferential adsorption around NPs only if the chain length and isotopic preference toward the surface combine their entropic and enthalpic driving forces. This structural evidence for interfacial polymer gradients will open the road for quantitative understanding of the dynamics of many-chain nanocomposite systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structural phase transition, s ± -wave pairing, and magnetic stripe order in bilayered superconductor La 3 Ni 2 O 7 under pressure

Motivated by the recently discovered high-T c superconductor La 3 Ni 2 O 7 , we comprehensively study this system using density functional theory and random phase approximation calculations. At low pressures, the Amam phase is stable, containing the Y 2- mode distortion from the Fmmm phase, while the Fmmm phase is unstable. Because of small differences in enthalpy and a considerable Y 2- mode amplitude, the two phases may coexist in the range between 10.6 and 14 GPa, beyond which the Fmmm phase dominates. In addition, the magnetic stripe-type spin order with wavevector (π, 0) was stable at the intermediate region. Pairing is induced in the s ± -wave channel due to partial nesting between the M = (π, π) centered pockets and portions of the Fermi surface centered at the X = (π, 0) and Y = (0, π) points. This resembles results for iron-based superconductors but has a fundamental difference with iron pnictides and selenides. Moreover, our present efforts also suggest La 3 Ni 2 O 7 is qualitatively different from infinite-layer nickelates and cuprate superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-energy magnetic excitations from heavy quasiparticles in CeCu 2 Si 2

Magnetic fluctuations is the leading candidate for pairing in cuprate, iron-based, and heavy fermion superconductors. This view is challenged by the recent discovery of nodeless superconductivity in CeCu 2 Si 2 , and calls for a detailed understanding of the corresponding magnetic fluctuations. Here, we mapped out the magnetic excitations in superconducting (S-type) CeCu 2 Si 2 using inelastic neutron scattering, finding a strongly asymmetric dispersion for E ≲ 1.5 meV, which at higher energies evolves into broad columnar magnetic excitations that extend to E ≳ 5 meV. While low-energy magnetic excitations exhibit marked three-dimensional characteristics, the high-energy magnetic excitations in CeCu 2 Si 2 are almost two-dimensional, reminiscent of paramagnons found in cuprate and iron-based superconductors. By comparing our experimental findings with calculations in the random-phase approximation, we find that the magnetic excitations in CeCu 2 Si 2 arise from quasiparticles associated with its heavy electron band, which are also responsible for superconductivity. Our results provide a basis for understanding magnetism and superconductivity in CeCu 2 Si 2 , and demonstrate the utility of neutron scattering in probing band renormalization in heavy fermion metals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Predicting nucleon-nucleus scattering observables using nuclear structure theory

Developing a predictive capability for inelastic scattering will find applications in multiple areas. Experimental data for neutron-nucleus inelastic scattering is limited and thus one needs a robust theoretical framework to complement it. Charged-particle inelastic scattering can be used as a surrogate for ( n, γ ) reactions to predict capture cross sections for unstable nuclei. Our work uses microscopic nuclear structure calculations for spherical nuclei to obtain nucleon-nucleus scattering potentials and calculate cross sections for these processes. We implement the Jeukenne, Lejeune, Mahaux (JLM) semi-microscopic folding approach, where the medium effects on nuclear interaction are parameterized in nuclear matter to obtain the nucleon-nucleon ( NN ) interaction in a medium at positive energies. We solve for the nuclear ground state using the Hartree-Fock-Bogliubov (HFB) many-body method, assuming the nucleons within the nucleus interact via the Gogny-D1M potential. The vibrational excited states of the target nucleus are calculated using the quasi-particle random phase approximation (QRPA). We demonstrate our approach for spherical nuclei in the medium-mass region, showing scattering results for the 90 Zr nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Microscopic calculations with noniterative finite amplitude methods and the application to neutron radiative captures and inelastic scatterings

We derive the fully self-consistent quasiparticle random-phase approximation (QRPA) equations with noniterative finite amplitude methods and calculate the transition strengths of giant resonances. Then, we apply the QRPA results to both neutron radiative capture calculations based on the statistical Hauser-Feshbach theory and inelastic scattering calculations based on distorted-wave Born approximation (DWBA). We compare the calculated results with available experimental data and demonstrate how our approach can reproduce giant resonances and various nuclear reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Electric dipole excitations near the neutron separation energies in 96 Mo

Electric dipole strength near the neutron separation energy significantly impacts nuclear structure properties and astrophysical scenarios. These excitations are complex in nature and may involve the so-called pygmy dipole resonance (PDR). Transition densities play a crucial role in understanding the nature of nuclear excited states, including collective excitations, as well as in constructing transition potentials in DWBA or coupled-channels equations. In this work, we focus on electric dipole excitations in spherical molybdenum isotopes, particularly 96 Mo, employing fully consistent Hartree-Fock-Bogoliubov (HFB) and Quasiparticle Random Phase Approximation (QRPA) methods. We analyze the dipole strength near the neutron separation energy, which represents the threshold for neutron capture processes, and examine the isospin characteristics of PDR states through transition density calculations. Examination of proton and neutron transition densities reveals distinctive features of each dipole state, indicating their isoscalar and isovector nature. We observe that the primary component in the enhanced low-energy region exhibits isovector character. The PDR displays a mixture of isoscalar and isovector nature, distinguishing it from the isovector giant dipole resonance (IVGDR). These findings lay the groundwork for future investigations into the role of transition densities in reaction models and for their application to inelastic scattering calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Cooper minima in high-Z atoms: effects of correlation and relativity on np photoionization

Abstract Photoionization dipole transition matrix elements pass through a zero or attain a minimum that leaves imprints on photoionization parameters like the cross-section, angular distribution asymmetry parameter, phase shift, and photoionization time delay. This minimum is commonly known as the ‘Cooper minimum’ (CM). The CM, in general, is strongly affected by relativistic and correlation effects. Previous works investigated CM in the 6pand 5psubshell photoionization up toZ= 100 using the single-particle Dirac-Slater (DS) method. The present work extends the earlier work toZup to 120 using more accurate methods; Dirac–Hartree–Fock (DHF) which includes the relativistic effects and exchange correlations, and the relativistic random phase approximation (RRPA) which includes both initial and final state electron-electron correlations along with relativistic effects. In addition to the study of photoionization from the 6pand 5psubshells, the 4psubshell has also been investigated in the present work. To demonstrate the prominent effects in the high-Z atoms, Rn (Z= 86), Ra (Z= 88), No (Z= 102), Cn (Z = 112), Og (Z= 118), and Ubn (Z= 120) are investigated.

Optics↗