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 793 records · Page 44

Improved accuracy and transferability of molecular-orbital-based machine learning: Organics, transition-metal complexes, non-covalent interactions, and transition states

We report molecular-orbital-based machine learning (MOB-ML) provides a general framework for the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. The application of Nesbet’s theorem makes it possible to recast a typical extrapolation task, training on correlation energies for small molecules and predicting correlation energies for large molecules, into an interpolation task based on the properties of orbital pairs. We demonstrate the importance of preserving physical constraints, including invariance conditions and size consistency, when generating the input for the machine learning model. Numerical improvements are demonstrated for different datasets covering total and relative energies for thermally accessible organic and transition-metal containing molecules, non-covalent interactions, and transition-state energies. MOB-ML requires training data from only 1% of the QM7b-T dataset (i.e., only 70 organic molecules with seven and fewer heavy atoms) to predict the total energy of the remaining 99% of this dataset with sub-kcal/mol accuracy. This MOB-ML model is significantly more accurate than other methods when transferred to a dataset comprising of 13 heavy atom molecules, exhibiting no loss of accuracy on a size intensive (i.e., per-electron) basis. It is shown that MOB-ML also works well for extrapolating to transition-state structures, predicting the barrier region for malonaldehyde intramolecular proton-transfer to within 0.35 kcal/mol when only trained on reactant/product-like structures. Finally, the use of the Gaussian process variance enables an active learning strategy for extending the MOB-ML model to new regions of chemical space with minimal effort. We demonstrate this active learning strategy by extending a QM7b-T model to describe non-covalent interactions in the protein backbone–backbone interaction dataset to an accuracy of 0.28 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Selectivity of nitrate and chloride ions in microporous carbons: the role of anisotropic hydration and applied potentials

Understanding ion transport in porous carbons is critical for a wide range of technologies, including supercapacitors and capacitive deionization for water desalination, yet many details remain poorly understood. For instance, an atomistic understanding of how ion selectivity is influenced by the molecular shape of ions, morphology of the micropores and applied voltages is largely lacking. In this work, we combined molecular dynamics simulations with enhanced sampling methods to elucidate the mechanism of nitrate and chloride selectivity in subnanometer graphene slit-pores. We demonstrate that nitrate is preferentially adsorbed over chloride in the slit-like micropores. This preferential adsorption was found to stem from the weaker hydration energy and unique anisotropy of the ion solvation of nitrate. Beside the effects of ion dehydration, we found that applied potential plays an important role in determining the ion selectivity, leading to a lower selectivity of nitrate over chloride at a high applied potential. We conclude that the measured ion selectivity results from a complex interplay between voltage, confinement, and specific ion effects-including ion shape and local hydration structure.

36 MATERIALS SCIENCE↗

Calculated Reduction Potentials of Electrolyte Species in Lithium–Sulfur Batteries

In this work, reduction potentials of electrolyte molecules in the lithium–sulfur (Li/S) battery and their variations in several solvent environments are studied using the density functional theory method with Dunning’s triple-ζ correlation consistent basis set. Reliable reduction potential values are key for electrolyte additive design needed for suppressing polysulfide dissolution shuttle mechanism, resulting in poor cycle performance and severe self-discharge of the Li/S battery. Although isolated electrolyte molecules have reduction potentials outside the operating voltage range of the Li/S battery, complexation with other electrolyte species enables the electrolyte molecules to be reduced within the operating voltage range. Among the electrolyte species considered in this study, bis(fluorosulfonyl)imide (FSI – ) and fluoroethylene carbonate (FEC) yield reduction potentials within the expected range, suggesting the development of fluorine-containing additives as a promising line of research.

25 ENERGY STORAGE↗

Hybrid quantum-classical simulations of magic angle spinning dynamic nuclear polarization in very large spin systems

