Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular Hamiltonian”

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 19 records

Feedback-based quantum algorithms for ground state preparation

The ground state properties of quantum many-body systems are a subject of interest across chemistry, materials science, and physics. Thus, algorithms for finding ground states can have broad impacts. Variational quantum algorithms are one class of ground state algorithms that has received significant attention in recent years. These algorithms utilize a hybrid quantum-classical computing framework to prepare ground states on quantum computers. However, this requires solving a classical optimization problem that can become prohibitively expensive in high dimensions. Here, we develop formulations of feedback-based quantum algorithms for ground state preparation that can be used to address this challenge for two broad classes of Hamiltonians: Fermi-Hubbard Hamiltonians, and molecular Hamiltonians represented in second quantization. Feedback-based quantum algorithms are optimization-free; in place of classical optimization, quantum circuit parameters are set according to a deterministic feedback law derived from quantum Lyapunov control principles. This feedback law guarantees a monotonic improvement in solution quality with respect to the depth of the quantum circuit. A variety of numerical illustrations are provided that analyze the convergence and robustness of feedback-based quantum algorithms for these problem classes. Published by the American Physical Society 2024

Larsen, James B. (ORCID:000000020777440X)↗

Fermi’s Golden Rule Rate Expression for Transitions Due to Nonadiabatic Derivative Couplings in the Adiabatic Basis

Starting from a general molecular Hamiltonian expressed in the basis of adiabatic electronic and nuclear position states, where a compact and complete expression for the nonadiabatic derivative coupling (NDC) Hamiltonian term is obtained, we provide a general analysis of the Fermi’s golden rule (FGR) rate expression for nonadiabatic transitions between adiabatic states. We then consider a quasi-adiabatic approximation that uses crude adiabatic states and NDC couplings, both evaluated at the minimum potential energy configuration of the initial adiabatic state, for the definition of the zeroth and first-order terms of the Hamiltonian. Although the application of this approximation is rather limited, it allows deriving a general FGR rate expression without further approximation while accounting for non-Condon contribution to the FGR rate arising from momentum operators of NDC terms and its coupling with vibronic displacements. For a generic and widely used model where all nuclear degrees of freedom and environmental effects are represented as linearly coupled harmonic oscillators, we derive a closed-form FGR rate expression that requires only Fourier transform. The resulting rate expression includes quadratic contributions of NDC terms and their couplings to Franck–Condon modes, which require evaluation of two additional bath spectral densities in addition to the conventional one that appears in a typical FGR rate theory based on the Condon approximation. Model calculations for the case where nuclear vibrations consist of both a sharp high-frequency mode and an Ohmic bath spectral density illustrate new features and implications of the rate expression. We then apply our theoretical expression to the nonradiative decay from the first excited singlet state of azulene, which illustrates the utility and implications of our theoretical results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Magnetoelectric Coupling in a Mn 4 Na–Organic Complex under Pulsed Magnetic Fields up to 73 T

Research on the magnetoelectric (ME) effect (or spin-electric coupling) in molecule-based magnetic materials is a relatively nascent but promising topic. Molecule-based magnetic materials have diverse magnetic functionalities that can be coupled to electrical properties. Here, in this study, we investigate a realization of ME coupling that is fundamental but not heavily studied─the coupling of magnetic spin level crossings to changes in electric polarization. A mixed-valence Mn 4 Na complex with a total ground-state spin S = 5/2 under zero magnetic field and S = 17/2 under high magnetic field undergoes a cascade of ground-state level crossings of the S z states with increasing magnetic field. Magnetization and electrical polarization measurements under pulsed magnetic fields up to 73 T show that each spin level crossing is accompanied by a significant change in electric polarization that is an even function of the applied magnetic field. A molecular Hamiltonian describing antiferromagnetic exchange in a distorted tetrahedron of three Mn III and one Mn II ions matches the data well. We conclude that the ME coupling is caused by magnetostriction within the polar molecule as it distorts to lower its magnetic exchange energy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantum Neural Networks: Issues, Training, and Applications

