Monte Carlo Wavefunction Approach to Singlet Fission Dynamics of Molecular Aggregates

We have developed a Monte Carlo wavefunction (MCWF) approach to the singlet fission (SF) dynamics of linear aggregate models composed of monomers with weak diradical character. As an example, the SF dynamics for a pentacene dimer model is investigated by considering the intermolecular electronic coupling and the vibronic coupling. By comparing with the results by the quantum master equation (QME) approach, we clarify the dependences of the MCWF results on the time step (Δt) and the number of MC trajectories (MC). The SF dynamics by the MCWF approach is found to quantitatively (within an error of 0.02% for SF rate and of 0.005% for double-triplet (TT) yield) reproduce that by the QME approach when using a sufficiently small Δt (~0.03 fs) and a sufficiently large MC (~105). The computational time (treq) in the MCWF approach also exhibits dramatic reduction with increasing the size of aggregates (N-mers) as compared to that in the QME approach, e.g., ~34 times faster at the 20-mer, and the size-dependence of treq shows significant reduction from N5.15 (QME) to N3.09 (MCWF). These results demonstrate the promising high performance of the MCWF approach to the SF dynamics in extended multiradical molecular aggregates including a large number of quantum dissipation, e.g., vibronic coupling, modes.

Calculating Trajectories for Atoms in Near-resonant Lightfields

We review several methods for calculating the time development of the internal state and the external motion of atoms in near-resonant light fields, with emphasis on studying the focussing of atomic beams into microscopic and potentially nanoscopic patterns. Three different approaches are considered: two-level semiclassical, multi-level semiclassical, and the Monte Carlo wavefunction method. The two-level semiclassical technique of McClelland and Scheinfein (1991) and McClelland (1995) is extended to three dimensions, and used to calculate the trajectories of atoms and the imaging properties of a simple lens formed from a near-resonant travelling TEM01 mode laser. The model is then extended to multi-level atoms, where we calculate the density matrix for the internal state of a sample of thermal atoms in a standing wave, and show how cooling processes can be simulated. Finally, we use the Monte Carlo wavefunction method to calculate the internal state of the atom, and compare the results and required computation time to those of the multi-level semiclassical technique.