Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quantum chemistry 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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Steady-state properties of multi-orbital systems using quantum Monte Carlo

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Chemistry↗
DOE OSTI2021

Comparing (stochastic-selection) ab initio multiple spawning with trajectory surface hopping for the photodynamics of cyclopropanone, fulvene, and dithiane

Ab Initio Multiple Spawning (AIMS) simulates the excited-state dynamics of molecular systems by representing nuclear wavepackets in a basis of coupled traveling Gaussian functions, called trajectory basis functions (TBFs). New TBFs are spawned when nuclear wavepackets enter regions of strong nonadiabaticity, permitting the description of non-Born–Oppenheimer processes. The spawning algorithm is simultaneously the blessing and the curse of the AIMS method: it allows for an accurate description of the transfer of nuclear amplitude between different electronic states, but it also dramatically increases the computational cost of the AIMS dynamics as all TBFs are coupled. Recently, a strategy coined stochastic-selection AIMS (SSAIMS) was devised to limit the ever-growing number of TBFs and tested on simple molecules. In this work, we use the photodynamics of three different molecules—cyclopropanone, fulvene, and 1,2-dithiane—to investigate (i) the potential of SSAIMS to reproduce reference AIMS results for challenging nonadiabatic dynamics, (ii) the compromise achieved by SSAIMS in obtaining accurate results while using the smallest average number of TBFs as possible, and (iii) the performance of SSAIMS in comparison to the mixed quantum/classical method trajectory surface hopping (TSH)—both in terms of its accuracy and computational cost. We show that SSAIMS can accurately reproduce the AIMS results for the three molecules considered at a much cheaper computational cost, often close to that of TSH. We deduce from these tests that an overlap-based criterion for the stochastic-selection process leads to the best agreement with the reference AIMS dynamics for the smallest average number of TBFs.

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

Computational Determination of the Refractive Index of TATB, LLM-105, β-HMX, & α-RDX

In this study, we use Materials Studio CASTEP, an ab initio quantum mechanical (QM) program employing density functional theory (DFT), to compute the refractive index (RI) of solid energetic materials, specifically, PETN (tetragonal phase), TATB, LLM-105, β-HMX, and α-RDX. The use of quantum mechanical DFT methods affords the investigation of the optical mate rial properties of a system without the need for any experimental input. The RI values of PETN and TATB have been experimentally determined previously and are compared here to the QM derived RI values as validation of the theoretical approach in determining material optical proper ties. To date, experimental measurements to determine the RI values for LLM-105 have yet to be determined. This fact obviates the need to determine LLM-105 RI values using quantum mechanical means.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Tunability and Long-Range Enhancement of Resonance Energy Transfer Facilitated by Plasmonic Nanorods

Resonance energy transfer (RET) between molecules or quantum dots is an important process in many energy-related applications. Different environmental structures have been studied and demonstrated to be able to enhance the RET rates between a nearby donor–acceptor pair. In particular, cylindrical silver nanorods and nanowires have shown an extraordinary ability to transfer energy along their longitudinal axes over large distances. However, the detailed mechanism of such transfer and the effects on the molecular RET process are as yet elusive. In this study, we use the recently developed computational tool based on the plasmon-coupled RET method to systematically study the effects of nanorods with different dimensions on the RET rates. We find that highly frequency-dependent coupling factor (CF) spectra, whose amplitudes determine RET rates, can be obtained due to the localized surface plasmon polariton modes of the rods with nanoscale dimensions. Simple phenomenological models can be derived for the wavelengths of CF peaks in relation to the length and width of the nanorods, providing easy tunability for enhancing the RET rate in specific wavelength ranges. When coupled to longer rods with mesoscale lengths, exponential decay of the CF over long donor–acceptor distances with a small decay constant is observed, leading to the possibility of long-range RET processes. Furthermore, drumhead resonance modes emerge on the flat ends of the rod when the rod’s diameter reaches 300 nm, resulting in extra enhancement to RET rate compared to certain thinner rods. Furthermore, these findings shed new light on the mechanism of plasmonic enhancement with silver nanorods and establish design principles for how to optimally utilize these structures to manipulate RET processes for various applications.

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

Electronic specific heat capacities and entropies from density matrix quantum Monte Carlo using Gaussian process regression to find gradients of noisy data

In this work, we present a machine learning approach to calculating electronic specific heat capacities for a variety of benchmark molecular systems. Our models are based on data from density matrix quantum Monte Carlo, which is a stochastic method that can calculate the electronic energy at finite temperature. As these energies typically have noise, numerical derivatives of the energy can be challenging to find reliably. In order to circumvent this problem, we use Gaussian process regression to model the energy and use analytical derivatives to produce the specific heat capacity. From there, we also calculate the entropy by numerical integration. We compare our results to cubic splines and finite differences in a variety of molecules in which Hamiltonians can be diagonalized exactly with full configuration interaction. We finally apply this method to look at larger molecules where exact diagonalization is not possible and make comparisons with more approximate ways to calculate the specific heat capacity and entropy.

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

