Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio methods”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Ab Initio Multiple Spawning Photochemical Dynamics of DMABN Using GPUs

The ultrafast decay dynamics of 4-( N,N -dimethylamino)benzonitrile (DMABN) following photoexcitation was studied with the ab initio multiple spawning (AIMS) method, combined with GPU-accelerated linear-response time-dependent density functional theory (LR-TDDFT). We validate the LR-TDDFT method for this case and then present a detailed analysis of the first ≈200 fs of DMABN excited-state dynamics. Almost complete nonadiabatic population transfer from S 2 (the initially populated bright state) to S 1 takes place in less than 50 fs, without significant torsion of the dimethylamino (DMA) group. Significant torsion of the DMA group is only observed after the nuclear wavepacket reaches S 1 and acquires locally excited electronic character. Here, our results show that torsion of the DMA group is not prerequisite for nonadiabatic transitions in DMABN, although such motion is indeed relevant on the lowest excited state (S 1 ).

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

Nonadiabatic Photodynamics of Retinal Protonated Schiff Base in Channelrhodopsin 2

Channelrhodopsin 2 (ChR2) is a light-gated ion channel and an important tool in optogenetics. Photoisomerization of retinal protonated Schiff base (RPSB) in ChR2 triggers channel activation. Despite the importance of ChR2 in optogenetics, the detailed mechanism for photoisomerization and channel activation is still not fully understood. Here, we report on computer simulations to investigate the photoisomerization mechanism and its effect on the activation of ChR2. Nonadiabatic dynamics simulation of ChR2 was carried out using the ab initio multiple spawning (AIMS) method and quantum mechanics/molecular mechanics (QM/MM) with a restricted ensemble Kohn–Sham (REKS) treatment of the QM region. Our results agree well with spectroscopic measurements and reveal that the RPSB isomerization is highly specific around the C 13 =C 14 bond and follows the “aborted bicycle-pedal” mechanism. In addition, RPSB photoisomerization facilitates its deprotonation and partially increases the hydration level in the channel, which could trigger subsequent channel opening and ion conduction.

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

Many-body enhancement of high-harmonic generation in monolayer MoS2

Abstract Many-body effects play an important role in enhancing and modifying optical absorption and other excited-state properties of solids in the perturbative regime, but their role in high harmonic generation (HHG) and other nonlinear response beyond the perturbative regime is not well-understood. We develop here an ab initio many-body method to study nonperturbative HHG based on the real-time propagation of the non-equilibrium Green’s function with the GW self energy. We calculate the HHG of monolayer MoS 2 and obtain good agreement with experiment, including the reproduction of characteristic patterns of monotonic and nonmonotonic harmonic yield in the parallel and perpendicular responses, respectively. Here, we show that many-body effects are especially important to accurately reproduce the spectral features in the perpendicular response, which reflect a complex interplay of electron-hole interactions (or exciton effects) in tandem with the many-body renormalization and Berry curvature of the independent quasiparticle bandstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Machine learning the electric field response of condensed phase systems using perturbed neural network potentials

Abstract The interaction of condensed phase systems with external electric fields is of major importance in a myriad of processes in nature and technology, ranging from the field-directed motion of cells (galvanotaxis), to geochemistry and the formation of ice phases on planets, to field-directed chemical catalysis and energy storage and conversion systems including supercapacitors, batteries and solar cells. Molecular simulation in the presence of electric fields would give important atomistic insight into these processes but applications of the most accurate methods such as ab-initio molecular dynamics (AIMD) are limited in scope by their computational expense. Here we introduce Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) to push back the accessible time and length scales of such simulations. We demonstrate that important dielectric properties of liquid water including the field-induced relaxation dynamics, the dielectric constant and the field-dependent IR spectrum can be machine learned up to surprisingly high field strengths of about 0.2 V Å −1 without loss in accuracy when compared to ab-initio molecular dynamics. This is remarkable because, in contrast to most previous approaches, the two neural networks on which PNNP MD is based are exclusively trained on molecular configurations sampled from zero-field MD simulations, demonstrating that the networks not only interpolate but also reliably extrapolate the field response. PNNP MD is based on rigorous theory yet it is simple, general, modular, and systematically improvable allowing us to obtain atomistic insight into the interaction of a wide range of condensed phase systems with external electric fields.

