SEARCH · Engineering Papers

Results for “Density functional perturbation theory”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Phonon-assisted optical absorption of SiC polytypes from first principles

Silicon carbide (SiC) is an indirect-gap semiconductor material widely used in electronic and optoelectronic applications. While experimental measurements of the phonon-assisted absorption coefficient of SiC across its indirect gap have existed for more than 50 years, theoretical investigations of phonon-assisted absorption have been hampered by their excessive computational cost. Here, in this work, we calculate the phonon-assisted temperature-dependent optical absorption spectra of the commonly occurring SiC polytypes (3C, 2H, 4H, 6H, and 15R), using first-principles approaches based on density functional theory and related techniques. We show that our results agree with experimentally determined absorption coefficients in the spectral region between the direct and indirect band gaps. The temperature dependence of the spectra can be well predicted with taking the temperature dependence of the band gaps into account. Lastly, we compare the spectra obtained with second-order perturbation theory to those determined by the special displacement method and we show that the full consideration of the electronic energy renormalization due to temperature is important to further improve the prediction of the phonon-assisted absorption in SiC. Our insights can be applied to predict the optical spectra of the less common SiC polytypes and other indirect-gap semiconductors in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1987

Nonlinear susceptibilities of finite conjugated organic polymers

Tight-binding calculations of the length dependence of the third-order molecular hyperpolarizability for polyenes and polyynes are reported. The pi-electron wave functions were determined by exploiting the limited translational symmetry of the molecules. Perturbation theory was used to calculate the longitudinal component of the electronic nonresonant hyperpolarizability. This is the first two-'band' calculation of third-order hyperpolarizabilities on finite pi-electron systems of varying length. In contrast to the results of the one-'band' models, the hyperpolarizability densities increase rapidly and then, after about 10-15 repeating units, approach an asymptotic value.

Beratan, David N.↗

DOE OSTI2025

Determination of Site Occupancy in the M–Pd–Zn (M = Cu, Ag, and Au) γ-Brass Phase by CALculation of PHAse Diagrams Modeling and Rietveld Refinement

The Pd–Zn γ-brass phase provides exciting opportunities for synthesizing site-isolated catalysts with precisely controlled Pd active site ensembles. Introducing a third metallic element into the γ-brass lattice further perturbs the catalytic active site ensembles. Here, in this work, we introduce coinage metallic elements M (M = Cu, Ag, and Au) into the Pd–Zn γ-brass phase and investigate the site occupation factors of each element in the γ-brass lattice. The CALculation of PHAse Diagrams (CALPHAD) modeling approach supported by energetics predicted by the density functional theory and X-ray and neutron diffraction with Rietveld refinement were used to identify the SOF on each Wyckoff site for various M amounts alloyed into the Pd–Zn γ-brass phase. The present analysis unveils the strong preference for Pd occupying the outer tetrahedral (OT) site in the γ-brass lattice, while the coinage metallic elements tend to substitute for Zn on the octahedral (OH) site. The determination of site occupancy in the bulk M–Pd–Zn γ-brass phase provides opportunities to investigate and tailor potential catalytically active site ensembles in the γ-brass phase materials.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Chiral effective field theory calculations of weak transitions in light nuclei

In this work, we report quantum Monte Carlo calculations of weak transitions in A≤10 nuclei, based on the Norfolk two- and three-nucleon chiral interactions, and associated one- and two-body axial currents. Furthermore, we find that the contribution from two-body currents is at the 2–3% level, with the exception of matrix elements entering the rates of 8Li, 8B, and 8He β decay. These matrix elements are suppressed in impulse approximation based on the (leading order) Gamow Teller transition operator alone; two-body currents provide a 20–30% correction, which is, however, insufficient to bring theory in agreement with experimental data. For the other transitions, the agreement with the data is satisfactory, and the results exhibit a negligible to mild model dependence when different combinations of Norfolk interactions are utilized to construct the nuclear wave functions. We report a complete study of two-body weak transition densities which reveals the expected universal behavior of two-body currents at short distances throughout the range of A=3 to A=10 systems considered here.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Machine learning dielectric screening for the simulation of excited state properties of molecules and materials

Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe–Salpeter equation (BSE) for electron–hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Small Molecule Inhibitor-Modulated Al 2 O 3 Atomic Layer Deposition on Monolayer MoS 2 for Controlled Nucleation

