Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HAMILTONIANS”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Simulating the Excited-State Dynamics of Polaritons with Ab Initio Multiple Spawning

Over the past decade, there has been a growth of interest in polaritonic chemistry, where the formation of hybrid light–matter states (polaritons) can alter the course of photochemical reactions. These hybrid states are created by strong coupling between molecules and photons in resonant optical cavities and can even occur in the absence of light when the molecule is strongly coupled with the electromagnetic fluctuations of the vacuum field. Here, we present a first-principles model to simulate nonadiabatic dynamics of such polaritonic states inside optical cavities by leveraging graphical processing units (GPUs). Our first implementation of this model is specialized for a single molecule coupled to a single-photon mode confined inside the optical cavity but with any number of excited states computed using complete active space configuration interaction (CASCI) and a Jaynes–Cummings-type Hamiltonian. Using this model, we have simulated the excited-state dynamics of a single salicylideneaniline (SA) molecule strongly coupled to a cavity photon with the ab initio multiple spawning (AIMS) method. We demonstrate how the branching ratios of the photodeactivation pathways for this molecule can be manipulated by coupling to the cavity. We also show how one can stop the photoreaction from happening inside of an optical cavity. Finally, we also investigate cavity-based control of the ordering of two excited states (one optically bright and the other optically dark) inside a cavity for a set of molecules, where the dark and bright states are close in energy.

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

Spatial Signatures of Electron Correlation in Least-Squares Tensor Hypercontraction

Least Squares Tensor Hypercontraction (LS-THC) has received some attention in recent years as an approach to reduce the significant computational costs of wavefunc- tion based methods in quantum chemistry. However, previous work has demonstrated that the LS-THC factorization performs disproportionately worse in the description of wavefunction components (e.g. cluster amplitudes T 2 ) than Hamiltonian compo- nents (e.g. electron repulsion integrals (pq|rs)). This work develops novel theoretical methods to study the source of these errors in the context of the real-space T 2 kernel, and reports, for the first time, the existence of a “correlation feature” in the errors of the LS-THC representation of the “exchange-like” correlation energy EX and T 2 that is remarkably consistent across ten molecular species, three correlated wavefunctions, and four basis sets. This correlation feature portends the existence of a “pair-point kernel” missing in the usual LS-THC representation of the wavefunction, which critically depends upon pairs of grid points situated close to atoms and with inter-pair distances between one and two Bohr radii. These findings point the way for future LS-THC developments to address these shortcomings.

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

High-Resolution Photoelectron Spectroscopy of NO3 – Vibrationally Excited Along Its ν3 Mode