Solid-state nuclear magnetic resonance can be enhanced using unpaired electron spins with a method known as dynamic nuclear polarization (DNP). Fundamentally, DNP involves ensembles of thousands of spins, a scale that is difficult to match computationally. This scale prevents us from gaining a complete understanding of the spin dynamics and applying simulations to design sample formulations. We recently developed an ab initio model capable of calculating DNP enhancements in systems of up to ~1000 nuclei; however, this scale is insufficient to accurately simulate the dependence of DNP enhancements on radical concentration or magic angle spinning (MAS) frequency. We build on this work by using ab initio simulations to train a hybrid model that makes use of a rate matrix to treat nuclear spin diffusion. We show that this model can reproduce the MAS rate and concentration dependence of DNP enhancements and build-up time constants. Finally, we then apply it to predict the DNP enhancements in core–shell metal-organic-framework nanoparticles and reveal new insights into the composition of the particles’ shells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Removal of Homogeneous Broadening from 1 H-Detected Multidimensional Solid-State NMR Spectra

1 H-detected magic-angle spinning (MAS) NMR experiments have revolutionized the NMR studies of biological and inorganic solids by providing unparalleled sensitivity and resolution. Despite these gains, homogeneous broadening, originating from the incomplete removal of homonuclear dipolar interactions under fast MAS, remains highly prevalent and limits the achievable resolution. In direct analogy to super-resolution microscopy methods, we show that resolution beyond that currently achievable by fast MAS alone can be obtained by experiment-driven deconvolution. Following the acquisition of a single 2D NMR spectrum to measure the frequency-dependent homogeneous lineshapes, any number of 1 H-detected spectra can be enhanced in resolution, yielding comparable spectra as obtained with twice the MAS frequency. In conclusion, the versatility of this approach is demonstrated in the enhancement of single- and double-quantum homonuclear correlation spectra, in addition to heteronuclear correlation spectra acquired on a surface organometallic complex and the protein GB1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Furthermore, fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pyroelectricity of Lead Sulfide (PbS) Quantum Dot Films Induced by Janus-Ligand Shells

Asymmetry is an essential property to control. To do that in nanocrystalline systems we have developed methods to produce Janus-ligand shells on otherwise symmetric PbS quantum dots (QDs). In this work, we demonstrate that control by constructing a system that exhibits pyroelectricity built from spherical PbS QDs. We observed a pyroelectric current in two different configurations. In one configuration, the QDs are self-assembled into close-packed arrays while in the second configuration, the QDs are dispersed into an electro-inactive polymer, polydimethylsiloxane. Both exhibit a pyroelectric response. In the first configuration we estimate a lower limit of the pyroelectric coefficient to be 1.97 × 10 -7 C/m 2 K, which is likely limited by the degree of QD alignment during film formation but is already on par with common pyroelectric systems. Compared with inorganic ceramic-like and polymeric pyroelectric materials, pyroelectric films self-assembled from polar QDs are easier to prepare, responsive to light with different energies based on QD exciton energy, and the polarization of each QD could be easily tuned by constructing different Janus-ligand shells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surrogate Hessian accelerated structural optimization for stochastic electronic structure theories

In this work, we present an efficient energy-based method for structural optimization with stochastic electronic structure theories, such as diffusion quantum Monte Carlo (DMC). This method is based on robust line-search energy minimization in reduced parameter space, exploiting approximate but accurate Hessian information from a surrogate theory, such as density functional theory. The surrogate theory is also used to characterize the potential energy surface, allowing for simple but reliable ways to maximize statistical efficiency while retaining controllable accuracy. We demonstrate the method by finding the minimum DMC energy structures of the selected flake-like aromatic molecules, such as benzene, coronene, and ovalene, represented by 2, 6, and 19 structural parameters, respectively. In each case, the energy minimum is found within two parallel line-search iterations. The method is near-optimal for a line-search technique and suitable for a broad range of applications. It is easily generalized to any electronic structure method where forces and stresses are still under active development and implementation, such as diffusion Monte Carlo, auxiliary-field Monte Carlo, and stochastic configuration interaction, as well as deterministic approaches such as the random-phase approximation. Accurate and efficient means of geometry optimization could shed light on a broad class of materials and molecules, showing high sensitivity of induced properties to structural variables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The kagomé metals RbTi 3 Bi 5 and CsTi 3 Bi 5