How Does Bending the Uranyl Unit Influence Its Spectroscopy and Luminescence?

We report bent uranyl complexes can be formed with chloride ligands and 1,10-phenanthroline (phen) ligands bound to the equatorial and axial planes of the uranyl(VI) moiety, as revealed by the crystal structures, IR and Raman spectroscopy, and quantum-chemical calculations. With the goal of probing the influence of chloride and phenanthroline coordination enforcing the bending on the absorption and emission spectra of this complex, spin–orbit time-dependent density functional theory calculations for the bare uranyl complexes as well as for the free UO 2 Cl 2 subunit and the UO 2 Cl 2 (phen) 2 complex were performed. The emission spectra have been fully simulated by ab initio methods and compared to experimental photoluminescence spectra, recorded for the first time for UO 2 Cl 2 (phen) 2 . Notably, the bending of uranyl in UO2Cl2 and UO 2 Cl 2 (phen) 2 triggers excitations of the uranyl bending mode, yielding a denser luminescence spectrum.

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

Advances in Singlet Fission Chromophore Design Enabled by Vibrational Spectroscopies

Singlet fission leads to the formation of two separate triplet T 1 excitons from an initial singlet S1 exciton through 1 (TT) and 1 (T...T), multiexcitonic intermediates that retain singlet character. Its ability to achieve external quantum efficiencies higher than 100% made it an attractive candidate for optoelectronic device applications. However, singlet fission has not been applied widely despite having been investigated by a myriad of spectroscopic methods, in part due to our poor understanding of how to optimize molecular structure and packing in chromophores well-suited to large-scale production. Vibrational spectroscopies provide a solution, because they directly probe nuclear motions, allowing us to monitor evolving structural changes in molecules undergoing singlet fission, thus providing us with roadmaps to design molecules suitable for optoelectronic applications. Furthermore, this Perspective reviews the contributions and analyzes the future directions of vibrational spectroscopies to the advancement in our knowledge about the mechanisms and rational designing of chromophores undergoing efficient singlet fission.

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

Ethylene oxidation on unpromoted silver catalysts: Reaction pathway and selectivity analysis using DFT calculations

Chemical conversions in catalytic partial oxidation processes of light hydrocarbons are responsible for the production of numerous industrial chemicals, plastics, and intermediates. These processes are relatively expensive to perform, and are typically operated at high thermodynamic inefficiency, so the development of novel, highly efficient catalysts would prove to be very cost effective. Herein, our study focused on surface catalytic mechanisms of the ethylene oxide (EO) formation process. Periodic plane-wave Density Functional Theory (DFT) methods were used to analyze related reaction mechanisms on the Ag(111) surface facet with low coverage. Energetic changes of related species and pathways were calculated. Key surface species are identified to suggest the factors for the observed selectivity during EO formation. Our results are consistent with previous kinetic modeling efforts in the literature which did not employ DFT analysis. Lastly, our study demonstrates how fundamental theoretical investigations and multi-scale modeling techniques are currently impacting the advancement of rational catalyst design and microkinetic modeling techniques in the light hydrocarbon processing industry.

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

Tests of the DFT Ladder for the Fulminic Acid Challenge

Properties of the historically pivotal fulminic acid (HCNO) molecule have been computed with a panoply of 473 density functionals of all varieties, providing a snapshot of the performance of contemporary density functional theory (DFT) for a challenging chemical system. Exhaustive tabulations and statistical analyses have been carried out for geometric parameters, vibrational frequencies, barriers to linearity, and the HCN–O dissociation energy. As the DFT ladder is climbed, confusion rather than consensus ensues regarding the details of the distinctive, extremely flat H–C–N bending potential of fulminic acid and whether the equilibrium structure is linear or bent. While high-ranking DFT functionals produce the smallest errors for the HCN + O( 3 P) → HCNO reaction energy, lower rungs emerge as the best performers for many of the bond distances and harmonic vibrational frequencies. This research shows that the current DFT zoo of approximations does not constitute a transparent ladder of increasingly accurate methods that consistently converges on definitive predictions for various properties of HCNO. Additional analyses are performed on the side effects of popular dispersion corrections on the covalently bonded properties and thermochemistry of HCNO.

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

Predictive Model of Charge Mobilities in Organic Semiconductor Small Molecules with Force-Matched Potentials

