For an ensemble of B850 rings of the light-harvesting system LH2 of purple bacteria the linear absorption spectrum is calculated. Using different Markovian and non-Markovian, time-dependent and time-independent methods based on second-order perturbation theory in the coupling between the excitonic system and its surrounding environment as well as the modified Redfield theory, the influence of the shape of the spectral density on the linear absorption spectrum is demonstrated for single samples and in the ensemble average. For long bath correlation times non-Markovian effects clearly show up in the static absorption line shapes. Among the different spectral densities studied is one of the purple bacterium Rhodospirillum molischianum obtained by a molecular-dynamics simulation earlier. The effect of static disorder on its line shapes in the ensemble average is analyzed and the results of the present calculations are compared to experimental data.

Full dimensional multi-configuration time-dependent Hartree calculations of the zero point energy and the tunneling splitting of malonaldehyde using a recently published potential energy surface [Y. Wang, B. J. Braams, J. M. Bowman, S. Carter, and D. P. Tew, J. Chem. Phys. 128, 224314 (2008)] are reported. The potential energy surface has been approximated by a modified version of the n-mode representation and careful convergence check has been performed to ensure accurate results. The obtained value for the splitting (23.4 cm(-1)) is in acceptable agreement with the experimental value of 21.583 cm(-1). The computed zero-point-energy is 14,670 cm(-1) which is lower than previous results of Wang et al., but likely to be about 4 cm(-1) too low because of shortcomings of the n-mode representation of the potential. The energies reported in this abstract contain a correction to account for neglected vibrational angular momentum terms.

A Monte Carlo method is proposed for transforming high-dimensional potential energy surfaces evaluated on discrete grid points into a sum-of-products form, more precisely into a Canonical Polyadic Decomposition form. To this end, a modified existing ansatz based on the alternating least squares method is used, in which numerically exact integrals are replaced with Monte Carlo integrals. This largely reduces the numerical cost by avoiding the evaluation of the potential on all grid points and allows the treatment of surfaces with many degrees of freedom. Calculations on the 15D potential of the protonated water dimer (Zundel cation) in a sum-of-products form are presented and compared to the results obtained in a previous work [M. Schröder and H.-D. Meyer, J. Chem. Phys. 147, 064105 (2017)], where a sum-of-products form of the potential was obtained in the Tucker format.

Several techniques to solve a hierarchical set of equations of motion for propagating a reduced density matrix coupled to a thermal bath have been developed in recent years. This is either done using the path integral technique as in the original proposal by Tanimura and Kubo [J. Phys. Soc. Jpn. 58, 101 (1998)] or by the use of stochastic fields as done by Yan et al. [Chem. Phys. Lett. 395, 216 (2004)]. Based on the latter ansatz a compact derivation of the hierarchy using a decomposition of the spectral density function is given in the present contribution. The method is applied to calculate the time evolution of the reduced density matrix describing the motion in a harmonic, an anharmonic, and two coupled oscillators where each system is coupled to a thermal bath. Calculations to several orders in the system-bath coupling with two different truncations of the hierarchy are performed. The respective density matrices are used to calculate the time evolution of various system properties and the results are compared and discussed with a special focus on the convergence with respect to the truncation scheme applied.

We report energies and tunneling splittings of vibrational excited states of malonaldehyde which have been obtained using full dimensional quantum mechanical calculations. To this end we employed the multi configuration time-dependent Hartree method. The results have been obtained using a recently published potential energy surface [Y. Wang, B. J. Braams, J. M. Bowman, S. Carter, and D. P. Tew, J. Chem. Phys. 128, 224314 (2008)] which has been brought into a suitable form by a modified version of the n-mode representation which was used with two different arrangements of coordinates. The relevant terms of the expansion have been identified with a Metropolis algorithm and a diffusion Monte-Carlo technique, respectively.