Search NASA⌕ Search

SEARCH · Search NASA

Results for “LIOUVILLE EQUATION”

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
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2020

Koopman–von Neumann approach to quantum simulation of nonlinear classical dynamics

Quantum computers can be used to simulate nonlinear non-Hamiltonian classical dynamics on phase space by using the generalized Koopman–von Neumann formulation of classical mechanics. The Koopman–von Neumann formulation implies that the conservation of the probability distribution function on phase space, as expressed by the Liouville equation, can be recast as an equivalent Schrödinger equation on Hilbert space with a Hermitian Hamiltonian operator and a unitary propagator. This Schrödinger equation is linear in the momenta because it derives from a constrained Hamiltonian system with twice the classical phase-space dimension. A quantum computer with finite resources can be used to simulate a finite-dimensional approximation of this unitary evolution operator. Quantum simulation of classical dynamics is exponentially more efficient than a deterministic Eulerian discretization of the Liouville equation if the Koopman–von Neumann Hamiltonian is sparse. Utilizing quantum walk techniques for state preparation and amplitude estimation for the calculation of observables leads to a quadratic improvement over classical probabilistic Monte Carlo algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Diffusion regime of electron–electron collisions in weakly ionized plasmas

We consider weakly ionized plasma where an elastic scattering of electrons on neutrals changes the individual acts and the rate of electron-electron collisions significantly. In this case, the kinetics of electron thermalization is very different from that in fully ionized plasma. Here, the colliding electrons do not move freely. They, instead, diffuse because of fast scattering on neutrals. We demonstrate how a proper account of this diffusion enables one to estimate the characteristic time of electron thermalization. We also present a rigorous derivation of the kinetic equation for electrons by using Bogolyubov’s method based on Liouville equations for multi-particle distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

An efficient explicit implementation of a near-optimal quantum algorithm for simulating linear dissipative differential equations

We propose an efficient block-encoding technique for the implementation of the Linear Combination of Hamiltonian Simulations (LCHS) for simulating dissipative initial-value problems. This algorithm approximates a target nonunitary operator as a weighted sum of Hamiltonian evolutions, thereby emulating a dissipative problem by mixing various time scales. We introduce an efficient encoding of the LCHS into a quantum circuit based on a simple coordinate transformation that turns the dependence on the summation index into a trigonometric function. Classically, this method is equivalent to the use of a highly accurate Fejér-Clenshaw-Curtis quadrature formula. Quantumly, this significantly simplifies block-encoding of a dissipative problem and allows one to perform an exponential number of Hamiltonian simulations by a single Quantum Signal Processing (QSP) circuit. The resulting LCHS circuit has high success probability and the selector scales logarithmically with the number of terms in the LCHS sum and linearly with time. Careful analysis of error convergence proves that this method is more efficient than other LCHS circuits that have recently appeared in the literature. We verify the quantum circuit and its scaling by simulating it on a digital emulator of fault-tolerant quantum computers and, as a test problem, solve the advection-diffusion equation. The proposed algorithm can be used for simulating a wide class of nonunitary initial-value problems including the Liouville equation with added dissipation and linear embeddings of nonlinear systems, such as the Koopman-von Neumann and Carleman embeddings.

