Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Pseudopotential”

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

Guiding centre motion for particles in a ponderomotive magnetostatic end plug

The Hamiltonian dynamics of a single particle in a rotating plasma column, interacting with an magnetic multipole is perturbatively solved for up to second order, using the method of Lie transformations. First, the exact Hamiltonian is expressed in terms of canonical action-angle variables, and then an approximate integrable Hamiltonian is introduced, using another set of actions and angles, which describe the centre of oscillation for the particle. The perturbation introduces an effective ponderomotive potential, which to leading order is positive. At the second order, the pseudopotential consists of a sum of terms of the Miller form, and can have either sign. Additionally, at second order, the ponderomotive interaction introduces a modification to the particle effective mass, when considering the motion along the column axis. It is found that particles can be axially confined by the ponderomotive potentials, but acquire radial excursions which scale as the confining potential. The radial excursions of the particle along its trajectory are investigated, and a condition for the minimal rotation frequency for which the particle remains radially confined is derived. Last, we comment on the changes to the aforementioned solution to the pseudopotentials and particle trajectory in the case of resonant motion, that is, a motion which has the same periodicity as the perturbation.

Rubin, T. (ORCID:0000000164850096)↗
DOE OSTI2024

Reproducibility of real-time time-dependent density functional theory calculations of electronic stopping power in warm dense matter

Real-time time-dependent density functional theory (TDDFT) is widely considered to be the most accurate available method for calculating electronic stopping powers from first principles, but there have been relatively few assessments of the consistency of its predictions across different implementations. This problem is particularly acute in the warm dense regime, where computational costs are high and experimental validation is rare and resource intensive. We report a comprehensive cross-verification of stopping power calculations in conditions relevant to inertial confinement fusion conducted using four different TDDFT implementations. We find excellent agreement among both the post-processed stopping powers and relevant time-resolved quantities for alpha particles in warm dense hydrogen. We also analyze sensitivities to a wide range of methodological details, including the exchange-correlation model, pseudopotentials, initial conditions, observable from which the stopping power is extracted, averaging procedures, projectile trajectory, and finite-size effects. We show that among these details, pseudopotentials, trajectory-dependence, and finite-size effects have the strongest influence, and we discuss different strategies for controlling the latter two considerations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Bilayer Ion Trap Design for 2D Arrays

Junctions are fundamental elements that support qubit locomotion in two-dimensional ion trap arrays and enhance connectivity in emerging trapped-ion quantum computers. In surface ion traps they have typically been implemented by shaping radio frequency (RF) electrodes in a single plane to minimize the disturbance to the pseudopotential. However, this method introduces issues related to RF lead routing that can increase power dissipation and the likelihood of voltage breakdown. Here, we propose and simulate a novel two-layer junction design incorporating two perpendicularly rotoreflected (rotated, then reflected) linear ion traps. The traps are vertically separated, and create a trapping potential between their respective planes. The orthogonal orientation of the RF electrodes of each trap relative to the other provides perpendicular axes of confinement that can be used to realize transport in two dimensions. While this design introduces manufacturing and operating challenges, as now two separate structures have to be precisely positioned relative to each other in the vertical direction and optical access from the top is obscured, it obviates the need to route RF leads below the top surface of the trap and eliminates the pseudopotential bumps that occur in typical junctions. Here in this paper the stability of idealized ion transfer in the new configuration is demonstrated, both by solving the Mathieu equation analytically to identify the stable regions and by numerically modeling ion dynamics. Our novel junction layout has the potential to enhance the flexibility of microfabricated ion trap control to enable large-scale trapped-ion quantum computing.

42 ENGINEERING↗
DOE OSTI2023

Computation of forces and stresses in solids: Towards accurate structural optimization with auxiliary-field quantum Monte Carlo

The accurate computation of forces and other energy derivatives has been a long-standing challenge for quantum Monte Carlo methods. A number of technical obstacles contribute to this challenge. We discuss how these obstacles can be removed with the auxiliary-field quantum Monte Carlo (AFQMC) approach. AFQMC is a general, high-accuracy, many-body total-energy method for molecules and solids. The implementation of back-propagation for pure estimators allows direct calculation of gradients of the energy via the Hellmann-Feynman theorem. A planewave basis with norm-conserving pseudopotentials is used for the study of periodic bulk materials. Completeness of the planewave basis minimizes the effect of so-called Pulay terms. The ionic pseudopotentials, which can be incorporated in AFQMC in exactly the same manner as in standard independent-electron methods, regulate the force and stress estimators and eliminate any potential divergence of the Monte Carlo variances. The resulting approach allows applications of full geometry optimizations in bulk materials. As a result, it also paves the way for many-body computations of the phonon spectrum in solids.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1978

Adhesion of a bimetallic interface