Charge mobility of crystalline organic semiconductors (OSC) is limited by local dynamic disorder. Recently, the charge mobility for several high mobility OSCs, including TIPS-pentacene, were accurately predicted from a density functional theory (DFT) simulation constrained by the crystal structure and the inelastic neutron scattering spectrum, which provide direct measures of the structure and the dynamic disorder in the length scale and energy range of interest. However, the computational expense required for calculating all of the atomic and molecular forces is prohibitive. In this work, we demonstrate the use of density functional tight binding (DFTB), a semiempirical quantum mechanical method that is 2 to 3 orders of magnitude more efficient than DFT. We show that force matching a many-body interaction potential to DFT derived forces yields highly accurate DFTB models capable of reproducing the low-frequency intricacies of experimental inelastic neutron scattering (INS) spectra and accurately predicting charge mobility. We subsequently predicted charge mobilities from our DFTB model of a number of previously unstudied structural analogues to TIPS-pentacene using dynamic disorder from DFTB and transient localization theory. The approach we establish here could provide a truly rapid simulation pathway for accurate materials properties prediction, in our vision applied to new OSCs with tailored properties.

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

2p x-ray absorption spectroscopy of 3d transition metal systems

This review provides an overview of the different methods and computer codes that are used to interpret 2p x-ray absorption spectra of 3d transition metal ions. We first introduce the basic parameters and give an overview of the methods used. We start with the semi-empirical multiplet codes and compare the different codes that are available. A special chapter is devoted to the user friendly interfaces that have been written on the basis of these codes. Next we discuss the first principle codes based on band structure, including a chapter on Density Functional theory based approaches. We also give an overview of the first-principle multiplet codes that start from a cluster calculation and we discuss the wavefunction based methods, including multi-reference methods. We end the review with a discussion of the link between theory and experiment and discuss the open issues in the spectral analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An initial investigation of accuracy required for the identification of small molecules in complex samples using quantum chemical calculated NMR chemical shifts

The majority of primary and secondary metabolites in nature have yet to be identified, representing a major challenge for metabolomics studies that currently require reference libraries from analyses of authentic compounds. Using currently available analytical methods, complete chemical characterization of metabolomes is infeasible for both technical and economic reasons. For example, unambiguous identification of metabolites is limited by the availability of authentic chemical standards, which, for the majority of molecules, do not exist. Computationally predicted or calculated data are a viable solution to expand the currently limited metabolite reference libraries, if such methods are shown to be sufficiently accurate. For example, determining nuclear magnetic resonance (NMR) spectroscopy spectra in silico has shown promise in the identification and delineation of metabolite structures. Many researchers have been taking advantage of density functional theory (DFT), a computationally inexpensive yet reputable method for the prediction of carbon and proton NMR spectra of metabolites. However, such methods are expected to have some error in predicted 13 >C and 1 H NMR spectra with respect to experimentally measured values. This leads us to the question–what accuracy is required in predicted 13 C and 1 H NMR chemical shifts for confident metabolite identification? Using the set of 11,716 small molecules found in the Human Metabolome Database (HMDB), we simulated both experimental and theoretical NMR chemical shift databases. We investigated the level of accuracy required for identification of metabolites in simulated pure and impure samples by matching predicted chemical shifts to experimental data. We found 90% or more of molecules in simulated pure samples can be successfully identified when errors of 1 H and 13 C chemical shifts in water are below 0.6 and 7.1 ppm, respectively, and below 0.5 and 4.6 ppm in chloroform solvation, respectively. In simulated complex mixtures, as the complexity of the mixture increased, greater accuracy of the calculated chemical shifts was required, as expected. However, if the number of molecules in the mixture is known, e.g., when NMR is combined with MS and sample complexity is low, the likelihood of confident molecular identification increased by 90%.

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

Non-Empirical and Self-Interaction Corrections for DFTB: Towards Accurate Quantum Simulations for Large Mesoscale Systems (Final Report)

This project was comprised of two complementary (but parallel) thrusts: (1) implementing massively-parallelized computing hardware (with new computational hardware that may replace GPUs) in the density functional tight binding (DFTB) approach and (2) developing new capabilities in DFTB to calculate the electronic structure and dynamics of large chemical systems. While classical molecular dynamics can handle hundreds of thousands of atoms, it cannot provide a first-principles based description of chemical systems at the quantum level. At the other extreme, conventional Kohn-Sham DFT methods can probe the true quantum mechanical nature of chemical systems; however, these methods cannot tackle the large sizes and length scales relevant to dynamics simulations of realistic systems. The DFTB formalism utilized in this project provides a viable approach for probing these large systems at a quantum mechanical level of detail. However, to utilize the DFTB approach for accurate calculations of electronic properties, it is crucial to incorporate quantum-based corrections in DFTB since exchange-correlation effects can still remain very strong in these large systems. At the same time, enhancing the computational efficiency of DFTB is also essential since optimal computational performance is required for addressing the large size scales associated with realistic chemical systems. As such, the new non-empirical corrections and computing hardware enhancements implemented in this project will enable accurate and computationally efficient approaches to directly probe electronic properties in these large, complex systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Proton-detected solid-state NMR spectroscopy of spin-1/2 nuclei with large chemical shift anisotropy

