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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Self-interaction-corrected Kohn–Sham effective potentials using the density-consistent effective potential method

Density functional theory (DFT) and beyond-DFT methods are often used in combination with photoelectron spectroscopy to obtain physical insights into the electronic structure of molecules and solids. The Kohn–Sham eigenvalues are not electron removal energies except for the highest occupied orbital. The eigenvalues of the highest occupied molecular orbitals often underestimate the electron removal or ionization energies due to the self-interaction (SI) errors in approximate density functionals. Here, we adapt and implement the density-consistent effective potential method of Kohut, Ryabinkin, and Staroverov (2014) to obtain SI-corrected local effective potentials from the SI-corrected Fermi–Löwdin orbitals and density in the Fermi–Löwdin orbital self-interaction correction scheme. The implementation is used to obtain the density of states (photoelectron spectra) and HOMO–LUMO gaps for a set of molecules and polyacenes. Good agreement with experimental values is obtained compared to a range of SI uncorrected density functional approximations.

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

Coupled charge and energy transfer dynamics in light harvesting complexes from a hybrid hierarchical equations of motion approach

Here, we describe a method for simulating exciton dynamics in protein–pigment complexes, including effects from charge transfer as well as fluorescence. The method combines the hierarchical equations of motion, which are used to describe quantum dynamics of excitons, and the Nakajima–Zwanzig quantum master equation, which is used to describe slower charge transfer processes. We study the charge transfer quenching in light harvesting complex II, a protein postulated to control non-photochemical quenching in many plant species. Using our hybrid approach, we find good agreement between our calculation and experimental measurements of the excitation lifetime. Furthermore, our calculations reveal that the exciton energy funnel plays an important role in determining quenching efficiency, a conclusion we expect to extend to other proteins that perform protective excitation quenching. This also highlights the need for simulation methods that properly account for the interplay of exciton dynamics and charge transfer processes.

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

Overcoming efficiency and cost barriers for large-area quantum dot photovoltaics through stable ink engineering

The bottom-up construction of electronics from colloidal quantum dots (CQDs) could innovate nanotechnology manufacturing through printing. However, the unstable and expensive semiconductive CQD inks make the scaling up of CQD electronics challenging. Here, in this study, we develop a strategy for engineering the solution chemistry of lead sulfide (PbS) CQD inks prepared from a low-cost direct synthesis method. By creating an iodine-rich environment in weakly coordinating solvents, we convert the iodoplumbates into functional anions, which condense into a robust surface shell. The fully charged electrostatic surface layer prevents aggregation and epitaxial fusion of CQDs, yielding stable inks. By eliminating the fusion-induced inter-band states, we print a compact CQD film with uniformity in three dimensions, flattened energy landscape and improved carrier transport. We achieved a certified efficiency of 13.40% on 0.04 cm 2 cells, with a 300-fold increase in active area, scaling up to a 12.60 cm 2 module with a certified efficiency of 10%.

14 SOLAR ENERGY↗
DOE OSTI2022

Final Technical Report for "Cloud-based Low-Scaling Quantum Chemistry Simulations for Materials"

The objective of phase I was to develop the basic infrastructure for performing efficient hybrid DFT calculations on solid materials with a turnaround time amenable to the use in large-scale data-driven approaches. To fulfill the goal we have developed a novel algorithm that significantly reduces the cost of exchange matrix formation. The algorithm does so by making use of three ingredients (a) interpolative decomposition of the electron integrals (b) robust pseudospectral method and (c) occ-RI approach. Our published work demonstrates that the algorithm is orders of magnitude faster than any other hybrid-DFT method and is on the order of only 3-4 times slower than pure DFT. All steps of the algorithm developed during phase I can be accelerated to the point that, for large enough systems, the diagonalization of the Fock matrix will become the most expensive step. To treat such large systems we have created an interface to the ASCR-funded PEXSI library. The PEXSI method is used to compute the density matrix directly from the Fock matrix in a manner that preserves the sparsity of the local representation.

Shiozaki, Toru↗
DOE OSTI2026

Integrative computational investigation of the spectroscopy, dynamics, and controlling of molecular excitons in complex environments (Final Technical Report)

The major goal of this research project was to advance theoretical and computational capability to describe the dynamics and the spectroscopy of molecular excitons as reliably and efficiently as possible. Outcomes of the project contributed to new understanding of the interaction between light and molecular materials and ensuing dynamics of excited electronic states. These can assist experimentalists in developing novel methods to utilize quantum properties of molecules and molecular excitons for new energy harvesting and quantum information processing.

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