The Hohenberg-Kohn and Kohn-Sham formalisms are used to examine binding (binding energy as a function of separation) for combinations of the simple metals Al(111), Zn(0001), Mg(0001), and Na(110) in contact. Similar metal contacts between Al, Zn, Mg, and Na are examined self-consistently in an ab initio calculation using the Kohn-Sham formalism. Crystallinity is included using the Aschroft pseudopotential via first order perturbation theory for the electron-ion interaction; and the ion-ion interaction is included exactly via a lattice sum. Binding energy was determined both in the local-density approximation and including gradient corrections to the exchange and correlation energy. Binding was found in all cases. In dissimilar metal contacts, interfacial bonding was greater than that in the weaker material predicting the possibility of metallic transfer. The nonzero position of the energy minimum in like metal contacts is explained in terms of consistency between the Ashcroft pseudopotential and the bulk charge density. Good agreement with experimental surface energies is obtained in the self-consistent calculation when nonlocal terms are included.

Ferrante, J.↗
NASA NTRS1985

Thermodynamic properties of liquid Si and Ge

Thermodynamic properties of liquid Si and Ge are calculated using standard variational techniques, with the hard-sphere system as a reference system. Third-order electron-ion pseudopotential contributions are included nonvariationally and are shown to have little effect on the free energy, although a larger one on the derivatives of the free energy. Two different model pseudopotentials are used and give similar results in good agreement with experiment.

Shih, W.-H.↗
DOE OSTI2023

Relation between the Hydrated Electron Solvation Structure and Its Partial Molar Volume

It is now generally accepted that the hydrated electron occupies a cavity in water, but the size of the cavity and the arrangements of the solvating water molecules are not fully characterized. Here, we use the Kirkwood-Buff (KB) approach to examine how the partial molar volume (V M ) provides insight into these issues. The KB method relates V M to an integral of the electron-water radial distribution function, a key measure of the hydrated electron structure. Here we have applied it to three widely-used pseudopotentials and the results show that V M is a sensitive measure of the fidelity of hydrated electron descriptions. Thus, the measured V M places constraints on the hydrated electron structure that are important in developing and evaluating model descriptions. Importantly, we find that V M does not reflect only the cavity size (and thus should not be used to infer the cavity radius), but is strongly dependent on the extended solvation structure.

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

The Effect of Monomer Size on Fusion and Coupling in Colloidal Quantum Dot Molecules

The fusion step in the formation of colloidal quantum dot molecules, constructed from two core/shell quantum dots, dictates the coupling strength and hence their properties and enriched functionalities compared to monomers. Herein, studying the monomer size effect on fusion and coupling, we observe a linear relation of the fusion temperature with the inverse nanocrystal radius. This trend, similar to that in nanocrystal melting, emphasizes the role of the surface energy. The suggested fusion mechanism involves intraparticle ripening where atoms diffuse to the reactive connecting neck region. Moreover, the effect of monomer size and neck filling on the degree of electronic coupling is studied by combined atomistic-pseudopotential calculations and optical measurements, uncovering strong coupling effects in small QD dimers, leading to significant optical changes. Understanding and controlling the fusion and hence coupling effect allows tailoring the optical properties of these nanoscale structures, with potential applications in photonic and quantum technologies.

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

Modelling Thomson Scattering in a Hydrogen Plasma at Stellar Interior Conditions Using the Hypernetted‐Chain Approach

Under the extreme conditions found in small stars, where electron degeneracy and Coulomb coupling are significant, accurate modeling of Thomson scattering is crucial for determining opacity, a primary quantity for stellar energy transport. We use hypernetted-chain calculations, incorporating quantum pseudopotentials and electron-exchange effects to obtain the electron–electron static structure factor to calculate the Thomson scattering transport cross-section for conditions prevailing in the interior of small stars. These results are compared to those from average-atom simulations and analytical calculations. Our findings support laboratory astrophysics experiments aimed at benchmarking opacity models for stellar interiors, particularly for red dwarf stars, and help to bridge theoretical models with observations.

ICF↗
DOE OSTI2023

Bridging adsorption behavior of confined CH 4 -CO 2 binary mixtures across scales

