Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum chemical computation”

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 631 records · Page 35

Coarse-Grained Density Functional Theory Predictions via Deep Kernel Learning

Scalable electronic predictions are critical for soft materials design. Recently, the Electronic Coarse-Graining (ECG) method was introduced to renormalize all-atom quantum chemical (QC) predictions to coarse-grained (CG) resolutions using deep neural networks (DNNs). While DNNs can learn complex representations that prove challenging for kernel-based methods, they are susceptible to overfitting and the overconfidence of uncertainty estimations. Here, we develop ECG within a GPU-accelerated Deep Kernel Learning (DKL) framework to enable CG QC predictions using range-separated hybrid density functional theory (DFT), obtaining a 107 speedup relative to naive all-atom QC. By treating the predicted electronic properties as random Gaussian Processes, DKL incorporates CG mapping degeneracy by learning the distribution of electronic energies as a function of CG configuration. DKL-ECG accurately reproduces molecular orbital energies from range-separated DFT while facilitating efficient training via active learning using the uncertainties provided by DKL. Further, we show that while active learning algorithms enable efficient sampling of a more diverse configurational space relative to random sampling, all explored query methods exhibit comparable performance for the examined system. We attribute this result to the significant overlap of the feature space and output property distributions across multiple temperatures.

97 MATHEMATICS AND COMPUTING↗

Zero-Field NMR and Millitesla-SLIC Spectra for >200 Molecules from Density Functional Theory and Spin Dynamics

NMR is usually performed at magnetic fields of 1 T and above to obtain sufficient sensitivity and spectral dispersion to identify chemicals based on chemical shifts and J couplings. At lower fields, the advent of hyperpolarization technologies and sensitive detectors can address sensitivity concerns. However, it remains disputed whether spectral signatures at zero and ultra-low fields are sufficient for chemical identification. Here, we report an all–electron DFT-based batch calculation of J-coupling constants, which are used to generate J coupling NMR spectra at zero field and 6.5 mT for over 200 small molecules. In the developed computational tool chain, we first used the all-electron FHI-aims code to calculate the molecular J couplings and chemical shifts. We then fed the calculated NMR parameters into the NMR simulation package SPINACH to simulate both heteronuclear J coupling spectra at zero-field, and homonuclear J coupling spectra as spin-lock induced crossing (SLIC) spectra at ultra-low field (6.5 mT). The resulting spectra demonstrate that zero and ultra-low field NMR spectra can represent unique identifiers of chemical structure for small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Origins of Conformational Heterogeneity in Peptoid Helices Formed by Chiral N -1-Phenylethyl Sidechains

N-substituted glycines (polypeptoids) containing chiral hydrophobic sidechains are known to fold into biomimetic alpha helices. These helix formers often produce conformationally heterogeneous structures and are difficult to characterize at a sub-nanometer resolution. Previously, peptoid N-1-phenylethyl (S)-enantiomer sidechains (Nspe) were inferred from various experiments to form right-handed helices and (R)-enantiomers (Nrpe), left-handed helices. Prior computational work for N(s/r)pe oligomers has struggled to reproduce this trend. Herein, quantum mechanics calculations and molecular dynamics simulations are used to understand the origins of this discrepancy. Results from DFT and molecular mechanics calculations on a variety of Nspe and Nrpe oligomers as a function of chain length are in agreement, showing that Nspe and Nrpe prefer left- and right-handed helices, respectively. Additional metadynamics simulations are used to study Nrpe and Nspe oligomers folding in water. These results show that the free-energy driving forces for assembly into a helical backbone configuration are very small (within ~k B T). Lastly, we compare DFT calculations for other experimentally characterized peptoid sidechains, N(r/s)sb, N(r/s)tbe, and N(r/s)npe. In this analysis, we show that peptoid sidechains determined to be more robust experimentally (tbe and npe) have helical preferences opposite the trend seen in less robust assemblies formed by N(r/s)pe and N(r/s)sb chemistries. Furthermore, the more robust tbe and nnpe favor the (S)-enantiomer to right-handed and the (R)-enantiomers to left-handed helices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scaling of silicon nanoparticle growth in low temperature flowing plasmas