Review of computational approaches to predict the thermodynamic stability of inorganic solids

Improvements in the efficiency and availability of quantum chemistry codes, supercomputing centers, and open materials databases have transformed the accessibility of computational materials design approaches. Thermodynamic stability predictions play a central role in the efficacy of these approaches and should be considered carefully. This review covers the fundamentals of calculating thermodynamic stability using first-principles methods. Stability is delineated into two main topics — stability with respect to decomposition into competing phases and stability with respect to phase transition into alternative structures at fixed composition. For each topic, a summary of the state-of-the-art is provided along with a tutorial overview of practical considerations. The application of machine learning to both kinds of stability predictions is also covered. Finally, the limitations of thermodynamic stability predictions are discussed within the context of predicting the synthesizability of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solvent structure and dynamics over Brønsted acid MWW zeolite nanosheets

In the liquid phase of heterogeneous catalysis, solvent plays an important role and governs the kinetics and thermodynamics of a reaction. Although it is often difficult to quantify the role of the solvent, it becomes particularly challenging when a zeolite is used as the catalyst. This difficulty arises from the complex nature of the liquid/zeolite interface and the different solvation environments around catalytically active sites. Here, we use ab initio molecular dynamics simulations to probe the local solvation structure and dynamics of methanol and water over MWW zeolite nanosheets with varying Brønsted acidity. We find that the zeolite framework and the number and location of the acid sites in the zeolite influence the structure and dynamics of the solvent. In particular, methanol is more likely to be in the vicinity of the aluminum (Al 3+ ) at the T4 site than at T1 due to easy accessibility. The methanol oxygen binds strongly to the Al at the T4 site, weakening the Al–O for the bridging acid site, which results in the formation of the silanol group, significantly reducing the acidity of the site. The behavior of methanol is in direct contrast to that of water, where protons can easily propagate from the zeolite to the solvent molecules regardless of the acid site location. Furthermore, our work provides molecular-level insights into how solvent interacts with zeolite surfaces, leading to an improved understanding of the catalytic site in the MWW zeolite nanosheet.

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

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗
DOE OSTI2021

Dimensionality reduction of the many-body problem using coupled-cluster subsystem flow equations: classical and quantum computing perspective

We discuss reduced-scaling strategies employing recently introduced sub-system embedding sub-algebras coupled-cluster formalism (SES-CC) to describe many-body systems. These strategies utilize properties of the SES-CC formulations where the equations describing certain classes of sub- systems can be integrated into a computational flows composed coupled eigenvalue problems of reduced dimensionality. Additionally, these flows can be defined at the level of the CC Ansatz defined by selected classes of cluster amplitudes, which define the wave function ”memory” of possible partitionings of the many-body system into constituent sub-systems. One of the possible ways of solving these coupled problems is through implementing procedures, where the information is passed between the sub-systems in a self-consistent manner. As a special case, we consider local flow formulations where the so-called local character of correlation effects can be closely related to properties of sub-system embedding sub-algebras employing localized molecular basis. We also generalize flow equations to the time domain and to downfolding methods utilizing double exponential unitary CC Ansatz (DUCC), where reduced dimensionality of constituent sub-problems offer a possibility of efficient utilization of limited quantum resources in modeling realistic systems.

Electron correlation, quantum chemistry, quantum c↗
DOE OSTI2026

Hydride and Seek: Comparing Crystallographic Hydride Placement Techniques with an Open-Shell Cobalt Complex

Locating hydrides is crucial in organometallic chemistry but difficult to do accurately using X-ray diffraction. Electron diffraction has been proposed as a way to overcome this problem but has not been systematically compared to neutron diffraction and to quantum crystallography (Hirshfeld atom refinement, HAR) to test this hypothesis. Here, we present a comparative analysis of methods for a terminal cobalt hydride complex by comparing a single-crystal neutron diffraction reference structure to results from single-crystal X-ray diffraction with and without Hirshfeld atom refinement (HAR, NoSpherA2), density functional theory (DFT), and electron diffraction (3D-ED/MicroED) refined under kinematical and dynamical formalisms. Conventional X-ray diffraction gives lower precision than neutron diffraction as expected. Despite expected improvements, HAR gives systematic deviation from the neutron benchmark. Interestingly, optimized DFT equilibrium geometries are closer to the neutron value than the value from HAR. On the other hand, electron diffraction with a high-quality data set coupled with dynamical refinement localizes the hydride in difference maps and gives excellent agreement with the neutron data. Dynamical refinement is crucial, as kinematical refinement does not allow assignment of a hydride peak. This cross-modal comparison defines the conditions under which 3D-ED/MicroED delivers high-precision metal–hydride distances for this open-shell cobalt hydride.