The nitrate (NO3) radical has long been the subject of both experimental and theoretical studies due to its complex electronic structure resulting from vibronic interactions between its X̃2A2' and B̃2E' states. In particular, the definite assignment of the fundamental of its degenerate stretching vibration (ν3) is still under debate. Here, we report high-resolution photoelectron spectra of vibrationally pre-excited NO3- using the recently developed IR-cryo-SEVI technique. The anions are excited through infrared (IR) excitation near 1350 cm-1, accessing the ν3 and 2ν3(e') vibrational levels with band centers at 1350.5 and ∼2700 cm-1, respectively. The IR-cryo-SEVI spectrum for 2ν3 pre-excitation shows clear evidence for an intense 321 transition. From the position of this feature (30031 cm-1), the electron affinity of NO3 also determined in this work (31680 cm-1), and the IR excitation energy, we obtain a fundamental frequency of 1051 cm-1 for the ν3 fundamental of the NO3 radical. This assignment and other features in the IR-cryo-SEVI spectra are supported by spectral simulations based on a vibronic Köppel-Domcke-Cederbaum Hamiltonian. The simulations also show that nearly all features in the IR-cryo-SEVI spectra arise because of pseudo-Jahn-Teller coupling between the X̃ and B̃ states of NO3. The results and analysis presented here settle a long-standing controversy regarding the ν3 frequency of NO3.

Anions↗
DOE OSTI2025

The Singlet–Triplet Gap of Cyclobutadiene: The CIPSI-Driven CC( P ; Q ) Study

An accurate determination of singlet−triplet gaps in biradicals, including cyclobutadiene in the automerization barrier region where one has to balance the substantial nondynamical many-electron correlation effects characterizing the singlet ground state with the predominantly dynamical correlations of the lowest-energy triplet, remains a challenge for many quantum chemistry methods. High-level coupled-cluster (CC) approaches, such as the CC method with a full treatment of singly, doubly, and triply excited clusters (CCSDT), are often capable of providing reliable results, but routine application of such methods is hindered by their high computational costs. We have recently proposed a practical alternative to converging the CCSDT energetics at small fractions of the computational effort, even when electron correlations become stronger and connected triply excited clusters are larger and nonperturbative, by merging the CC(P;Q) moment expansions with the selected configuration interaction methodology abbreviated as CIPSI. We demonstrate that one can accurately approximate the highly accurate CCSDT potential surfaces characterizing the lowest singlet and triplet states of cyclobutadiene along the automerization coordinate and the gap between them using tiny fractions of triply excited cluster amplitudes identified with the help of relatively inexpensive CIPSI Hamiltonian diagonalizations.

Basis sets↗
DOE OSTI2026

Improved Accuracy in Semi-Experimental Structure Determination by Resolving Problems Associated with Rotation of Principal Inertial Axes of Isotopologues: Structures of 1,3-Oxazole ( c -C 3 H 3 NO)

The rotational spectrum of the normal isotopologue of 1,3-oxazole (c-C 3 H 3 NO) was observed from 43 to 750 GHz. Over 3900 transitions for the ground vibrational state are measured, assigned, and least-squares fit to sextic centrifugally distorted-rotor Hamiltonians. The measured frequencies and resulting spectroscopic constants from this extended spectral range, combined with previous measurements of the nuclear quadrupole coupling constants, will facilitate astronomical searches for oxazole across the majority of the range of modern radiotelescopes. Spectra for a set of 30 oxazole isotopologues, which include multiple isotopic substitutions of each atom, are used to determine the first semi-experimental equilibrium ($r$$^{SE}_{e}$) structure and semi-experimental substitution structure ($r$$^{SE}_{e}$), each using CCSD(T) computed values for the vibration–rotation interaction and electron-mass corrections. The large number of isotopologues, including 21 isotopologues observed for the first time, and the redundant substitutions of each atom provide sufficient spectroscopic information to determine the $r$$^{SE}_{e}$ structure with the expected high level of accuracy and precision (0.0001 or 0.0002 Å in bond distances and 0.013 to 0.025° in bond angles). In the course of this study, we analyzed a known issue for some $r$$^{SE}_{e}$ structure determinations of near-oblate asymmetric tops in which inclusion of individual isotopologues degrades the structure determination. We demonstrate that this problem primarily arises from the difference in the values of the computed vibration–rotation interaction corrections as evaluated at the computed re geometry vs the $r$$^{SE}_{e}$ geometry of the “real” molecule. Our solution to this problem substantially improves the $r$$^{SE}_{e}$ structure of oxazole and likely can be generalized to many other molecules.

Chemical structure↗
DOE OSTI2026

Restricted Partition Function Semiclassical Transition State Theory (RPF-SCTST): Applications to Reactions Involving Hydrogen Cyanide, Hydrogen Peroxide, and Formaldehyde Oxide

We demonstrate the efficiency of the numerical calculation of thermal semiclassical transition state theory (SCTST) rates across several representative chemical reactions using the restricted partition function (RPF)─viz a function that depends only on the imaginary action associated with a reaction. Here, we treat the potential energy surfaces (PESs) through fourth order expansions around the saddle point and integrate the Hamiltonian up to second order in employing vibrational perturbation theory. We apply this formalism to uni- and bimolecular reactions with different─though relatively high─barrier heights to uncover the influence of the barrier properties, minimal energies, and separability of the rotational motion on rate constants and tunneling corrections. Although all modes are coupled within the RPF, we found in our numerical examples that the rotational component in the absence of significant rotational distortions can be separated from the vibrational contributions. Moreover, a classical treatment of the rotational contribution is adequate over the temperature range from approximately 100 to 1000 K. We also found that the choice of DFT basis set can lead to variations in rate constants of up to an order of magnitude. Lastly, different schemes for counting eligible energy levels result in rate constants that differ by no more than a factor of 2, with the remaining discrepancies attributed to the treatment of near-convergent levels in perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Exciton Delocalization in Indolenine Squaraine Aggregates Templated by DNA Holliday Junction Scaffolds

Exciton delocalization plays a prominent role in the photophysics of molecular aggregates, ultimately governing their particular function or application. DNA is a compelling scaffold in which to template molecular aggregates and promote exciton delocalization. As individual dye molecules are the basis of exciton delocalization in molecular aggregates, their judicious selection is important. Motivated by their excellent photostability and spectral properties, here we examine the ability of squaraine dyes to undergo exciton delocalization when aggregated via a DNA Holliday junction (HJ) template. A commercially available indolenine squaraine dye was chosen for the study given its strong structural resemblance to Cy5, a commercially available cyanine dye previously shown to undergo exciton delocalization in DNA HJs. Three types of DNA-dye aggregate configurations—transverse dimer, adjacent dimer, and tetramer—were investigated. Signatures of exciton delocalization were observed in all squaraine-DNA aggregates. Specifically, strong blue shift and Davydov splitting were observed in steady-state absorption spectroscopy and exciton-induced features were evident in circular dichroism spectroscopy. Strongly suppressed fluorescence emission provided additional, indirect evidence for exciton delocalization in the DNA-templated squaraine dye aggregates. To quantitatively evaluate and directly compare the excitonic Coulombic coupling responsible for exciton delocalization, the strength of excitonic hopping interactions between the dyes were obtained by simultaneous fitting the experimental steady-state absorption and CD spectra via a Holstein-like Hamiltonian in which, following the theoretical approach of Kühn, Renger, and May, the dominant vibrational mode is explicitly considered. The excitonic hopping strength within indolenine squaraines was found to be comparable to that of the analogous Cy5 DNA-templated aggregate. The squaraine aggregates adopted primarily an H-type (dyes oriented parallel to each other) spatial arrangement. Extracted geometric details of dye mutual orientation in the aggregates enabled close comparison of aggregate configurations and the elucidation of the influence of dye angular relationship on excitonic hopping interactions in squaraine aggregates. Furthermore, these results encourage the application of squaraine-based aggregates in next generation systems driven by molecular excitons.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular Dynamical and Quantum Mechanical Exploration of the Site-Specific Dynamics of Cy3 Dimers Internally Linked to dsDNA

Performing spectroscopic measurements on biomolecules labeled with fluorescent probes is a powerful approach to locating the molecular behavior and dynamics of large systems at specific sites within their local environments. The indocarbocyanine dye Cy3 has emerged as one of the most commonly used chromophores. The incorporation of Cy3 dimers into DNA enhances experimental resolution owing to the spectral characteristics influenced by the geometric orientation of excitonically coupled monomeric units. Various theoretical models and simulations have been utilized to aid in the interpretation of the experimental spectra. In this study, we employ all-atom molecular dynamics simulations to study the structural dynamics of Cy3 dimers internally linked to the dsDNA backbone. We used quantum mechanical calculations to derive insights from both the linear absorption spectra and the circular dichroism data. Furthermore, we explore potential limitations within a commonly used force field for cyanine dyes. The molecular dynamics simulations suggest the presence of four possible Cy3 dimeric populations. The spectral simulations on the four populations show one of them to agree better with the experimental signatures, suggesting it to be the dominant population. Furthermore, the relative orientation of Cy3 in this population compares very well with previous predictions from the Holstein–Frenkel Hamiltonian model.

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

Inverse Design of a Graphene-Based Quantum Transducer via Neuroevolution

Here we introduce an inverse design framework based on artificial neural networks, genetic algorithms, and tight-binding calculations, capable to optimize the very large configuration space of nanoelectronic devices. Our non-linear optimization procedure operates on trial Hamiltonians through superoperators controlling growth policies of regions of distinct doping. We demonstrate that our algorithm optimizes the doping of graphene-based three-terminal devices for valleytronics applications, monotonously converging to synthesizable devices with high merit functions in a few thousand evaluations (out of similar or equal to ≃ 2 3800 possible configurations). The best-performing device allowed for a terminal-specific separation of valley currents with similar or equal to ≃ 96% (≃ 94%) K(K') valley purity. Importantly, the devices found through our non-linear optimization procedure have both a higher merit function and higher robustness to defects than the ones obtained through geometry optimization.

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

Voltage-Dependent Barrier Height of Electron Transport through Iron Porphyrin Molecular Junctions

Electron transport through iron porphyrin (FeP) molecules self-assembled on a gold (Au) substrate was investigated using conductive atomic force microscopy (AFM) to measure current–voltage (I–V) characteristics. In the direct tunneling region (|V| ≤ 0.1 V), the Simmons model was used to characterize the electron transport. Furthermore, the energy barrier between the Fermi energy level of Au and the highest occupied molecular orbital (HOMO) level of the FeP molecule was determined to be between 0.3 and 0.6 eV; the range of the electron attenuation coefficient was 0.6–0.8 Å–1. Instead of a constant barrier height, a voltage-dependent barrier height was adapted to simulate the experimental I–V curves over the entire voltage range (|V| ≤ 2 V) using the Simmons model for the intermediate case. The voltage-dependent barrier height is supported by a previously predicted response of molecular-projected self-consistent Hamiltonian orbitals. The dependence showed that the HOMO level relative to the Fermi energy level of the Au electrode decreased as the bias voltage increased. To verify the deposition of the FeP on the Au substrate, Raman spectroscopy and AFM analysis were performed.

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

Experimental and Theoretical Study of the Electronic Structures of Lanthanide Indium Perovskites LnInO 3

Ternary lanthanide indium oxides LnInO 3 (Ln = La, Pr, Nd, Sm) were synthesized by high-temperature solid-state reaction and characterized by Xray powder diffraction. Rietveld refinement of the powder patterns showed the LnInO 3 materials to be orthorhombic perovskites belonging to the space group Pnma, based on almost-regular InO 6 octahedra and highly distorted LnO 12 polyhedra. Experimental structural data were compared with results from density functional theory (DFT) calculations employing a hybrid Hamiltonian. Valence region X-ray photoelectron and K-shell X-ray emission and absorption spectra of the LnInO 3 compounds were simulated with the aid of the DFT calculations. Photoionization of lanthanide 4f orbitals gives rise to a complex final-state multiplet structure in the valence region for the 4f n compounds PrInO 3 , NdInO 3 , and SmInO 3 , and the overall photoemission spectral profiles were shown to be a superposition of final-state 4f n-1 terms onto the cross-section weighted partial densities of states from the other orbitals. The occupied 4f states are stabilized in moving across the series Pr-Nd-Sm. Band gaps were measured using diffuse reflectance spectroscopy. These results demonstrated that the band gap of LaInO 3 is 4.32 eV, in agreement with DFT calculations. This is significantly larger than a band gap of 2.2 eV first proposed in 1967 and based on the idea that In 4d states lie above the top of the O 2p valence band. However, both DFT and X-ray spectroscopy show that In 4d is a shallow core level located well below the bottom of the valence band. Band gaps greater than 4 eV were observed for NdInO 3 and SmInO 3 , but a lower gap of 3.6 eV for PrInO 3 was shown to arise from the occupied Pr 4f states lying above the main O 2p valence band.

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

Inelastic Neutron Scattering Study of Magnetic Exchange Pathways in MnS

We report an investigation of the magnetic structure and magnetic exchange pathways in MnS via neutron scattering methods, aided by density functional theory (DFT) modeling. The material has been confirmed to undergo antiferromagnetic (AFM) ordering at 152 K, with the magnetic structure representing AFM stacking of ferromagnetic (FM) (111) planes of Mn magnetic moments. Correspondingly, the magnetic structure is described by a propagation vector k = (1/2, 1/2, 1/2), with the volume of the magnetic unit cell being 8 times larger than the volume of the nuclear unit cell. Analysis of inelastic neutron scattering (INS) data collected on a powder sample of MnS revealed that the next-nearest-neighbor magnetic exchange constant (J 2 ) exceeds the nearest-neighbor exchange constant (J 1 ) by more than 3 times, while in the case of MnO, which exhibits the same nuclear and magnetic structures as MnS, the J 2 /J 1 ratio was reported to be below 1.5. Although for MnO the signs of both J 1 and J 2 indicated AFM exchange interactions, machine-learning INS data analysis in combination with DFT calculations suggests that the INS data collected on MnS are best described with J 1 < 0 and J 2 > 0, corresponding to FM and AFM exchange couplings, respectively. To achieve a satisfactory fit to the experimentally observed data, the Hamiltonian used to model the INS spectra also included the next-next-nearest-neighbor magnetic exchange constant (J 3 ). The best-fit model has been obtained with the values of the exchange constants J 1 = -0.27 meV, J 2 = 1.05 meV, and J 3 = -0.19 meV.

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

Unraveling Electronic Trends in O* and OH* Surface Adsorption in the MO 2 Transition-Metal Oxide Series

Understanding the bond strength of O* and OH* intermediates to metal-oxide surfaces is key to predicting the catalytic activity in oxygen-based electrochemistry. Here, we uncover highly non-linear trends in O* and OH* adsorption energies across the 3d, 4d, and 5d series of MO 2 transition-metal (TM) oxide surfaces computed within Hubbard- U corrected density functional theory (DFT + U ). Investigating the electronic structure with crystal orbital Hamiltonian populations (COHP) of the relevant metal–oxygen bonds reveals that the spin-dependent coupling strength between metal-d and oxygen-2p atomic orbitals together with the extent of filling of bonding and anti-bonding orbitals are the primary contributors to the adsorption energy. Importantly, we show that the integrated COHP obtained purely from bulk calculations is a highly accurate descriptor for surface adsorption energetics that captures trends across the group 5–12 TM oxide series within 0.19–0.36 eV. Our results suggest a pathway to prediction of adsorption energies for an arbitrary metal–ligand catalyst system.

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

Optical Properties of Neutral F Centers in Bulk MgO with Density Matrix Embedding

The optical spectra of neutral oxygen vacancies (F 0 centers) in the bulk MgO lattice are investigated using density matrix embedding theory. The impurity Hamiltonian is solved with the complete active space self-consistent field and second-order n-electron valence state perturbation theory (NEVPT2-DMET) multireference methods. Here, to estimate defect-localized vertical excitation energies at the nonembedding and thermodynamic limits, a double extrapolation scheme is employed. The extrapolated NEVPT2-DMET vertical excitation energy value of 5.24 eV agrees well with the experimental absorption maxima at 5.03 eV, whereas the excitation energy value of 2.89 eV at the relaxed triplet defect-localized state geometry overestimates the experimental emission at 2.4 eV by only nearly 0.5 eV, indicating the involvement of the triplet–singlet decay pathway.

embedding↗
DOE OSTI2025

Visualizing Stereodynamics in Cold Collisions through Shape Resonance Wavefunctions

Shape resonances in cold collisions are often strongly affected by stereodynamics. These resonances can sometimes be assigned by single-channel quantization along the scattering coordinate. However, sensitive steric control of collision implies strong anisotropy of the interaction potential energy surface, which usually leads to coupling among multiple scattering channels in the strongly interacting region. Hence, the resonances might be the result of quantization in a multidimensional space. Using the cold rotationally inelastic collision between para-H 2 (v 1 = 1, j 1 = 2) and HF (v 2 = 0, j 2 = 0) as an example, we analyze four low-lying shape resonances via diagonalization of the full-dimensional Hamiltonian with a stabilization method. While some resonances can indeed be assigned with a single partial wave, others apparently involve more than one scattering channel. Furthermore, a new model based on these resonance wavefunctions is developed to better understand the stereodynamics through shape resonances in cold collisions.

Collisions↗
DOE OSTI2021

Signature of Many-Body Localization of Phonons in Strongly Disordered Superlattices

Many-body localization (MBL) has attracted significant research attention due to its immunity to thermalization, role in logarithmic entanglement entropy growth, and opportunities to reach exotic quantum orders. However, experimental realization of MBL in solid-state systems has remained challenging. Here we report evidence of a possible phonon MBL phase in disordered GaAs/AlAs superlattices. Through grazing-incidence inelastic X-ray scattering, we observe a strong deviation of the phonon population from equilibrium in disordered samples doped with ErAs nanodots at low temperature, signaling a departure from thermalization. This behavior occurs within finite phonon energy and wavevector windows, suggesting a localization-thermalization crossover. We support our observation by proposing a theoretical model for the effective phonon Hamiltonian in disordered superlattices, and showing that it can be mapped exactly to a disordered 1D Bose-Hubbard model with a known MBL phase. The large on-site repulsion qualitatively explains the experimental findings. Furthermore, our work provides mode-resolved experimental evidence of phonon localization, potentially extending the scope of MBL to disordered solid-state systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Finite-Temperature Structures of Supported Subnanometer Catalysts Inferred via Statistical Learning and Genetic Algorithm-Based Optimization

Single-atom catalysts (SACs) minimize noble metal utilization and can alter the activity and selectivity of supported metal nanoparticles. However, the morphology of active centers, including single atoms and subnanometer clusters of a few atoms, remains elusive due to experimental challenges. The computational cost to describe numerous cluster shapes and sizes makes direct first-principles calculations impractical. We present a computational framework to enable structure determination for single-atom and subnanometer cluster catalysts. As a case study, we obtained the low energy structures of Pd n (n = 1-21) clusters supported on CeO 2 (111), which are critical components of automobile three-way catalysts. Trained on density functional theory data, a three-dimensional cluster expansion is established using statistical learning to describe the Hamiltonian and predict energies of supported Pdn clusters of any structure. Low energy stable and metastable structures are identified using a Metropolis Monte Carlo-based genetic algorithm in the canonical ensemble at 300 K. We observe that supported single atoms sinter to form bilayer clusters and large cluster isomers share similarities in both shape and energy, and elucidate the significance of the support and microstructure on cluster stability. We discovered a simple surrogate structure-energy model, where the energy per atom scales with the square root of the average first coordination number, which can be used to estimate energies and compare the stability of clusters. Our framework, applicable to any metal/support system, fills an important methodological gap to predict the stability of supported metal catalysts in the subnanometer regime.

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

Artificial Graphene Nanoribbons: A Test Bed for Topology and Low-Dimensional Dirac Physics

Here, we synthesize artificial graphene nanoribbons by positioning carbon monoxide molecules on a copper surface to confine its surface state electrons into artificial atoms positioned to emulate the low-energy electronic structure of graphene derivatives. We demonstrate that the dimensionality of artificial graphene can be reduced to one dimension with proper "edge" passivation, with the emergence of an effectively gapped one-dimensional nanoribbon structure. These one-dimensional structures show evidence of topological effects analogous to graphene nanoribbons. Guided by first-principles calculations, we spatially explore robust, zero-dimensional topological states by altering the topological invariants of quasi-one-dimensional artificial graphene nanostructures. The robustness and flexibility of our platform allow us to toggle the topological invariants between trivial and nontrivial on the same nanostructure. Ultimately, we spatially manipulate the states to understand fundamental coupling between adjacent topological states that are finely engineered and simulate complex Hamiltonians.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