Hybrid quantum-classical systems combine both classical and quantum degrees of freedom. Typically, in chemistry, molecular physics, or saterials science, the classical degrees of freedom describe atomic nuclei (or cations with frozen core electrons), whereas the quantum particles are the electrons. Although many possible hybrid dynamical models exist, the essential one is that the so-called Ehrenfest dynamics that results from the simple partial classical limit applied to the complete quantum Schrodinger equation. Few numerical methods are developed specifically for the mixing of this sort of systems. Here we present a preliminary study of the performance of a family of recently developed propagators: the (quasi) commutator-free Magnus expansions. These methods, however, were initially designed for nonautonomous linear equations. We employ them for the nonlinear Ehrenfest system, by approximating the state value at whenever step within the propagation, using an extrapolation from previous time steps.
Hybrid quantum-classical systems combine both classical and quantum degrees of freedom. Typically, in chemistry, molecular physics, or saterials science, the classical degrees of freedom describe atomic nuclei (or cations with frozen core electrons), whereas the quantum particles are the electrons. Although many possible hybrid dynamical models exist, the essential one is that the so-called Ehrenfest dynamics that results from the simple partial classical limit applied to the complete quantum Schrodinger equation. Few numerical methods are developed specifically for the mixing of this sort of systems. Here we present a preliminary study of the performance of a family of recently developed propagators: the (quasi) commutator-free Magnus expansions. These methods, however, were initially designed for nonautonomous linear equations. We employ them for the nonlinear Ehrenfest system, by approximating the state value at whenever step within the propagation, using an extrapolation from previous time steps.
Hybrid quantum-classical systems combine both classical and quantum degrees of freedom. Typically, in Chemistry, Molecular Physics, or Materials Science, the classical degrees of freedom describe atomic nuclei (or cations with frozen core electrons), whereas the quantum particles are the electrons. Although many possible hybrid dynamical models exist, the basic one is the so-called Ehrenfest dynamics that results from the straightforward partial classical limit applied to the full quantum Schrödinger equation. Few numerical methods have been developed specifically for the integration of this type of systems. Here we present a preliminary study of the performance of a family of recently developed propagators: the (quasi) commutator-free Magnus expansions. These methods, however, were initially designed for nonautonomous linear equations. We employ them for the nonlinear Ehrenfest system, by approximating the state value at each time step in the propagation, using an extrapolation from previous time steps.
We consider the numerical propagation of models that combine both quantum and classical degrees of freedom, usually, electrons and nuclei, respectively. We focus, in our computational examples, on the case in which the quantum electrons are modeled with time-dependent density-functional theory, although the methods discussed below can be used with any other level of theory. Often, for these so-called quantum-classical molecular dynamics models, one uses some propagation technique to deal with the quantum part and a different one for the classical equations. While the resulting procedure may, in principle, be consistent, it can however spoil some of the properties of the methods, such as the accuracy order with respect to the time step or the preservation of the geometrical structure of the equations. Few methods have been developed specifically for hybrid quantum-classical models. We propose using the same method for both the quantum and classical particles, in particular, one family of techniques that proves to be very efficient for the propagation of Schrödinger-like equations: the (quasi)-commutator free Magnus expansions. These have been developed, however, for linear systems, yet our problem is nonlinear: formally, the full quantum-classical system can be rewritten as a nonlinear Schrödinger equation, i.e., one in which the Hamiltonian depends on the system itself. The Magnus expansion algorithms for linear systems require the application of the Hamiltonian at intermediate points in a given propagating interval. For nonlinear systems, this poses a problem as this Hamiltonian is unknown due to its dependence on the state. We approximate it by employing a higher order extrapolation using previous steps as input. The resulting technique can then be regarded as a multistep technique or, alternatively, as a predictor corrector formula.
A classical limit of quantum dynamics can be defined by compensation of the quantum potential in the time-dependent Schrödinger equation. The quantum potential is a non-local quantity, defined in the trajectory-based form of the Schrödinger equation, due to Madelung, de Broglie, and Bohm, which formally generates the quantum-mechanical features in dynamics. Selective inclusion of the quantum potential for the degrees of freedom deemed "quantum," defines a hybrid quantum/classical dynamics, appropriate for molecular systems comprised of light and heavy nuclei. The wavefunction is associated with all of the nuclei, and the Ehrenfest, or mean-field, averaging of the force acting on the classical degrees of freedom, typical of the mixed quantum/classical methods, is avoided. The hybrid approach is used to examine evolution of light/heavy systems in the harmonic and double-well potentials, using conventional grid-based and approximate quantum-trajectory time propagation. The approximate quantum force is defined on spatial domains, which removes unphysical coupling of the wavefunction fragments corresponding to distinct classical channels or configurations. The quantum potential, associated with the quantum particle, generates forces acting on both quantum and classical particles to describe the backreaction.
Chemical dynamics, in principle, should be understood by solving the time-dependent Schrödinger equation for a molecular system, describing motion of the nuclei and electrons. However, the computational efforts to solve this partial second-order differential equation scales exponentially with the system size, which prevents us from getting exact numerical solutions for systems larger than 4-5 atoms. Thus, approximations simplifying the picture are necessary. The so-called Born-Oppenheimer approximation, separating motion of the electrons and nuclei is the central one: solution to the electronic Schrödinger equation defines the potential energy surface on which the nuclear motion unfolds, and there are standard quantum chemistry software packages for solving the electronic Schrödinger equation. For the nuclear Schrödinger equation, however, there are no widely applicable quantum-mechanical approaches, and most simulations are performed using classical Newtonian mechanics which is often adequate due to large nuclear masses. However, the nuclear quantum effects are significant for chemical processes involving light nuclei at low energies, and including these effects into simulation, even approximately, is highly desirable. In this dissertation, an approximate methodology of including quantum-mechanical effects within the quantum trajectory or the de Broglie-Bohm formulation of the Schrödinger equations is developed. Use of the trajectory framework makes the approach scalable to hundreds of degrees of freedom. The methodology is applied to study high-dimensional systems (solid He4 and others) relevant to chemistry.
We apply the MQCB (mixed quantum/classical mechanics based on Bohmian trajectories) [E. Gindensperger, C. Meier, and J. A. Beswick, J. Chem. Phys. 113, 9369 (2000)] to the case where the quantum degree of freedom contains both bound and continuum states. The quantum degree of freedom is treated by a wave packet that is propagated on a spatial grid with a Hamiltonian that depends parametrically on the classical degrees of freedom, while the classical degrees of freedom themselves are coupled to the wave function via quantum trajectories. The method is applied to a simple model of a light particle colliding with a heavy (classical) particle absorbed on a surface, which has been used by other authors to develop methods to combine quantum degrees of freedom including continuum states with classical mechanics. The results are compared to the Ehrenfest mean-field method as well as to full quantum results.
Introduction We now move to considering the dynamics of a system of nuclei and electrons. Of course, both electrons and nuclei are subject to the laws of quantum mechanics, but since nuclei are 2000 to 200,000 times heavier than electrons, we expect that classical mechanics will be a much better approximation for the motion of nuclei than for the motion of electrons. This means that we expect a level of approximation to be valid, where some of the degrees of freedom (d.o.f.) of a system behave essentially classically and others behave essentially quantum-mechanically. The system then is of a mixed quantum/classical nature. Most often the quantum subsystem consists of system of electrons in a dynamical field of classical nuclei, but the quantum subsystem may also be a selection of generalized nuclear coordinates (e.g., corresponding to high-frequency vibrations) while other generalized coordinates are supposed to behave classically, or describe the motion of a proton in a classical environment. So, in this chapter we consider the dynamics of a quantum system in a non-stationary potential. In Section 5.2 we consider the time-dependent potential as externally given, without taking notice of the fact that the sources of the time-dependent potential are moving nuclei, which are quantum particles themselves, feeling the interaction with the quantum d.o.f. Thus we consider the time evolution of the quantum system, which now involves mixing-in of excited states, but we completely ignore the back reaction of the quantum system onto the d.o.f. that cause the time-dependent potential, i.e., the moving nuclei.