Our work in the field aims at explaining the limitations and expressive power of Quantum Machine Learning models, as well as finding feasible training algorithms that could be implemented in near-term Quantum Computers. The promise of Quantum Machine Learning is that by incorporating quantum effects, such as entanglement, into machine learning models researchers can improve model performance and understand more complex datasets. This pledge is particularly pronounced in the design of Quantum neural networks (QNNs), a promising framework for creating quantum algorithms, that promise to outperform classical models by combining the speedups of quantum computation with the widespread successes of deep learning. We show that applying this approach alone to quantum deep learning is problematic given that an excess of entanglement between the hidden and visible layers can destroy the predictive power of our QNN models. We address the barren plateau problem by suggesting the use of a generative, unbounded, nonlinear loss function with simple gradients. The loss function quantifies how much the quantum states generated by the QNNs differ from the data and the goal during training is to minimize it. Finally, we showcase how to use generative training to construct a "classical-quantum" neural network to accurately interpolate between the ground states of a Molecular Hamiltonian, a central question in Quantum Chemistry.

97 MATHEMATICS AND COMPUTING↗

Diabatization with Electrostatic Embedding for Studying Photophysics in Organic Molecular Crystals

Highly emissive organic molecular crystals find applications in several areas, such as organic electronics, solar cells, and sensors. Understanding the excited-state mechanisms underlying these applications is essential for optimizing and controlling them effectively. Exciton models coupled with nonadiabatic dynamics, particularly quantum dynamics, provide crucial insights into photochemical and photophysical processes in molecular crystals. Nevertheless, there remains a lack of general tools and automated workflows to facilitate such simulations. In this paper, we present a computational strategy to investigate the photoactivated dynamics of organic molecular crystals, bridging methodologies traditionally used for molecular systems and materials science, with a particular focus on the interplay between local excitations and charge transfer (CT) processes. We have implemented an interface between the fromage and Overdia programs, enabling the construction of vibronic Hamiltonians for molecular crystals within an excited-state ONIOM(QM:QM′) framework, incorporating long-range electrostatics through a RESP-based Ewald summation. Fragment-based diabatization provides a route to quantum dynamics simulations in weak-to-intermediate coupling regimes. The method was applied to the photophysics of dibenzo[g,p]chrysene (DBC) crystals using time-dependent DFT. The fromage/ Overdia interface was employed to compute the couplings of local excitations and CT states for 18 unique DBC dimers in the crystal and to quantify the influence of electrostatic embedding, which was found to be modest (10−20%). Simulations on π-stacked dimers reproduced the small red shift observed experimentally from solution to crystal, attributed to electronic interactions among fixed monomers rather than crystal electrostatics. Quantum dynamics simulations revealed ultrafast population transfer from bright local excitations to CT states. This approach establishes a robust framework linking molecular and solid-state excited-state dynamics, with potential applications for studying excitations, defects, and impurities in molecular crystals.

Crystals↗

A parallel, distributed memory implementation of the adaptive sampling configuration interaction method

The many-body simulation of quantum systems is an active field of research that involves several different methods targeting various computing platforms. Many methods commonly employed, particularly coupled cluster methods, have been adapted to leverage the latest advances in modern high-performance computing. Selected configuration interaction (sCI) methods have seen extensive usage and development in recent years. However, the development of sCI methods targeting massively parallel resources has been explored only in a few research works. Here, we present a parallel, distributed memory implementation of the adaptive sampling configuration interaction approach (ASCI) for sCI. In particular, we will address the key concerns pertaining to the parallelization of the determinant search and selection, Hamiltonian formation, and the variational eigenvalue calculation for the ASCI method. Load balancing in the search step is achieved through the application of memory-efficient determinant constraints originally developed for the ASCI-PT2 method. The presented benchmarks demonstrate near optimal speedup for ASCI calculations of Cr 2 (24e, 30o) with 10 6 , 10 7 , and 3 × 10 8 variational determinants on up to 16 384 CPUs. Importantly, to the best of the authors’ knowledge, this is the largest variational ASCI calculation to date.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scattering solution of an interacting Hamiltonian for the electronic control of molecular spin qubits

We theoretically study how a scattered electron can entangle molecular spin qubits (MSQs). This requires solving the inelastic transport of a single electron through a scattering region described by a tight-binding interacting Hamiltonian. We accomplish this using a Green's-function solution. We can model realistic physical implementations of MSQs by parametrizing the tight-binding Hamiltonian with first-principles descriptions of magnetic anisotropy and exchange interactions. As a result, we find that, for two-MSQ systems with inversion symmetry, projective measurement of the spin degree of freedom of the scattered electron offers probabilistic control of the degree of entanglement between the MSQs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effective many-body interactions in reduced-dimensionality spaces through neural network models