Low temperature plasmas are an emerging method to synthesize high quality nanoparticles (NPs). An established and successful technique to produce NPs is using a capacitively coupled plasma (CCP) in cylindrical geometry. Although a robust synthesis technique, optimizing or specifying NP properties using CCPs, is challenging. In this paper, results from a computational investigation for the growth of silicon NPs in flowing inductively coupled plasmas (ICPs) using Ar/SiH 4 gas mixtures of up to a few Torr are discussed. ICPs produce more locally constrained and quiescent plasma potentials. These positive plasma potentials produce an electrostatic trap for negatively charged NPs, which can significantly extend the residence time of NPs in the plasma, which in turn provides a controllable period for particle growth. The computational platforms used in this study consist of a two-dimensional plasma hydrodynamics model, a three-dimensional nanoparticle growth and trajectory tracking model, and a molecular dynamics simulation for deriving reactive sticking coefficients of silane radicals on Si NPs. Trends for the nanoparticle growth as a function of SiH 4 inlet fraction, gas residence time, energy deposition per particle, pressure, and reactor diameter are discussed. The general path for particle synthesis is the trapping of small NPs in the positive electrostatic potential, followed by entrainment in the gas flow upon reaching a critical particle size. Optimizing or controlling NP synthesis then depends on the spatial distribution of plasma potential, the density of growth species, and the relative time that particles spend in the electrostatic trap and flowing through higher densities of growth species upon leaving the trap.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mechanistic insights into rubidium ion adsorption at the quartz (101) surface from quantum chemical metadynamics

Ion sorption extent and mechanism depend in part on the mineral surface termination, which can be highly complex. Variations in surface functional groups, particularly with defect density and surface roughness, influence mineral reactivity towards solutes. In this work, we investigate the adsorption of a rubidium cation (Rb + ) at pristine and defect quartz (101) surface sites using well-tempered metadynamics, based on simulations with the quantum chemical density-functional tight-binding (DFTB) method. We compare the relative energetics of Rb + adsorption across selected sites for each surface, with nanosecond-level sampling, highlighting similarities between vicinal and geminal silanol sites. We find that the positive Rb + partial atomic charge can increase by as much as 0.5e as it approaches the surface, with implications for modulation of ion adsorption strength and extent at the quartz (101) surface with surface vacancies, silanol coverage, and charge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Automated Development of Molten Salt Machine Learning Potentials: Application to LiCl

The in silico modeling of molten salts is critical for emerging "carbon-free" energy applications but is inhibited by the cost of quantum mechanically treating the high polarizabilities of molten salts. Here, we integrate configurational sampling using classical force fields with active learning to automate and accelerate the generation of Gaussian approximation potentials (GAP) for molten salts. This methodology reduces the number of expensive ab initio evaluations required for training set generation to O(100), enabling the facile parametrization of a molten LiCl GAP model that exhibits a 19 000-fold speedup relative to AIMD. The developed molten LiCl GAP model is applied to sample extended spatiotemporal scales, permitting new physical insights into molten LiCl's coordination structure as well as experimentally validated predictions of structures, densities, self-diffusion constants, and ionic conductivities. The developed methodology significantly lowers the barrier to the in silico understanding and design of molten salts across the periodic table.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Realistic Ion Dynamics through Charge Renormalization in Nonaqueous Electrolytes