The integration of ultrathin dielectrics on two-dimensional (2D) semiconductors is essential for advancing beyond-Si electronics. However, the intrinsic inertness of van der Waals 2D basal planes remains a primary bottleneck to achieving uniform dielectric nucleation and growth. Here, in this study, we introduce a small molecule inhibitor (SMI)-modulated thermal atomic layer deposition (ALD) strategy, exemplified by aluminum oxide (Al 2 O 3 ) ALD on monolayer molybdenum disulfide (1L MoS 2 ) with acetic acid (HAc) SMI. The ABC-type sequence comprises HAc inhibitor (A), trimethylaluminum (TMA) precursor (B), and deionized H 2 O coreactant (C). In situ quartz crystal microbalance (QCM) studies reveal robust HAc adsorption on Al 2 O 3 and suppression of subsequent oxide growth on HAc-passivated surfaces. When applied to 1L MoS 2 , this inhibitory pathway enables HAc to selectively passivate nascent Al 2 O 3 nuclei formed on the MoS 2 surface, limiting their three-dimensional (3D) island coarsening and redirecting precursor adsorption toward the uncovered basal plane. Consequently, nearly continuous ultrathin (∼1.5 nm) Al 2 O 3 films are achieved on 1L MoS 2 with markedly improved uniformity compared to standard Al 2 O 3 ALD using TMA and H 2 O, as validated by atomic force microscopy (AFM), cross-sectional scanning transmission electron microscopy (STEM), and energy-dispersive X-ray spectroscopy (EDS). Density functional theory (DFT) calculations further provide atomistic insight into HAc-modulated Al 2 O 3 nucleation, corroborating the energetic preference of HAc for Al 2 O 3 over MoS 2 and attenuated TMA adsorption on HAc-passivated surfaces. Spatially resolved Raman spectroscopy also confirms that the HAc-modulated process preserves the structural integrity of 1L MoS 2 , with only minimal strain and doping perturbations observed after dielectric deposition. This SMI-modulated approach offers a broadly applicable framework for controlling ALD nucleation across various inhibitors, ALD chemistries, and 2D materials, opening opportunities for reliable dielectric integration in next-generation nanoelectronics.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Nuclei with Up to $$\varvec{A=6}$$ Nucleons with Artificial Neural Network Wave Functions

The ground-breaking works of Weinberg have opened the way to calculations of atomic nuclei that are based on systematically improvable Hamiltonians. Solving the associated many-body Schrödinger equation involves non-trivial difficulties, due to the non-perturbative nature and strong spin-isospin dependence of nuclear interactions. Artificial neural networks have proven to be able to compactly represent the wave functions of nuclei with up to $A=4$ nucleons. In this work, we extend this approach to $^6$Li and $^6$He nuclei, using as input a leading-order pionless effective field theory Hamiltonian. We successfully benchmark their binding energies, point-nucleon densities, and radii with the highly-accurate hyperspherical harmonics method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Halting Oxygen Evolution to Achieve Long Cycle Life in Sodium Layered Cathodes

Oxygen redox chemistries at high voltage have materialized as a revolutionary paradigm for cathodes with high-energy density; however, they are plagued by the challenges of labile oxygen loss and rapid degradations upon cycling, even after concerted endeavors from the research community. Here we propose a multi-concentration stratagem propelled by entropy reinforcement to enhance the electronic structure disorder (ESD) at high desodiation states for impeding undesired oxygen mobility and ensuring controlled oxygen activity, elucidated by density functional theory calculations. The increased disorder strengthens the reversible electrochemistry of lattice oxygen redox, leading to effectively suppressed P−O structural evolution and highly stable localized TMO 6 octahedral environments, as demonstrated by soft/hard X-ray absorption spectroscopy. Furthermore, through a comparative analysis of sodium-layered cathodes with different configuration entropy, we reveal that a high-entropy state induced by cationic disordering has the capacity to perturb cationic redox boundaries, significantly restraining the formation of detrimental O′3 phases. As a consequence, the high-voltage cycling stability has been greatly upgraded, up to 4.4 V versus Na + /Na, with an impressive 90.1 % capacity retention at 1 C over 100 cycles and 76.1 % capacity retention at 2 C over 300 cycles. In conclusion, the resilient oxygen redox, enabled through the control of ESD, broadens the horizons for entropy engineering and lays the foundation for advancements in high-energy, long-cycling, and safe batteries.