Science & Technology - Other Topics↗
DOE OSTI2021

Machine learning for deep elastic strain engineering of semiconductor electronic band structure and effective mass

Abstract The controlled introduction of elastic strains is an appealing strategy for modulating the physical properties of semiconductor materials. With the recent discovery of large elastic deformation in nanoscale specimens as diverse as silicon and diamond, employing this strategy to improve device performance necessitates first-principles computations of the fundamental electronic band structure and target figures-of-merit, through the design of an optimal straining pathway. Such simulations, however, call for approaches that combine deep learning algorithms and physics of deformation with band structure calculations to custom-design electronic and optical properties. Motivated by this challenge, we present here details of a machine learning framework involving convolutional neural networks to represent the topology and curvature of band structures in k -space. These calculations enable us to identify ways in which the physical properties can be altered through “deep” elastic strain engineering up to a large fraction of the ideal strain. Algorithms capable of active learning and informed by the underlying physics were presented here for predicting the bandgap and the band structure. By training a surrogate model with ab initio computational data, our method can identify the most efficient strain energy pathway to realize physical property changes. The power of this method is further demonstrated with results from the prediction of strain states that influence the effective electron mass. We illustrate the applications of the method with specific results for diamonds, although the general deep learning technique presented here is potentially useful for optimizing the physical properties of a wide variety of semiconductor materials.

Tsymbalov, Evgenii↗
DOE OSTI2022

Dissociative electron attachment to 5-bromo-uracil: non-adiabatic dynamics on complex-valued potential energy surfaces

Electron induced dissociation reactions are relevant to many fields, ranging from prebiotic chemistry to cancer treatments. However, the simulation of dissociation electron attachment (DEA) dynamics is very challenging because the auto-ionization widths of the transient negative ions must be accounted for. We propose an adaptation of the ab initio multiple spawning (AIMS) method for complex-valued potential energy surfaces, along the lines of recent developments based on surface hopping dynamics. Our approach combines models for the energy dependence of the auto-ionization widths, obtained from scattering calculations, with survival probabilities computed for the trajectory basis functions employed in the AIMS dynamics. The method is applied to simulate the DEA dynamics of 5-bromouracil in full dimensionality, i.e., taking all the vibrational modes into consideration. The propagation starts on the $π_{2}^{*}$ resonance state and describes the formation of Br - anions mediated by non-adiabatic couplings. The potential energies, gradients and non-adiabatic couplings were computed with the fractional-occupancy molecular orbital complete-active-space configuration-interaction method, and the calculated DEA cross section are consistent with the observed DEA intensities.

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

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗
DOE OSTI2015

Ab initio molecular orbital-configuration interaction based quantum master equation (MOQME) approach to the dynamic first hyperpolarizabilities of asymmetric π-conjugated systems

In this study, we apply the ab initio molecular orbital - configuration interaction based quantum master equation (MOQME) approach to the calculation and analysis of the dynamic first hyperpolarizabilities (β) of asymmetric π-conjugated molecules. In this approach, we construct the excited state models by the ab initio configuration interaction singles method. Then, time evolutions of system reduced density matrix ρ(t) and system polarization p(t) are calculated by the QME approach. Dynamic β in the second harmonic generation is calculated based on the nonperturbative definition of nonlinear optical susceptibility, using the frequency domain system polarization p(ω). Spatial contributions of electrons to β are analyzed based on the dynamic hyperpolarizability density map, which visualizes the second-order response of charge density oscillating with a frequency of 2ω. We apply the present method to the calculation of the dynamic β of a series of donor/acceptor substituted polyene oligomers, and then discuss the applicability of the MOQME method to the calculation and analysis of dynamic NLO properties of molecular systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Specific interactions between DNA and regulatory protein controlled by ligand-binding: Ab initio molecular simulation