While many practically important electrolytes contain lithium ions, interactions of these ions are particularly difficult to probe experimentally because of their small X-ray and neutron scattering cross sections and large neutron absorption cross sections. Molecular dynamics (MD) is a powerful tool for understanding the properties of nonaqueous electrolyte solutions from the atomic level, but the accuracy of this computational method crucially depends on the physics built into the classical force field. In this study, we demonstrate that several force fields for lithium bistriflimide (LiTFSI) in acetonitrile yield a solution structure that is consistent with the neutron scattering experiments, yet these models produce dramatically different ion dynamics in solution. Such glaring discrepancies indicate that inadequate representation of long-range interactions leads to excessive ionic association and ion-pair clustering. We show that reasonable agreement with the experimental observations can be achieved by renormalization of the ion charges using a “titration” method suggested herewith. This simple modification produces realistic concentration dependencies for ionic diffusion and conductivity in <2 M solutions, without loss in quality for simulation of the structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum Davidson algorithm for excited states

Abstract Excited state properties play a pivotal role in various chemical and physical phenomena, such as charge separation and light emission. However, the primary focus of most existing quantum algorithms has been the ground state, as seen in quantum phase estimation and the variational quantum eigensolver (VQE). Although VQE-type methods have been extended to explore excited states, these methods grapple with optimization challenges. In contrast, the quantum Krylov subspace (QKS) method has been introduced to address both ground and excited states, positioning itself as a cost-effective alternative to quantum phase estimation. However, conventional QKS methodologies depend on a pre-generated subspace through real or imaginary-time evolutions. This subspace is inherently expansive and can be plagued with issues like slow convergence or numerical instabilities, often leading to relatively deep circuits. Our research presents an economic QKS algorithm, which we term the quantum Davidson (QDavidson) algorithm. This innovation hinges on the iterative expansion of the Krylov subspace and the incorporation of a pre-conditioner within the Davidson framework. By using the residues of eigenstates to expand the Krylov subspace, we manage to formulate a compact subspace that aligns closely with the exact solutions. This iterative subspace expansion paves the way for a more rapid convergence in comparison to other QKS techniques, such as the quantum Lanczos. Using quantum simulators, we employ the novel QDavidson algorithm to delve into the excited state properties of various systems, spanning from the Heisenberg spin model to real molecules. Compared to the existing QKS methods, the QDavidson algorithm not only converges swiftly but also demands a significantly shallower circuit. This efficiency establishes the QDavidson method as a pragmatic tool for elucidating both ground and excited state properties on quantum computing platforms.

97 MATHEMATICS AND COMPUTING↗

Physical patterning of high-Q superconducting niobium resonators via ion beam etching

The development of superconducting quantum circuits increasingly involves the exploration of chemically distinct materials and complex multilayered structures. Accelerating this trend may benefit from low-damage, materials-agnostic patterning techniques that are compatible with a broad range of materials. Here, in this work, we investigate the utility of low-energy ion beam etching (IBE), a physical patterning technique, as an alternative to reactive ion etching for fabricating low-loss superconducting resonators. We use niobium (Nb) resonators as a test platform, leveraging their well-characterized performance metrics for benchmarking. To address IBE-induced surface redeposition, we introduce an in situ aluminum capping layer combined with targeted post-fabrication chemical treatment. This strategy yields resonators with internal quality factors as high as 6 × 10 5 in the single-photon regime at 50 mK. These results establish low-energy IBE as a promising patterning technique for superconducting devices, with the potential to accelerate development across chemically diverse and multilayered material platforms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pulsed-laser epitaxy of metallic delafossite PdCrO 2 films