anions↗
DOE OSTI2025

Synthesis of nanodiamonds encapsulated by zeolitic imidazole framework-8 for quantum sensing applications

Nitrogen vacancy (NV)-containing nanodiamonds are widely used in quantum sensing applications due to their high sensitivity to magnetic fields, relatively low cost, and ability to be initialized, manipulated, and read out at room temperature. Quantum sensing techniques such as optically detected magnetic resonance (ODMR) and spin relaxometry have exploited the sensitivity of the NV nanodiamonds to magnetic fields to detect a range of analytes, such as pH, metal ions, and biomolecules. However, diversifying the sensing targets accessible by NV diamond quantum sensors typically requires careful engineering of the diamond surface chemistry with stimuli-responsive functional groups. Here, a simple protocol for coating NV nanodiamonds with the zeolitic imidazole framework 8 (ZIF-8), a widely used metal-organic framework, is presented. ZIF-8 is a highly porous material that has been used as a selective sensor for gasses, metal ions, and other analytes. The material is well-characterized by x-ray diffraction, transmission electron microscopy, scanning electron microscopy, x-ray photoelectron spectroscopy, and luminescence spectroscopy. Encapsulation of NV nanodiamonds with a porous scaffold such as ZIF-8 provides a promising method for improving the selectivity for the quantum sensing of various analytes. Importantly, the ZIF-8 coating does not impact the luminescence properties of the NV diamond, which is a key readout in ODMR and spin relaxometry sensing approaches. Indeed, the ODMR spectra with and without the ZIF-8 shell is nearly identical. Moreover, the ZIF-8 coating increases the longitudinal spin relaxation time of the NV nanodiamond by a factor of 4 relative to aggregated diamond, a desirable outcome for spin relaxation-based quantum sensing. Metal-organic framework composites with nanodiamonds thus are an exciting strategy for enhancing NV nanodiamond performance in applications such as quantum sensing and quantum-enhanced nuclear magnetic resonance spectroscopy.

nitrogen vacancy nanodiamond↗
DOE OSTI2023

Light-Induced Spin Trapping in Transition Metal Compounds (DOE-DE-SC0017889-Final-Report)

In solar cells, devices that convert the energy of light directly into electricity, dye sensitized solar cells (DSSC) are some of the most capable devices for energy conversion. Ruthenium is a key component for these high-power conversion efficiencies in DSSC. Ruthenium, unfortunately, is among one of the rarest elements, and, thus, more earth-abundant replacements are needed. Though the properties of ruthenium are unique, these properties may be possible to mimic with other elements. To do this, the ability to control how other more earth-abundant elements like iron behave is essential. In this work, a method was developed to enable better control over iron species, for example, and their interaction with light.

14 SOLAR ENERGY↗
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Nitrogen Out-of-Plane (NOOP) Mechanism for Imine-Based Light-Driven Molecular Motors

Light-driven molecular motors have generated considerable interest due to their potential applications in material and biological systems. Recently, Greb and Lehn reported a new class of molecular motors, chiral N-alkyl imines, which undergo unidirectional rotation induced by light and heat. The mechanism of unidirectional motion in molecular motors containing a C═N group has been assumed to consist of photoinduced torsion about the double bond. In this work, we present a computational study of the photoisomerization dynamics of a chiral N-alkyl imine motor. We find that the location and energetics of minimal energy conical intersections (MECIs) alone are insufficient to understand the mechanism of the motor. Furthermore, a key part of the mechanism consists of out-of-plane distortions of the N atom (followed by isomerization about the double bond). Dynamic effects and out-of-plane distortions are critical to understand the observed (rather low) quantum yield for photoisomerization. Finally, our results provide hints as to how the photoisomerization quantum yield might be increased, improving the efficiency of this class of molecular motors.

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

Operator learning for predicting multiscale bubble growth dynamics