Abstract The kagomé metals RbTi 3 Bi 5 and CsTi 3 Bi 5 were synthesized both as polycrystalline powders by heating the elements in an argon atmosphere and as single crystals grown using a self-flux method. The compounds crystallize in the hexagonal crystal system isotypically to KV 3 Sb 5 ( P 6/ mmm , Z = 1, CsTi 3 Bi 5 : a = 5.7873(1), c = 9.2062(1) Å; RbTi 3 Bi 5 : a = 5.773(1), c = 9.065(1) Å). The titanium atoms form a kagomé net with bismuth atoms in the hexagons as well as above and below the triangles. The alkali metal atoms are coordinated by 12 bismuth atoms and form AlB 2 -like slabs between the kagomé layers. Magnetic susceptibility measurements with CsTi 3 Bi 5 and RbTi 3 Bi 5 single crystals reveal Pauli-paramagnetism and traces of superconductivity caused by CsBi 2 /RbBi 2 impurities. Magnetotransport measurements reveal conventional Fermi liquid behavior and quantum oscillations indicative of a single dominant orbit at low temperature. DFT calculations show the characteristic metallic kagomé band structure similar to that of CsV 3 Sb 5 with reduced band filling. A symmetry analysis of the band structure does not reveal an obvious and unique signature of a nontrivial topology.

Chemistry↗

Quantum Time Dynamics Mediated by the Yang–Baxter Equation and Artificial Neural Networks

Quantum computing shows great potential, but errors pose a significant challenge. This study explores new strategies for mitigating quantum errors using artificial neural networks (ANNs) and the Yang–Baxter equation (YBE). Unlike traditional error mitigation methods, which are computationally intensive, we investigate artificial error mitigation. We developed a novel method that combines ANNs for noise mitigation combined with the YBE to generate noisy data. This approach effectively reduces noise in quantum simulations, enhancing the accuracy of the results. The YBE rigorously preserves quantum correlations and symmetries in spin chain simulations in certain classes of integrable lattice models, enabling effective compression of quantum circuits while retaining linear scalability with the number of qubits. This compression facilitates both full and partial implementations, allowing the generation of noisy quantum data on hardware alongside noiseless simulations using classical platforms. By introducing controlled noise through the YBE, we enhance the data set for error mitigation. We train an ANN model on partial data from quantum simulations, demonstrating its effectiveness in mitigating errors in time-evolving quantum states, providing a scalable framework to enhance quantum computation fidelity, particularly in noisy intermediate-scale quantum (NISQ) systems. We demonstrate the efficacy of this approach by performing quantum time dynamics simulations using the Heisenberg XY Hamiltonian on real quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Substituent effect on napthodithiophene-fused porphyrins: Understanding the unusual trend of fluorescence quantum yield quenching

