Koopmon trajectories in nonadiabatic quantum-classical dynamics

Werner Bauer,Paul Bergold,François Gay-Balmaz,Cesare Tronci

2024-06-18

Abstract:In order to alleviate the computational costs of fully quantum nonadiabatic dynamics, we present a mixed quantum-classical (MQC) particle method based on the theory of Koopman wavefunctions. Although conventional MQC models often suffer from consistency issues such as the violation of Heisenberg's principle, we overcame these difficulties by blending Koopman's classical mechanics on Hilbert spaces with methods in symplectic geometry. The resulting continuum model enjoys both a variational and a Hamiltonian structure, while its nonlinear character calls for suitable closures. Benefiting from the underlying action principle, here we apply a regularization technique previously developed within our team. This step allows for a singular solution ansatz which introduces the trajectories of computational particles - the koopmons - sampling the Lagrangian classical paths in phase space. In the case of Tully's nonadiabatic problems, the method reproduces the results of fully quantum simulations with levels of accuracy that are not achieved by standard MQC Ehrenfest simulations. In addition, the koopmon method is computationally advantageous over similar fully quantum approaches, which are also considered in our study. As a further step, we probe the limits of the method by considering the Rabi problem in both the ultrastrong and the deep strong coupling regimes, where MQC treatments appear hardly applicable. In this case, the method succeeds in reproducing parts of the fully quantum results.

Numerical Analysis,Mathematical Physics,Chemical Physics,Quantum Physics

What problem does this paper attempt to address?

This paper attempts to address the computational cost issue in non - adiabatic quantum - classical dynamics simulations. Specifically, the authors propose a hybrid quantum - classical (MQC) particle method based on Koopman wave function theory to overcome the consistency problems in traditional MQC models, such as the violation of the Heisenberg uncertainty principle.

### Main problems

1. **Computational cost**: The computational cost of fully quantum non - adiabatic dynamics is extremely high, which limits its feasibility in practical applications.
2. **Consistency problems**: Traditional MQC models (such as the Ehrenfest model) often have consistency problems. For example, they cannot maintain the positivity of the classical density and the quantum density matrix, resulting in the failure of the uncertainty principle.

### Solutions

To solve the above problems, the authors propose the following methods:

1. **Introducing Koopman wave functions**: By combining Koopman's classical mechanics with methods on Hilbert space and using methods in symplectic geometry, a continuous model with variational and Hamiltonian structures is constructed.
2. **Regularization techniques**: Applying the previously developed regularization techniques, the singular solution hypothesis, namely "koopmons", is introduced. These are computational particles that sample the Lagrangian classical paths in the phase space.
3. **Multi - trajectory scheme**: For Tully's non - adiabatic problem, this method can reproduce the results of fully quantum simulations and is more computationally efficient than the standard MQC Ehrenfest simulation.
4. **Extension to the strong - coupling regime**: This method also successfully partially reproduces the fully quantum results in the ultrastrong - coupling and deep - strong - coupling regimes of the Rabi problem.

### Specific steps

1. **From Koopman wave functions to MQC dynamics**:
- Use the tensor product space of Koopman wave functions and Schrödinger wave functions to describe MQC dynamics.
- Make the classical phase unobservable to ensure that the classical Liouville density is always positive.

2. **Regularization and particle scheme**:
- Introduce a smooth and strictly positive convolution kernel \(K(\alpha)\) such that it converges to the delta distribution in the limit \( \alpha \to 0\).
- Use the singular solution hypothesis \(f(z, t)=\sum_{a = 1}^N w_a\delta(z-\zeta_a(t))\) and \(bP(z, t)=\sum_{a = 1}^N w_a\hat{\varrho}_a\delta(z-\zeta_a(t))\), where each koopmon carries a positive weight \(w_a\) and a quantum state \(\hat{\varrho}_a\).

### Results

1. **Tully model**: For single - avoided - crossing, double - avoided - crossing, and extended - coupling - region problems, the koopmon method can accurately reproduce the results of fully quantum simulations.
2. **Rabi model**: In the ultrastrong - coupling and deep - strong - coupling regimes, the koopmon method successfully partially reproduces the fully quantum results, demonstrating its applicability under extreme conditions.

### Summary

By introducing Koopman wave functions and regularization techniques, this method not only solves the consistency problems of traditional MQC models but also shows significant advantages in computational efficiency. This provides a new approach for the efficient simulation of non - adiabatic quantum - classical dynamics.