Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

NaCo 2 (SeO 3 ) 2 (OH): competing magnetic ground states of a new sawtooth structure with 3d 7 Co 2+ ions

While certain magnetic sublattices have long been known theoretically to give rise to emergent physics via competing magnetic interactions and quantum effects, finding such configurations in real materials is often deeply challenging. Here we report the synthesis and characterization of a new such material, NaCo 2 (SeO 3 ) 2 (OH) which crystallizes with a highly frustrated sublattice of sawtooth Co 2+ chains. Single crystals of NaCo 2 (SeO 3 ) 2 (OH) were synthesized using a low-temperature hydrothermal method. X-ray single crystal structure analysis reveals that the material crystallizes in orthorhombic space group of Pnma (no. 62). Its crystal structure exhibits one-dimensional chains of corner-sharing isosceles triangles that are made of two crystallographically distinct 3d 7 Co 2+ sites (Co(1) and Co(2)). The chains run along the b-axis and are interconnected via [SeO 3 ] groups to form a three-dimensional structure mediating super-exchange interactions. The temperature dependent magnetization data show a ferromagnetic-like (FM) transition at 11 K (T 1 ) followed by an antiferromagnetic (AFM) transition at about 6 K (T 2 ). Neutron-powder diffraction measurements reveal that at T 1 = 11 K only the Co(2) site orders magnetically, forming ferromagnetic zigzag chains along the b-axis. Below T 2 = 6 K, both Co(1) and Co(2) sites order in a nearly orthogonal configuration, with Co(1) moments lying inside the plane of the sawtooth chain while Co(2) moments cant out of the plane. The canting of the magnetic moments leads to a net ferromagnetic component along b-axis, parallel to the chain direction. The ordered moments are fully compensated in the ac-plane. Inelastic neutron scattering measurements reveal crystal field excitations that are consistent with the presence of a spin–orbital entangled pseudo-spin state J eff = 1/2 for the Co 2+ ions. Low-energy spin-wave excitations are observed below the second magnetic transition. The analysis of powder excitation spectrum suggests complex exchange interactions that go beyond a Heisenberg Hamiltonian model with nearest neighbor couplings. Our results demonstrate the richness of the magnetic properties of sawtooth-type structures and encourage the study of similar structures with different oxyanion groups.

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

Spin dynamics of [1,2- 13 C 2 ]pyruvate hyperpolarization by parahydrogen in reversible exchange at micro Tesla fields

Hyperpolarization of 13 C-pyruvate via Signal Amplification By Reversible Exchange (SABRE) is an important recent discovery because of both the relative simplicity of hyperpolarization and the central biological relevance of pyruvate as a biomolecular probe for in vitro or in vivo studies. Here, we analyze the [1,2- 13 C 2 ]pyruvate-SABRE spin system and its field dependence theoretically and experimentally. We provide first-principles analysis of the governing 4-spin dihydride- 13 C 2 Hamiltonian and numerical spin dynamics simulations of the 7-spin dihydride- 13 C 2 –CH 3 system. Furthermore, the analytical and the numerical results are compared to matching systematic experiments. With these methods we unravel the observed spin state mixing of singlet states and triplet states at microTesla fields and we also analyze the dynamics during transfer from micro-Tesla field to high field for detection to understand the resulting spectra from the [1,2- 13 C 2 ]pyruvate-SABRE system.

13C pyruvate↗
DOE OSTI2023

Polaron absorption in aligned conjugated polymer films: breakdown of adiabatic treatments and going beyond the conventional mid-gap state model

