We present a new partially linearized mapping-based approach for approximating real-time quantum correlation functions in condensedphase nonadiabatic systems, called the spin partially linearized density matrix (spin-PLDM) approach. Within a classical trajectory picture, partially linearized methods treat the electronic dynamics along forward and backward paths separately by explicitly evolving two sets of mapping variables. Unlike previously derived partially linearized methods based on the Meyer-Miller-Stock-Thoss mapping, spin-PLDM uses the Stratonovich-Weyl transform to describe the electronic dynamics for each path within the spin-mapping space; this automatically restricts the Cartesian mapping variables to lie on a hypersphere and means that the classical equations of motion can no longer propagate the mapping variables out of the physical subspace. The presence of a rigorously derived zero-point energy parameter also distinguishes spin-PLDM from other partially linearized approaches. These new features appear to give the method superior accuracy for computing dynamical observables of interest when compared with other methods within the same class. The superior accuracy of spin-PLDM is demonstrated in this paper through application of the method to a wide range of spin-boson models as well as to the Fenna-Matthews-Olsen complex.

The exciton relaxation dynamics of photoexcited electronic states in poly(p-phenylenevinylene) are theoretically investigated within a coarse-grained model, in which both the exciton and nuclear degrees of freedom are treated quantum mechanically. The Frenkel-Holstein Hamiltonian is used to describe the strong exciton-phonon coupling present in the system, while external damping of the internal nuclear degrees of freedom is accounted for by a Lindblad master equation. Numerically, the dynamics are computed using the time evolving block decimation and quantum jump trajectory techniques. The values of the model parameters physically relevant to polymer systems naturally lead to a separation of time scales, with the ultra-fast dynamics corresponding to energy transfer from the exciton to the internal phonon modes (i.e., the C-C bond oscillations), while the longer time dynamics correspond to damping of these phonon modes by the external dissipation. Associated with these time scales, we investigate the following processes that are indicative of the system relaxing onto the emissive chromophores of the polymer: (1) Exciton-polaron formation occurs on an ultra-fast time scale, with the associated exciton-phonon correlations present within half a vibrational time period of the C-C bond oscillations. (2) Exciton decoherence is driven by the decay in the vibrational overlaps associated with exciton-polaron formation, occurring on the same time scale. (3) Exciton density localization is driven by the external dissipation, arising from "wavefunction collapse" occurring as a result of the system-environment interactions. Finally, we show how fluorescence anisotropy measurements can be used to investigate the exciton decoherence process during the relaxation dynamics.

In a previous paper [J. R. Mannouch and J. O. Richardson, J. Chem. Phys. 153, 194109 (2020)], we derived a new partially linearized mappingbased classical-trajectory technique called the spin partially linearized density matrix (spin-PLDM) approach. This method describes the dynamics associated with the forward and backward electronic path integrals using a Stratonovich-Weyl approach within the spin-mapping space. While this is the first example of a partially linearized spin-mapping method, fully linearized spin-mapping is already known to be capable of reproducing dynamical observables for a range of nonadiabatic model systems reasonably accurately. Here, we present a thorough comparison of the terms in the underlying expressions for the real-time quantum correlation functions for spin-PLDM and fully linearized spin mapping in order to ascertain the relative accuracy of the two methods. In particular, we show that spin-PLDM contains an additional term within the definition of its real-time correlation function, which diminishes many of the known errors that are ubiquitous for fully linearized approaches. One advantage of partially linearized methods over their fully linearized counterparts is that the results can be systematically improved by re-sampling the mapping variables at intermediate times. We derive such a scheme for spin-PLDM and show that for systems for which the approximation of classical nuclei is valid, numerically exact results can be obtained using only a few "jumps." Additionally, we implement focused initial conditions for the spin-PLDM method, which reduces the number of classical trajectories that are needed in order to reach convergence of dynamical quantities, with seemingly little difference to the accuracy of the result.