Accurately describing properties of challenging problems in physical sciences often requires complex mathematical models that are unmanageable to tackle head on. Therefore, developing reduced-dimensionality representations that encapsulate complex correlation effects in many-body systems is crucial to advance the understanding of these complicated problems. However, a numerical evaluation of these predictive models can still be associated with a significant computational overhead. To address this challenge, in this paper we discuss a combined framework that integrates recent advances in the development of active-space representations of coupled cluster (CC) downfolded Hamiltonians with neural network approaches. The primary objective of this effort is to train neural networks to eliminate the computationally expensive steps required for evaluating hundreds or thousands of Hugenholtz diagrams, which correspond to multidimensional tensor contractions necessary for evaluating a many-body form of downfolded effective Hamiltonians. Using small molecular systems (the H 2 O and HF molecules) as examples, we demonstrate that training neural networks employing effective Hamiltonians for a few nuclear geometries of molecules can accurately interpolate or extrapolate their forms to other geometrical configurations characterized by different intensities of correlation effects. We also discuss differences between effective interactions that define CC downfolded Hamiltonians with those of bare Hamiltonians defined by Coulomb interactions in the active spaces. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗

Cartesian equivariant representations for learning and understanding molecular orbitals

Qualitative and quantitative orbital properties such as bonding/antibonding character, localization, and orbital energies are critical to how chemists understand reactivity, catalysis, and excited-state behavior. Despite this, representations of orbitals in deep learning models have been very underdeveloped relative to representations of molecular geometries and Hamiltonians. Here, we apply state-of-the-art equivariant deep learning architectures to the task of assigning global labels to orbitals, namely energies characterizations, given the molecular coefficients from Hartree–Fock or density functional theory. The architecture we have developed, the Cartesian Equivariant Orbital Network (CEONET), shows how molecular orbital coefficients are readily featurized as equivariant node features common to all graph-based machine-learned potentials. We find that CEONET performs well at predicting difficult quantitative labels such as the orbital energy and orbital entropy. Furthermore, we find that the CEONET representation provides an intuitive latent space for differentiating orbital character for the qualitative assignment of e.g. bonding or antibonding character. In addition to providing a useful representation for further integrating deep learning with electronic structure theory, we expect CEONET to be useful for automatizing and interpreting the results of advanced electronic structure methods such as complete active space self-consistent field theory. In particular, the ability of CEONET to infer multireference character via the orbital entropy paves the way toward the machine-learned selection of active spaces.

chemical reactions↗

Bromine Dioxide, OBrO: The Rotational Spectrum and Molecular Properties

The spectra are well descrie by Hamiltonian which included centrifugal distortion efects for fine and hyperfine terms. The molecular structure and the harmonic force field have been derived, and they as well as fine and hyperfine structure constants, are compared with data of related molecules and ESR data from OBrO isolated in cryogenic salt matrices.

Bromine Dioxide Rotational Spectrum↗

The far-infrared laser magnetic resonance spectrum of the SiH radical and determination of ground state parameters

The far-infrared laser magnetic resonance spectrum of the SiH radical in the v = O level of its X2Pi state has been recorded. The signals are rather weak. The molecules were generated in the reaction between fluorine atoms and SiH4. Rotational transitions have been detected in both 2Pi1/2 and 2Pi3/2 spin components but no fine structure transitions between the spin components were observed. Proton hyperfine splittings were resolved on some lines. The measurements have been analyzed, subjected to a least-squares fit using an effective Hamiltonian, and the appropriate molecular parameters determined. The weakness of the spectrum and the failure of attempts to power saturate favorable lines are both consistent with a small value for the electric dipole moment for SiH.

Brown, J. M.↗

Mesospheric Observations and Modeling of the Zeeman Split 233.9 GHz O-18O-16 line

Observations made from Kitt Peak, AZ, of the 233.9 GHz emission line of O-18O-16 in the upper stratosphere and lower mesosphere are reported. A good model fit to the line is obtained by incorporation of pressure and temperature broadening effects, as well as a Zeeman splitting algorithm that uses a standard geomagnetic field model and a paramagnetic Hamiltonian description of the molecular energy states. These observations are used, along with the well known O-18O-16 mixing ratio, to establish absolute calibration for observations of other chemical species from Kitt Peak. Repeated measurements show no change in this absolute calibration between observation dates. The wide magnetic splitting (+/-1.8 MHz) exhibited by this line with only six Zeeman components provides a unique test of middle atmosphere Zeeman effect model calculations, supporting the use Of 02 lines by microwave atmospheric sounders to measure pressure and temperature.