This study provides the first experimental polarized intermolecular and intramolecular optical absorption components of field-induced polarons in regioregular poly(3-hexylthiophene-2,5-diyl), rr-P3HT, a polymer semiconductor. Highly aligned rr-P3HT thin films were prepared by a high temperature shear-alignment process that orients polymer backbones along the shearing direction. rr-P3HT in-plane molecular orientation was measured by electron diffraction, and out-of-plane orientation was measured through series of synchrotron X-ray scattering techniques. Then, with molecular orientation quantified, polarized charge modulation spectroscopy was used to probe mid-IR polaron absorption in the ℏω = 0.075 — 0.75 eV range and unambiguously assign intermolecular and intramolecular optical absorption components of hole polarons in rr-P3HT. This data represents the first experimental quantification of these polarized components and allowed long-standing theoretical predictions to be compared to experimental results. The experimental data is discrepant with predictions of polaron absorption based on an adiabatic framework that works under the Born–Oppenheimer approximation, but the data is entirely consistent with a more recent nonadiabatic treatment of absorption based on a modified Holstein Hamiltonian. This nonadiabatic treatment was used to show that both intermolecular and intramolecular polaron coherence break down at length scales significantly smaller than estimated structural coherence in either direction. Importantly, this strongly suggests that polaron delocalization is fundamentally limited by energetic disorder in rr-P3HT.

36 MATERIALS SCIENCE↗
DOE OSTI2025

SANE: strategic autonomous non-smooth exploration for multiple optima discovery in multi-modal and non-differentiable black-box functions

Both computational and experimental material discovery bring forth the challenge of exploring multidimensional and multimodal parameter spaces, such as phase diagrams of Hamiltonians with multiple interactions, composition spaces of combinatorial libraries, material structure image spaces, and molecular embedding spaces. Often these systems are black-boxes and time-consuming to evaluate, which resulted in strong interest towards active learning methods such as Bayesian optimization (BO). However, these systems are often noisy which make the black box function severely multi-modal and non-differentiable, where a vanilla BO can get overly focused near a single or faux optimum, deviating from the broader goal of scientific discovery. To address these limitations, here we developed Strategic Autonomous Non-Smooth Exploration (SANE) to facilitate an intelligent Bayesian optimized navigation with a proposed cost-driven probabilistic acquisition function to find multiple global and local optimal regions, avoiding the tendency to becoming trapped in a single optimum. To distinguish between a true and false optimal region due to noisy experimental measurements, a human (domain) knowledge driven dynamic surrogate gate is integrated with SANE. We implemented the gate-SANE into pre-acquired piezoresponse spectroscopy data of a ferroelectric combinatorial library with high noise levels in specific regions, and piezoresponse force microscopy (PFM) hyperspectral data. SANE demonstrated better performance than classical BO to facilitate the exploration of multiple optimal regions and thereby prioritized learning with higher coverage of scientific values in autonomous experiments. Our work showcases the potential application of this method to real-world experiments, where such combined strategic and human intervening approaches can be critical to unlocking new discoveries in autonomous research.

