Optimization of energy transport in the Fenna-Matthews-Olson complex via site-varying pigment-protein interactions.

Energy transport in photosynthetic systems can be tremendously efficient. In particular, we study exciton transport in the Fenna-Mathews-Olson (FMO) complex found in green sulphur bacteria. The exciton dynamics and energy transfer efficiency depend on the interaction of excited chromophores with their environment. Based upon realistic, site-dependent models of the system-bath coupling, we present results that suggest that this interaction may be optimized in the case of FMO. Furthermore we verify two transport pathways and note that one is dominated by coherent dynamics and the other by incoherent energy dissipation. In particular, we note a significant correlation between energy transport efficiency and coherence for exciton transfer from bacteriochlorophyll (BChl) 8 to BChl 4.

Oxygen defects in phosphorene.

Surface reactions with oxygen are a fundamental cause of the degradation of phosphorene. Using first-principles calculations, we show that for each oxygen atom adsorbed onto phosphorene there is an energy release of about 2 eV. Although the most stable oxygen adsorbed forms are electrically inactive and lead only to minor distortions of the lattice, there are low energy metastable forms which introduce deep donor and/or acceptor levels in the gap. We also propose a mechanism for phosphorene oxidation involving reactive dangling oxygen atoms and we suggest that dangling oxygen atoms increase the hydrophilicity of phosphorene.

Iterative linearized density matrix propagation for modeling coherent excitation energy transfer in photosynthetic light harvesting.

Rather than incoherent hopping between chromophores, experimental evidence suggests that the excitation energy transfer in some biological light harvesting systems initially occurs coherently, and involves coherent superposition states in which excitation spreads over multiple chromophores separated by several nanometers. Treating such delocalized coherent superposition states in the presence of decoherence and dissipation arising from coupling to an environment is a significant challenge for conventional theoretical tools that either use a perturbative approach or make the Markovian approximation. In this paper, we extend the recently developed iterative linearized density matrix (ILDM) propagation scheme [E. R. Dunkel et al., J. Chem. Phys. 129, 114106 (2008)] to study coherent excitation energy transfer in a model of the Fenna-Matthews-Olsen light harvesting complex from green sulfur bacteria. This approach is nonperturbative and uses a discrete path integral description employing a short time approximation to the density matrix propagator that accounts for interference between forward and backward paths of the quantum excitonic system while linearizing the phase in the difference between the forward and backward paths of the environmental degrees of freedom resulting in a classical-like treatment of these variables. The approach avoids making the Markovian approximation and we demonstrate that it successfully describes the coherent beating of the site populations on different chromophores and gives good agreement with other methods that have been developed recently for going beyond the usual approximations, thus providing a new reliable theoretical tool to study coherent exciton transfer in light harvesting systems. We conclude with a discussion of decoherence in independent bilinearly coupled harmonic chromophore baths. The ILDM propagation approach in principle can be applied to more general descriptions of the environment.

Ultrafast H2 and D2 rotational Raman responses in near critical CO2: an experimental and theoretical study of anisotropic solvation dynamics.

The optical heterodyne detected anisotropic rotational Raman responses of H(2) and D(2) (22 mol %) in a near critical CO(2) (rho(*) = rho/rho(c) = 0.8, T = 308 K) solution are reported. J-specific rotational Raman correlation functions (RCFs) for the S(J) transitions of H(2) (J = 0,1,2) and D(2) (J = 0,1,2,3) in this CO(2) solution are determined from these measurements. A mixed classical-quantum simulation methodology results in RCFs that are in excellent agreement with the experimentally derived J-specific responses. The observed S(J) coherence decay time scales, J-dependence, rotor mass dependence, and solvent-induced transition frequency shifts are well captured by these simulations. Pure dephasing of these rotational Raman transitions is shown to be close to the homogeneous limit of the standard Kubo line shape analysis and attributable to the rotor center-of-mass translation in an anisotropic solvent cage. Rotor translational motion in the vicinity of a single CO(2) appears to dominate this dephasing mechanism. Mixed classical-quantum simulations, incorporating the effects of solution fluctuation driven nonadiabatic coupling of instantaneous adiabatic states, including full J-mixing, are required for the agreement between theory and experiment obtained here. Simulations of the classically excited angular kinetic energy of D(2) rotors are used as an estimate of T(1) relaxation rates and are found to be negligible compared to the D(2) rotational Raman coherence time scale. These results are discussed in the context of previous mixed classical-quantum and rotational friction calculations of the dephasing and energy relaxation contributions to H(2) rotational Raman coherence decays. Advantages of time domain acquisition of these rotational Raman responses as compared to spontaneous Raman measurements are illustrated here.

Quantum initial condition sampling for linearized density matrix dynamics: Vibrational pure dephasing of iodine in krypton matrices.

This paper reviews the linearized path integral approach for computing time dependent properties of systems that can be approximated using a mixed quantum-classical description. This approach is applied to studying vibrational pure dephasing of ground state molecular iodine in a rare gas matrix. The Feynman-Kleinert optimized harmonic approximation for the full system density operator is used to sample initial conditions for the bath degrees of freedom. This extremely efficient approach is compared to alternative initial condition sampling techniques at low temperatures where classical initial condition sampling yields dephasing rates that are nearly an order of magnitude too slow compared to quantum initial condition sampling and experimental results.

Iterative linearized approach to nonadiabatic dynamics.

This paper presents a new approach to propagating the density matrix based on a time stepping procedure arising from a Trotter factorization and combining the forward and backward incremental propagators. The sums over intermediate states of the discrete quantum subsystem are implemented by a Monte Carlo surface hopping-like procedure, while the integrals over the continuous variables are performed using a linearization in the difference between the forward and backward paths of these variables leading to classical-like equations of motion with forces determined by the quantum subsystem states. The approach is tested on several models and numerical convergence is explored.

Creating a Stable Oxide at the Surface of Black Phosphorus.

The stability of the surface of in situ cleaved black phosphorus crystals upon exposure to atmosphere is investigated with synchrotron-based photoelectron spectroscopy. After 2 days atmosphere exposure a stable subnanometer layer of primarily P2O5 forms at the surface. The work function increases by 0.1 eV from 3.9 eV for as-cleaved black phosphorus to 4.0 eV after formation of the 0.4 nm thick oxide, with phosphorus core levels shifting by <0.1 eV. The results indicate minimal charge transfer, suggesting that the oxide layer is suitable for passivation or as an interface layer for further dielectric deposition.

LAND-map, a linearized approach to nonadiabatic dynamics using the mapping formalism.

We present a new approach for calculating quantum time correlation functions for systems whose dynamics exhibits relevant nonadiabatic effects. The method involves partial linearization of the full quantum path-integral expression for the time correlation function written in the nonadiabatic mapping Hamiltonian formalism. Our analysis gives an algorithm which is both numerically efficient and accurate as we demonstrate in test calculations on the spin-boson model where we find results in good agreement with exact calculations. The accuracy of our new approach is comparable to that of calculations performed using other approximate methods over a relatively broad range of model parameters. However, our method converges relatively quickly when compared with most alternative schemes. These findings are very encouraging in view of the application of the new method for studying realistic nonadiabatic model problems in the condensed phase.

Ultrafast nonadiabatic dynamics: quasiclassical calculation of the transient photoelectron spectrum of I2(-).(CO2)8.

In this paper we investigate the transient photoelectron spectrum of I2(-) in CO2 clusters recently measured by Neumark and co-workers. This work reveals a rich excited state dynamics with various competing electronic output channels. We find good agreement with experiments and we are able to relate the transient signal to different dynamical events that occur during the evolution of the cluster and its fragmentation products.