Search

Fewest-switches surface hopping (FSSH) has emerged as one of the leading methods for modeling the quantum dynamics of molecular systems. While its original formulation was limited to adiabatic populations, the growing interest in the application of FSSH to coherent phenomena prompts the question of how one should construct a complete density matrix based on FSSH trajectories. A straightforward solution is to define adiabatic coherences based on wavefunction coefficients. In this paper, we demonstrate that inconsistencies introduced in the density matrix through such treatment may lead to a violation of positivity. We furthermore show that a recently proposed coherent generalization of FSSH results in density matrices that satisfy positivity while yielding improved accuracy throughout much (but not all) of parameter space. [ABSTRACT FROM AUTHOR]

We study the exciton localization and resulting optical response for disordered tubular aggregates of optically active molecules. It has previously been shown that such tubular structures allow for excitons delocalized over more than a thousand molecules, owing to the combined effects of long-range dipole–dipole interactions and the higher-dimensional (not truly one-dimensional) nature of the aggregate. Such large delocalization sizes prompt the question to what extent in experimental systems the delocalization may still be determined by the aggregate size (diameter and length) and how this affects the aggregate's optical response and dynamics. We perform a systematic study of the size effects on the localization properties using numerical simulations of the exciton states in a cylindrical model structure inspired by the previously derived geometry of a cylindrical aggregate of cyanine dye molecules (C8S3). To characterize the exciton localization, we calculate the participation ratio and the autocorrelation function of the exciton wave function. We also calculate the density of states and absorption spectrum. We find strong effects of the tube's radius on the localization and optical properties in the range of parameters relevant to the experiment. In addition, surprisingly, we find that even for tubes as long as 750 nm, the localization size is limited by the tube's length for disorder values that are relevant to experimental circumstances, while observable effects of the tube's length in the absorption spectrum still occur for tube lengths up to about 150 nm. The latter may explain the changes in the optical spectra observed during the aging process of bromine-substituted C8S3 aggregates. For weak disorder, the exciton wave functions exhibit a scattered, fractal-like nature, similar to the quasi-particles in two-dimensional disordered systems. [ABSTRACT FROM AUTHOR]

In this paper, we present a novel benchmarking method for validating the modelling of vibrational spectra for the amide I region of proteins. We use the linear absorption spectra and two-dimensional infrared spectra of four experimentally well-studied proteins as a reference and test nine combinations of molecular dynamics force fields, vibrational frequency mappings, and coupling models. We find that two-dimensional infrared spectra provide a much stronger test of the models than linear absorption does. The best modelling approach in the present study still leaves significant room for future improvement. The presented benchmarking scheme, thus, provides a way of validating future protocols for modelling the amide I band in proteins. [ABSTRACT FROM AUTHOR]

We compare theoretical methods for calculating excitation energy transfer rates in multichromophoric systems. The employed methods are the multichromophoric Förster resonance energy transfer (MC-FRET), the numerical integration of the Schrödinger equation (NISE), and the Haken-Strobl-Reineker (HSR) model. As a reference, we use the numerically exact Hierarchy of Equations of Motion (HEOM). We examine these methods in various system and bath parameter regimes including the regime relevant to biological light-harvesting systems. We apply them to a model system of a monomeric donor coupled to a multichromophoric acceptor ring of varying size in two limiting configurations, namely symmetric and asymmetric. We find that the symmetric case is more sensitive to the approximations of the methods studied. The NISE method gives the most reasonable results throughout the parameter regimes tested, while still being computationally tractable. [ABSTRACT FROM AUTHOR]