Release of TheoDORE 2.0 (beta)

Version 2.0 of the TheoDORE wavefunction analysis package has been released, download below. The two main features of TheoDORE 2.0 are the computation of conditional electron densities and compatibility with python3.

Conditional electron densities can be used for the visualisation of excited-state electron correlation, see ChemPhotoChem (2019). Below, the application of this method to a PPV oligomer is shown. Here, the probe hole (red) is always fixed on the terminal phenyl ring and the different shapes for the conditional electron density (blue) for the first six excited states is observed. One can see that for the different states the electron is either repelled, attracted or unaffected by the hole.

  • S1
  • S2
  • S3
  • S4
  • S5
  • S6

Continue reading

Paper: Dynamics within an Exciton Model

Another paper working on improving the efficiency of surface hopping dynamics just appeared, this time in JCTC: “Surface hopping within an exciton picture – An electrostatic embedding scheme.” authored by M. F. S. J. Menger, F. Plasser, B. Mennucci, and L. González. In this paper, we explored the possibility of running nonadiabatic dynamics simulations within an exciton model. The main challenge in this endeavour was to derive a consistent energy expression for combining QM/MM electrostatic embedding calculations of the different chromophores.

Surface Hopping within an Exciton Scheme

To test the implementation, we ran simulations on a molecular dyad, where full TDDFT nonadiabatic dynamics simulations were available. Good agreement was found.

The method was implemented in the SHARC molecular dynamics package.

Paper: Vibronic coupling constants

You can find our new paper “Interstate vibronic coupling constants between electronic excited states for complex molecules” that recently appeared in JCP. The purpose of this paper was the development of a method that allows to determine interstate vibronic coupling constants, which are a decisive ingredient for model Hamiltonians used in quantum dynamics. Our idea was to start with a method based on wavefunction overlaps that is commonly used for trajectory dynamics simulations and adapt it for the case of quantum dynamics.