We report simulating and predicting multiscale problems that couple multiple physics and dynamics across many orders of spatiotemporal scales is a great challenge that has not been investigated systematically by deep neural networks (DNNs). Herein, we develop a framework based on operator regression, the so-called deep operator network (DeepONet), with the long-term objective to simplify multiscale modeling by avoiding the fragile and time-consuming “hand-shaking” interface algorithms for stitching together heterogeneous descriptions of multiscale phenomena. To this end, as a first step, we investigate if a DeepONet can learn the dynamics of different scale regimes, one at the deterministic macroscale and the other at the stochastic microscale regime with inherent thermal fluctuations. Specifically, we test the effectiveness and accuracy of the DeepONet in predicting multirate bubble growth dynamics, which is described by a Rayleigh–Plesset (R–P) equation at the macroscale and modeled as a stochastic nucleation and cavitation process at the microscale by dissipative particle dynamics (DPD). First, we generate data using the R–P equation for multirate bubble growth dynamics caused by randomly time-varying liquid pressures drawn from Gaussian random fields (GRFs). Our results show that properly trained DeepONets can accurately predict the macroscale bubble growth dynamics and can outperform long short-term memory networks. We also demonstrate that the DeepONet can extrapolate accurately outside the input distribution using only very few new measurements. Subsequently, we train the DeepONet with DPD data corresponding to stochastic bubble growth dynamics. Although the DPD data are noisy and we only collect sparse data points on the trajectories, the trained DeepONet model is able to predict accurately the mean bubble dynamics for time-varying GRF pressures. Taken together, our findings demonstrate that DeepONets can be employed to unify the macroscale and microscale models of the multirate bubble growth problem, hence providing new insight into the role of operator regression via DNNs in tackling realistic multiscale problems and in simplifying modeling with heterogeneous descriptions.

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

Broadband 2DES detection of vibrational coherence in the S x state of canthaxanthin

The nonadiabatic mechanism that mediates nonradiative decay of the bright S 2 state to the dark S 1 state of carotenoids involves population of a bridging intermediate state, S x , in several examples. The nature of S x remains to be determined definitively, but it has been recently suggested that S x corresponds to conformationally distorted molecules evolving along out-of-plane coordinates of the isoprenoid backbone near a low barrier between planar and distorted conformations on the S 2 potential surface. Here, the electronic and vibrational dynamics accompanying the formation of S x in toluene solutions of the ketocarotenoid canthaxanthin (CAN) are characterized with broadband two-dimensional electronic spectroscopy (2DES) with 7.8 fs excitation pulses and detection of the linear polarization components of the third-order nonlinear optical signal. A stimulated-emission cross peak in the 2DES spectrum accompanies the formation of S x in <20 fs following excitation of the main absorption band. S x is prepared instantaneously, however, with excitation of hot-band transitions associated with distorted conformations of CAN’s isoprenoid backbone in the low frequency onset of the main absorption band. Vibrational coherence oscillation maps and modulated anisotropy transients show that S x undergoes displacements from the Franck–Condon S 2 state along out-of-plane coordinates as it passes to the S 1 state. The results are consistent with the conclusion that CAN’s carbonyl-substituted β-ionone rings impart an intramolecular charge-transfer character that frictionally slows the passage from S x to S 1 compared to carotenoids lacking carbonyl substitution. Despite the longer lifetime, the S 1 state of CAN is formed with retention of vibrational coherence after passing through a conical intersection seam with the S x state.

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

Excitons: Energetics and spatiotemporal dynamics

The concept of an exciton as a quasiparticle that represents collective excited states was originally adapted from solid-state physics and has been successfully applied to molecular aggregates by relying on the well-established limits of the Wannier exciton and the Frenkel exciton. However, the study of excitons in more complex chemical systems and solid materials over the past two decades has made it clear that simple concepts based on Wannier or Frenkel excitons are not sufficient to describe detailed excitonic behavior, especially in nano-structured solid materials, multichromophoric macromolecules, and complex molecular aggregates. In addition, important effects such as vibronic coupling, the influence of charge-transfer (CT) components, spin-state interconversion, and electronic correlation, which had long been studied but not fully understood, have turned out to play a central role in many systems. This has motivated new experimental approaches and theoretical studies of increasing sophistication. This article provides an overview of works addressing these issues that were published for A Special Topic of the Journal of Chemical Physics on ``Excitons: Energetics and spatio-temporal dynamics" and discusses their implications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