Search papers, labs, and topics across Lattice.
This paper introduces an efficient on-the-fly nonadiabatic coupling (NAC) framework implemented directly within the CP2K software package. By integrating the NAC module with CP2K's Gaussian and Plane Waves (GPW) method and incorporating a phase correction scheme, the framework enables efficient and stable NAC calculations for large-scale systems, overcoming limitations of traditional approaches. Benchmark results on crystalline pentacene and larger supercells validate the framework's accuracy and scalability for nonadiabatic molecular dynamics simulations.
Calculating nonadiabatic couplings for large molecular dynamics simulations just got a whole lot faster thanks to a new on-the-fly implementation directly within CP2K.
Nonadiabatic molecular dynamics (NAMD) is widely used to describe hot electron relaxation and nonradiative recombination processes, but high computational costs limit its application to large supercells. Here, we implement a nonadiabatic coupling (NAC) module directly into CP2K, enabling on-the-fly NAC calculation during ab-initio molecular dynamics. Unlike conventional approaches relying on interfaces between external NAMD programs and electronic structure codes, this integration streamlines the workflow. Combined with CP2K's inherent Gaussian and Plane Waves (GPW) method, it achieves efficient NAC calculations for large-scale systems. To ensure numerical stability, a phase correction scheme is introduced to remove inconsistencies of wavefunction phases. Benchmark calculations of hot electron relaxation in crystalline pentacene agree with previous studies, while the results of simulations with and without phase correction show significant difference, confirming the necessity of the phase correction. Simulations on larger pentacene supercells further demonstrate the capability of this implementation for large-scale NAMD simulations.