A series of π-extended porphyrins containing thiophene units was prepared to investigate an unusual trend observed in the fluorescence quantum yield of naphtho[2,1-b:3,4-b']dithiophene-fused porphyrins. Monobenzoporphyrins carrying 2-thiophenyl (vinyl thiophene porphyrin) groups (2VTP, 2VTBr and Br2VTP) (numbering 2 and 3 refers to the position of sulfur on the thiophene ring) were synthesized through a Heck-based cascade reaction followed by a newly developed bromination method. These were converted to naphtho[2,1-b:3,4-b']dithiophene-fused porphyrin derivatives (FBr2VTP and F2VTBr) via ring closure with an intramolecular Scholl reaction. Increasingly larger Stokes shifts with a higher number of bromo groups were observed in the unfused π-extended molecular systems, reflecting the heavy atom effect. Fluorescence spectroscopy further confirmed the unusual trend seen in the previous work: structural rigidification in naphtho[2,1-b:3,4-b']dithiophene-fused porphyrins leads to longer fluorescence lifetimes but unexpectedly lowers the quantum yield. Adding one bromo group to the naphtho[2,1-b:3,4-b']dithiophene unit does not change this trend. Furthermore, the presence of two bromo groups prevents the quantum yield from dropping. DFT, TDDFT, and NICS analyses suggest a drastic change in aromaticity in the pyrrole ring where the naphtho[2,1-b:3,4-b']dithiophene is fused in FBr2VTP, which might contribute to the quantum yield reductions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Compositionally Tuning Electron Transfer from Photoexcited Core/Shell Quantum Dots via Cation Exchange

It is critical to find methods to control the thermodynamic driving force for photoexcited charge transfer from quantum dots (QDs) and explore how this affects charge transfer rates, since the efficiency of QD-based photovoltaic and photocatalysis technologies depends on both this rate and the associated energetic losses. In this work, we introduce a single-pot shell growth and Cu-catalyzed cation exchange method to synthesize Cd x Zn 1-x Se/Cd y Zn 1-y S QDs with tunable driving forces for electron transfer. Functionalizing them with two molecular electron acceptors—naphthalenediimide (NDI) and anthraquinone (AQ)—allowed us to probe nearly 1 eV of driving forces. For AQ, at lower driving forces, we find that higher Zn content results in a 130-fold increase of electron transfer rate constants. However, at higher driving forces electron transfer dynamics are unaltered. Here, the data are understood using an Auger-assisted electron transfer model and analyzed with computational work to determine approximate binding geometries of these electron acceptors. Our work provides a method to tune QD reducing power and produces useful metrics for optimizing QD charge transfer systems that maximize rates of electron transfer while minimizing energetic losses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PySpawn: Software for Nonadiabatic Quantum Molecular Dynamics

The ab initio multiple spawning (AIMS) method enables nonadiabatic quantum molecular dynamics simulations in an arbitrary number of dimensions, with potential energy surfaces provided by electronic structure calculations performed on-the-fly. However, the intricacy of the AIMS algorithm complicates software development, deployment on modern shared computer resources, and post-simulation data analysis. PySpawn is a nonadiabatic molecular dynamics software package that addresses these issues. Here, the program is designed to be easily interfaced with electronic structure software, and an interface to the TeraChem software package is described here. PySpawn introduces a task-based reorganization of the AIMS algorithm, allowing fine-grained restart capability and setting the stage for efficient parallelization in a future release. PySpawn includes a user-friendly and interactive Python analysis module that will enable novice users to painlessly adopt AIMS. As a demonstration of PySpawn’s simulation capability and analysis module, we report complete active space self-consistent field–based AIMS simulations of the 1,2- dithienyl-1,2-dicyanoethene molecule, a promising molecular photoswitch.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phase transitions and dimensional cross-over in layered confined solids

The nature of solid phases and cross-over of order–disorder phase transitions from two-dimensional (2D) layers to three-dimensional (3D) bulk in confined atomic systems remain largely unexplained. To this end, we consider noble gases and aluminum confined between graphene sheets at different pressures and temperatures. Using crystal structure search methods and molecular dynamics based on machine-learned potentials with quantum-mechanical accuracy, we identify structures of multilayer confined solids that deviate from simple close packing. Upon heating, we find that confined 2D monolayers melt according to the two-step continuous Kosterlitz–Thouless–Halperin–Nelson–Young theory. However, multilayer solids transition continuously into an intermediate layered-hexatic phase before melting discontinuously into an isotropic liquid. This intermediate phase persists at least up to 12 layers studied here. This change can be qualitatively understood based on the cross-over from 2D topological defects toward 3D ones during melting as the number of layers increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advances in Molecular Beam Epitaxy Growth of Ultra-Wide Bandgap Ga2O3 Based Alloys

Gallium oxide (Ga2O3) is an emerging ultra-wide bandgap semiconductor material that has attracted attention for its potential to outperform existing SiC and GaN based devices operating at high breakdown voltages and high temperature. Isovalent alloying of In and Al in Ga2O3 provides the ability to engineer bandgap energy and strain of the material. Alloying with Al increases the bandgap energy and the theoretically achievable Baliga's figure of merit, a key measure of a material's ultimate performance limits for high power switching devices. Alloying with In introduces compressive strain and can be used to counteract the tensile strain of Al incorporation. The resulting (AlxGa1-x-yIny)2O3 alloy can be lattice-matched to commercially available Ga2O3 wafers and has a tunable bandgap energy greater than that of Ga2O3, 4.76 eV. Such lattice-matched material can be grown arbitrarily thick without the detrimental effects of elastic strain and relaxation, making it suitable for high voltage diodes and transistors. However, efforts to synthesize isovalent alloys are complicated by their tendency to phase separate into corundum Al2O3 or bixbyite In2O3. Literature reports of the quaternary (AlxGa1-x-yIny)2O3 are limited to <1% unintentional indium incorporation in In-catalyzed (AlxGa1-x)2O3. The primary limitation to quaternary growth is the limited incorporation of indium at elevated growth temperatures. This limited incorporation is due to both the volatility of indium oxide and Al and Ga cation exchange reactions which replace indium in In2O3. We report on the development of a novel high-throughput molecular beam epitaxy (MBE) technique to screen the growth conditions for the ternary alloy (InyGa1-y)2O3, and the application of these findings to the first successful synthesis of phase pure monoclinic (AlxGa1-x-yIny)2O3 by MBE. By leveraging the unique sub-oxide chemistry of Ga2O3 and in-situ monitoring of crystal properties by reflection high-energy electron diffraction (RHEED), a cyclical growth and etch-back method is developed and applied to rapidly characterize the (InyGa1-y)2O3 growth space. This cyclical method provides approximately 10x increase in experimental throughput and up to 46x improvement in Ga2O3 substrate utilization. Appropriate growth conditions for monoclinic (InyGa1- y)2O3 are identified by machine learning analysis of RHEED patterns and targeted growths are characterized ex-situ to confirm improved In incorporation. These growth conditions are then combined with established (AlxGa1-x)2O3 growth conditions to grow quaternary (AlxGa1-x-yIny)2O3 with Al mole fractions ranging from 1.4% - 24.4% and In mole fractions ranging from 3.1% to 15.5%. The chemical and optical properties of the alloys are investigated by XRD, XPS, and spectroscopic ellipsometry. A lattice-matched (AlxGa1-x-yIny)2O3 alloy is examined by 4D-STEM and the chemical and physical uniformity of Al and In incorporation are discussed.

alloy↗

Raman Spectroscopy Analysis of the Structure and Surface Chemistry of Ti3C2Tx MXene

Raman spectroscopy is one of the most useful tools for the analysis of two-dimensional (2D) materials. While MXenes are a very large family of 2D transition metal carbides and nitrides, there have been just a few Raman studies of materials from this family. Here, we report on a systematic study of the most widely used and most important MXene to date: Ti 3 C 2 T x . By synthesizing material using different methods, we show that Raman spectra of Ti 3 C 2 T x are affected not only by the composition and surface groups but also by intercalated species and stacking. Due to a plasmonic peak of Ti 3 C 2 T x around 785 nm, resonant conditions are achieved, enabling us to observe an extra peak at ~120 cm –1 , when excited with a red diode laser. We report differences in Raman spectra collected from single flakes of Ti 3 C 2 T x , colloidal solutions, and multilayer films. Lastly, we show how an undesirable photoluminescent background could serve as evidence of material degradation, which leads to the formation of defective titania and amorphous carbon. Overall, this study shows how Raman spectroscopy can be used for the characterization of important emerging 2D materials: MXenes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. Here, a vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