1
|
Allan L, Zuehlsdorff TJ. Taming the third order cumulant approximation to linear optical spectroscopy. J Chem Phys 2024; 160:074108. [PMID: 38380749 DOI: 10.1063/5.0182745] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Accepted: 01/23/2024] [Indexed: 02/22/2024] Open
Abstract
The second order cumulant method offers a promising pathway to predicting optical properties in condensed phase systems. It allows for the computation of linear absorption spectra from excitation energy fluctuations sampled along molecular dynamics (MD) trajectories, fully accounting for vibronic effects, direct solute-solvent interactions, and environmental polarization effects. However, the second order cumulant approximation only guarantees accurate line shapes for energy gap fluctuations obeying Gaussian statistics. A third order correction has recently been derived but often yields unphysical spectra or divergent line shapes for moderately non-Gaussian fluctuations due to the neglect of higher order terms in the cumulant expansion. In this work, we develop a corrected cumulant approach, where the collective effect of neglected higher order contributions is approximately accounted for through a dampening factor applied to the third order cumulant term. We show that this dampening factor can be expressed as a function of the skewness and kurtosis of energy gap fluctuations and can be parameterized from a large set of randomly sampled model Hamiltonians for which exact spectral line shapes are known. This approach is shown to systematically remove unphysical contributions in the form of negative absorbances from cumulant spectra in both model Hamiltonians and condensed phase systems sampled from MD and dramatically improves over the second order cumulant method in describing systems exhibiting Duschinsky mode mixing effects. We successfully apply the approach to the coumarin-153 dye in toluene, obtaining excellent agreement with experiment.
Collapse
Affiliation(s)
- Lucas Allan
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| |
Collapse
|
2
|
Mandal A, Taylor MA, Weight BM, Koessler ER, Li X, Huo P. Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics. Chem Rev 2023; 123:9786-9879. [PMID: 37552606 PMCID: PMC10450711 DOI: 10.1021/acs.chemrev.2c00855] [Citation(s) in RCA: 31] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Indexed: 08/10/2023]
Abstract
When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.
Collapse
Affiliation(s)
- Arkajit Mandal
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- Department
of Chemistry, Columbia University, New York, New York 10027, United States
| | - Michael A.D. Taylor
- The
Institute of Optics, Hajim School of Engineering, University of Rochester, Rochester, New York 14627, United States
| | - Braden M. Weight
- Department
of Physics and Astronomy, University of
Rochester, Rochester, New York 14627, United
States
| | - Eric R. Koessler
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Xinyang Li
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Pengfei Huo
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- The
Institute of Optics, Hajim School of Engineering, University of Rochester, Rochester, New York 14627, United States
| |
Collapse
|
3
|
Dawson W, Degomme A, Stella M, Nakajima T, Ratcliff LE, Genovese L. Density functional theory calculations of large systems: Interplay between fragments, observables, and computational complexity. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1574] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
| | | | - Martina Stella
- Department of Materials Imperial College London London UK
| | | | | | - Luigi Genovese
- Université Grenoble Alpes, INAC‐MEM, L_Sim Grenoble France
| |
Collapse
|
4
|
Shedge SV, Zuehlsdorff TJ, Khanna A, Conley S, Isborn CM. Explicit environmental and vibronic effects in simulations of linear and nonlinear optical spectroscopy. J Chem Phys 2021; 154:084116. [PMID: 33639769 DOI: 10.1063/5.0038196] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurately simulating the linear and nonlinear electronic spectra of condensed phase systems and accounting for all physical phenomena contributing to spectral line shapes presents a significant challenge. Vibronic transitions can be captured through a harmonic model generated from the normal modes of a chromophore, but it is challenging to also include the effects of specific chromophore-environment interactions within such a model. We work to overcome this limitation by combining approaches to account for both explicit environment interactions and vibronic couplings for simulating both linear and nonlinear optical spectra. We present and show results for three approaches of varying computational cost for combining ensemble sampling of chromophore-environment configurations with Franck-Condon line shapes for simulating linear spectra. We present two analogous approaches for nonlinear spectra. Simulated absorption spectra and two-dimensional electronic spectra (2DES) are presented for the Nile red chromophore in different solvent environments. Employing an average Franck-Condon or 2DES line shape appears to be a promising method for simulating linear and nonlinear spectroscopy for a chromophore in the condensed phase.
Collapse
Affiliation(s)
- Sapana V Shedge
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Ajay Khanna
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Stacey Conley
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
5
|
Zuehlsdorff TJ, Hong H, Shi L, Isborn CM. Nonlinear spectroscopy in the condensed phase: The role of Duschinsky rotations and third order cumulant contributions. J Chem Phys 2020; 153:044127. [PMID: 32752702 DOI: 10.1063/5.0013739] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
First-principles modeling of nonlinear optical spectra in the condensed phase is highly challenging because both environment and vibronic interactions can play a large role in determining spectral shapes and excited state dynamics. Here, we compute two dimensional electronic spectroscopy (2DES) signals based on a cumulant expansion of the energy gap fluctuation operator, with specific focus on analyzing mode mixing effects introduced by the Duschinsky rotation and the role of the third order term in the cumulant expansion for both model and realistic condensed phase systems. We show that for a harmonic model system, the third order cumulant correction captures effects introduced by a mismatch in curvatures of ground and excited state potential energy surfaces, as well as effects of mode mixing. We also demonstrate that 2DES signals can be accurately reconstructed from purely classical correlation functions using quantum correction factors. We then compute nonlinear optical spectra for the Nile red and methylene blue chromophores in solution, assessing the third order cumulant contribution for realistic systems. We show that the third order cumulant correction is strongly dependent on the treatment of the solvent environment, revealing the interplay between environmental polarization and the electronic-vibrational coupling.
Collapse
Affiliation(s)
- Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Hanbo Hong
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Liang Shi
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
6
|
Prentice JCA, Aarons J, Womack JC, Allen AEA, Andrinopoulos L, Anton L, Bell RA, Bhandari A, Bramley GA, Charlton RJ, Clements RJ, Cole DJ, Constantinescu G, Corsetti F, Dubois SMM, Duff KKB, Escartín JM, Greco A, Hill Q, Lee LP, Linscott E, O'Regan DD, Phipps MJS, Ratcliff LE, Serrano ÁR, Tait EW, Teobaldi G, Vitale V, Yeung N, Zuehlsdorff TJ, Dziedzic J, Haynes PD, Hine NDM, Mostofi AA, Payne MC, Skylaris CK. The ONETEP linear-scaling density functional theory program. J Chem Phys 2020; 152:174111. [PMID: 32384832 DOI: 10.1063/5.0004445] [Citation(s) in RCA: 51] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints. By taking advantage of localization, onetep is able to perform calculations including thousands of atoms with computational effort, which scales linearly with the number or atoms. The code has a large and diverse range of capabilities, explored in this paper, including different boundary conditions, various exchange-correlation functionals (with and without exact exchange), finite electronic temperature methods for metallic systems, methods for strongly correlated systems, molecular dynamics, vibrational calculations, time-dependent DFT, electronic transport, core loss spectroscopy, implicit solvation, quantum mechanical (QM)/molecular mechanical and QM-in-QM embedding, density of states calculations, distributed multipole analysis, and methods for partitioning charges and interactions between fragments. Calculations with onetep provide unique insights into large and complex systems that require an accurate atomic-level description, ranging from biomolecular to chemical, to materials, and to physical problems, as we show with a small selection of illustrative examples. onetep has always aimed to be at the cutting edge of method and software developments, and it serves as a platform for developing new methods of electronic structure simulation. We therefore conclude by describing some of the challenges and directions for its future developments and applications.
Collapse
Affiliation(s)
- Joseph C A Prentice
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Jolyon Aarons
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - James C Womack
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Alice E A Allen
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Lampros Andrinopoulos
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Lucian Anton
- UKAEA, Culham Science Centre, Abingdon OX14 3DB, United Kingdom
| | - Robert A Bell
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Arihant Bhandari
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Gabriel A Bramley
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Robert J Charlton
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Rebecca J Clements
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Daniel J Cole
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU, United Kingdom
| | - Gabriel Constantinescu
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Fabiano Corsetti
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Simon M-M Dubois
- Institute of Condensed Matter and Nanosciences, Université Catholique de Louvain, Louvain-la-Neuve, Belgium
| | - Kevin K B Duff
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - José María Escartín
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Andrea Greco
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Quintin Hill
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Louis P Lee
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Edward Linscott
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - David D O'Regan
- School of Physics, AMBER, and CRANN Institute, Trinity College Dublin, The University of Dublin, Dublin 2, Ireland
| | - Maximillian J S Phipps
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Laura E Ratcliff
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Álvaro Ruiz Serrano
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Edward W Tait
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Gilberto Teobaldi
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Valerio Vitale
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Nelson Yeung
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Peter D Haynes
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Nicholas D M Hine
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Arash A Mostofi
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Mike C Payne
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
7
|
Zuehlsdorff TJ, Hong H, Shi L, Isborn CM. Influence of Electronic Polarization on the Spectral Density. J Phys Chem B 2019; 124:531-543. [DOI: 10.1021/acs.jpcb.9b10250] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Tim J. Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Hanbo Hong
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Christine M. Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
8
|
Shedge SV, Zuehlsdorff TJ, Servis MJ, Clark AE, Isborn CM. Effect of Ions on the Optical Absorption Spectra of Aqueously Solvated Chromophores. J Phys Chem A 2019; 123:6175-6184. [DOI: 10.1021/acs.jpca.9b03163] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Sapana V. Shedge
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Tim J. Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Michael J. Servis
- Department of Chemistry and the Material Science and Engineering Program, Washington State University, Pullman, Washington 99164, United States
| | - Aurora E. Clark
- Department of Chemistry and the Material Science and Engineering Program, Washington State University, Pullman, Washington 99164, United States
- Voiland School of Chemical Engineering and Bioengineering, Washington State University, Pullman, Washington 99164, United States
| | - Christine M. Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
9
|
Carli M, Turelli M, Faccioli P. Microscopic calculation of absorption spectra of macromolecules: An analytic approach. J Chem Phys 2019; 150:144103. [PMID: 30981270 DOI: 10.1063/1.5084120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We develop a cross-disciplinary approach to analytically compute optical response functions of open macromolecular systems by exploiting the mathematical formalism of quantum field theory (QFT). Indeed, the entries of the density matrix for the electronic excitations interacting with their open dissipative environment are mapped into vacuum-to-vacuum Green's functions in a fictitious relativistic closed quantum system. We show that by re-summing appropriate self-energy diagrams in this dual QFT, it is possible to obtain analytic expressions for the response functions in Mukamel's theory. This yields physical insight into the structure and dynamics of vibronic resonances, since their frequency and width is related to fundamental physical constants and microscopic model parameters. For illustration, we apply this scheme to compute the linear absorption spectrum of the Fenna-Matthews-Olson light harvesting complex, comparing analytic calculations, numerical Monte Carlo simulations, and experimental data.
Collapse
Affiliation(s)
- Matteo Carli
- Physics Department of Trento University, Via Sommarive 14, Povo, Trento 38123, Italy
| | - Michele Turelli
- Physics Department of Trento University, Via Sommarive 14, Povo, Trento 38123, Italy
| | - Pietro Faccioli
- Physics Department of Trento University, Via Sommarive 14, Povo, Trento 38123, Italy
| |
Collapse
|
10
|
|
11
|
Zuehlsdorff TJ, Napoli JA, Milanese JM, Markland TE, Isborn CM. Unraveling electronic absorption spectra using nuclear quantum effects: Photoactive yellow protein and green fluorescent protein chromophores in water. J Chem Phys 2018; 149:024107. [PMID: 30007372 DOI: 10.1063/1.5025517] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Many physical phenomena must be accounted for to accurately model solution-phase optical spectral line shapes, from the sampling of chromophore-solvent configurations to the electronic-vibrational transitions leading to vibronic fine structure. Here we thoroughly explore the role of nuclear quantum effects, direct and indirect solvent effects, and vibronic effects in the computation of the optical spectrum of the aqueously solvated anionic chromophores of green fluorescent protein and photoactive yellow protein. By analyzing the chromophore and solvent configurations, the distributions of vertical excitation energies, the absorption spectra computed within the ensemble approach, and the absorption spectra computed within the ensemble plus zero-temperature Franck-Condon approach, we show how solvent, nuclear quantum effects, and vibronic transitions alter the optical absorption spectra. We find that including nuclear quantum effects in the sampling of chromophore-solvent configurations using ab initio path integral molecular dynamics simulations leads to improved spectral shapes through three mechanisms. The three mechanisms that lead to line shape broadening and a better description of the high-energy tail are softening of heavy atom bonds in the chromophore that couple to the optically bright state, widening the distribution of vertical excitation energies from more diverse solvation environments, and redistributing spectral weight from the 0-0 vibronic transition to higher energy vibronic transitions when computing the Franck-Condon spectrum in a frozen solvent pocket. The absorption spectra computed using the combined ensemble plus zero-temperature Franck-Condon approach yield significant improvements in spectral shape and width compared to the spectra computed with the ensemble approach. Using the combined approach with configurations sampled from path integral molecular dynamics trajectories presents a significant step forward in accurately modeling the absorption spectra of aqueously solvated chromophores.
Collapse
Affiliation(s)
- Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Joseph A Napoli
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Joel M Milanese
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
12
|
Khmelnitskiy A, Reinot T, Jankowiak R. Impact of Single-Point Mutations on the Excitonic Structure and Dynamics in a Fenna-Matthews-Olson Complex. J Phys Chem Lett 2018; 9:3378-3386. [PMID: 29863366 DOI: 10.1021/acs.jpclett.8b01396] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Hole burning (HB) spectroscopy and modeling studies reveal significant changes in the excitonic structure and dynamics in several mutants of the FMO trimer from the Chlorobaculum tepidum. The excited-state decay times ( T1) of the high-energy excitons are significantly modified when mutation occurs near bacteriochlorophyll (BChl) 1 (V152N mutant) or BChl 6 (W184F). Longer (averaged) T1 times of highest-energy excitons in V152N and W184F mutants suggest that site energies of BChls 1 and 6, believed to play an important role in receiving excitation from the baseplate BChls, likely play a critical role to ensure the femtosecond (fs) energy relaxation observed in wild-type FMO. HB spectroscopy reveals preferentially slower T1 times (about 1 ps on average) because fs times prohibit HB due to an extremely low HB quantum yield. Uncorrelated (incoherent) excitation energy transfer times between monomers, the composition of exciton states, and average, frequency-dependent, excited-state decay times ( T1) are discussed.
Collapse
|
13
|
Zuehlsdorff TJ, Isborn CM. Combining the ensemble and Franck-Condon approaches for calculating spectral shapes of molecules in solution. J Chem Phys 2018; 148:024110. [PMID: 29331131 DOI: 10.1063/1.5006043] [Citation(s) in RCA: 55] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The correct treatment of vibronic effects is vital for the modeling of absorption spectra of many solvated dyes. Vibronic spectra for small dyes in solution can be easily computed within the Franck-Condon approximation using an implicit solvent model. However, implicit solvent models neglect specific solute-solvent interactions on the electronic excited state. On the other hand, a straightforward way to account for solute-solvent interactions and temperature-dependent broadening is by computing vertical excitation energies obtained from an ensemble of solute-solvent conformations. Ensemble approaches usually do not account for vibronic transitions and thus often produce spectral shapes in poor agreement with experiment. We address these shortcomings by combining zero-temperature vibronic fine structure with vertical excitations computed for a room-temperature ensemble of solute-solvent configurations. In this combined approach, all temperature-dependent broadening is treated classically through the sampling of configurations and quantum mechanical vibronic contributions are included as a zero-temperature correction to each vertical transition. In our calculation of the vertical excitations, significant regions of the solvent environment are treated fully quantum mechanically to account for solute-solvent polarization and charge-transfer. For the Franck-Condon calculations, a small amount of frozen explicit solvent is considered in order to capture solvent effects on the vibronic shape function. We test the proposed method by comparing calculated and experimental absorption spectra of Nile red and the green fluorescent protein chromophore in polar and non-polar solvents. For systems with strong solute-solvent interactions, the combined approach yields significant improvements over the ensemble approach. For systems with weak to moderate solute-solvent interactions, both the high-energy vibronic tail and the width of the spectra are in excellent agreement with experiments.
Collapse
Affiliation(s)
- T J Zuehlsdorff
- School of Natural Sciences, University of California Merced, N. Lake Road, Merced, California 95344, USA
| | - C M Isborn
- School of Natural Sciences, University of California Merced, N. Lake Road, Merced, California 95344, USA
| |
Collapse
|
14
|
Khmelnitskiy A, Saer RG, Blankenship RE, Jankowiak R. Excitonic Energy Landscape of the Y16F Mutant of the Chlorobium tepidum Fenna-Matthews-Olson (FMO) Complex: High Resolution Spectroscopic and Modeling Studies. J Phys Chem B 2018; 122:3734-3743. [PMID: 29554425 DOI: 10.1021/acs.jpcb.7b11763] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
We report high-resolution (low-temperature) absorption, emission, and nonresonant/resonant hole-burned (HB) spectra and results of excitonic calculations using a non-Markovian reduced density matrix theory (with an improved algorithm for parameter optimization in heterogeneous samples) obtained for the Y16F mutant of the Fenna-Matthews-Olson (FMO) trimer from the green sulfur bacterium Chlorobium tepidum. We show that the Y16F mutant is a mixture of FMO complexes with three independent low-energy traps (located near 817, 821, and 826 nm), in agreement with measured composite emission and HB spectra. Two of these traps belong to mutated FMO subpopulations characterized by significantly modified low-energy excitonic states. Hamiltonians for the two major subpopulations (Sub821 and Sub817) provide new insight into extensive changes induced by the single-point mutation in the vicinity of BChl 3 (where tyrosine Y16 was replaced with phenylalanine F16). The average decay time(s) from the higher exciton state(s) in the Y16F mutant depends on frequency and occurs on a picosecond time scale.
Collapse
Affiliation(s)
| | - Rafael G Saer
- Departments of Biology and Chemistry , Washington University in St. Louis , Saint Louis , Missouri 63130 , United States
| | - Robert E Blankenship
- Departments of Biology and Chemistry , Washington University in St. Louis , Saint Louis , Missouri 63130 , United States
| | | |
Collapse
|
15
|
Khmelnitskiy A, Kell A, Reinot T, Saer RG, Blankenship RE, Jankowiak R. Energy landscape of the intact and destabilized FMO antennas from C. tepidum and the L122Q mutant: Low temperature spectroscopy and modeling study. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2018; 1859:165-173. [DOI: 10.1016/j.bbabio.2017.11.008] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/24/2017] [Revised: 11/23/2017] [Accepted: 11/27/2017] [Indexed: 12/21/2022]
|
16
|
Baker LA, Habershon S. Photosynthesis, pigment-protein complexes and electronic energy transport: simple models for complicated processes. Sci Prog 2017; 100:313-330. [PMID: 28779762 PMCID: PMC10365183 DOI: 10.3184/003685017x14967574639964] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
In this review, we discuss our recent work on modelling biological pigment-protein complexes, such as the Fenna-Matthews-Olson complex and light-harvesting complex-II, to explain their electronic energy transport properties. In particular, we highlight how a network-based analysis approach, where the light-absorbing pigments are treated as a network of interconnected nodes, can provide a qualitative picture of quantum dynamic energy transport. With this in mind, we demonstrate how other properties such as robustness to environmental changes can be assessed in a simple and computationally tractable manner. Such analyses could prove useful for the design of artificial energy transport networks such as those which might find application in solar cells.
Collapse
|
17
|
Lee MK, Bravaya KB, Coker DF. First-Principles Models for Biological Light-Harvesting: Phycobiliprotein Complexes from Cryptophyte Algae. J Am Chem Soc 2017; 139:7803-7814. [DOI: 10.1021/jacs.7b01780] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Mi Kyung Lee
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Ksenia B. Bravaya
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - David F. Coker
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
18
|
Fokas AS, Cole DJ, Hine NDM, Wells SA, Payne MC, Chin AW. Evidence of Correlated Static Disorder in the Fenna-Matthews-Olson Complex. J Phys Chem Lett 2017; 8:2350-2356. [PMID: 28485971 DOI: 10.1021/acs.jpclett.7b00669] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Observation of excitonic quantum beats in photosynthetic antennae has prompted wide debate regarding the function of excitonic coherence in pigment-protein complexes. Much of this work focuses on the interactions of excitons with the femto-to-picosecond dynamical fluctuations of their environment. However, in experiments these effects can be masked by static disorder of the excited-state energies across ensembles, whose microscopic origins are challenging to predict. Here the excited-state properties of ∼2000 atom clusters of the Fenna-Matthews-Olson complex are simulated using a unique combination of linear-scaling density functional theory and constrained geometric dynamics. While slow, large amplitude protein motion leads to large variations in the Qy transitions of two pigments, we identify pigment-protein correlations that greatly reduce variations in the energy gap across the ensemble, which is consistent with experimental observations of suppressed inhomogeneous dephasing of quantum beats.
Collapse
Affiliation(s)
- Alexander S Fokas
- TCM Group, Cavendish Laboratory , 19 J J Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Daniel J Cole
- TCM Group, Cavendish Laboratory , 19 J J Thomson Avenue, Cambridge CB3 0HE, United Kingdom
- School of Chemistry, Newcastle University , Newcastle upon Tyne NE1 7RU, United Kingdom
| | - Nicholas D M Hine
- Department of Physics, University of Warwick , Coventry CV4 7AL, United Kingdom
| | - Stephen A Wells
- Department of Chemistry, University of Bath , Claverton Down BA2 7AY, United Kingdom
| | - Michael C Payne
- TCM Group, Cavendish Laboratory , 19 J J Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Alex W Chin
- TCM Group, Cavendish Laboratory , 19 J J Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
19
|
Baker LA, Habershon S. Photosynthetic pigment-protein complexes as highly connected networks: implications for robust energy transport. Proc Math Phys Eng Sci 2017; 473:20170112. [PMID: 28588417 PMCID: PMC5454362 DOI: 10.1098/rspa.2017.0112] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2017] [Accepted: 05/04/2017] [Indexed: 02/01/2023] Open
Abstract
Photosynthetic pigment-protein complexes (PPCs) are a vital component of the light-harvesting machinery of all plants and photosynthesizing bacteria, enabling efficient transport of the energy of absorbed light towards the reaction centre, where chemical energy storage is initiated. PPCs comprise a set of chromophore molecules, typically bacteriochlorophyll species, held in a well-defined arrangement by a protein scaffold; this relatively rigid distribution leads to a viewpoint in which the chromophore subsystem is treated as a network, where chromophores represent vertices and inter-chromophore electronic couplings represent edges. This graph-based view can then be used as a framework within which to interrogate the role of structural and electronic organization in PPCs. Here, we use this network-based viewpoint to compare excitation energy transfer (EET) dynamics in the light-harvesting complex II (LHC-II) system commonly found in higher plants and the Fenna-Matthews-Olson (FMO) complex found in green sulfur bacteria. The results of our simple network-based investigations clearly demonstrate the role of network connectivity and multiple EET pathways on the efficient and robust EET dynamics in these PPCs, and highlight a role for such considerations in the development of new artificial light-harvesting systems.
Collapse
Affiliation(s)
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, UK
| |
Collapse
|
20
|
Jassas M, Reinot T, Kell A, Jankowiak R. Toward an Understanding of the Excitonic Structure of the CP47 Antenna Protein Complex of Photosystem II Revealed via Circularly Polarized Luminescence. J Phys Chem B 2017; 121:4364-4378. [PMID: 28394609 DOI: 10.1021/acs.jpcb.7b00362] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Identification of the lowest energy pigments in the photosynthetic CP47 antenna protein complex of Photosystem II (PSII) is essential for understanding its excitonic structure, as well as excitation energy pathways in the PSII core complex. Unfortunately, there is no consensus concerning the nature of the low-energy state(s), nor chlorophyll (Chl) site energies in this important photosynthetic antenna. Although we raised concerns regarding the estimations of Chl site energies obtained from modeling studies of various types of CP47 optical spectra [Reinot, T; et al., Anal. Chem. Insights 2016, 11, 35-48] recent new assignments imposed by the shape of the circularly polarized luminescence (CPL) spectrum [Hall, J.; et al., Biochim. Biophys. Acta 2016, 1857, 1580-1593] necessitate our comments. We demonstrate that other combinations of low-energy Chls provide equally good or improved simultaneous fits of various optical spectra (absorption, emission, CPL, circular dichroism, and nonresonant hole-burned spectra), but more importantly, we expose the heterogeneous nature of the recently studied complexes and argue that the published composite nature of the CPL (contributed to by CPL685, CPL691, and CPL695) does not represent an intact CP47 protein. A positive CPL695 is extracted for the intact protein, which, when simultaneously fitted with multiple other optical spectra, provides new information on the excitonic structure of intact and destabilized CP47 complexes and their lowest energy state(s).
Collapse
Affiliation(s)
- Mahboobe Jassas
- Department of Chemistry and ‡Department of Physics, Kansas State University , Manhattan, Kansas 66506, United States
| | - Tonu Reinot
- Department of Chemistry and ‡Department of Physics, Kansas State University , Manhattan, Kansas 66506, United States
| | - Adam Kell
- Department of Chemistry and ‡Department of Physics, Kansas State University , Manhattan, Kansas 66506, United States
| | - Ryszard Jankowiak
- Department of Chemistry and ‡Department of Physics, Kansas State University , Manhattan, Kansas 66506, United States
| |
Collapse
|
21
|
Sisto A, Stross C, van der Kamp MW, O’Connor M, McIntosh-Smith S, Johnson GT, Hohenstein EG, Manby FR, Glowacki DR, Martinez TJ. Atomistic non-adiabatic dynamics of the LH2 complex with a GPU-accelerated ab initio exciton model. Phys Chem Chem Phys 2017; 19:14924-14936. [DOI: 10.1039/c7cp00492c] [Citation(s) in RCA: 50] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
We present GPU-accelerated ab initio molecular dynamics simulations of nonadiabatic dynamics in the LH2 complex in full atomistic detail.
Collapse
Affiliation(s)
- Aaron Sisto
- PULSE Institute and Department of Chemistry
- Stanford University
- Stanford
- USA
- SLAC National Accelerator Laboratory
| | - Clem Stross
- School of Chemistry
- University of Bristol
- Bristol
- UK
| | | | - Michael O’Connor
- School of Chemistry
- University of Bristol
- Bristol
- UK
- Department of Computer Science
| | | | - Graham T. Johnson
- California Institute for Quantitative Biosciences (QB3)
- University of California
- San Francisco
- USA
- Department of Bioengineering and Therapeutic Sciences
| | | | | | - David R. Glowacki
- School of Chemistry
- University of Bristol
- Bristol
- UK
- Department of Computer Science
| | - Todd J. Martinez
- PULSE Institute and Department of Chemistry
- Stanford University
- Stanford
- USA
- SLAC National Accelerator Laboratory
| |
Collapse
|
22
|
Morgan SE, Cole DJ, Chin AW. Nonlinear network model analysis of vibrational energy transfer and localisation in the Fenna-Matthews-Olson complex. Sci Rep 2016; 6:36703. [PMID: 27827409 PMCID: PMC5101523 DOI: 10.1038/srep36703] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2016] [Accepted: 10/19/2016] [Indexed: 11/10/2022] Open
Abstract
Collective protein modes are expected to be important for facilitating energy transfer in the Fenna-Matthews-Olson (FMO) complex of photosynthetic green sulphur bacteria, however to date little work has focussed on the microscopic details of these vibrations. The nonlinear network model (NNM) provides a computationally inexpensive approach to studying vibrational modes at the microscopic level in large protein structures, whilst incorporating anharmonicity in the inter-residue interactions which can influence protein dynamics. We apply the NNM to the entire trimeric FMO complex and find evidence for the existence of nonlinear discrete breather modes. These modes tend to transfer energy to the highly connected core pigments, potentially opening up alternative excitation energy transfer routes through their influence on pigment properties. Incorporating localised modes based on these discrete breathers in the optical spectra calculations for FMO using ab initio site energies and excitonic couplings can substantially improve their agreement with experimental results.
Collapse
Affiliation(s)
- Sarah E Morgan
- Theory of Condensed Matter Group, Physics Department, University of Cambridge, CB3 0HE, United Kingdom
| | - Daniel J Cole
- Theory of Condensed Matter Group, Physics Department, University of Cambridge, CB3 0HE, United Kingdom.,School of Chemistry, Newcastle University, Newcastle upon Tyne, NE1 7RU, United Kingdom
| | - Alex W Chin
- Theory of Condensed Matter Group, Physics Department, University of Cambridge, CB3 0HE, United Kingdom
| |
Collapse
|
23
|
Cole DJ, Hine NDM. Applications of large-scale density functional theory in biology. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:393001. [PMID: 27494095 DOI: 10.1088/0953-8984/28/39/393001] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Density functional theory (DFT) has become a routine tool for the computation of electronic structure in the physics, materials and chemistry fields. Yet the application of traditional DFT to problems in the biological sciences is hindered, to a large extent, by the unfavourable scaling of the computational effort with system size. Here, we review some of the major software and functionality advances that enable insightful electronic structure calculations to be performed on systems comprising many thousands of atoms. We describe some of the early applications of large-scale DFT to the computation of the electronic properties and structure of biomolecules, as well as to paradigmatic problems in enzymology, metalloproteins, photosynthesis and computer-aided drug design. With this review, we hope to demonstrate that first principles modelling of biological structure-function relationships are approaching a reality.
Collapse
Affiliation(s)
- Daniel J Cole
- Theory of Condensed Matter group, Cavendish Laboratory, 19 JJ Thomson Ave, Cambridge CB3 0HE, UK. School of Chemistry, Newcastle University, Newcastle upon Tyne NE1 7RU, UK
| | | |
Collapse
|
24
|
Dziedzic J, Mao Y, Shao Y, Ponder J, Head-Gordon T, Head-Gordon M, Skylaris CK. TINKTEP: A fully self-consistent, mutually polarizable QM/MM approach based on the AMOEBA force field. J Chem Phys 2016; 145:124106. [DOI: 10.1063/1.4962909] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
- Faculty of Applied Physics and Mathematics, Gdańsk University of Technology, Gdańsk, Poland
| | - Yuezhi Mao
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yihan Shao
- Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Jay Ponder
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, USA
| | - Teresa Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA
- Department of Bioengineering, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
25
|
Reinot T, Chen J, Kell A, Jassas M, Robben KC, Zazubovich V, Jankowiak R. On the Conflicting Estimations of Pigment Site Energies in Photosynthetic Complexes: A Case Study of the CP47 Complex. ANALYTICAL CHEMISTRY INSIGHTS 2016; 11:35-48. [PMID: 27279733 PMCID: PMC4892206 DOI: 10.4137/aci.s32151] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/15/2016] [Revised: 04/10/2016] [Accepted: 04/26/2016] [Indexed: 01/14/2023]
Abstract
We focus on problems with elucidation of site energies
(E0n) for photosynthetic complexes (PSCs) in order to raise some genuine concern regarding the conflicting estimations propagating in the literature. As an example, we provide a stern assessment of the site energies extracted from fits to optical spectra of the widely studied CP47 antenna complex of photosystem II from spinach, though many general comments apply to other PSCs as well. Correct values of
E0n for chlorophyll (Chl) a in CP47 are essential for understanding its excitonic structure, population dynamics, and excitation energy pathway(s). To demonstrate this, we present a case study where simultaneous fits of multiple spectra (absorption, emission, circular dichroism, and nonresonant hole-burned spectra) show that several sets of parameters can fit the spectra very well. Importantly, we show that variable emission maxima (690–695 nm) and sample-dependent bleaching in nonresonant hole-burning spectra reported in literature could be explained, assuming that many previously studied CP47 samples were a mixture of intact and destabilized proteins. It appears that the destabilized subpopulation of CP47 complexes could feature a weakened hydrogen bond between the 131-keto group of Chl29 and the PsbH protein subunit, though other possibilities cannot be entirely excluded, as discussed in this work. Possible implications of our findings are briefly discussed.
Collapse
Affiliation(s)
- Tonu Reinot
- Department of Chemistry, Kansas State University, Manhattan, KS, USA
| | - Jinhai Chen
- Department of Chemistry, Kansas State University, Manhattan, KS, USA
| | - Adam Kell
- Department of Chemistry, Kansas State University, Manhattan, KS, USA
| | - Mahboobe Jassas
- Department of Chemistry, Kansas State University, Manhattan, KS, USA
| | - Kevin C Robben
- Department of Chemistry, Kansas State University, Manhattan, KS, USA
| | | | - Ryszard Jankowiak
- Department of Chemistry, Kansas State University, Manhattan, KS, USA.; Department of Physics, Kansas State University, Manhattan, KS, USA
| |
Collapse
|
26
|
Lee MK, Huo P, Coker DF. Semiclassical Path Integral Dynamics: Photosynthetic Energy Transfer with Realistic Environment Interactions. Annu Rev Phys Chem 2016; 67:639-68. [DOI: 10.1146/annurev-physchem-040215-112252] [Citation(s) in RCA: 96] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Mi Kyung Lee
- Department of Chemistry, Boston University, Boston, Massachusetts 02215;
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, Rochester, New York 14627;
| | - David F. Coker
- Department of Chemistry, Boston University, Boston, Massachusetts 02215;
| |
Collapse
|
27
|
Cole DJ, Vilseck JZ, Tirado-Rives J, Payne MC, Jorgensen WL. Biomolecular Force Field Parameterization via Atoms-in-Molecule Electron Density Partitioning. J Chem Theory Comput 2016; 12:2312-23. [PMID: 27057643 PMCID: PMC4864407 DOI: 10.1021/acs.jctc.6b00027] [Citation(s) in RCA: 91] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
![]()
Molecular mechanics
force fields, which are commonly used in biomolecular
modeling and computer-aided drug design, typically treat nonbonded
interactions using a limited library of empirical parameters that
are developed for small molecules. This approach does not account
for polarization in larger molecules or proteins, and the parametrization
process is labor-intensive. Using linear-scaling density functional
theory and atoms-in-molecule electron density partitioning, environment-specific
charges and Lennard-Jones parameters are derived directly from quantum
mechanical calculations for use in biomolecular modeling of organic
and biomolecular systems. The proposed methods significantly reduce
the number of empirical parameters needed to construct molecular mechanics
force fields, naturally include polarization effects in charge and
Lennard-Jones parameters, and scale well to systems comprised of thousands
of atoms, including proteins. The feasibility and benefits of this
approach are demonstrated by computing free energies of hydration,
properties of pure liquids, and the relative binding free energies
of indole and benzofuran to the L99A mutant of T4 lysozyme.
Collapse
Affiliation(s)
- Daniel J Cole
- Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States.,TCM Group, Cavendish Laboratory, 19 JJ Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Jonah Z Vilseck
- Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States
| | - Julian Tirado-Rives
- Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States
| | - Mike C Payne
- TCM Group, Cavendish Laboratory, 19 JJ Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - William L Jorgensen
- Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States
| |
Collapse
|
28
|
Zuehlsdorff TJ, Haynes PD, Hanke F, Payne MC, Hine NDM. Solvent Effects on Electronic Excitations of an Organic Chromophore. J Chem Theory Comput 2016; 12:1853-61. [DOI: 10.1021/acs.jctc.5b01014] [Citation(s) in RCA: 52] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- T. J. Zuehlsdorff
- Cavendish
Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, U.K
| | | | - F. Hanke
- Dassault Systèmes
BIOVIA, 334 Science Park, Cambridge CB4 0WN, U.K
| | - M. C. Payne
- Cavendish
Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, U.K
| | - N. D. M. Hine
- Department
of Physics, University of Warwick, Coventry CV4 7AL, U.K
| |
Collapse
|
29
|
Abstract
The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.
Collapse
Affiliation(s)
- Carles Curutchet
- Departament de Fisicoquímica, Facultat de Farmàcia, Universitat de Barcelona , Av. Joan XXIII s/n, 08028 Barcelona, Spain
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa , via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
30
|
Zuehlsdorff TJ, Hine NDM, Payne MC, Haynes PD. Linear-scaling time-dependent density-functional theory beyond the Tamm-Dancoff approximation: Obtaining efficiency and accuracy with in situ optimised local orbitals. J Chem Phys 2015; 143:204107. [PMID: 26627950 DOI: 10.1063/1.4936280] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present a solution of the full time-dependent density-functional theory (TDDFT) eigenvalue equation in the linear response formalism exhibiting a linear-scaling computational complexity with system size, without relying on the simplifying Tamm-Dancoff approximation (TDA). The implementation relies on representing the occupied and unoccupied subspaces with two different sets of in situ optimised localised functions, yielding a very compact and efficient representation of the transition density matrix of the excitation with the accuracy associated with a systematic basis set. The TDDFT eigenvalue equation is solved using a preconditioned conjugate gradient algorithm that is very memory-efficient. The algorithm is validated on a small test molecule and a good agreement with results obtained from standard quantum chemistry packages is found, with the preconditioner yielding a significant improvement in convergence rates. The method developed in this work is then used to reproduce experimental results of the absorption spectrum of bacteriochlorophyll in an organic solvent, where it is demonstrated that the TDA fails to reproduce the main features of the low energy spectrum, while the full TDDFT equation yields results in good qualitative agreement with experimental data. Furthermore, the need for explicitly including parts of the solvent into the TDDFT calculations is highlighted, making the treatment of large system sizes necessary that are well within reach of the capabilities of the algorithm introduced here. Finally, the linear-scaling properties of the algorithm are demonstrated by computing the lowest excitation energy of bacteriochlorophyll in solution. The largest systems considered in this work are of the same order of magnitude as a variety of widely studied pigment-protein complexes, opening up the possibility of studying their properties without having to resort to any semiclassical approximations to parts of the protein environment.
Collapse
Affiliation(s)
- T J Zuehlsdorff
- Cavendish Laboratory, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - N D M Hine
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - M C Payne
- Cavendish Laboratory, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - P D Haynes
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| |
Collapse
|
31
|
Steinmann C, Kongsted J. Electronic Energy Transfer in Polarizable Heterogeneous Environments: A Systematic Investigation of Different Quantum Chemical Approaches. J Chem Theory Comput 2015; 11:4283-93. [DOI: 10.1021/acs.jctc.5b00470] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Casper Steinmann
- Department of Physics, Chemistry,
and Pharmacy, University of Southern Denmark, DK-5230 Odense
M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry,
and Pharmacy, University of Southern Denmark, DK-5230 Odense
M, Denmark
| |
Collapse
|
32
|
Jurinovich S, Viani L, Curutchet C, Mennucci B. Limits and potentials of quantum chemical methods in modelling photosynthetic antennae. Phys Chem Chem Phys 2015; 17:30783-92. [DOI: 10.1039/c5cp00986c] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
A critical overview of quantum chemical approaches to simulate the light-harvesting process in photosynthetic antennae is presented together with a perspective on the developments that need to be introduced to reach a quantitative predictive power.
Collapse
Affiliation(s)
- Sandro Jurinovich
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
| | - Lucas Viani
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
- Institute for Fluid Dynamics
| | - Carles Curutchet
- Departament de Fisicoquímica
- Facultat de Farmàcia
- Universitat de Barcelona
- 08028 Barcelona
- Spain
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
| |
Collapse
|
33
|
Li JH, Zuehlsdorff TJ, Payne MC, Hine NDM. Identifying and tracing potential energy surfaces of electronic excitations with specific character via their transition origins: application to oxirane. Phys Chem Chem Phys 2015; 17:12065-79. [DOI: 10.1039/c5cp01018g] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We show that the transition origins of electronic excitations identified by quantified natural transition orbital (QNTO) analysis can be employed to connect potential energy surfaces (PESs) according to their character across a wide range of molecular geometries.
Collapse
Affiliation(s)
- Jian-Hao Li
- TCM Group
- Cavendish Laboratory
- Cambridge CB3 0HE
- UK
| | | | - M. C. Payne
- TCM Group
- Cavendish Laboratory
- Cambridge CB3 0HE
- UK
| | | |
Collapse
|
34
|
Andreussi O, Caprasecca S, Cupellini L, Guarnetti-Prandi I, Guido CA, Jurinovich S, Viani L, Mennucci B. Plasmon Enhanced Light Harvesting: Multiscale Modeling of the FMO Protein Coupled with Gold Nanoparticles. J Phys Chem A 2014; 119:5197-206. [DOI: 10.1021/jp509870b] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Oliviero Andreussi
- Department of Chemistry, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Stefano Caprasecca
- Department of Chemistry, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Lorenzo Cupellini
- Department of Chemistry, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | | | - Ciro A. Guido
- Department of Chemistry, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Sandro Jurinovich
- Department of Chemistry, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Lucas Viani
- Department of Chemistry, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Department of Chemistry, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| |
Collapse
|
35
|
Fokas AS, Cole DJ, Chin AW. Constrained geometric dynamics of the Fenna-Matthews-Olson complex: the role of correlated motion in reducing uncertainty in excitation energy transfer. PHOTOSYNTHESIS RESEARCH 2014; 122:275-292. [PMID: 25034014 DOI: 10.1007/s11120-014-0027-3] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/21/2014] [Accepted: 06/30/2014] [Indexed: 06/03/2023]
Abstract
The trimeric Fenna-Mathews-Olson (FMO) complex of green sulphur bacteria is a well-studied example of a photosynthetic pigment-protein complex, in which the electronic properties of the pigments are modified by the protein environment to promote efficient excitonic energy transfer from antenna complexes to the reaction centres. By a range of simulation methods, many of the electronic properties of the FMO complex can be extracted from knowledge of the static crystal structure. However, the recent observation and analysis of long-lasting quantum dynamics in the FMO complex point to protein dynamics as a key factor in protecting and generating quantum coherence under laboratory conditions. While fast inter- and intra-molecular vibrations have been investigated extensively, the slow, conformational dynamics which effectively determine the optical inhomogeneous broadening of experimental ensembles has received less attention. The following study employs constrained geometric dynamics to study the flexibility in the protein network by efficiently generating the accessible conformational states from the published crystal structure. Statistical and principle component analyses reveal highly correlated low frequency motions between functionally relevant elements, including strong correlations between pigments that are excitonically coupled. Our analysis reveals a hierarchy of structural interactions which enforce these correlated motions, from the level of monomer-monomer interfaces right down to the α-helices, β-sheets and pigments. In addition to inducing strong spatial correlations across the conformational ensemble, we find that the overall rigidity of the FMO complex is exceptionally high. We suggest that these observations support the idea of highly correlated inhomogeneous disorder of the electronic excited states, which is further supported by the remarkably low variance (typically <5%) of the excitonic couplings of the conformational ensemble.
Collapse
Affiliation(s)
- Alexander S Fokas
- The Theory of Condensed Matter Group Group, Cavendish Laboratory, 19 JJ Thomson Avenue, Cambridge, CB3 0HE, UK,
| | | | | |
Collapse
|
36
|
Lever G, Cole DJ, Lonsdale R, Ranaghan KE, Wales DJ, Mulholland AJ, Skylaris CK, Payne MC. Large-Scale Density Functional Theory Transition State Searching in Enzymes. J Phys Chem Lett 2014; 5:3614-9. [PMID: 26278727 DOI: 10.1021/jz5018703] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Linear-scaling quantum mechanical density functional theory calculations have been applied to study the rearrangement of chorismate to prephenate in large-scale models of the Bacillus subtilis chorismate mutase enzyme. By treating up to 2000 atoms at a consistent quantum mechanical level of theory, we obtain an unbiased, almost parameter-free description of the transition state geometry and energetics. The activation energy barrier is calculated to be lowered by 10.5 kcal mol(-1) in the enzyme, compared with the equivalent reaction in water, which is in good agreement with experiment. Natural bond orbital analysis identifies a number of active site residues that are important for transition state stabilization in chorismate mutase. This benchmark study demonstrates that linear-scaling density functional theory techniques are capable of simulating entire enzymes at the ab initio quantum mechanical level of accuracy.
Collapse
Affiliation(s)
- Greg Lever
- †Theory of Condensed Matter Group, Cavendish Laboratory, 19 JJ Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Daniel J Cole
- †Theory of Condensed Matter Group, Cavendish Laboratory, 19 JJ Thomson Avenue, Cambridge CB3 0HE, United Kingdom
- ‡Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - Richard Lonsdale
- ¶Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Kara E Ranaghan
- ¶Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - David J Wales
- §University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Adrian J Mulholland
- ¶Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Chris-Kriton Skylaris
- ∥School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Mike C Payne
- †Theory of Condensed Matter Group, Cavendish Laboratory, 19 JJ Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
37
|
Jurinovich S, Curutchet C, Mennucci B. The Fenna-Matthews-Olson protein revisited: a fully polarizable (TD)DFT/MM description. Chemphyschem 2014; 15:3194-204. [PMID: 25080315 DOI: 10.1002/cphc.201402244] [Citation(s) in RCA: 61] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2014] [Indexed: 01/21/2023]
Abstract
We report a combined molecular dynamics and quantum mechanics (QM)/molecular mechanics (MM) analysis of the excitonic properties of the Fenna-Matthews-Olson (FMO) protein by using a polarizable MM model combined with a time-dependent density functional theory description. Overall, our results indicate that structural fluctuations, electrostatic interactions, and short-range quantum effects can significantly modulate the model Hamiltonian parameters (site energies and couplings). We find that the specific interactions with the axial ligand and the hydrogen-bonded residues are responsible for the energy ladder, with their effects being mainly due to electrostatic interactions, but with short-range quantum contributions that are not negligible. In addition, a striking modulation of the screening effects experienced by the BChl pairs, due to the heterogeneous polarizability of the FMO and solvent environment, is observed. Finally, we find that the exciton model gives a reliable description of the delocalized excited states in the complex.
Collapse
Affiliation(s)
- Sandro Jurinovich
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via Risorgimento 35, I-56126 Pisa (Italy)
| | | | | |
Collapse
|