The catabolite activator protein (CAP) is one of the regulatory proteins controlling the transcription mechanism of gene. Biochemical experiments elucidated that the complex of CAP with cyclic AMP (cAMP) is indispensable for controlling the mechanism, while previous molecular simulations for the monomer of CAP+cAMP complex revealed the specific interactions between CAP and cAMP. However, the effect of cAMP-binding to CAP on the specific interactions between CAP and DNA is not elucidated at atomic and electronic levels. We here considered the ternary complex of CAP, cAMP and DNA in solvating water molecules and investigated the specific interactions between them at atomic and electronic levels using ab initio molecular simulations based on classical molecular dynamics and ab initio fragment molecular orbital methods. The results highlight the important amino acid residues of CAP for the interactions between CAP and cAMP and between CAP and DNA.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Communication: XFAIMS—eXternal Field Ab Initio Multiple Spawning for electron-nuclear dynamics triggered by short laser pulses

Attoscience is an emerging field where attosecond pulses or few cycle IR pulses are used to pump and probe the correlated electron-nuclear motion of molecules. We present the trajectory-guided eXternal Field Ab Initio Multiple Spawning (XFAIMS) method that models such experiments “on-the-fly,” from laser pulse excitation to fragmentation or nonadiabatic relaxation to the ground electronic state. For the photoexcitation of the LiH molecule, we show that XFAIMS gives results in close agreement with numerically exact quantum dynamics simulations, both for atto- and femtosecond laser pulses. As a result, we then show the ability of XFAIMS to model the dynamics in polyatomic molecules by studying the effect of nuclear motion on the photoexcitation of a sulfine (H 2 CSO).

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

Development of reactive force fields using ab initio molecular dynamics simulation minimally biased to experimental data

Incorporation of quantum mechanical electronic structure data is necessary to properly capture the physics of many chemical processes. Proton hopping in water, which involves rearrangement of chemical and hydrogen bonds, is one such example of an inherently quantum mechanical process. Standard ab initio molecular dynamics (AIMD) methods, however, do not yet accurately predict the structure of water and are therefore less than optimal for developing force fields. We have instead utilized here a recently developed method which minimally biases AIMD simulations to match limited experimental data to develop novel multiscale reactive molecular dynamics (MS-RMD) force fields by using relative entropy minimization. In this paper, we present two new MS-RMD models using such a parameterization: one which employs water with harmonic internal vibrations and another which uses anharmonic water. We show that the newly developed MS-RMD models very closely reproduce the solvation structure of the hydrated excess proton in the target AIMD data. We also find that the use of anharmonic water increases proton hopping, thereby increasing the proton diffusion constant.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Ultrafast electron diffraction of photoexcited gas-phase cyclobutanone predicted by ab initio multiple cloning simulations

We present the result of our calculations of ultrafast electron diffraction (UED) for cyclobutanone excited into the S2 electronic state, which is based on the non-adiabatic dynamics simulations with the Ab Initio Multiple Cloning (AIMC) method with the electronic structure calculated at the SA(3)-CASSCF(12,12)/aug-cc-pVDZ level of theory. The key features in the UED pattern were identified, which can be used to distinguish between the reaction pathways observed in the AIMC dynamics, although there is a significant overlap between representative signals due to the structural similarity of the products. The calculated UED pattern can be compared with the experiment.

Chemistry↗
DOE OSTI2020

New ideas in constraining nuclear forces

In recent years, nuclear physics has benefited greatly from the development of powerful ab initio many-body methods and their combination with interactions from chiral effective field theory. With increasing computational power and continuous development of these methods, we are entering an era of precision nuclear physics. Indeed, uncertainties from nuclear Hamiltonians now dominate over uncertainties from many-body methods. This review summarizes the current status of, and future directions in, deriving and constraining nuclear Hamiltonians.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Ab initio calculation of atomic solid hydrogen phases based on Gutzwiller many-body wave functions

We apply two ab initio many-body methods based on Gutzwiller wave functions, i.e., correlation matrix renormalization theory (CMRT) and Gutzwiller conjugate gradient minimization (GCGM), to the study of crystalline phases of atomic hydrogen. Both methods avoid empirical Hubbard U parameters and are free from double-counting issues. CMRT employs a Gutzwiller-type approximation that enables efficient calculations, while GCGM goes beyond this approximation to achieve higher accuracy at higher computational cost. By benchmarking against available quantum Monte Carlo (QMC) results, we demonstrate that while both methods are more accurate than the widely used density-functional theory, GCGM systematically captures additional correlation energy missing in CMRT, leading to significantly improved total energy predictions. We also show that by including the correlation energy Ec from local density approximation in the CMRT calculation, CMRT + E c produces energy in better agreement with the QMC results in these hydrogen lattice systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

