Surface hopping simulations on photoexcitation dynamics of conjugated polymer
Introduction
In recent decades, surface hopping method aroused much attention in scientific community and was in rapid developing in the field of nonadiabatic molecular dynamics [1]. Surface hopping method is the so-called semiclassical (mixed quantum-classical) approach, where the nuclear motions are described by classical mechanics and the electrons quantum mechanically. The nuclei evolve on a certain potential energy surface (PES), but the transitions between different adiabatic states are allowed. Ehrenfest method is another nonadiabatic molecular dynamics methods, where the nuclei evolve on an effective PES which is an average of all adiabatic states involved weighted by their populations (also called mean-field method).
The relaxation process of a molecular system from a high-lying excited state to lower excited state is controlled by the strength of nonadiabatic coupling terms (NACTs) between different excited states. Some excited states have strong nonadiabatic couplings, some are noninteracting, depending on the characters of their wave functions and the energy level difference between them. When two noninteracting excited states are approaching degeneracy, the nonadiabatic coupling between them presents sharp peaks that are strongly localized in time. They actually cross each other, called “trivial" crossings. These trivial crossings may lead to “unphysical" hoppings and long range energy transfer in numerical simulations, and they must be avoided.
Different methods have been proposed to deal with the trivial crossing problem [2], [3], [4], [5], [6], [7], [8], [9]. Our group proposed a method in which the excited states are classified into two groups according to their transition density matrices and different group of excited states are noninteracting [9]. Because of this, the system can separately evolve within one group of excited states, and the trivial crossings between noninteracting states are completely eliminated. However, this method only applies to those systems whose excited states can be classified into noninteracting groups, such as a single polymer chain. For those systems composed by interacting polymer chains, it is not applicable. Another method is the so-called Min-Cost algorithm, which is used in the present work [4]. The key point of this algorithm is that the identities of excited states must be followed during dynamics based on the characters of their wave functions, instead of the energy-ordering criterion. In other words, at each time step the excited states are tracked according to a maximum overlap principle so that the old states can be mapped to the new states. However, there is a problem here that can very easily be ignored by people. At time zero, we are not able to track the excited states like at any other time and have to assign them based on energy-ordering criterion.This may lead to fake populations of an excited state during dynamical simulations. We will discuss this problem and give a solution to it in the present work.
People often have to consider how to reduce the computational load because the computational cost is a bottleneck for most surface hopping simulations. In surface hopping simulations, the most time-consuming parts are calculating the NACTs and the forces exerting on atoms. The higher the level of excited-state calculations, the more expensive the calculating NACTs and forces. In general, configuration interaction singles (CIS) is thought to be the zeroth-order approximation for excited-state calculations. It is generally used in surface hopping approaches [10], [11]. The time-dependent density functional theory (TDDFT) is also an efficient approach for calculating excited states of complex molecular systems, which has similar computational cost with CIS. TDDFT is also widely used in surface hopping approaches [12], [13]. However, even for CIS and TDDFT, surface hopping simulations on conjugated polymers are still expensive.
In order to speed up the surface hopping simulations, people have developed many approximation methods, especially on calculation the NACTs since they are the most time-consuming step. Pittner et al. proposed a scheme to neglect all NACTs between other states than the current state.[14] For a surface hopping simulation involving Ns states, taking into account that NACT σJJ = 0 and σIJ = σJI, Ns(Ns − 1)∕2 coupling terms should in principle be evaluated. Using this scheme, the number of coupling terms that should be computed is reduced to Ns − 1. Fabiano et al. neglected the second term of the nonadiabatic coupling vector expression, since they thought the first term dominate the expression in the region of a conical intersection [2]. Nelson et al. thought that the contribution of MO response to the NACTs, reflected by the Z matrix, is small. They truncated the summation up to the number of involved excited states when calculating the Z matrix [11]. Lobaugh and Rossky used approximate CIS wave functions to calculate NACTs and forces. These wave functions are generated using only 50 single excitation configurations [10].
In this work, we fully introduce the surface hopping method which is based on one-dimensional Pariser-Parr-Pople (PPP) Hamiltonian and CIS formalism. This method is easy to implement and has low computational cost [9], [15], [16]. In this method, we deduced analytic expressions for the forces and the NACTs. Using them, the computation speed is greatly accelerated. In order to further reduce the computational cost, we discussed the contribution of Z matrix to see if it is possible to neglect it.
In Section II, we fully describe the implementation of the surface hopping method. In Section III, the results of numerical simulations of photoexcitation dynamics in cis-polyacetylene are presented. Finally, in Section IV we conclude our investigation.
Section snippets
Total Hamiltonian
Before proceeding, we summarize our notation conventions: {i, j, k, . . . } index occupied molecular orbitals (MOs); {a, b, c, . . . } virtual MOs; and {μ, ν, λ, . . . } atomic sites. For simplicity, we consider a cis-polyacetylene chain with 20 sites. All π electrons of this polymer chain are described by the PPP Hamiltonian [17]. The sites are treated as classical nuclei cores. The total Hamiltonian H = He + Hn, where He and Hn are written as
Results and discussion
We start with analysis of excited state properties of the cis-polyacetylene chain at it’s ground-state optimized geometry. Table 1 shows the energy levels and oscillator strengths for the lowest 6 excited states. We can see that states S1, S5 and S6 are “bright” states (f0K ≠ 0). These excited states can be directly photoexcited upon absorbing a photon. States S2, S3 and S4 are “dark” states (f0K = 0). They can not be directly photoexcited upon absorbing a photon due to the forbidden transition
Summary
In summary, we have derived all necessary equations for the surface hopping method based on PPP Hamiltonian and CIS excited-state approximation, and then use this method to simulate the excited-state dynamics in a cis-polyacetylene chain. We find that an excited state is falsely populated when the system is relaxed from a high-lying excited state. It happens because it is impossible to track excited states at time zero using the traditional Min-Cost method. To solve this problem, we need to
CRediT authorship contribution statement
Zhen Sun: Conceptualization, Methodology, Software, Writing - original draft. Sheng Li: Resources, Formal analysis. Shijie Xie: Resources, Formal analysis, Writing - review & editing. Zhong An: Formal analysis, Validation.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgements
This work was supported by Natural Science Foundation of Zhejiang Province (LY19A040007), and National Natural Science Foundation of China (NSF) (11974212).
References (32)
- et al.
Implementation of surface hopping molecular dynamics using semiempirical methods
Chem. Phys.
(2008) - et al.
Artifacts due to trivial unavoided crossings in the modeling of photoinduced energy transfer dynamics in extended conjugated molecules
Chem. Phys. Lett.
(2013) - et al.
Optimization of mixed quantum-classical dynamics: time-derivative coupling terms and selected couplings
Chem. Phys.
(2009) - et al.
Constant temperature simulations using the langevin equation with velocity verlet integration
Chem. Phys.
(1998) - et al.
A study on excited-state properties of conjugated polymers using the pariser-parr-pople-peierls model combined with configuration-interaction-singles
Org. Electron.
(2018) - et al.
Exploring interchain polaron pair formation in neat conjugated polymers
Eur. Phys. J. B
(2020) - et al.
Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications
J. Phys.-Condens. Matter
(2019) - et al.
Direct semiclassical simulation of photochemical processes with semiempirical wave functions
J. Chem. Phys.
(2001) - et al.
Identification of unavoided crossings in nonadiabatic photoexcited dynamics involving multiple electronic states in polyatomic conjugated molecules
J. Chem. Phys.
(2012) - et al.
Surface hopping dynamics using a locally diabatic formalism: charge transfer in the ethylene dimer cation and excited state dynamics in the 2-pyridone dimer
J. Chem. Phys.
(2012)
Identification of unavoided crossings in nonadiabatic photoexcited dynamics involving multiple electronic states in polyatomic conjugated molecules
J. Chem. Phys.
A simple solution to the trivial crossing problem in surface hopping
J. Phys. Chem. Lett.
Solution for the trivial crossing problem in surface hopping simulations by the classification on excited states
J. Phys. Chem. C
Computer simulation of the excited state dynamics of betaine-30 in acetonitrile
J. Phys. Chem. A
Nonadiabatic excited-state molecular dynamics modeling of photoinduced dynamics in conjugated molecules
J. Phys. Chem. B
Nonadiabatic dynamics and simulation of time resolved photoelectron spectra within time-dependent density functional theory: ultrafast photoswitching in benzylideneaniline
J. Chem. Phys.
Cited by (6)
Surface hopping simulations on charge photogeneration in conjugated polymers
2024, Physical Chemistry Chemical PhysicsStrain-Tuneable Bipolaron Stability on Ultranarrow Bilayer Graphene Nanoribbon
2024, Journal of Physical Chemistry CThe effect of infrared push pulse on the relaxed exciton in single-component organic solar cells
2023, Applied Physics LettersSurface hopping simulations on charge separation in an organic donor-acceptor system
2023, Physical Chemistry Chemical PhysicsA reciprocal-space formulation of surface hopping
2022, Journal of Chemical Physics