Constant-time (CT) dipolar heteronuclear multiple quantum coherence (D-HMQC) has previously been demonstrated as a method for proton detection of high-resolution wideline NMR spectra of spin-1/2 nuclei with large chemical shift anisotropy (CSA). However, 1 H transverse relaxation and t 1 -noise often reduce the sensitivity of D-HMQC experiments, preventing the theoretical gains in sensitivity provided by 1 H detection from being realized. In this paper, we demonstrate a series of improved pulse sequences for 1 H detection of spin-1/2 nuclei under fast MAS, with 195 Pt SSNMR experiments on cisplatin as an example. First, a t 1 -incrementation protocol for D-HMQC dubbed Arbitrary Indirect Dwell (AID) is demonstrated. AID allows the use of arbitrary, rotor asynchronous t 1 -increments, but removes the constant time period from CT D-HMQC, resulting in improved sensitivity by reducing transverse relaxation losses. Next, we show that short high-power adiabatic pulses (SHAPs), which efficiently invert broad MAS sideband manifolds, can be effectively incorporated into 1 H detected symmetry-based resonance echo double resonance (S-REDOR) and t 1 -noise eliminated (TONE) D-HMQC experiments. The S-REDOR experiments with SHAPs provide approximately double the dipolar dephasing, as compared to experiments with rectangular inversion pulses. We lastly show that sensitivity and resolution can be further enhanced with the use of swept excitation pulses as well as adiabatic magic angle turning (aMAT).

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

2D HSQC Chemical Shifts of Impurities from Biomass Pretreatment

Two dimensional (2D) heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) is a powerful analytical method which can be used to elucidate the structure of biomass. During processing, biomass is typically subjected to some form of chemical treatment which can be performed with a variety of compounds. The presence of these compounds, even in trace amounts, has the potential to contaminate the sample and lead to misinterpretation of the HSQC spectra. Here we report the chemical shifts of 29 compounds commonly used in biomass processing which have the potential to contaminate the biomass samples and lead to the misinterpretation of peaks associated with biomass (Populus trichocarpa ) pretreated via autohydrolysis. The identification of these chemical shifts could serve as a valuable tool in preventing errors in characterizing biomass via HSQC.

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

Operation of Cs–Sb–O activated GaAs in a high voltage DC electron gun at high average current

Negative Electron Affinity (NEA) activated GaAs photocathodes are the most popular option for generating a high current (>1 mA) spin-polarized electron beam. Despite its popularity, a short operational lifetime is the main drawback of this material. Recent works have shown that the lifetime can be improved by using a robust Cs–Sb–O NEA layer with minimal adverse effects. In this work, we operate GaAs photocathodes with this new activation method in a high voltage environment to extract a high current. We demonstrate that improved chemical resistance of Cs–Sb–O activated GaAs photocathodes allowed them to survive a day-long transport process from a separate vacuum system using a vacuum suitcase. During beam running, we observed spectral dependence on lifetime improvement. In particular, we saw a 45% increase in the lifetime at 780 nm on average for Cs–Sb–O activated GaAs compared to Cs–O activated GaAs.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Large Single Crystals of V 1– x Mo x O 2 from a Two-Step Chemical Vapor Transport Synthesis

This article details the synthesis of single crystals with the chemical formulas V 1-x Mo x O 2 ( 0 ≤ x ≤ 0.6) using the chemical vapor transport method at high temperatures (850 to 950 °C). These compounds are of interest due to the special relationship between their structure and electronic properties. Understanding this relationship will help guide future efforts in designing energy-saving devices. While some of these compounds have been synthesized before, this is the first time their synthesis has been reported for such a wide range of compounds in the form of bulk single crystals. The single crystal form is required for detailed characterization of the compounds. The materials were characterized chemically and structurally in order to validate the precise chemical compositions that were reported. Characterization methods included x-ray diffraction, inductively-coupled plasma atomic emission spectroscopy, and electron dispersive spectroscopy. The synthesis method is tunable, and the introduction of multiple steps both overcome the previous chemical challenges precluding synthesis of high molybdenum content crystals and allow for the creation of larger single crystals that were otherwise possible. The crystals synthesized this way are suitable for specialized measurements for determining the precise physical properties and atomic structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