25 ENERGY STORAGE↗

DOE OSTI2025

Multiple Slater Determinants and Strong Spin-Fluctuations as Key Ingredients of the Electronic Structure of Electron- and Hole-Doped Pb 10− x Cu x (PO 4 ) 6 O

LK-99, with chemical formula Pb 10−x Cu x (PO 4 ) 6 O, was recently reported to be a room-temperature superconductor. While this claim has met with little support in a flurry of ensuing work, a variety of calculations (mostly based on density-functional theory) have demonstrated that the system possesses some unusual characteristics in the electronic structure, in particular flat bands. We have established previously that within DFT, the system is insulating with many characteristics resembling the classic cuprates, provided the structure is not constrained to the P3(143) symmetry nominally assigned to it. Here we describe the basic electronic structure of LK-99 within self-consistent many-body perturbative approach, quasiparticle self-consistent GW (QSGW) approximation and their diagrammatic extensions. QSGW predicts that pristine LK-99 is indeed a Mott/charge transfer insulator, with a bandgap gap in excess of 3 eV, whether or not constrained to the P3(143) symmetry. When Pb 9 Cu(PO 4 ) 6 O is hole-doped, the valence bands modify only slightly, and a hole pocket appears. However, two solutions emerge: a high-moment solution with the Cu local moment aligned parallel to neighbors, and a low-moment solution with Cu aligned antiparallel to its environment. In the electron-doped case the conduction band structure changes significantly: states of mostly Pb character merge with the formerly dispersionless Cu d state, and high-spin and low spin solutions once again appear. Thus we conclude that with suitable doping, the ground state of the system is not adequately described by a band picture, and that strong correlations are likely. Irrespective of whether this system class hosts superconductivity or not, the transition of Pb 10 (PO 4 ) 6 O from being a band insulator to Pb 9 Cu(PO 4 ) 6 O, a Mott insulator, and multi-determinantal nature of doped Mott physics make this an extremely interesting case-study for strongly correlated many-body physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1987

A simple theory of bimodal star formation

A model of bimodal star formation is presented, wherein massive stars form in giant molecular clouds (GNC), at a rate regulated by supernovae energy feedback through the interstellar medium, the heat input also ensuring that the initial mass function (IMF) remains skewed towards massive stars. The low mass stars form at a constant rate. The formation of the GMC is governed by the dynamics of the host galaxy through the rotation curve and potential perturbations such as a spiral density wave. The characteristic masses, relative normalizations, and rates of formation of the massive and low mass modes of star formation may be tightly constrained by the requirements of the chemical evolution in the Solar Neighborhood. Good fits were obtained for the age metallicity relation and the metallicity structure of thin disk and spheroid stars only for a narrow range of these parameters.

Wyse, Rosemary F. G.↗

DOE OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗

DOE OSTI2025

Real-time scattering in Ising field theory using matrix product states

We study scattering in Ising field theory (IFT) using matrix product states and the time-dependent variational principle. IFT is a one-parameter family of strongly coupled nonintegrable quantum field theories in 1+1 dimensions, interpolating between massive free fermion theory and Zamolodchikov's integrable massive 𝐸 8 theory. Particles in IFT may scatter either elastically or inelastically. In the postcollision wave function, particle tracks from all final-state channels occur in superposition; processes of interest can be isolated by projecting the wave function onto definite particle sectors, or by evaluating energy density correlation functions. Using numerical simulations we determine the time delay of elastic scattering and the probability of inelastic particle production as a function of collision energy. We also study the mass and width of the lightest resonance near the 𝐸 8 point in detail. Close to both the free fermion and 𝐸 8 theories, our results for both elastic and inelastic scattering are in good agreement with expectations from form-factor perturbation theory. Using numerical computations to go beyond the regime accessible by perturbation theory, we find that the high-energy behavior of the two-to-two particle scattering probability in IFT is consistent with a conjecture of Zamolodchikov. Our results demonstrate the efficacy of tensor-network methods for simulating the real-time dynamics of strongly coupled quantum field theories in 1+1 dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2001

Physical and Electronic Isolation of Carbon Nanotube Conductors

