The point is that the picture of isolated electronic states breaks down, i.e. the Born-Oppenheimer approximation. In solid state chemistry that is what they call electron-phonon-coupling.

If you consider the nuclei in terms of an external perturbation, you can take at look at the adiabatic theorem: If nuclear motion is fast and there are close lying electronic states, the states will mix. You get the fast nuclear motion from the excitation energy. About close lying states: Typically all the excited states are close and the molecule quickly reaches S

_{1}if it was initially excited to a higher state (Kasha's rule). If there are suitable intersections to the ground state it can even relax completely.

For doing molecular dynamics you need the energy of your electronic state Ψ

_{i}as always. But you also need couplings to the other states:

They are called the non-adiabatic coupling vectors and give the major post Born-Oppenheimer contributions.

Here is a little dynamics movie. This is a test MCSCF run on cytosine. The electronic structure is from the Columbus program and the non-adiabatic surface-hopping dynamics are with Newton-X (both from our group actually).

The idea of surface-hopping is that we still want to have a classical trajectory as a basis, but post Born-Oppenheimer corrections are introduced through jumping in between states. Here you can see the cytosine molecule with color coded electronic states: green S

_{1}, orange S

_{2}. You see how the molecule jumps into the second excited state intermediately and then stays in the first excited state. There is a lot of motion because of the high excess energy.

There is no decay to the ground state in these 341 fs simulated. It should happen soon after or tanning would be extremely dangerous. But the time simulated is still rather short, it is for example a ten-thousandth of the typical fluorescence life-time.

## 2 comments:

Cool post! How long did the calculation for the movie take?

And what do you mean by "The electronic structure is from the Columbus program" ? Don't you need a PES for the dynamics?

Greetings

Till

thanks.

this trajectory took one month on a single processor. one hour per time step, this is actually kind of too slow if you run more trajectories to get some statistics

these are semi-classical on-the-fly dynamics. they are classical except for a post Born-Oppenheimer correction.

you compute the gradient at every geometry and propagate it accordingly just like in regular classical dynamics. but you also compute non-adiabatic couplings as the quantum part.

this is called suface hopping dynamics (http://dx.doi.org/10.1039/a801824c). for the implementation: http://dx.doi.org/10.1016/j.jphotochem.2006.12.008

Post a Comment