Alternate stacking of a highly conducting metallic layer with a magnetic triangular layer found in delafossite PdCrO 2 provides an excellent platform for discovering intriguing correlated quantum phenomena. Thin film growth of delafossites may enable not only the tuning of the basic physical properties beyond what bulk materials can exhibit, but also the development of novel hybrid materials by interfacing with dissimilar materials, yet this has proven to be extremely challenging. Here, we report the epitaxial growth of metallic delafossite PdCrO 2 films by pulsed laser epitaxy (PLE). The fundamental role of the PLE growth conditions, epitaxial strain, and chemical and structural characteristics of the substrate is investigated by growing under various growth conditions and on various types of substrates. While strain plays a large role in improving the crystallinity, the direct growth of epitaxial PdCrO 2 films without impurity phases was not successful. We attribute this difficulty to both the chemical and structural dissimilarities with the substrate and volatile nature of the PdO sublayer, which make nucleation of the right phase difficult. This difficulty was overcome by growing CuCrO 2 buffer layers before PdCrO 2 films were grown. Unlike PdCrO 2 , CuCrO 2 films were readily grown with a relatively wide growth window. Only a monolayer thick buffer layer was sufficient to grow the correct PdCrO 2 phase. This result indicates that the epitaxy of Pd-based delafossites is extremely sensitive to the chemistry and structure of the interface, necessitating near perfect substrate materials. The resulting films are commensurately strained and show an antiferromagnetic transition at 40 K that persists down to as thin as 3.6 nm in thickness. This work provides key insights into advancing the epitaxial growth of the broader class of metallic delafossites for both studying the basic physical properties and developing new spintronic and computing devices.

36 MATERIALS SCIENCE↗

Room Temperature Lasing from Semiconducting Single-Walled Carbon Nanotubes

Miniaturized near-infrared semiconductor lasers that are able to generate coherent light with low energy consumption have widespread applications in fields such as optical interconnects, neuromorphic computing, and deep-tissue optogenetics. With optical transitions at near-infrared wavelengths, diameter-tunable electronic structures, and superlative optoelectronic properties, semiconducting single-walled carbon nanotubes (SWCNTs) are promising candidates for nanolaser applications. However, despite significant efforts in this direction and recent progress toward enhancing spontaneous emission from SWCNTs through Purcell effects, SWCNT-based excitonic lasers have not yet been demonstrated. Leveraging an optimized cavity-emitter integration scheme enabled by a self-assembly process, here we couple SWCNT emission to the whispering gallery modes supported by polymer microspheres, resulting in room temperature excitonic lasing with an average lasing threshold of 4.5 kW/cm 2 . The high photostability of SWCNTs allows stable lasing for prolonged duration with minimal degradation. Finally, this experimental realization of excitonic lasing from SWCNTs, combined with their versatile electronic and optical properties that can be further controlled by chemical modification, offers far-reaching opportunities for tunable near-infrared nanolasers that are applicable for optical signal processing, in vivo biosensing, and optoelectronic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stoichiometrically-informed symbolic regression for extracting chemical reaction mechanisms from data

A data-driven computational method is introduced to extract chemical reaction mechanisms from time series chemical concentration data. It is realized through the use of dynamic symbolic regression in which a sparse analytical form for a dynamical system is discoverable from the underlying data. We specifically develop the stoichiometrically-informed symbolic regression (SISR) method to address a standing challenge in complex chemical reaction networks: given a time-series dataset of concentrations of several components, what is the mechanism and the associated rate constants? SISR finds the optimal mechanism, kinetic equations and rate constants by combining differential optimization with a genetic optimization approach that searches a symbolic space of possible reaction mechanisms. Use of SISR in several paradigmatic examples spanning linear and nonlinear reaction schemes results in excellent agreement between true and predicted mechanisms, including when the method is applied to noisy data. The advantages of a stoichiometrically-informed approach such as SISR to address reaction discovery is illustrated through comparison with the use of generic state-of-the-art data-driven approaches.

36 MATERIALS SCIENCE↗

Study of Self-Interaction Errors in Density Functional Calculations of Magnetic Exchange Coupling Constants Using Three Self-Interaction Correction Methods