Multi-walled nanotubes are proposed as a method to electrically and physically isolate nanoscale conductors from their surroundings. We use tight binding (TB) and density functional theory (DFT) to simulate the effects of an external electric field on multi-wall nanotubes. Two categories of multi-wall nanotube are investigated, those with metallic and semiconducting outer shells. In the metallic case, simulations show that the outer wall effectively screens the inner core from an applied electric field. This offers the ability to reduce crosstalk between nanotube conductors. A semiconducting outer shell is found not to perturb an electric field incident on the inner core, thereby providing physical isolation while allowing the tube to remain electrically coupled to its surroundings.

OKeeffe, James↗

DOE OSTI2022

Electron dynamics in extended systems within real-time time-dependent density-functional theory

Due to a beneficial balance of computational cost and accuracy, real-time time-dependent density-functional theory has emerged as a promising first-principles framework to describe electron real-time dynamics. Here we discuss recent implementations around this approach, in particular in the context of complex, extended systems. Results include an analysis of the computational cost associated with numerical propagation and when using absorbing boundary conditions. Here we extensively explore the shortcomings for describing electron–electron scattering in real time and compare to many-body perturbation theory. Modern improvements of the description of exchange and correlation are reviewed. In this work, we specifically focus on the Qb@ll code, which we have mainly used for these types of simulations over the last years, and we conclude by pointing to further progress needed going forward.

2D materials↗

DOE OSTI2020

Biased tracers in redshift space in the EFTofLSS with exact time dependence

We study the effect of the Einstein-de Sitter (EdS) approximation on the one-loop power spectrum of galaxies in redshift space in the Effective Field Theory of Large-Scale Structure. The dark matter density perturbations and velocity divergence are treated with exact time dependence. Splitting the density perturbation into its different temporal evolutions naturally gives rise to an irreducible basis of biases. While, as in the EdS approximation, at each time this basis spans a seven-dimensional space, this space is a slightly different one, and the difference is captured by a single calculable time- and vec k-dependent function. We then compute the redshift-space galaxy one-loop power spectrum with the EdS approximation (P EdS-approx ) and without (P Exact ). For the monopole we find P 0 Exact /P 0 Ed S-approx ~ 1.003 and for the quadrupole P 2 Exact /P 2 EdS-approx ~ 1.007 at z=0.57, and sharply increasing at lower redshifts. Finally, we show that a substantial fraction of the effect remains even after allowing the bias coefficients to shift within a physically allowed range. This suggests that the EdS approximation can only fit the data to a level of precision that is roughly comparable to the precision of the next generation of cosmological surveys. Furthermore, we find that implementing the exact time dependence formalism is not demanding and is easily applicable to data. Both of these points motivate a direct study of this effect on the cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1972

Dynamics of stratified liquids in the presence of space charge.

The dispersion characteristics of small signal electrohydrodynamic space charge and polarization waves in continuously and discretely stratified liquids are developed. The liquids, modeled as incompressible, inviscid, and perfectly insulating have distributions of mass density, charge density, permittivity, and electric field. A general set of relations for perturbation field and flow variables defined for a prototype layer are used to describe the dynamics of many superposed planar layers. It is shown that the general relations can be used to represent systems of coupled layers which approximate continuous distributions by a series of step functions. Specific examples of weak gradient and exponential distributions are presented showing that the solution found directly from the distributed theory is approached by the system of coupled layers, if the limit is taken in which the number of layers become large while each layer thickness approaches zero.

Zahn, M.↗

DOE OSTI2023

Shape polarization in the tin isotopes near N = 60 from precision g-factor measurements on short-lived 11/2 - isomers

The g factors of 11/2 - isomers in semimagic 109 Sn and 111 Sn (isomeric lifetimes $\tau$ = 2.9(3) ns and $\tau$ = 14.4(7) ns, respectively) were measured by an extension of the Time Differential Perturbed Angular Distribution technique, which uses LaBr3 detectors and the hyperfine fields of a gadolinium host to achieve precise measurements in a new regime of short-lived isomers. The results, g(11/2 - ; 109 Sn) = -0.186(8) and g(11/2 - ; 111 Sn) = -0.214(4), are significantly lower in magnitude than those of the 11/2 - isomers in the heavier isotopes and depart from the value expected for a near pure neutron h 11/2 configuration. Broken-symmetry density functional theory calculations applied to the sequence of 11/2 - states reproduce the magnitude and location of this deviation. The g(11/2 - ) values are affected by shape core polarization; the odd Oh 11/2 neutron couples to J π = 2 + , 4 + , 6 + ... configurations in the weakly-deformed effective core, causing a decrease in the g-factor magnitudes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