CTMQC, a module for excited-state nonadiabatic dynamics
CTMQC is a module for excited-state nonadiabatic dynamics. It is used to simulate the coupled dynamics of electrons and nuclei (ideally in gas phase molecular systems) in response to, for instance, an initial electronic excitation.
The CTMQC module is based on the coupled-trajectory mixed quantum-classical (CT-MQC) algorithm [1,2] that has been derived starting from the evolution equations in the framework the exact factorization of the electron-nuclear wavefunction [3,4,5]. The CTMQC algorithm belongs to the family of quantum-classical methods, as the time evolution of the nuclear degrees of freedom is treated within the classical approximation, whereas electronic dynamics is treated fully quantum mechanically. Basically, the nuclei evolve as point particles, following classical trajectories, while the electrons generate the potential inducing such time evolution.
In its current implementation (used in Refs. [6,7]), the module cannot deal with arbitrary nuclear dimensions, but it is restricted to treat up to 3-dimensional problems, which gives the possibility to compare quantum-classical results easily and directly with quantum wavepacket dynamics. CTMQC has been analyzed and benchmarked against exact propagation results on typical low-dimensional model systems [1,2,6,7], and applied for the simulation of the photo-initiated ring-opening process of Oxirane . For this study, CTMQC has been implemented in a developer version of the CPMD electronic structure package based on time-dependent density functional theory. Concerning electronic input properties, the CTMQC module requires a grid representation of the adiabatic potential energy surfaces and of the nonadiabatic coupling vectors, since the electronic dynamics is represented and solved in the adiabatic basis.
This feature allows the algorithm to be easily adaptable, in the current form, to any quantum chemistry electronic structure package. The number of electronic states to be included is not limited and can be specified as input.
Practical application and exploitation of the code
Software documentation can be found in our E-CAM software Library here.