An accurate understanding of the competitive adsorption of CH 4 -CO 2 binary mixtures in nano-confined systems is critical for engineering CO 2 storage in shale gas reservoirs. Due to difficulties in making reliable experimental observations in nano-scale, atomistic simulations (ASs), such as the Grand Canonical Monte Carlo (GCMC) method, provide a viable approach to studying the adsorption behavior of confined fluids. ASs are, however, limited in the size of the compositional domain due to the high computational cost. This work proposes a framework that combines AS and the lattice Boltzmann (LB) method to bridge the physics of confined fluids across scales. The Peng–Robinson equation of state (PR-EoS) produces fugacity coefficients, which serve as input for conducting multi-component GCMC simulations. These GCMC simulations explore the competitive adsorption behavior of CH 4 -CO 2 in nano-slits at various composition, pressure, and channel-width conditions. Both components generate adsorption layers with high densities near the walls with CO 2 preferentially adsorbing compared to CH 4 on the organic walls of carbon sheets. At the mesoscale, a pseudopotential model represents the intermolecular forces in multi-component, multiple-relaxation-time LB simulations. The LB simulations are in good agreement with the GCMC results, allowing us to obtain values for tunable LB parameters. We then extend the use of LB to simulate adsorption behavior in complex networks with nano-sized channels. The phase behavior and fluid properties in the complex geometries of nano-channels differ from nano-slits and bulk systems. Furthermore, the bridging of physics from GCMC (microscale) to LB (mesoscale) via the macroscale PR-EoS connects the adsorption behavior of binary systems across scales.

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

Bridging confined phase behavior of CH 4 -CO 2 binary systems across scales

Phase behavior of confined fluids may deviate significantly from that of the bulk fluid due to the fluid-wall interactions being a significant portion of all intermolecular interactions under confinement. Despite recent advancements in understanding confined phase behavior of pure fluids, confined phase behavior of mixtures remains an understudied topic. In this work, we examine the confined phase behavior of a CH 4 -CO 2 binary system by combining Monte Carlo (MC) simulations, a cubic equation of state (EoS), and the lattice Boltzmann method (LBM). First, the effects of confinement on density and phase distribution in nano-size pores are established using Gibbs Ensemble MC calculations, which produce precise results of liquid and vapor confined pressures and account for the modification of the phase change location. By comparing the phase envelopes of bulk and confined mixtures at a fixed temperature, here it is observed that the phase envelopes shrink with reductions in pore size. Based on this observation, we extend a modified Peng-Robinson EoS, which was originally developed for pure fluids under confinement, to mixtures via van-der-Waals-type mixing rules and by accounting for shifts in the critical properties of confined CH 4 -CO 2 . The resulting phase envelopes are in good agreement with the MC data. In addition, a local density model is used in combination with the confined EoS to calculate adsorption isotherms of CH 4 -CO 2 mixtures and to characterize the behavior of confined matter in nanopores. Finally, we incorporate this EoS in a multicomponent multiphase LBM that uses a pseudopotential model to represent intermolecular forces. This workflow utilizes multiscale simulation techniques to bridge the behavior of multicomponent systems across scales and to shed light on the confined phase behavior of CH 4 -CO 2 binary systems.

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

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2024

Prediction of Redox Potentials for U, Np, Pu, and Am in Aqueous Solution

The redox properties of the actinides in aqueous solution are important for fuel production/reprocessing and understanding the environmental impact of nuclear waste. The redox potentials for U, Np, Pu, and Am in oxidation states from 0 up to VII (as appropriate) in aqueous solutions have been predicted at the density functional theory level with the B3LYP functional, Stuttgart small core pseudopotential basis sets for the actinides, and explicit (30H 2 O molecules)/implicit treatment of the aqueous solvent using the self-consistent reaction field COSMO and SMD approaches for the implicit solvation. The predictions of the structural parameters of clusters incorporating first and second solvation shells are consistent with the available experimental data., Our results are typically within 0.2 V of the available experimental data using two explicit solvation shells with an implicit solvent model. The use of the PW91 functional substantially improved the prediction of the Pu(VI/V) redox couple. The redox couples for An(VI/IV) and An(V/IV) which involve the addition of protons and removal of the actinyl oxygens led to slightly larger differences from experiment. Here, the An(IV/0) and An(III/0) couples were reliably predicted with our approach. Predictions of the unknown An(II/I) redox potentials were negative, consistent with expectations, and predictions for unknown An(VII/VI), An(III/II), and An(II/0) redox couples improve prior estimates.

Actinides↗
DOE OSTI2024

Prediction of Redox Potentials for Ac, Th, and Pa in Aqueous Solution

Density functional theory in conjunction with small core pseudopotentials and the associated basis sets was used to calculate potentials for multiple redox couples, covering a range of oxidation states for Ac (0 to III), Th (0 to IV), and Pa (0 to V) in aqueous solution. Solvation effects were incorporated using a supermolecule-continuum approach, with 30 water molecules representing two solvation shells, and the COSMO and SMD implicit solvation models. The calculated geometries for Ac(III), Th(IV), and Pa(V) were in reasonable agreement with the available experimental data. Using the COSMO model with the B3LYP functional, the calculated redox potentials were within ± 0.2 V from experiment for most redox couples. Several pathways were explored for the Pa(V/IV) redox couple for different forms of Pa(V) and Pa(IV). Most Pa(V/IV) redox couples have very similar potentials, ranging from 0 to -0.4 V up to a pH of 1.4. At pH = 1.4, the potentials shift to values that are more negative than -0.7 V, reflecting the growing unfavorable nature of the redox process at higher pH levels. The calculated values for An(III/II) potentials were consistent with prior estimates and the available experimental data. The predicted redox potentials for An(II/I) were highly negative, as expected. For An(I/0) potentials, Th and Pa exhibited positive values, contrasting with the negative values calculated for Ac. Furthermore, the An +m /An(0) potentials agreed better with the experimental data when using the COSMO solvation model as compared to the SMD model.