Novikau, I [Lawrence Livermore National Laboratory↗
DOE OSTI2020

Machine learning a molecular Hamiltonian for predicting electron dynamics

We develop a computational method to learn a molecular Hamiltonian matrix from matrix-valued time series of the electron density. As we demonstrate for three small molecules, the resulting Hamiltonians can be used for electron density evolution, producing highly accurate results even when propagating 1,000 time steps beyond the training data. As a more rigorous test, we use the learned Hamiltonians to simulate electron dynamics in the presence of an applied electric field, extrapolating to a problem that is beyond the field-free training data. We find that the resulting electron dynamics predicted by our learned Hamiltonian are in close quantitative agreement with the ground truth. Our method relies on combining a reduced-dimensional, linear statistical model of the Hamiltonian with a time-discretization of the quantum Liouville equation within time-dependent Hartree Fock theory. In conclusion, we train the model using a least-squares solver, avoiding numerous, CPU-intensive optimization steps. For both field-free and field-on problems, we quantify training and propagation errors, highlighting areas for future development.

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

On the nature of two-photon transitions for a collection of molecules in a Fabry–Perot cavity

Here, we investigate the effect of a cavity on nonlinear two-photon transitions of a molecular system and we analyze how such an effect depends on the cavity quality factor, the field enhancement, and the possibility of dephasing. We find that the molecular response to strong light fields in a cavity with a variable quality factor can be understood as arising from a balance between (i) the ability of the cavity to enhance the field of an external probe and promote multiphoton transitions more easily and (ii) the fact that the strict selection rules on multiphoton transitions in a cavity support only one resonant frequency within the excitation range. Although our simulations use a classical level description of the radiation field (i.e., we solve Maxwell–Bloch or Maxwell–Liouville equations within the Ehrenfest approximation for the field–molecule interaction), based on experience with this level of approximation in the past studies of plasmonic and polaritonic systems, we believe that our results are valid over a wide range of external probing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Strong ultrafast demagnetization due to the intraband transitions

Abstract Demagnetization in ferromagnetic transition metals driven by a femtosecond laser pulse is a fundamental problem in solid state physics, and its understanding is essential to the development of spintronic devices.Ab initiocalculation of time-dependent magnetic moment in the velocity gauge so far has not been successful in reproducing the large amount of demagnetization observed in experiments. In this work, we propose a method to incorporate intraband transitions within the velocity gauge through a convective derivative in the crystal momentum space. Our results for transition-element bulk crystals (bcc Fe, hcp Co and fcc Ni) based on the time-dependent quantum Liouville equation show a dramatic enhancement in the amount of demagnetization after the inclusion of an intraband term, in agreement with experiments. We also find that the effect of intraband transitions on each ferromagnetic material is distinctly different because of their band structure and spin property differences. Our finding has a far-reaching impact on understanding of ultrafast demagnetization.

Physics↗
DOE OSTI2020

Liouville theory, AdS 2 string, and three-point functions

This is a write-up of the lectures given in Young Researchers Integrability School 2017. The main goal is to explain the connection between the ODE/IM correspondence and the classical integrability of strings in AdS. As a warm up, we first discuss the classical three-point function of the Liouville theory. The starting point is the well-known fact that the classical solutions to the Liouville equation can be constructed by solving a Schrodinger-like differential equation. We then convert it into a set of functional equations using a method similar to the ODE/IM correspondence. The classical three-point functions can be computed directly from these functional equations, and the result matches with the classical limit of the celebrated DOZZ formula. We then discuss the semi-classical three-point function of strings in AdS 2 and show that one can apply a similar idea by making use of the classical integrability of the string sigma model on AdS 2 . Furthermore, the result is given in terms of the "massive" generalization of Gamma functions, which show up also in string theory on pp-wave backgrounds and the twistorial generalization of topological string.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

High-order harmonic generation in solid C 60

High-order harmonic generation (HHG) has unleashed the power of strong laser physics in solids. Here we investigate HHG from a large system, solid C 60 , with 240 valence electrons engaging harmonic generation at each crystal momentum. We employ the density functional theory and the time-dependent Liouville equation of the density matrix to compute HHG signals. We find that under a moderately strong laser pulse, HHG signals reach 15th order, consistent with the experimental results from C 60 plasma. The helicity dependence in solid C 60 is weak, due to the high symmetry. In contrast to the general belief, HHG is unsuitable for band structure mapping in C 60 . However, we find a window of opportunity using a long wavelength, where harmonics are generated through multiple-photon excitation. In particular, the fifth-order harmonic energies closely follow the transition energy dispersion between the valence and conduction bands. Finally, this finding is expected to motivate future experimental investigations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Probing magnetic configuration-mediated topological phases via high harmonic generation in MnBi 2 Te 4

Compared to a space group, a magnetic space group is much more complex, as both magnetic structure and magnetic moment direction can change the original symmetries of materials. The interplay between space symmetry and magnetism can generate versatile novel quantum states. However, detecting these topological phases experimentally, achieved through manipulating the magnetic configuration, has been restricted. It is mainly because the intrinsic link between the theory and the available experimental technique remains elusive. Here, we show that high harmonic generation (HHG) can identify these topological quantum states. In this work, we use rhombohedral MnBi 2 Te 4 film as an example, and analyze the symmetry-dependent harmonic order and signal by combining first-principles calculations and time-dependent Liouville equation numerical computations. Our results provide a fundamental basis for using HHG to study the topological quantum phases mediated by various magnetic configurations that can be easily realized by applying an external magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Improved Precision Scaling for Simulating Coupled Quantum-Classical Dynamics

In this paper, we present a superpolynomial improvement in the precision scaling of quantum simulations for coupled quantum-classical systems. Such systems are found in, e.g., molecular-dynamics simulations within the Born-Oppenheimer approximation. By employing a framework based on the Koopman–von Neumann formulation of classical mechanics, we express the Liouville equation of motion as unitary dynamics and utilize phase kickback from a dynamical quantum simulation to calculate the quantum forces acting on classical particles. This approach allows us to simulate the dynamics of these classical particles without the overheads associated with measuring gradients and solving the equations of motion on a classical computer, resulting in a superpolynomial advantage at the price of increased space complexity. We demonstrate that these simulations can be performed in both microcanonical and canonical ensembles, enabling the estimation of thermodynamic properties from the prepared probability density. Published by the American Physical Society 2024

Simon, Sophia (ORCID:0000000332612338)↗
DOE OSTI2021

Hyperbolic three-string vertex

We begin developing tools to compute off-shell string amplitudes with the recently proposed hyperbolic string vertices of Costello and Zwiebach. Exploiting the relation between a boundary value problem for Liouville’s equation and a monodromy problem for a Fuchsian equation, we construct the local coordinates around the punctures for the generalized hyperbolic three-string vertex and investigate their various limits. This vertex corresponds to the general pants diagram with three boundary geodesics of unequal lengths. We derive the conservation laws associated with such vertex and perform sample computations. We note the relevance of our construction to the calculations of the higher-order string vertices using the pants decomposition of hyperbolic Riemann surfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Inner Belt Electron Decay Timescales: A Comparison of Van Allen Probes and DREAM3D Losses Following the June 2015 Storm

Abstract NASA's Van Allen Probes observed significant, long‐lived fluxes of inner belt electrons up to ∼1 MeV after geomagnetic storms in March and June 2015. Reanalysis of Magnetic Electron Ion Spectrometer (MagEIS) data with improved background correction showed a clearer picture of the relativistic electron population that persisted through 2016 and into 2017 above the Fennell et al. (2015, https://doi.org/10.1002/2014gl062874 ) limit. The intensity and duration of these enhancements allow estimation of decay timescales for comparison with simulated decay rates and theoretical lifetimes. We compare decay timescales from these data and DREAM3D simulations based on them using geomagnetic activity‐dependent pitch angle diffusion coefficients derived from plasmapause‐indexed wave data (Malaspina et al., 2016, https://doi.org/10.1002/2016gl069982 , 2018, https://doi.org/10.1029/2018gl078564 ) and phase space densities derived from MagEIS observations. Simulated decay rates match observed decay rates more closely than the theoretical lifetime due to significantly nonequilibrium pitch angle distributions in simulation and data. We conclude that nonequilibrium effects, rather than a missing diffusion or loss process, account for observed short decay rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗
DOE OSTI2023

Density Matrix Implementation of the Fermi–Löwdin Orbital Self-Interaction Correction Method

The Fermi–Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi–Löwdin orbitals which are used to remove the self-interaction error in the Perdew–Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi–Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn–Sham with one additional term in the Kohn–Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn–Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. Here, we analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville–von Neumann equation to calculate excitation energies.

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

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

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

BPZ equations for higher degenerate fields and nonperturbative Dyson-Schwinger equations

In the two-dimensional Liouville conformal field theory, correlation functions involving a degenerate field satisfy partial differential equations due to the decoupling of the null descendant field. On the other hand, the instanton partition function of a four-dimensional N = 2 supersymmetric theory in the Ω -background at a special point of the parameter space also satisfies a partial differential equation resulting from the constraints of the gauge field configurations. This partial differential equation can be proved using the nonperturbative Dyson-Schwinger equations. We show for the next-to-simplest case that the partial differential equations obtained from two different perspectives can be identified, thereby confirming an assertion of the Bogomol’nyi-Prasad-Sommerfield/conformal field theory correspondence. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2021

Liouville operators over the Hardy space

The role of Liouville operators in the study of dynamical systems through the use of occupation measures has been an active area of research in control theory over the past decade. This article investigates Liouville operators over the Hardy space, which encode complex ordinary differential equations in an operator over a reproducing kernel Hilbert space.

97 MATHEMATICS AND COMPUTING↗