We examine the role of self-interaction errors (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular we analyze the effect of scaling down the self-interaction-correction (SIC) for three non-empirical density functional approximations (DFAs) namely, the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods: Perdew-Zunger [Perdew, J. P.; Zunger, A. Phys. Rev. B, 1981, 23, 5048] SIC, the orbitalwise scaled SIC method [Vydrov, O. A. et al., J. Chem. Phys. 2006, 124, 094108], and the recent local scaling method [Zope, R. R. et al., J. Chem. Phys. 2019, 151, 214108]. We compute the magnetic exchange coupling constants using the spin projection and non projection approaches for sets of molecules composed of dinuclear and polynuclear H-He models, organic radical molecules, and chlorocuprate, and compare these results against accurate theories and experiment. Our results show that for the systems that mainly consist of single electron regions, PZSIC performs well but for more complex organic systems and the chlorcuprates, an overcorrecting tendency of PZSIC combined with the DFAs utilized in this work is more pronounced, and in such cases LSIC with kinetic energy density ratio performs better than PZSIC. Furthermore, analysis of the results in terms of SIC corrections to the density and to the total energy shows that both density and energy correction are required to obtain an improved prediction of magnetic exchange couplings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computer-aided design of a proton pump

The use of transmembrane proton gradients in energy transduction is an almost universal feature of life on earth. These proton gradients are established and maintained by specialized assemblies of proteins which actively pump protons across membranes. One broad class of proton pumps uses captured light energy to drive the proton pumping. Our goal is to elucidate the minimum structural requirements of a light-driven proton-pump. There are two basic components to a simple light-driven proton pump: a source of photo-generated protons and a "gate-keeper" which prevents these protons from reattaching themselves to their source. A wide variety of molecules in the membrane, even as simple as polycyclic aromatic hydrocarbons, are capable of releasing protons when illuminated. Our work is therefore focused on the design of the "gate-keeper." Our initial model involves a pair of proton acceptors, coupled to each other by a transient water bridge, and supported in the membrane by a small bundle of peptide helices. Upon illumination, the proton source transfers its proton to the:- first acceptor of the gate-keeper. While the reverse reaction is highly probable, all that is needed to ensure irreversibility is a nonvanishing probability that the proton will be transferred to the second acceptor across a transient water bridge. Back transfer of the proton to the first acceptor, and thence to the proton source, is impeded by the free energy required to move the proton uphill towards the. proton source and by the disruption of the transient water bridge. As a prototypical water-bridged proton transfer system, we are studying the transfer of a proton across a water bridge from a formic acid to a formate anion. With a pK(sub alpha), of 3.7. formic acid is a good model for the acidic amino acids glutamate and aspartate which are good candidates for gate-keeper proton acceptors. Simulations of proton transfer reactions in a membrane are complicated by the quantum mechanical nature of the process breaking and forming chemical bonds.

New, Michael H.↗

Formation and Microwave Losses of Hydrides in Superconducting Niobium Thin Films Resulting from Fluoride Chemical Processing

Abstract Superconducting niobium (Nb) thin films have recently attracted significant attention due to their utility for quantum information technologies. In the processing of Nb thin films, fluoride‐based chemical etchants are commonly used to remove surface oxides that are known to affect superconducting quantum devices adversely. However, these same etchants can also introduce hydrogen to form Nb hydrides, potentially negatively impacting microwave loss performance. Here, comprehensive materials characterization of Nb hydrides formed in Nb thin films as a function of fluoride chemical treatments is presented. In particular, secondary‐ion mass spectrometry, X‐ray scattering, and transmission electron microscopy reveal the spatial distribution and phase transformation of Nb hydrides. The rate of hydride formation is determined by the fluoride solution acidity and the etch rate of Nb 2 O 5, which acts as a diffusion barrier for hydrogen into Nb. The resulting Nb hydrides are detrimental to Nb superconducting properties and lead to increased power‐independent microwave loss in coplanar waveguide resonators. However, Nb hydrides do not correlate with two‐level system loss or device aging mechanisms. Overall, this work provides insight into the formation of Nb hydrides and their role in microwave loss, thus guiding ongoing efforts to maximize coherence time in superconducting quantum devices.

36 MATERIALS SCIENCE↗

GeSn-on-GaAs with photoconductive carrier lifetime >400 ns: role of substrate orientation and atomistic simulation

Group IV GeSn quantum material finds application in electronics and silicon-compatible photonics. Synthesizing these materials with low defect density and high carrier lifetime is a potential challenge due to lattice mismatch induced defects and tin segregation at higher growth temperature. Recent advancements in the growth of these GeSn materials on Si, Ge, GaAs, and with substrate orientations, demonstrated different properties using epitaxial and chemical deposition methods. This article addresses the effect of GaAs substrate orientation and misorientation on the materials’ properties and carrier lifetimes in epitaxial Ge 0.94 Sn 0.06 layers. With starting GaAs substrates of (100)/2°, (100)/6°, (110) and (111)A orientations, Ge 0.94 Sn 0.06 epitaxial layers were grown with an intermediate Ge buffer layer by molecular beam epitaxy and analyzed by several analytical tools. X-ray analysis displayed good crystalline quality, and Raman spectroscopy measurements showed blue shifts in phonon wavenumber due to biaxial compressive strain in Ge 0.94 Sn 0.06 epilayers. Cross-sectional transmission electron microscopy analysis confirmed the defect-free heterointerface of Ge 0.94 Sn 0.06 /Ge/GaAs heterostructure. Minority carrier lifetimes of the unintentionally doped n-type Ge 0.94 Sn 0.06 epilayers displayed photoconductive carrier lifetimes of >400 ns on (100)/6°, 319 ns on (100)/2°, and 434 ns on (110) GaAs substrate at 1500 nm excitation wavelength. On the other hand, Ge 0.94 Sn 0.06 layer showed poor carrier lifetime on (111)A GaAs substrate. The observed differences in carrier lifetimes were correlated with the formation energy of the Ge on (100)/6° and (100)/2° GaAs heterointerface using Stillinger–Weber interatomic potential model-based atomistic simulation with different heterointerfacial bonding by Synopsys QuantumATK tool. Total energy computation of 6280-atom Ge/GaAs supercell on (100)/6° leads to lower formation energy than (100)/2°, making it more thermodynamically stable. Hence, the growth of the GeSn/III–V material system using misoriented (100) substrates that are more thermodynamically stable will enhance the performances of optoelectronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spatial Interactions in Hydrogenated Perovskite Nickelate Synaptic Networks

A key aspect of how the brain learns and enables decision-making processes is through synaptic interactions. Electrical transmission and communication in a network of synapses are modulated by extracellular fields generated by ionic chemical gradients. Emulating such spatial interactions in synthetic networks can be of potential use for neuromorphic learning and the hardware implementation of artificial intelligence. Here, in this work, we demonstrate that in a network of hydrogen-doped perovskite nickelate devices, electric bias across a single junction can tune the coupling strength between the neighboring cells. Electrical transport measurements and spatially resolved diffraction and nanoprobe X-ray and scanning microwave impedance spectroscopic studies suggest that graded proton distribution in the inhomogeneous medium of hydrogen-doped nickelate film enables this behavior. We further demonstrate signal integration through the coupling of various junctions.

36 MATERIALS SCIENCE↗

Spin–orbit configuration interaction study of spectral properties of PbO

Relativistic calculations of the structural and spectral properties of the PbO molecule can provide fundamental information about the importance of a proper treatment of angular momentum coupling among electrons in order to achieve accurate computational results for spectral properties. Specifically, the nature of these couplings in PbO is expected to be intermediate between the LS- and jj-coupling limits because of its light/heavy element composition. This article reports potential energy curves, transition energies, electric dipole transition moments, permanent dipole moments and spectroscopic constants of PbO calculated using a multireference single plus double excitations spin–orbit configuration interaction approach in the context of relativistic effective core potentials and their concomitant spin–orbit coupling operators. The calculated results are in general agreement with both available experimental results as well as earlier calculations. New values for properties of excited states are also reported. It is noteworthy that certain properties show larger deviations from previous calculations. Furthermore, these deviations are attributed to direct and indirect relativistic effects resulting from diatomic electron–electron angular momentum coupling effects, which are included consistently in the calculations reported herein.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