Chemical calculations↗
DOE OSTI2025

Prediction of Redox Potentials for Different Oxidation States of U, Np, Pu, and Am in Alkaline Aqueous Solution

The redox potentials for U, Np, Pu, and Am for oxidation states +III up to +VIII in alkaline aqueous solutions were predicted using density functional theory (DFT) and small-core pseudopotentials and their basis sets, with a hybrid explicit/implicit solvent model using SHE = 4.28 V. For each oxidation state, various oxo/hydroxo complexes were evaluated resulting in a variety of one-electron redox pathways. For An(VIII/VII) couples, the predicted redox potentials for the [An(VIII)O 5 (OH)] -3 /[An(VII)O 4 (OH) 2 ] -3 , or [An(VIII)O 4 (OH) 2 ] -2 /[An(VII)O 4 (OH) 2 ] -3 couples are in good agreement with existing estimates. For An(VII/VI) redox couples, all couples, particularly [An(VII)O 4 (OH) 2 ] -3 /[An(VI)O 2 (OH) 4 ] -2 , were in agreement with experimental values for U, Np, and Pu, but the results for Am showed larger differences from the estimated potentials. The An(VI/V) couples were consistent with experiment for dioxo/tetrahydroxo couples, and the An(V/IV) couples showed acceptable agreement based on actinide-specific couples, with neutral hydroxides often favored in the +IV state. The An(IV/III) couples were consistent with literature values when modeled as soluble neutral hydroxides. The use of our approach yielded calculated redox potentials that were within ±0.2 V from experimental or estimated values consistent with our prior calculations on redox potentials of actinides from Ac to Am in acidic aqueous solutions. Finally, as a result, this supports the robustness of our DFT-based methodology for predicting actinide redox potentials, offering valuable insights for actinide chemistry in aqueous solutions.

Actinides↗
DOE OSTI2025

Prediction of Redox Potentials for the Late Actinides Cm to Lr Using Electronic Structure Methods

Our previously developed computational method for calculating the aqueous redox potentials of the early actinides has been extended to the later elements in the actinide series: Cm, Bk, Cf, Es, Fm, Md, No, and Lr in multiple oxidation states. These calculations were performed using density functional theory with small-core pseudopotentials and their associated basis sets. Solvation effects were considered via a supermolecule-continuum approach, with 30 water molecules representing two solvation shells. Both the COSMO and SMD implicit solvation models were utilized. The structural parameters and hydration numbers for Cm(III), Bk(III), Bk(IV), and Cf(III) are in reasonable agreement with the available experimental data. For redox processes involving atomic cations in solution, the B3LYP/COSMO approach predicted redox potentials to within ±0.2 V of experiment for most redox couples, consistent with our prior work. Inclusion of spin-orbit corrections in specific redox pairs, especially those with the later actinides in high oxidation states, yields improved results relative to calculations including only scalar-relativistic corrections. The An +m /An(0) redox potentials were calculated using a Born-Haber cycle incorporating sublimation, ionization, and hydration energies. Due to a lack of experimental data, three sets of ionization energies were used for the Born-Haber cycle. The calculated An(III/0) potentials showed better agreement with experimental data when using the COSMO solvation model and the test set comprising the NIST recommended ionization energies. Furthermore, the Md(II/0) potential was better described with the SMD model, whereas No(II/0) was not well described by all methods. Finally, the computational approach was able to predict redox potentials that for most cases agreed with the current available experimental or estimated data.

Actinides↗
DOE OSTI2024

Ligand-Induced Size-Dependent Circular Dichroism in Quantum Dots

Recent experiments have probed the chiral properties of semiconductor nanocrystal (NC) quantum dots (QDs), but understanding the circular dichroism line shape, excitonic features, and chirality induction mechanism remains a challenge. We propose an atomistic pseudopotential method to model chiral ligand passivated QDs, computing circular dichroism (CD) spectra for CdSe QDs (2.6-3.8 nm). We find strong agreement between calculated and measured line shapes, predicting consistent bisignate line shapes with decreasing CD magnitude as size increases. Further, our analysis reveals the origin of bisignate line shapes, arising from nondegenerate excitons with opposing angular momenta. We also explore the impact of chiral ligand orientation on QD surfaces, observing changes in the optical activity magnitude and sign. This orientation sensitivity offers the means to distinguish ordered from disordered ligand configurations, facilitating the study of order-disorder transitions at ligand-QD interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