We present a nonadiabatic classical-trajectory approach that offers the best of both worlds between fewest-switches surface hopping (FSSH) and quasiclassical mapping dynamics. This mapping approach to surface hopping (MASH) propagates the nuclei on the active adiabatic potential-energy surface, like in FSSH. However, unlike in FSSH, transitions between active surfaces are deterministic and occur when the electronic mapping variables evolve between specified regions of the electronic phase space. This guarantees internal consistency between the active surface and the electronic degrees of freedom throughout the dynamics. MASH is rigorously derivable from exact quantum mechanics as a limit of the quantum-classical Liouville equation (QCLE), leading to a unique prescription for momentum rescaling and frustrated hops. Hence, a quantum-jump procedure can in principle be used to systematically converge the accuracy of the results to that of the QCLE. This jump procedure also provides a rigorous framework for deriving approximate decoherence corrections similar to those proposed for FSSH. We apply MASH to simulate the nonadiabatic dynamics in various model systems and show that it consistently produces more accurate results than FSSH at a comparable computational cost.

We develop a model of excitons coupled to the rotational motion of monomers to study the torsionally induced relaxation and decoherence of excitons in π-conjugated polymers. The model assumes that the monomer units are described by elastically uncoupled harmonic oscillators and that there is a linear exciton-roton coupling. Although the rotational degrees of freedom are much slower than the exciton, so that the adiabatic approximation is generally expected to be valid, we also investigate possible quantized roton corrections via coupled time evolving block decimation-Ehrenfest equations of motion. For the relaxation of the lowest-excited exciton, we find that (1) for a polymer chain with a ground state spiral torsional conformation, the equilibrium angular displacement of each monomer is proportional to the difference of the exciton bond-orders on the neighboring bridging bonds. Consequently, this displacement vanishes in the long chain limit and a classical (Landau) exciton-polaron is not formed.(2) For a polymer chain with a ground state staggered torsional conformation, the equilibrium angular displacement of each monomer is proportional to the sum of the exciton bond-orders on the neighboring bridging bonds. Consequently, there is significant angular displacement and local planarization causing exciton density localization. A classical (Landau) exciton-polaron is formed where the staggered angular displacement is proportional to the exciton density. (3) Generally, in the adiabatic limit, the decay of off-diagonal long-range order (i.e., exciton decoherence) mirrors the localization of the exciton density. However, quantum corrections to the rotational motion alter this adiabatic prediction because of correlated exciton-roton dynamics within the first rotational half-period. In particular, exciton-polaron quasiparticle formation causes more rapid and oscillatory exciton decoherence and slower exciton density localization.

Photosynthetic organisms are known to use a mechanism of vibrationally assisted exciton energy transfer to efficiently harvest energy from light. The importance of quantum effects in this mechanism is a long-standing topic of debate, which has traditionally focused on the role of excitonic coherences. Here, we address another recent claim: that the efficient energy transfer in the Fenna–Matthews–Olson complex relies on nuclear quantum uncertainty and would not function if the vibrations were classical. We present a counter-example to this claim, showing by trajectory-based simulations that a description in terms of quantum electrons and classical nuclei is indeed sufficient to describe the funneling of energy to the reaction center. We analyze and compare these findings to previous classical-nuclear approximations that predicted the absence of an energy funnel and conclude that the key difference and the reason for the discrepancy is the ability of the trajectories to properly account for Newton’s third law.

The periodic Anderson model (PAM) is studied within the framework of dynamical mean-field theory, with particular emphasis on the interaction-driven Mott transition it contains, and on resultant Mott insulators of both Mott-Hubbard and charge-transfer type. The form of the PAM phase diagram is first deduced on general grounds using two exact results, over the full range of model parameters and including metallic, Mott, Kondo and band insulator phases. The effective low-energy model which describes the PAM in the vicinity of a Mott transition is then shown to be a one-band Hubbard model, with effective hoppings that are not in general solely nearest neighbour, but decay exponentially with distance. This mapping is shown to have a range of implications for the physics of the problem, from phase boundaries to single-particle dynamics; all of which are confirmed and supplemented by NRG calculations. Finally we consider the locally degenerate, non-Fermi liquid Mott insulator, to describe which requires a two-self-energy description. This is shown to yield a number of exact results for the associated local moment, charge, and interaction-renormalised levels, together with a generalisation of Luttinger's theorem to the Mott insulator.

No abstract