Sandor, Brad J.↗

2D Quantum Simulation of MOSFET Using the Non Equilibrium Green's Function Method

The objectives this viewgraph presentation summarizes include: (1) the development of a quantum mechanical simulator for ultra short channel MOSFET simulation, including theory, physical approximations, and computer code; (2) explore physics that is not accessible by semiclassical methods; (3) benchmarking of semiclassical and classical methods; and (4) study other two-dimensional devices and molecular structure, from discretized Hamiltonian to tight-binding Hamiltonian.

Svizhenko, Alexel↗

Molecular Static Third-Order Polarizabilities of Carbon-Cage Fullerene and Their Correlation with Three Geometric Properties: Symmetry, Aromaticity, and Size

The static third-order polarizabilities (gamma) of C60, C70, five isomers of C78 and two isomers of C84 were analyzed in terms of three properties, from a geometric point of view: symmetry, aromaticity and size. The polarizability values were based on the finite field approximation using a semiempirical Hamiltonian (AM1) and applied to molecular structures obtained from density functional theory calculations. Symmetry was characterized by the molecular group order. The selection of 6-member rings as aromatic was determined from an analysis of bond lengths. Maximum interatomic distance and surface area were the parameters considered with respect to size. Based on triple linear regression analysis, it was found that the static linear polarizability (alpha) and gamma in these molecules respond differently to geometrical properties: alpha depends almost exclusively on surface area while gamma is affected by a combination of number of aromatic rings, length and group order, in decreasing importance. In the case of alpha, valence electron contributions provide the same information as all-electron estimates. For gamma, the best correlation coefficients are obtained when all-electron estimates are used and when the dependent parameter is ln(gamma) instead of gamma.

Moore, C. E.↗

Dissipative ground state preparation in ab initio electronic structure theory

Dissipative engineering is a powerful tool for quantum state preparation, and has drawn significant attention in quantum algorithms and quantum many-body physics in recent years. In this work, we introduce a novel approach using the Lindblad dynamics to efficiently prepare the ground state for general ab initio electronic structure problems on quantum computers, without variational parameters. These problems often involve Hamiltonians that lack geometric locality or sparsity structures, which we address by proposing two generic types of jump operators for the Lindblad dynamics. Type-I jump operators break the particle number symmetry and should be simulated in the Fock space. Type-II jump operators preserves the particle number symmetry and can be simulated more efficiently in the full configuration interaction space. For both types of jump operators, we prove that in a simplified Hartree-Fock framework, the spectral gap of our Lindbladian is lower bounded by a universal constant. For physical observables such as energy and reduced density matrices, the convergence rate of our Lindblad dynamics with Type-I jump operators remains universal, while the convergence rate with Type-II jump operators only depends on coarse grained information such as the number of orbitals and the number of electrons. To validate our approach, we employ a Monte Carlo trajectory-based algorithm for simulating the Lindblad dynamics for full ab initio Hamiltonians, demonstrating its effectiveness on molecular systems amenable to exact wavefunction treatment.

Quantum chemistry↗

Benchmarking machine learning interatomic potentials via phonon anharmonicity

Abstract Machine learning approaches have recently emerged as powerful tools to probe structure-property relationships in crystals and molecules. Specifically, machine learning interatomic potentials (MLIPs) can accurately reproduce first-principles data at a cost similar to that of conventional interatomic potential approaches. While MLIPs have been extensively tested across various classes of materials and molecules, a clear characterization of the anharmonic terms encoded in the MLIPs is lacking. Here, we benchmark popular MLIPs using the anharmonic vibrational Hamiltonian of ThO 2 in the fluorite crystal structure, which was constructed from density functional theory (DFT) using our highly accurate and efficient irreducible derivative methods. The anharmonic Hamiltonian was used to generate molecular dynamics (MD) trajectories, which were used to train three classes of MLIPs: Gaussian approximation potentials, artificial neural networks (ANN), and graph neural networks (GNN). The results were assessed by directly comparing phonons and their interactions, as well as phonon linewidths, phonon lineshifts, and thermal conductivity. The models were also trained on a DFT MD dataset, demonstrating good agreement up to fifth-order for the ANN and GNN. Our analysis demonstrates that MLIPs have great potential for accurately characterizing anharmonicity in materials systems at a fraction of the cost of conventional first principles-based approaches.

interatomic potentials↗