Biswas, Arpan [University of Tennessee, Knoxville,↗
DOE OSTI2024

Rapidly convergent quantum Monte Carlo using a Chebyshev projector

The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.

Zhao, Zijun↗
DOE OSTI2026

Towards a theoretical understanding of excitonic properties of phthalocyanine thin films. I. Low-temperature exciton absorption spectra

Phthalocyanine (Pc)-based molecular thin films have emerged in recent years as a promising class of organic semiconductor materials for optoelectronic applications owing to their long exciton coherence length and fast exciton diffusion. However, the dependence of their exciton properties on the dimensionality and thermodynamic conditions, presence of metal ions, and effects of chemical modifications to the PC systems is not yet fully understood. As a first step towards a more comprehensive theoretical understanding of the excitonic properties of Pc thin films, we model their low-temperature exciton absorption spectra by employing the Frenkel Hamiltonian. The latter is derived from quantum-chemical estimates of site energies and exciton–exciton couplings. The predicted exciton absorption spectra of octabutoxy phthalocyanine (H 2 OBPc) is found to be strongly dependent on the dimensionality of the model as well as the distance cutoff for including the monomer–monomer exciton coupling. We also caution that the widely used dipole–dipole approximation could substantially overestimate the excitonic coupling between different monomers compared to a more accurate evaluation using the respective transition densities.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Advanced spectroscopic studies of (PPh 4 ) 2 [Co(N 3 ) 4 ], a field-induced single-ion magnet

The high-spin Co II complex (PPh 4 ) 2 [Co(N 3 ) 4 ] (Co-N 3 )has been investigated using advanced spectroscopic techniques [far-IR magneto-spectroscopy (FIRMS), high-frequency and high-field EPR (HFEPR), and inelastic neutron scattering (INS)] to study its zero-field-splitting (ZFS), giving spin-Hamiltonian (SH) parameters. The analysis of multi-frequency HFEPR reveals the easy-axis anisotropy with a D value of −10.39(5) cm −1 and a rhombic ratio (E/D) of 0.21(1). The magnetic properties have also been probed by direct-current (DC) magnetometry, suggesting minor differences in anisotropy from the previously reported polymorph (Co-N 3 ′). Ligand-field theory (LFT) analysis indicates that the structures of Co-N 3 and Co-N 3 ′ are closer to D 2d symmetry than other symmetries considered. Alternate-current (AC) susceptibility reveals slow magnetic relaxation under an applied field, indicating that Co-N 3 is a field-induced single-ion magnet (SIM). Here, while both Co-N 3 and Co-N 3 ′ were studied by DC magnetometry, one unusual aspect of the current work on Co-N 3 is that advanced spectroscopies HFEPR, FIRMS, and INS were used to directly observe transitions between ZFS split states, giving accurate SH parameters.

Hand, Adam T. [Univ. of Tennessee, Knoxville, TN (↗
DOE OSTI2026

A new layered kagome strip structure Na 2 Co 3 (AsO 4 ) 2 (OH) 2 : static and dynamic magnetic properties

One-dimensional kagome strip chains share much of the same frustrated structural motif as two-dimensional kagome antiferromagnets, making them valuable for deepening our understanding of kagome lattice magnetism. In this paper, we report the hydrothermal synthesis and detailed structural and property characterization of Na 2 Co 3 (AsO 4 ) 2 (OH) 2 , a striped kagome system. The crystal structure was characterized using single crystal X-ray diffraction which reveals that Na 2 Co 3 (AsO 4 ) 2 (OH) 2 crystallizes in monoclinic crystal system C2/m. The structure features a one-dimensional kagome strip lattice built from Co 2+ ions and undergoes an antiferromagnetic transition at T N = 14 K. The magnetic ground state at zero field was characterized using neutron powder diffraction. Below the magnetic transition, Na 2 Co 3 (AsO 4 ) 2 (OH) 2 orders into an antiferromagnetic structure with a k-vector (0.5, 0.5, 0.5). In the proposed model, the Co1 moment is predominantly confined to the ac-plane while the Co2 moment is primarily aligned along the b-axis. Two flat bands were observed in the inelastic neutron spectra below the magnetic transition, 5 and 10 meV. Inelastic neutron spectra were modeled with a Heisenberg Hamiltonian including three nearest-neighbor exchange interactions (J 1 , J 2 , J 3 ) and strong single-ion anisotropy to stabilize the observed magnetic structure. Our study highlights the complexity of Co 2+ -based kagome strip magnetic lattice compound Na 2 Co 3 (AsO 4 ) 2 (OH) 2 which provides an excellent platform to broaden our understanding of the frustrated kagome magnetic lattice space.

Liurukara, Duminda S. [Oak Ridge National Laborato↗
DOE OSTI2026

Gaunt and Breit two-electron contributions to mean-field transformations and fine structure splitting

Materials utilized by novel energy systems are often studied using weakly correlated mean-field theories. However, if these systems incorporate heavy elements, relativistic effects must be included. Therefore, a Kramers unrestricted coupled cluster with singles and doubles excitation formalism within a molecular mean-field exact two-component framework (X2C mmf ) using a four-component Dirac–Hartree–Fock (DHF) reference state is presented. The exact X2C mmf transformed normal-order Hamiltonian incorporates all one-electron and two-electron (2e) contributions from the Coulomb, Gaunt, and Breit operators and is used with the equation of motion method to calculate the excitation energies of the alkali group of elements. Using this framework, the effects of 2e Gaunt and Breit integrals are studied. Results demonstrate growing contributions from these integrals to the generated X2C mmf mean-fields and electronic fine structure calculations with increasing atomic number. Overall, this paper outlines the method, its effect within the X2C mmf approach, and lays the foundation for future theoretical development of relativistic calculations within this framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Solving differential‐algebraic equations in power system dynamic analysis with quantum computing

Abstract Power system dynamics are generally modeled by high dimensional non‐linear differential‐algebraic equations (DAEs) given a large number of components forming the network. These DAEs' complexity can grow exponentially due to the increasing penetration of distributed energy resources, whereas their computation time becomes sensitive due to the increasing interconnection of the power grid with other energy systems. This paper demonstrates the use of quantum computing algorithms to solve DAEs for power system dynamic analysis. We leverage a symbolic programming framework to equivalently convert the power system's DAEs into ordinary differential equations (ODEs) using index reduction methods and then encode their data into qubits using amplitude encoding. The system non‐linearity is captured by Hamiltonian simulation with truncated Taylor expansion so that state variables can be updated by a quantum linear equation solver. Our results show that quantum computing can solve the power system's DAEs accurately with a computational complexity polynomial in the logarithm of the system dimension. We also illustrate the use of recent advanced tools in scientific machine learning for implementing complex computing concepts, that is, Taylor expansion, DAEs/ODEs transformation, and quantum computing solver with abstract representation for power engineering applications.

computational complexity↗
DOE OSTI2020

The role of rotation-vibration coupling in symmetric and asymmetric isotopomers of ozone

A theoretical framework and a computer code (SpectrumSDT) are developed for accurate calculations of coupled rotational–vibrational states in triatomic molecules using hyper-spherical coordinates and taking into account the Coriolis coupling effect. Concise final formulas are derived for the construction of the Hamiltonian matrix using an efficient combination of the variational basis representation and discrete variable representation methods with locally optimized basis sets and grids. First, the new code is tested by comparing its results with those of the APH3D program of Kendrick et al. [Kendrick, Pack, Walker, and Hayes, J. Chem. Phys. 110, 6673 (1999)]. Then, accurate calculations of the rovibrational spectra are carried out for doubly substituted symmetric ( 18 O 16 O 18 O) and asymmetric ( 18 O 18 O 16 O) ozone isotopomers for the total angular momentum up to J = 5. Together with similar data recently reported for the singly substituted symmetric ( 16 O 18 O 16 O) and asymmetric ( 16 O 16 O 18 O) ozone isotopomers, these calculations quantify the role of the Coriolis coupling effect in the large mass-independent isotopic enrichment of ozone, observed in both laboratory experiments and the atmosphere of the Earth. It is found that the Coriolis effect in ozone is relatively small, as evidenced by deviations of its rotational constants from the symmetric-top-rotor behavior, magnitudes of parity splittings (Λ-doubling), and ratios of rovibrational partition functions for asymmetric vs symmetric ozone molecules. It is concluded that all of these characteristics are influenced by the isotopic masses as much as they are influenced by the overall symmetry of the molecule. It is therefore unlikely that the Coriolis coupling effect could be responsible for symmetry-driven mass-independent fractionation of oxygen isotopes in ozone.

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

27 AlNMR chemical shift of Al(OH)$^{-}_{4}$ calculated from first principles: Assessment of error cancellation in chemically distinct reference and target systems

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O)$^{3+}_{6}$, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)$^{-}_{4}$ requires an understanding of how these errors cancel, in order to define the limits of accurately predicting 27 Al chemical shielding in Al(OH)$^{-}_{4}$. In this work, we estimate the absolute shielding of the 27 Al nucleus in Al(OH)$^{-}_{4}$ at the coupled cluster level (515.1± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab-initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O)$^{3+}_{6}$ and Al(OH)$^{-}_{4}$ ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the 27 Al chemical shift of Al(OH)$^{-}_{4}$ from first principles.

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

Molecular structure determination: Equilibrium structure of pyrimidine (m-C 4 H 4 N 2 ) from rotational spectroscopy (reSE) and high-level ab initio calculation (re) agree within the uncertainty of experimental measurement

The pure rotational spectrum of pyrimidine (m-C 4 H 4 N 2 ), the meta-substituted dinitrogen analog of benzene, has been studied in the millimeter-wave region from 235 GHz to 360 GHz. The rotational spectrum of the ground vibrational state has been assigned and fit to yield accurate rotational and distortion constants. Over 1700 distinct transitions were identified for the normal isotopologue in its ground vibrational state and least-squares fit to a partial sextic S-reduced Hamiltonian. Transitions for all four singly substituted 13 C and 15 N isotopologues were observed at natural abundance and were likewise fit. Deuterium-enriched samples of pyrimidine were synthesized, giving access to all eleven possible deuterium-substituted isotopologues, ten of which were previously unreported. Experimental values of rotational constants and computed values of vibration–rotation interaction constants and electron-mass corrections were used to determine semi-experimental equilibrium structures (r e SE ) of pyrimidine. The r e SE structure obtained using coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] corrections shows exceptional agreement with the re structure computed at the CCSD(T)/cc-pCV5Z level (≤0.0002 Å in bond distance and ≤0.03° in bond angle). Of the various computational methods examined, CCSD(T)/cc-pCV5Z is the only method for which the computed value of each geometric parameter lies within the statistical experimental uncertainty (2σ) of the corresponding semi-experimental coordinate. The exceptionally high accuracy and precision of the structure determination is a consequence of the large number of isotopologues measured, the precision and extent of the experimental frequency measurements, and the sophisticated theoretical treatment of the effects of vibration–rotation coupling and electron mass. Taken together, these demanding experimental and computational studies establish the capabilities of modern structural analysis for a prototypical monocyclic aromatic compound.

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

High-resolution infrared spectroscopy of jet cooled CH 2 Br radicals: The symmetric CH stretch manifold and absence of nuclear spin cooling

Direct laser absorption of a slit supersonic discharge expansion provides the first high-resolution spectroscopic results on the symmetric CH stretch excitation (ν 1 ) of the bromomethyl (CH 2 Br) radical in the ground electronic state. Narrowband (<1 MHz) mid-infrared radiation is produced by difference-frequency generation of two visible laser beams, with the open shell halohydrocarbon radical generated by electron dissociative attachment of CH 2 Br 2 in a discharge and rapidly cooled to T rot = 18 ± 1 K in the subsequent slit-jet supersonic expansion. A rovibrational structure in the radical spectrum is fully resolved, as well as additional splittings due to spin-rotation effects and 79 Br/ 81 Br isotopologues in natural abundance. Spectroscopic constants and band origins are determined by fitting the transition frequencies to a non-rigid Watson Hamiltonian, yielding results consistent with a vibrationally averaged planar radical and an unpaired electron in the out-of-plane p π orbital. Additionally, extensive satellite band structure from a vibrational hot band is observed and analyzed. The hot band data is compared to CFOUR/VPT2 (CCSD(T)cc-pVQZ) ab initio anharmonic predictions of the vibration rotation alpha matrix, which permits unambiguous assignment to CH 2 symmetric-stretch excitation built on the singly excited CH 2 out-of-plane bending mode (ν 1 + ν 4 ← ν 4 ). Longitudinal cooling of the Doppler width in the slit-jet expansion geometry also reveals partially resolved hyperfine structure on transitions out of the lowest angular momentum states in excellent agreement with predictions based on microwave studies. High level ab initio MOLPRO calculations [CCSD(T)-f12b/VnZ-f12 (n = 3, 4, CBS)] are also performed with explicitly correlated f12 electron methods for the out-of-plane CH 2 bending mode over the halogen series CH 2 X (X = F, Cl, Br, I), which clearly reveals a non-planar geometry for X = F (with a ΔE ≈ 0.3 kcal/mol barrier) and yet planar equilibrium geometries for X = Cl, Br, and I. Finally, a detailed Boltzmann analysis of the transition intensities provides support for negligible collisional equilibration of the entangled H atom nuclear spin states on the few hundred microsecond time scale and high collision densities of a slit supersonic expansion.

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

Reduced scaling extended multi-state CASPT2 (XMS-CASPT2) using supporting subspaces and tensor hyper-contraction

We present a reduced scaling formulation of the extended multi-state CASPT2 (XMS-CASPT2) method, which is based on our recently developed state-specific CASPT2 (SS-CASPT2) formulation using supporting subspaces and tensor hyper-contraction. By using these two techniques, the off-diagonal elements of the effective Hamiltonian can be computed with only O(N 3 ) operations and O(N 2 ) memory, where N is the number of basis functions. Furthermore, this limits the overall computational scaling to O(N 4 ) operations and O(N 2 ) memory. Thus, excited states can now be obtained at the same reduced (relative to previous algorithms) scaling we achieved for SS-CASPT2.

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

Multichannel dynamics in the OH+ n -butane reaction revealed by crossed-beam slice imaging and quasiclassical trajectory calculations

Multidimensional reactions present various channels that can exhibit very different dynamics and give products of varying subsequent reactivity. Here, we present a combination of experiment and theory to reveal the dynamics of hydrogen abstraction by OH radical at primary and secondary sites in n-butane at a collision energy of 8 kcal/mol. Crossed molecular beam slice imaging experiments unequivocally probe the secondary abstraction channel showing backward angular distributions with mild energy release to product translation, which are accurately captured by trajectory calculations using a specific-reaction-parameter Hamiltonian. Additionally, experiments containing both reaction channels indicate a less marked backward character in the angular distribution, whose origin is shown by trajectory calculations to appear as an evolution toward more sideways scattering from the secondary to primary channel. While the two channels have markedly different angular distributions, their energy release is largely comparable, showing ample energy release into the water product. The synergistic combination of crossed-beam imaging and trajectories opens the door to detailed reaction-dynamics studies of chemical reactions with ever-increasing complexity.

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

On the proper derivation of the Floquet-based quantum classical Liouville equation and surface hopping describing a molecule or material subject to an external field

We investigate different approaches to derive the proper Floquet-based quantum–classical Liouville equation (F-QCLE) for laser-driven electron-nuclear dynamics. The first approach projects the operator form of the standard QCLE onto the diabatic Floquet basis and then transforms to the adiabatic representation. The second approach directly projects the QCLE onto the Floquet adiabatic basis. Both approaches yield a form that is similar to the usual QCLE with two modifications: (1) The electronic degrees of freedom are expanded to infinite dimension and (2) the nuclear motion follows Floquet quasi-energy surfaces. However, the second approach includes an additional cross derivative force due to the dual dependence on time and nuclear motion of the Floquet adiabatic states. Our analysis and numerical tests indicate that this cross derivative force is a fictitious artifact, suggesting that one cannot safely exchange the order of Floquet state projection with adiabatic transformation. Our results are in accord with similar findings by Izmaylov et al., [J. Chem. Phys. 140, 084104 (2014)] who found that transforming to the adiabatic representation must always be the last operation applied, although now we have extended this result to a time-dependent Hamiltonian. Finally, this paper and the proper derivation of the F-QCLE should lay the basis for further improvements of Floquet surface hopping.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