Iterative Linearized Density Matrix Propagation for Modeling Coherent Energy Transfer in Photosynthetic Light Harvesting

We present results of calculations [1] that employ a new mixed quantum classical iterative density matrix propagation approach (ILDM , or so called Is‐Landmap) [2] to explore the survival of coherence in different photo synthetic models. Our model studies confirm the long lived quantum coherence , while conventional theoretical tools (such as Redfield equation) fail to describe these phenomenon [3,4]. Our ILDM method is a numerical exactly propagation scheme and can be served as a bench mark calculation tools[2]. Result get from ILDM and from other recent methods have been compared and show agreement with each other[4,5]. Long lived coherence plateau has been attribute to the shift of harmonic potential due to the system bath interaction, and the harvesting efficiency is a balance between the coherence and dissipation[1]. We use this approach to investigate the excitation energy transfer dynamics in various light harvesting complex include Fenna‐Matthews‐Olsen light harvesting complex[1] and Cryptophyte Phycocyanin 645 [6]. [1] P.Huo and D.F.Coker ,J. Chem. Phys. 133, 184108 (2010) . [2] E.R. Dunkel, S. Bonella, and D.F. Coker, J. Chem. Phys. 129, 114106 (2008). [3] A. Ishizaki and G.R. Fleming, J. Chem. Phys. 130, 234111 (2009). [4] A. Ishizaki and G.R. Fleming, Proc. Natl. Acad. Sci. 106, 17255 (2009). [5] G. Tao and W.H. Miller, J. Phys. Chem. Lett. 1, 891 (2010). [6] P.Huo and D.F.Coker in preparationNational Science Foundation (CHE-0911635

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 with 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 with quantum initial condition sampling and experimental results.Comment: 20 pages and 8 figure

Single electron states in polyethylene

We report computer simulations of an excess electron in various structural motifs of polyethylene at room temperature, including lamellar and interfacial regions between amorphous and lamellae, as well as nanometre-sized voids. Electronic properties such as density of states, mobility edges, and mobilities are computed on the different phases using a block Lanczos algorithm. Our results suggest that the electronic density of states for a heterogeneous material can be approximated by summing the single phase density of states weighted by their corresponding volume fractions. Additionally, a quantitative connection between the localized states of the excess electron and the local atomic structure is presented.The US National Science Foundation under grant CHE-0911635 and from his Stokes Professorship in Nano Biophysics from Science Foundation Ireland thanks the Irish Centre for High End Computing (ICHEC) for computer resources and Science Foundation Ireland for support from grant 08-IN.1-I1869

Computational Study Of Molecular Hydrogen In Zeolite Na-A. II. Density Of Rotational States And Inelastic Neutron Scattering Spectra

Part I of this series [J. Chem. Phys. 111, 7599 (1999)] describes a simulation of H(2) adsorbed within zeolite Na-A in which a block Lanczos procedure is used to generate the first several (9) rotational eigenstates of H(2), modeled as a rigid rotor, and equilibrated at a given temperature via Monte Carlo sampling. Here, we show that rotational states are strongly perturbed by the electrostatic fields in the solid. Wave functions and densities of rotational energy states are presented. Simulated neutron spectra are compared with inelastic neutron scattering data. Comparisons are made with IR spectra in which rotational levels may appear due to rovibrational coupling. (C) 2001 American Institute of Physics

Computational Study Of Molecular Hydrogen In Zeolite Na-A. I. Potential Energy Surfaces And Thermodynamic Separation Factors For Ortho And Para Hydrogen

We simulate H-2 adsorbed within zeolite Na-A. We use a block Lanczos procedure to generate the first several (9) rotational eigenstates of the molecule, which is modeled as a rigid, quantum rotor with an anisotropic polarizability and quadrupole moment. The rotor interacts with Na cations and O anions; interaction parameters are chosen semiempirically and the truncation of electrostatic fields is handled with a switching function. A Monte Carlo proceedure is used to sample a set of states based on the canonical distribution. Potential energy surfaces, favorable adsorbtion sites, and distributions of barriers to rotation are analyzed. Separation factors for ortho-parahydrogen are calculated; at low temperatures, these are controlled by the ease of rotational tunneling through barriers. (C) 1999 American Institute of Physics

Low angular momentum flow model of Sgr A* activity

Sgr A* is the closest massive black hole and can be observed with the highest angular resolution. Nevertheless, our current understanding of the accretion process in this source is very poor. The inflow is almost certainly of low radiative efficiency and it is accompanied by a strong outflow and the flow is strongly variable but the details of the dynamics are unknown. Even the amount of angular momentum in the flow is an open question. Here we argue that low angular momentum scenario is better suited to explain the flow variability. We present a new hybrid model which describes such a flow and consists of an outer spherically symmetric Bondi flow and an inner axially symmetric flow described through MHD simulations. The assumed angular momentum of the matter is low, i.e. the corresponding circularization radius in the equatorial plane of the flow is just above the innermost stable circular orbit in pseudo-Newtonian potential. We compare the radiation spectrum from such a flow to the broad band observational data for Sgr A*.Comment: Proceedings of the AHAR 2008 Conference: The Universe under the Microscope; Astrophysics at High Angular Resolution, Bad Honef