The sea of quarks and antiquarks in the nucleon

The quark and gluon structure of the proton has been under intense experimental and theoretical investigation for five decades. Even for the distributions of the well-studied valence quarks, challenges such as the value of the down quark to up quark ratio at high fractional momenta remain. Much of the sea of quark-antiquark pairs emerges from the splitting of gluons and is well described by perturbative evolution in quantum chromodynamics. However, experiments confirm that there is a non-perturbative component to the sea that is not well understood and hitherto has been difficult to calculate with ab initio non-perturbative methods. This non-perturbative structure shows up, perhaps most directly, in the flavor dependence of the sea antiquark distributions. Even though some of the general trends can be reproduced by models, there are features of the data that do not seem to be well described. Conclusively, this article discusses the experimental situation, the status of calculations and models, and the directions where these studies will progress in the near future.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

First-principles studies on the equation-of-state, thermal-conductivity, and opacity of deuterium-tritium and polystyrene (CH) for inertial confinement fusion applications

Using first-principles (FP) methods, we have performed ab initio compute for the equation of state (EOS), thermal conductivity, and opacity of deuterium-tritium (DT) in a wide range of densities and temperatures for inertial confinement fusion (ICF) applications. These systematic investigations have recently been expanded to accurately compute the plasma properties of CH ablators under extreme conditions. In particular, the first-principles EOS and thermal-conductivity tables of CH are self-consistently built from such FP calculations, which are benchmarked by experimental measurements. When compared with the traditional models used for these plasma properties in hydrocodes, significant differences have been identified in the warm dense plasma regime. When these FP-calculated properties of DT and CH were used in our hydrodynamic simulations of ICF implosions, we found that the target performance in terms of neutron yield and energy gain can vary by a factor of 2 to 3, relative to traditional model simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

First-principles investigation of Sc-III/IV under high pressure

Using an ab initio evolutionary structure prediction method in conjunction with density functional theory, we performed a systematic investigation of the structural transition of elemental scandium under pressure up to 250 GPa. Our prediction successfully reproduced several allotropes which have been reported in the literature, including the Sc-I, Sc-II, and Sc-V. Moreover, we observed a series of energetically degenerate and geometrically similar structures at 110–195 GPa, which can partly explain the experimental observations regarding the unsolved phases III and IV reported by Y. Akahama et al. [Phys. Rev. Lett. 94, 195503 (2005)]. A detailed comparison of the powder x-ray diffraction pattern suggested that the Ccca–20 phase is a likely candidate for the observed Sc-III, while Sc-IV may be explained by a structure with random stacking of two different structural units. We also used the candidate Sc-III structure as the model system to explore its superconducting behavior under pressures between 80 and 130 GPa. As a result, the predicted superconducting transition temperature Tc values are in satisfactory agreement with previous experimental results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Towards Highly Scalable Ab Initio Molecular Dynamics (AIMD) Simulations on the Intel Knights Landing Manycore Processor

The Ab Initio Molecular Dynamics (AIMD) method allows scientists to treat the dynamics of molecular and condensed phase systems while retaining a first-principles-based description of their interactions. This extremely important method has tremendous computational requirements, because the electronic Schr¨odinger equation, approximated using Kohn-Sham Density Functional Theory (DFT), is solved at every time step. With the advent of manycore architectures, application developers have a significant amount of processing power within each compute node that can only be exploited through massive parallelism. A compute intensive application such as AIMD forms a good candidate to leverage this processing power. In this paper, we focus on adding thread level parallelism to the plane wave DFT methodology implemented in NWChem. Through a careful optimization of tall-skinny matrix products, which are at the heart of the Lagrange multiplier and nonlocal pseudopotential kernels, as well as 3D FFTs, our OpenMP implementation delivers excellent strong scaling on the latest Intel Knights Landing (KNL) processor. We assess the efficiency of our Lagrange multiplier kernels by building a Roofline model of the platform, and verify that our implementation is close to the roofline for various problem sizes. Finally, we present strong scaling results on the complete AIMD simulation for a 64 water molecules test case, that scales up to all 68 cores of the Knights Landing processor.

Environmental Molecular Sciences Laboratory↗