1
|
Zhang Q, Shao X, Li W, Mi W, Pavanello M, Akimov AV. Nonadiabatic molecular dynamics with subsystem density functional theory: application to crystalline pentacene. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2024; 36:385901. [PMID: 38866023 DOI: 10.1088/1361-648x/ad577d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2024] [Accepted: 06/12/2024] [Indexed: 06/14/2024]
Abstract
In this work, we report the development and assessment of the nonadiabatic molecular dynamics approach with the electronic structure calculations based on the linearly scaling subsystem density functional method. The approach is implemented in an open-source embedded Quantum Espresso/Libra software specially designed for nonadiabatic dynamics simulations in extended systems. As proof of the applicability of this method to large condensed-matter systems, we examine the dynamics of nonradiative relaxation of excess excitation energy in pentacene crystals with the simulation supercells containing more than 600 atoms. We find that increased structural disorder observed in larger supercell models induces larger nonadiabatic couplings of electronic states and accelerates the relaxation dynamics of excited states. We conduct a comparative analysis of several quantum-classical trajectory surface hopping schemes, including two new methods proposed in this work (revised decoherence-induced surface hopping and instantaneous decoherence at frustrated hops). Most of the tested schemes suggest fast energy relaxation occurring with the timescales in the 0.7-2.0 ps range, but they significantly overestimate the ground state recovery rates. Only the modified simplified decay of mixing approach yields a notably slower relaxation timescales of 8-14 ps, with a significantly inhibited ground state recovery.
Collapse
Affiliation(s)
- Qingxin Zhang
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, NY 14260, United States of America
| | - Xuecheng Shao
- Department of Physics, Rutgers University, The State University of New Jersey, Newark, NJ 07102, United States of America
| | - Wei Li
- School of Chemistry and Materials Science, Hunan Agricultural University, Changsha 410128, People's Republic of China
| | - Wenhui Mi
- Key Laboratory of Material Simulation Methods & Software of Ministry of Education, College of Physics, Jilin University, Changchun 130012, People's Republic of China
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, People's Republic of China
| | - Michele Pavanello
- Department of Physics, Rutgers University, The State University of New Jersey, Newark, NJ 07102, United States of America
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, NY 14260, United States of America
| |
Collapse
|
2
|
Käfer S, Niemeyer N, Tölle J, Neugebauer J. Triplet Excitation-Energy Transfer Couplings from Subsystem Time-Dependent Density-Functional Theory. J Chem Theory Comput 2024; 20:2475-2490. [PMID: 38450637 DOI: 10.1021/acs.jctc.3c01365] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/08/2024]
Abstract
We present an implementation of triplet excitation-energy transfer (TEET) couplings based on subsystem-based time-dependent density-functional theory (sTDDFT). TEET couplings are systematically investigated by comparing "exact" and approximate variants of sTDDFT. We demonstrate that, while sTDDFT utilizing explicit approximate non-additive kinetic energy (NAKE) density functionals is well-suited for describing singlet EET processes, it is inadequate for characterizing TEET. However, we show that projection-based embedding (PbE)-based sTDDFT addresses the challenges faced by NAKE-sTDDFT and emerges as a promising method for accurately describing electronic couplings in TEET processes. We also introduce the mixed PbE-/NAKE-embedding procedure to investigate the TEET effects in solvated pairs of chromophores. This approach offers a good balance between accuracy and efficiency, enabling comprehensive studies of TEET processes in complex environments.
Collapse
Affiliation(s)
- Sabine Käfer
- Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, University of Münster, Corrensstraße 36, Münster 48149, Germany
| | - Niklas Niemeyer
- Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, University of Münster, Corrensstraße 36, Münster 48149, Germany
| | - Johannes Tölle
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Johannes Neugebauer
- Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, University of Münster, Corrensstraße 36, Münster 48149, Germany
| |
Collapse
|
3
|
Cherepanov DA, Milanovsky GE, Neverov KV, Obukhov YN, Maleeva YV, Aybush AV, Kritsky MS, Nadtochenko VA. Exciton interactions of chlorophyll tetramer in water-soluble chlorophyll-binding protein BoWSCP. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 309:123847. [PMID: 38217986 DOI: 10.1016/j.saa.2024.123847] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Revised: 01/02/2024] [Accepted: 01/03/2024] [Indexed: 01/15/2024]
Abstract
The exciton interaction of four chlorophyll a (Chl a) molecules in a symmetrical tetrameric complex of the water-soluble chlorophyll-binding protein BoWSCP was analyzed in the pH range of 3-11. Exciton splitting ΔE = 232 ± 2 cm-1 of the Qy band of Chl a into two subcomponents with relative intensities of 78.1 ± 0.7 % and 21.9 ± 0.7 % was determined by a joint decomposition of the absorption and circular dichroism spectra into Gaussian functions. The exciton coupling parameters were calculated based on the BoWSCP atomic structure in three approximations: the point dipole model, the distributed atomic monopoles, and direct ab initio calculations in the TDDFT/PCM approximation. The Coulomb interactions of monomers were calculated within the continuum model using three values of optical permittivity. The models based on the properties of free Chl a in solution suffer from significant errors both in estimating the absolute value of the exciton interaction and in the relative intensity of exciton transitions. Calculations within the TDDFT/PCM approximation reproduce the experimentally determined parameters of the exciton splitting and the relative intensities of the exciton bands. The following factors of pigment-protein and pigment-pigment interactions were examined: deviation of the macrocycle geometry from the planar conformation of free Chl; the formation of hydrogen bonds between the macrocycle and water molecules; the overlap of wave functions of monomers at close distances. The most significant factor is the geometrical deformation of the porphyrin macrocycle, which leads to an increase in the dipole moment of Chl monomer from 5.5 to 6.9 D and to a rotation of the dipole moment by 15° towards the cyclopentane ring. The contributions of resonant charge-transfer states to the wave functions of the Chl dimer were determined and the transition dipole moments of the symmetric and antisymmetric charge-transfer states were estimated.
Collapse
Affiliation(s)
- D A Cherepanov
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 119991 Moscow, Kosygina str., 4, Russian Federation; A.N. Belozersky Institute Of Physico-Chemical Biology, Moscow State University, 119992 Moscow, Leninskye gory, 1b.40, Russian Federation.
| | - G E Milanovsky
- A.N. Belozersky Institute Of Physico-Chemical Biology, Moscow State University, 119992 Moscow, Leninskye gory, 1b.40, Russian Federation
| | - K V Neverov
- A.N. Bach Institute of Biochemistry, Federal Research Center "Fundamentals of Biotechnology", Russian Academy of Sciences", 119071 Moscow, Leninsky prospect, 33b.2, Russian Federation; Faculty of Biology, Moscow State University, 119234 Moscow, Leninskye gory, 1b.12, Russian Federation
| | - Yu N Obukhov
- A.N. Bach Institute of Biochemistry, Federal Research Center "Fundamentals of Biotechnology", Russian Academy of Sciences", 119071 Moscow, Leninsky prospect, 33b.2, Russian Federation
| | - Yu V Maleeva
- Faculty of Biology, Moscow State University, 119234 Moscow, Leninskye gory, 1b.12, Russian Federation
| | - A V Aybush
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 119991 Moscow, Kosygina str., 4, Russian Federation
| | - M S Kritsky
- A.N. Bach Institute of Biochemistry, Federal Research Center "Fundamentals of Biotechnology", Russian Academy of Sciences", 119071 Moscow, Leninsky prospect, 33b.2, Russian Federation
| | - V A Nadtochenko
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 119991 Moscow, Kosygina str., 4, Russian Federation; Department of Chemistry, Moscow State University, 119991 Moscow, Leninskye gory, 1b.3, Russian Federation
| |
Collapse
|
4
|
Götze JP, Maity S, Kleinekathöfer U. Incoherent Energy Transfer between the Baseplate and FMO Protein Explored at Ideal Geometries. J Phys Chem B 2023; 127:7829-7838. [PMID: 37691433 DOI: 10.1021/acs.jpcb.3c02568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/12/2023]
Abstract
The Förster resonance energy transfer (FRET) between the Fenna-Matthews-Olson (FMO) protein complex and the chlorosomal baseplate (CBP) is investigated by using an idealized model. This simplified model is based on crystal structure and molecular dynamics conformations. Some of the further input, such as the transition dipole moments, was extracted from earlier molecular-level simulations. The resulting model mimics the effects of the relative position between the CBP and the FMO complex on the corresponding FRET efficiency under ideal conditions, involving about 1.3 billion FRET calculations per investigated model. In this idealized model and employing some approximations, it is found that FRET efficiency is almost completely independent of the FMO trimer orientation (displacement, distance, and rotation), despite FMO and CBP being highly structured complexes. Even removing individual FMO BChl triples will only reduce the FRET efficiency by up to 8.6%. An FMO containing only the least efficient BChl triple will retain about 25% of the FRET efficiency of a full FMO complex. In addition to its proposed function as an energetic funnel, FMO is thus identified to act as a highly robust spatial funnel for CBP excitation harvesting, independent of the mutual CBP-FMO orientation.
Collapse
Affiliation(s)
- Jan P Götze
- Institut für Chemie und Biochemie, Physikalische und Theoretische Chemie, Freie Universität Berlin, Arnimallee 22, 14195 Berlin, Germany
| | - Sayan Maity
- School of Science, Constructor University, Campusring 1, 28759 Bremen, Germany
| | | |
Collapse
|
5
|
Cignoni E, Slama V, Cupellini L, Mennucci B. The atomistic modeling of light-harvesting complexes from the physical models to the computational protocol. J Chem Phys 2022; 156:120901. [DOI: 10.1063/5.0086275] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
The function of light-harvesting complexes is determined by a complex network of dynamic interactions among all the different components: the aggregate of pigments, the protein, and the surrounding environment. Complete and reliable predictions on these types of composite systems can be only achieved with an atomistic description. In the last few decades, there have been important advances in the atomistic modeling of light-harvesting complexes. These advances have involved both the completeness of the physical models and the accuracy and effectiveness of the computational protocols. In this Perspective, we present an overview of the main theoretical and computational breakthroughs attained so far in the field, with particular focus on the important role played by the protein and its dynamics. We then discuss the open problems in their accurate modeling that still need to be addressed. To illustrate an effective computational workflow for the modeling of light harvesting complexes, we take as an example the plant antenna complex CP29 and its H111N mutant.
Collapse
Affiliation(s)
- Edoardo Cignoni
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Vladislav Slama
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
6
|
ZHU ZHE, Higashi M, Saito S. Excited states of chlorophyll a and b in solution by time-dependent density functional theory. J Chem Phys 2022; 156:124111. [DOI: 10.1063/5.0083395] [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
The ground state and excited state electronic properties of chlorophyll (Chl) a and Chl b in diethyl ether, acetone, and ethanol solutions are investigated using quantum mechanical and molecular mechanical calculations with density functional theory (DFT) and time-dependent DFT (TDDFT). Although the DFT/TDDFT methods are widely used, the electronic structures of molecules, especially large molecules, calculated with these methods are known to be strongly dependent on the functionals and the parameters used in functionals. Here, we optimize the range-separated parameter, µ, of the CAM-B3LYP functional of Chl a and Chl b to reproduce the experimental excitation energy differences of these Chl molecules in solution. The optimal values of µ for Chl a and Chl b are smaller than the default value of µ and that for bacteriochlorophyll a, indicating the change in electronic distribution, i.e., an increase in electron delocalization, within the molecule. We find that the electronic distribution of Chl b with an extra formyl group is different from that of Chl a. We also find that the polarity of solution and hydrogen bond cause the decrease in the excitation energies and the increase in the widths of excitation energy distributions of Chl a and Chl b. The present results are expected to be useful for understanding the electronic properties of each pigment molecule in a local heterogeneous environment, which will play an important role in the excitation energy transfer in light-harvesting complex II.
Collapse
Affiliation(s)
| | - Masahiro Higashi
- Department of Molecular Engineering, Kyoto University - Katsura Campus, Japan
| | - Shinji Saito
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Japan
| |
Collapse
|
7
|
Tölle J, Neugebauer J. The Seamless Connection of Local and Collective Excited States in Subsystem Time-Dependent Density Functional Theory. J Phys Chem Lett 2022; 13:1003-1018. [PMID: 35061387 DOI: 10.1021/acs.jpclett.1c04023] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
The theoretical understanding of photoinduced processes in multichromophoric systems requires, as an essential ingredient, the possibility of accurately describing their electronically excited states. However, the size of these systems often prohibits the usage of conventional electronic-structure methods, so that often multiscale approaches based on phenomenologically motivated models are employed. In contrast, subsystem time-dependent density functional theory (sTDDFT) allows for a subsystem-based ab initio description of multichromophoric systems and therefore allows for, in principle, an exact description of photoinduced processes. This Perspective aims to outline the theoretical foundations and commonly used practical realizations as well as to illustrate benefits of recent developments and open issues in the field of sTDDFT. Prospective, potential future applications and possible methodological developments are discussed.
Collapse
Affiliation(s)
- Johannes Tölle
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
8
|
Abstract
Numerous linear and non-linear spectroscopic techniques have been developed to elucidate structural and functional information of complex systems ranging from natural systems, such as proteins and light-harvesting systems, to synthetic systems, such as solar cell materials and light-emitting diodes. The obtained experimental data can be challenging to interpret due to the complexity and potential overlapping spectral signatures. Therefore, computational spectroscopy plays a crucial role in the interpretation and understanding of spectral observables of complex systems. Computational modeling of various spectroscopic techniques has seen significant developments in the past decade, when it comes to the systems that can be addressed, the size and complexity of the sample types, the accuracy of the methods, and the spectroscopic techniques that can be addressed. In this Perspective, I will review the computational spectroscopy methods that have been developed and applied for infrared and visible spectroscopies in the condensed phase. I will discuss some of the questions that this has allowed answering. Finally, I will discuss current and future challenges and how these may be addressed.
Collapse
Affiliation(s)
- Thomas L C Jansen
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
9
|
Scholz L, Neugebauer J. Protein Response Effects on Cofactor Excitation Energies from First Principles: Augmenting Subsystem Time-Dependent Density-Functional Theory with Many-Body Expansion Techniques. J Chem Theory Comput 2021; 17:6105-6121. [PMID: 34524815 DOI: 10.1021/acs.jctc.1c00551] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
We investigate the possibility of describing protein response effects on a chromophore excitation by means of subsystem time-dependent density-functional theory (sTDDFT) in combination with a many-body expansion (MBE) approach. While sTDDFT is in principle intrinsically able to include such contributions, addressing cofactor excitations in protein models or entire proteins with full environment-response treatments is currently out of reach. Taking different model structures of the green fluorescent protein (GFP) and bovine rhodopsin as examples, we demonstrate that an embedded-MBE approach based on sTDDFT in its simplest version leads to a good agreement of the predicted protein response effect already at second order. To reproduce reference response effects from nonsubsystem TDDFT calculations quantitatively (error ≤ 5%), however, a third- or even fourth-order MBE may be required. For the latter case, we explore a selective inclusion of fourth-order terms that drastically reduces the computational burden. In addition, we demonstrate how this sTDDFT-MBE treatment can be utilized as an analysis tool to identify residues with dominant response contributions. This, in turn, can be employed to arrive at smaller structural models for light-absorbing proteins, which still feature all of the main characteristics in terms of photoresponse properties.
Collapse
Affiliation(s)
- Linus Scholz
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
10
|
Tölle J, Deilmann T, Rohlfing M, Neugebauer J. Subsystem-Based GW/Bethe-Salpeter Equation. J Chem Theory Comput 2021; 17:2186-2199. [PMID: 33683119 DOI: 10.1021/acs.jctc.0c01307] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Subsystem Density-Functional Theory and its extension to excited states, namely, subsystem Time-Dependent Density-Functional Theory, have been proven to be efficient and accurate fragmentation approaches for ground and excited states. In the present study we extend this approach to the subsystem-based description of total systems by means of GW and the Bethe-Salpeter equation (BSE). For this, we derive the working equations starting from a subsystem-based partitioning of the screened-Coulomb interaction for an arbitrary number of subsystems. Making use of certain approximations, we develop a parameter-free approach in which environmental screening contributions are effectively included for each subsystem. We demonstrate the applicability of these approximations by comparing quasi-particle energies and excitation energies from subsystem-based GW/BSE calculations to the supermolecular reference. Furthermore, we demonstrate the computational efficiency and the usefulness of this method for the description of photoinduced processes in complex chemical environments.
Collapse
Affiliation(s)
- Johannes Tölle
- Theoretische Organische Chemie Organisch-Chemisches Institut, Westfälische Wilhelms-Universität, Corrensstraße 40, Münster, 48149, Germany
| | | | | | - Johannes Neugebauer
- Theoretische Organische Chemie Organisch-Chemisches Institut, Westfälische Wilhelms-Universität, Corrensstraße 40, Münster, 48149, Germany
| |
Collapse
|
11
|
Wolter M, von Looz M, Meyerhenke H, Jacob CR. Systematic Partitioning of Proteins for Quantum-Chemical Fragmentation Methods Using Graph Algorithms. J Chem Theory Comput 2021; 17:1355-1367. [PMID: 33591754 DOI: 10.1021/acs.jctc.0c01054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Quantum-chemical fragmentation methods offer an efficient approach for the treatment of large proteins, in particular if local target quantities such as protein-ligand interaction energies, enzymatic reaction energies, or spectroscopic properties of embedded chromophores are sought. However, the accuracy that is achievable for such local target quantities intricately depends on how the protein is partitioned into smaller fragments. While the commonly employed naı̈ve approach of using fragments with a fixed size is widely used, it can result in large and unpredictable errors when varying the fragment size. Here, we present a systematic partitioning scheme that aims at minimizing the fragmentation error of a local target quantity for a given maximum fragment size. To this end, we construct a weighted graph representation of the protein, in which the amino acids constitute the nodes. These nodes are connected by edges weighted with an estimate for the fragmentation error that is expected when cutting this edge. This allows us to employ graph partitioning algorithms provided by computer science to determine near-optimal partitions of the protein. We apply this scheme to a test set of six proteins representing various prototypical applications of quantum-chemical fragmentation methods using a simplified molecular fractionation with conjugate caps (MFCC) approach with hydrogen caps. We show that our graph-based scheme consistently improves upon the naı̈ve approach.
Collapse
Affiliation(s)
- Mario Wolter
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstrasse 17, 38106 Braunschweig, Germany
| | - Moritz von Looz
- Department of Computer Science, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany
| | - Henning Meyerhenke
- Department of Computer Science, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany
| | - Christoph R Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstrasse 17, 38106 Braunschweig, Germany
| |
Collapse
|
12
|
Mališ M, Luber S. ΔSCF with Subsystem Density Embedding for Efficient Nonadiabatic Molecular Dynamics in Condensed-Phase Systems. J Chem Theory Comput 2021; 17:1653-1661. [DOI: 10.1021/acs.jctc.0c01200] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Momir Mališ
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
13
|
Osella S. Artificial Photosynthesis: Is Computation Ready for the Challenge Ahead? NANOMATERIALS 2021; 11:nano11020299. [PMID: 33498961 PMCID: PMC7911014 DOI: 10.3390/nano11020299] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/07/2021] [Revised: 01/19/2021] [Accepted: 01/20/2021] [Indexed: 12/13/2022]
Abstract
A tremendous effort is currently devoted to the generation of novel hybrid materials with enhanced electronic properties for the creation of artificial photosynthetic systems. This compelling and challenging problem is well-defined from an experimental point of view, as the design of such materials relies on combining organic materials or metals with biological systems like light harvesting and redox-active proteins. Such hybrid systems can be used, e.g., as bio-sensors, bio-fuel cells, biohybrid photoelectrochemical cells, and nanostructured photoelectronic devices. Despite these efforts, the main bottleneck is the formation of efficient interfaces between the biological and the organic/metal counterparts for efficient electron transfer (ET). It is within this aspect that computation can make the difference and improve the current understanding of the mechanisms underneath the interface formation and the charge transfer efficiency. Yet, the systems considered (i.e., light harvesting protein, self-assembly monolayer and surface assembly) are more and more complex, reaching (and often passing) the limit of current computation power. In this review, recent developments in computational methods for studying complex interfaces for artificial photosynthesis will be provided and selected cases discussed, to assess the inherent ability of computation to leave a mark in this field of research.
Collapse
Affiliation(s)
- Silvio Osella
- Chemical and Biological Systems Simulation Lab, Center of New Technologies, University of Warsaw, Banacha 2C, 02-097 Warsaw, Poland
| |
Collapse
|
14
|
Tölle J, Cupellini L, Mennucci B, Neugebauer J. Electronic couplings for photo-induced processes from subsystem time-dependent density-functional theory: The role of the diabatization. J Chem Phys 2020; 153:184113. [PMID: 33187428 DOI: 10.1063/5.0022677] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Subsystem time-dependent density-functional theory (sTDDFT) making use of approximate non-additive kinetic energy (NAKE) functionals is known to be capable of describing excitation energy transfer processes in a variety of applications. Here, we show that sTDDFT, especially when combined with projection-based embedding (PbE), can be employed for the entire range of photo-induced electronic couplings essential for modeling photophysical properties of complex chemical and biological systems and therefore represents a complete toolbox for this class of problems. This means that it is capable of capturing the interaction/coupling associated with local- and charge-transfer (CT) excitons. However, this requires the choice of a reasonable diabatic basis. We therefore propose different diabatization strategies of the virtual orbital space in PbE-sTDDFT and show how CT excitations can be included in sTDDFT using NAKE functionals via a phenomenological approach. Finally, these electronic couplings are compared to couplings from a multistate fragment excitation difference (FED)-fragment charge difference (FCD) diabatization procedure. We show that both procedures, multistate FED-FCD and sTDDFT (with the right diabatization procedure chosen), lead to an overall good agreement for the electronic couplings, despite differences in their general diabatization strategy. We conclude that the entire range of photo-induced electronic couplings can be obtained using sTDDFT (with the right diabatization procedure chosen) in a black-box manner.
Collapse
Affiliation(s)
- Johannes Tölle
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster Corrensstraße 40, 48149 Münster, Germany
| | - Lorenzo Cupellini
- Università di Pisa, Dipartimento di Chimica e Chimica Industriale, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Università di Pisa, Dipartimento di Chimica e Chimica Industriale, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
15
|
Bondarenko AS, Patmanidis I, Alessandri R, Souza PCT, Jansen TLC, de Vries AH, Marrink SJ, Knoester J. Multiscale modeling of molecular structure and optical properties of complex supramolecular aggregates. Chem Sci 2020; 11:11514-11524. [PMID: 34094396 PMCID: PMC8162738 DOI: 10.1039/d0sc03110k] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022] Open
Abstract
Supramolecular aggregates of synthetic dye molecules offer great perspectives to prepare biomimetic functional materials for light-harvesting and energy transport. The design is complicated by the fact that structure–property relationships are hard to establish, because the molecular packing results from a delicate balance of interactions and the excitonic properties that dictate the optics and excited state dynamics, in turn sensitively depend on this packing. Here we show how an iterative multiscale approach combining molecular dynamics and quantum mechanical exciton modeling can be used to obtain accurate insight into the packing of thousands of cyanine dye molecules in a complex double-walled tubular aggregate in close interaction with its solvent environment. Our approach allows us to answer open questions not only on the structure of these prototypical aggregates, but also about their molecular-scale structural and energetic heterogeneity, as well as on the microscopic origin of their photophysical properties. This opens the route to accurate predictions of energy transport and other functional properties. Multiscale modeling resolves the molecular structure of a synthetic light-harvesting complex, unraveling the microscopic origin of its photophysical properties.![]()
Collapse
Affiliation(s)
- Anna S Bondarenko
- University of Groningen, Zernike Institute for Advanced Materials Groningen The Netherlands
| | - Ilias Patmanidis
- University of Groningen, Zernike Institute for Advanced Materials Groningen The Netherlands .,University of Groningen, Groningen Biomolecular Sciences and Biotechnology Institute Groningen The Netherlands
| | - Riccardo Alessandri
- University of Groningen, Zernike Institute for Advanced Materials Groningen The Netherlands .,University of Groningen, Groningen Biomolecular Sciences and Biotechnology Institute Groningen The Netherlands
| | - Paulo C T Souza
- University of Groningen, Zernike Institute for Advanced Materials Groningen The Netherlands .,University of Groningen, Groningen Biomolecular Sciences and Biotechnology Institute Groningen The Netherlands
| | - Thomas L C Jansen
- University of Groningen, Zernike Institute for Advanced Materials Groningen The Netherlands
| | - Alex H de Vries
- University of Groningen, Zernike Institute for Advanced Materials Groningen The Netherlands .,University of Groningen, Groningen Biomolecular Sciences and Biotechnology Institute Groningen The Netherlands
| | - Siewert J Marrink
- University of Groningen, Zernike Institute for Advanced Materials Groningen The Netherlands .,University of Groningen, Groningen Biomolecular Sciences and Biotechnology Institute Groningen The Netherlands
| | - Jasper Knoester
- University of Groningen, Zernike Institute for Advanced Materials Groningen The Netherlands
| |
Collapse
|
16
|
Analysis of Photosynthetic Systems and Their Applications with Mathematical and Computational Models. APPLIED SCIENCES-BASEL 2020. [DOI: 10.3390/app10196821] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
In biological and life science applications, photosynthesis is an important process that involves the absorption and transformation of sunlight into chemical energy. During the photosynthesis process, the light photons are captured by the green chlorophyll pigments in their photosynthetic antennae and further funneled to the reaction center. One of the most important light harvesting complexes that are highly important in the study of photosynthesis is the membrane-attached Fenna–Matthews–Olson (FMO) complex found in the green sulfur bacteria. In this review, we discuss the mathematical formulations and computational modeling of some of the light harvesting complexes including FMO. The most recent research developments in the photosynthetic light harvesting complexes are thoroughly discussed. The theoretical background related to the spectral density, quantum coherence and density functional theory has been elaborated. Furthermore, details about the transfer and excitation of energy in different sites of the FMO complex along with other vital photosynthetic light harvesting complexes have also been provided. Finally, we conclude this review by providing the current and potential applications in environmental science, energy, health and medicine, where such mathematical and computational studies of the photosynthesis and the light harvesting complexes can be readily integrated.
Collapse
|
17
|
Jones LO, Mosquera MA, Schatz GC, Ratner MA. Embedding Methods for Quantum Chemistry: Applications from Materials to Life Sciences. J Am Chem Soc 2020; 142:3281-3295. [PMID: 31986877 DOI: 10.1021/jacs.9b10780] [Citation(s) in RCA: 58] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Quantum mechanical embedding methods hold the promise to transform not just the way calculations are performed, but to significantly reduce computational costs and improve scaling for macro-molecular systems containing hundreds if not thousands of atoms. The field of embedding has grown increasingly broad with many approaches of different intersecting flavors. In this perspective, we lay out the methods into two streams: QM:MM and QM:QM, showcasing the advantages and disadvantages of both. We provide a review of the literature, the underpinning theories including our contributions, and we highlight current applications with select examples spanning both materials and life sciences. We conclude with prospects and future outlook on embedding, and our view on the use of universal test case scenarios for cross-comparisons of the many available (and future) embedding theories.
Collapse
Affiliation(s)
- Leighton O Jones
- Department of Chemistry , Northwestern University , Evanston , Illinois 60208 , United States
| | - Martín A Mosquera
- Department of Chemistry , Northwestern University , Evanston , Illinois 60208 , United States
| | - George C Schatz
- Department of Chemistry , Northwestern University , Evanston , Illinois 60208 , United States
| | - Mark A Ratner
- Department of Chemistry , Northwestern University , Evanston , Illinois 60208 , United States
| |
Collapse
|
18
|
Tölle J, Böckers M, Niemeyer N, Neugebauer J. Inter-subsystem charge-transfer excitations in exact subsystem time-dependent density-functional theory. J Chem Phys 2019; 151:174109. [DOI: 10.1063/1.5121908] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Johannes Tölle
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Michael Böckers
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Niklas Niemeyer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
19
|
Saito S, Higashi M, Fleming GR. Site-Dependent Fluctuations Optimize Electronic Energy Transfer in the Fenna–Matthews–Olson Protein. J Phys Chem B 2019; 123:9762-9772. [DOI: 10.1021/acs.jpcb.9b07456] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Shinji Saito
- Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan
- The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | - Masahiro Higashi
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto daigaku-katsura, Nishikyo-ku, Kyoto 615-8510, Kyoto, Japan
- JST, PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Graham R. Fleming
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Lab, Berkeley, California 94720, United States
- Kavli Energy Nanosciences Institute at Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
20
|
Schelter I, Foerster JM, Gardiner AT, Roszak AW, Cogdell RJ, Ullmann GM, de Queiroz TB, Kümmel S. Assessing density functional theory in real-time and real-space as a tool for studying bacteriochlorophylls and the light-harvesting complex 2. J Chem Phys 2019; 151:134114. [DOI: 10.1063/1.5116779] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Affiliation(s)
- Ingo Schelter
- Theoretical Physics IV, University of Bayreuth, Bayreuth, Germany
| | - Johannes M. Foerster
- Theoretical Physics IV and Computational Biochemistry, University of Bayreuth, Bayreuth, Germany
| | | | - Aleksander W. Roszak
- Institute of Molecular, Cell and Systems Biology, University of Glasgow, Glasgow, United Kingdom
| | - Richard J. Cogdell
- Institute of Molecular, Cell and Systems Biology, University of Glasgow, Glasgow, United Kingdom
| | | | | | - Stephan Kümmel
- Theoretical Physics IV, University of Bayreuth, Bayreuth, Germany
| |
Collapse
|
21
|
Cupellini L, Bondanza M, Nottoli M, Mennucci B. Successes & challenges in the atomistic modeling of light-harvesting and its photoregulation. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2019; 1861:148049. [PMID: 31386831 DOI: 10.1016/j.bbabio.2019.07.004] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/10/2019] [Revised: 07/10/2019] [Accepted: 07/15/2019] [Indexed: 12/21/2022]
Abstract
Light-harvesting is a crucial step of photosynthesis. Its mechanisms and related energetics have been revealed by a combination of experimental investigations and theoretical modeling. The success of theoretical modeling is largely due to the application of atomistic descriptions combining quantum chemistry, classical models and molecular dynamics techniques. Besides the important achievements obtained so far, a complete and quantitative understanding of how the many different light-harvesting complexes exploit their structural specificity is still missing. Moreover, many questions remain unanswered regarding the mechanisms through which light-harvesting is regulated in response to variable light conditions. Here we show that, in both fields, a major role will be played once more by atomistic descriptions, possibly generalized to tackle the numerous time and space scales on which the regulation takes place: going from the ultrafast electronic excitation of the multichromophoric aggregate, through the subsequent conformational changes in the embedding protein, up to the interaction between proteins.
Collapse
Affiliation(s)
- Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy.
| |
Collapse
|
22
|
Macrocycle ring deformation as the secondary design principle for light-harvesting complexes. Proc Natl Acad Sci U S A 2018; 115:E9051-E9057. [PMID: 30194231 DOI: 10.1073/pnas.1719355115] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Natural light-harvesting is performed by pigment-protein complexes, which collect and funnel the solar energy at the start of photosynthesis. The identity and arrangement of pigments largely define the absorption spectrum of the antenna complex, which is further regulated by a palette of structural factors. Small alterations are induced by pigment-protein interactions. In light-harvesting systems 2 and 3 from Rhodoblastus acidophilus, the pigments are arranged identically, yet the former has an absorption peak at 850 nm that is blue-shifted to 820 nm in the latter. While the shift has previously been attributed to the removal of hydrogen bonds, which brings changes in the acetyl moiety of the bacteriochlorophyll, recent work has shown that other mechanisms are also present. Using computational and modeling tools on the corresponding crystal structures, we reach a different conclusion: The most critical factor for the shift is the curvature of the macrocycle ring. The bending of the planar part of the pigment is identified as the second-most important design principle for the function of pigment-protein complexes-a finding that can inspire the design of novel artificial systems.
Collapse
|
23
|
Jansen TLC. Simple Quantum Dynamics with Thermalization. J Phys Chem A 2018; 122:172-183. [PMID: 29199829 PMCID: PMC5770886 DOI: 10.1021/acs.jpca.7b10380] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2017] [Revised: 12/04/2017] [Indexed: 02/05/2023]
Abstract
In this paper, we introduce two simple quantum dynamics methods. One is based on the popular surface-hopping method, and the other is based on rescaling of the propagation on the bath ground-state potential surface. The first method is special, as it avoids specific feedback from the simulated quantum system to the bath and can be applied for precalculated classical trajectories. It is based on the equipartition theorem to determine if hops between different potential energy surfaces are allowed. By comparing with the formally exact Hierarchical Equations Of Motion approach for four model systems we find that the method generally approximates the quantum dynamics toward thermal equilibrium very well. The second method is based on rescaling of the nonadiabatic coupling and also neglect the effect of the state of the quantum system on the bath. By the nature of the approximations, they cannot reproduce the effect of bath relaxation following excitation. However, the methods are both computationally more tractable than the conventional fewest switches surface hopping, and we foresee that the methods will be powerful for simulations of quantum dynamics in systems with complex bath dynamics, where the system-bath coupling is not too strong compared to the thermal energy.
Collapse
Affiliation(s)
- Thomas L. C. Jansen
- Zernike Institute for Advanced
Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
24
|
Anda A, De Vico L, Hansen T. Intermolecular Modes between LH2 Bacteriochlorophylls and Protein Residues: The Effect on the Excitation Energies. J Phys Chem B 2017; 121:5499-5508. [PMID: 28485594 DOI: 10.1021/acs.jpcb.7b02071] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Light-harvesting system 2 (LH2) executes the primary processes of photosynthesis in purple bacteria; photon absorption, and energy transportation to the reaction center. A detailed mechanistic insight into these operations is obscured by the complexity of the light-harvesting systems, particularly by the chromophore-environment interaction. In this work, we focus on the effects of the protein residues that are ligated to the bacteriochlorophylls (BChls) and construct potential energy surfaces of the ground and first optically excited state for the various BChl-residue systems where we in each case consider two degrees of freedom in the intermolecular region. We find that the excitation energies are only slightly affected by the considered modes. In addition, we see that axial ligands and hydrogen-bonded residues have opposite effects on both excitation energies and oscillator strengths by comparing to the isolated BChls. Our results indicate that only a small part of the chromophore-environment interaction can be associated with the intermolecular region between a BChl and an adjacent residue, but that it may be possible to selectively raise or lower the excitation energy at the axial and planar residue positions, respectively.
Collapse
Affiliation(s)
- André Anda
- Department of Chemistry, H. C. Ørsted Institute, University of Copenhagen , Universitetsparken 5, DK-2100 Copenhagen, Denmark
| | - Luca De Vico
- Department of Chemistry, H. C. Ørsted Institute, University of Copenhagen , Universitetsparken 5, DK-2100 Copenhagen, Denmark.,Department of Biotechnology, Chemistry and Pharmacy, University of Siena , via Aldo Moro 2, 53100 Siena, Italy
| | - Thorsten Hansen
- Department of Chemistry, H. C. Ørsted Institute, University of Copenhagen , Universitetsparken 5, DK-2100 Copenhagen, Denmark
| |
Collapse
|
25
|
Varsano D, Caprasecca S, Coccia E. Theoretical description of protein field effects on electronic excitations of biological chromophores. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:013002. [PMID: 27830666 DOI: 10.1088/0953-8984/29/1/013002] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.
Collapse
Affiliation(s)
- Daniele Varsano
- S3 Center, CNR Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| | | | | |
Collapse
|
26
|
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
|
27
|
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
|
28
|
Stross C, Van der Kamp MW, Oliver TAA, Harvey JN, Linden N, Manby FR. How Static Disorder Mimics Decoherence in Anisotropy Pump-Probe Experiments on Purple-Bacteria Light Harvesting Complexes. J Phys Chem B 2016; 120:11449-11463. [PMID: 27723973 DOI: 10.1021/acs.jpcb.6b09916] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Anisotropy pump-probe experiments have provided insights into the character of excitons formed in photosynthetic complexes. Rapid decay in the observed anisotropy is cited as evidence of the strength of coupling of the excitonic degrees of freedom to their environment. Here we show that ensemble averaging over realistic model Hamiltonians leads to a rapid decay of anisotropy to a value close to the observed asymptote, and at a rate comparable to observed decay rates, even in the absence of coupling to the environment. While coupling to the environment will clearly play a role in the dynamics of such systems, our calculations suggest that caution is needed in deducing the strength of this coupling from anisotropy experiments. We also set out to clarify the nature of the quantum states and processes involved in the dynamics of such systems and the associated terminology.
Collapse
Affiliation(s)
- Clement Stross
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K.,School of Mathematics, University of Bristol , Bristol BS8 1TW, U.K
| | - Marc W Van der Kamp
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | | | - Jeremy N Harvey
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K.,Department of Chemistry, KU Leuven , Celestijnenlaan 200F, B-3001 Leuven, Belgium
| | - Noah Linden
- School of Mathematics, University of Bristol , Bristol BS8 1TW, U.K
| | - Frederick R Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| |
Collapse
|
29
|
Kim CW, Rhee YM. Constructing an Interpolated Potential Energy Surface of a Large Molecule: A Case Study with Bacteriochlorophyll a Model in the Fenna-Matthews-Olson Complex. J Chem Theory Comput 2016; 12:5235-5246. [PMID: 27760297 DOI: 10.1021/acs.jctc.6b00647] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Constructing a reliable potential energy surface (PES) is a key step toward computationally studying the chemical dynamics of any molecular system. The interpolation scheme is a useful tool that can closely follow the accuracy of quantum chemical means at a dramatically reduced computational cost. However, applying interpolation to building a PES of a large molecule is not a straightforward black-box approach, as it frequently encounters practical difficulties associated with its large dimensionality. Here, we present detailed courses of applying interpolation toward building a PES of a large chromophore molecule. We take the example of S0 and S1 electronic states of bacteriochlorophyll a (BChla) molecules in the Fenna-Matthews-Olson light harvesting complex. With a reduced model molecule that bears BChla's main π-conjugated ring, various practical approaches are designed for improving the PES quality in a stable manner and for fine-tuning the final surface such that the surface can be adopted for long time molecular dynamics simulations. Combined with parallel implementation, we show that interpolated mechanics/molecular mechanics (IM/MM) simulations of the entire complex in the nanosecond time scale can be conducted readily without any practical issues. With 1500 interpolation data points for each chromophore unit, the PES error relative to the reference quantum chemical calculation is found to be ∼0.15 eV in the thermally accessible region of the conformational space, together with ∼0.01 eV error in S0 - S1 transition energies. The performance issue related to the use of a large interpolation database within the framework of our parallel routines is also discussed.
Collapse
Affiliation(s)
- Chang Woo Kim
- Center for Self-assembly and Complexity, Institute for Basic Science (IBS) , Pohang 37673, Korea.,Department of Chemistry, Pohang University of Science and Technology (POSTECH) , Pohang 37673, Korea
| | - Young Min Rhee
- Center for Self-assembly and Complexity, Institute for Basic Science (IBS) , Pohang 37673, Korea.,Department of Chemistry, Pohang University of Science and Technology (POSTECH) , Pohang 37673, Korea
| |
Collapse
|
30
|
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: 30] [Impact Index Per Article: 3.8] [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
|
31
|
Higashi M, Saito S. Quantitative Evaluation of Site Energies and Their Fluctuations of Pigments in the Fenna–Matthews–Olson Complex with an Efficient Method for Generating a Potential Energy Surface. J Chem Theory Comput 2016; 12:4128-37. [DOI: 10.1021/acs.jctc.6b00516] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Masahiro Higashi
- Department
of Chemistry, Biology and Marine Science, University of the Ryukyus, 1 Senbaru, Nishihara, Okinawa 903-0213, Japan
| | - Shinji Saito
- Department
of Theoretical and Computational Molecular Science, Institute for Molecular Science, 38 Nishigo-Naka, Myodaiji, Okazaki, Aichi 444-8585, Japan
- The Graduate University for Advanced Studies (SOKENDAI), 38 Nishigo-Naka, Myodaiji, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
32
|
Challenges facing an understanding of the nature of low-energy excited states in photosynthesis. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2016; 1857:1627-1640. [PMID: 27372198 DOI: 10.1016/j.bbabio.2016.06.010] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/18/2016] [Revised: 06/27/2016] [Accepted: 06/28/2016] [Indexed: 01/09/2023]
Abstract
While the majority of the photochemical states and pathways related to the biological capture of solar energy are now well understood and provide paradigms for artificial device design, additional low-energy states have been discovered in many systems with obscure origins and significance. However, as low-energy states are naively expected to be critical to function, these observations pose important challenges. A review of known properties of low energy states covering eight photochemical systems, and options for their interpretation, are presented. A concerted experimental and theoretical research strategy is suggested and outlined, this being aimed at providing a fully comprehensive understanding.
Collapse
|
33
|
Andermatt S, Cha J, Schiffmann F, VandeVondele J. Combining Linear-Scaling DFT with Subsystem DFT in Born–Oppenheimer and Ehrenfest Molecular Dynamics Simulations: From Molecules to a Virus in Solution. J Chem Theory Comput 2016; 12:3214-27. [DOI: 10.1021/acs.jctc.6b00398] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
| | - Jinwoong Cha
- Department
of Materials, ETH Zürich, Zürich, Switzerland
| | - Florian Schiffmann
- Department
of Materials, ETH Zürich, Zürich, Switzerland
- Centre
of Policy Studies, Victoria University, Melbourne, Australia
| | | |
Collapse
|
34
|
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
|
35
|
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
|
36
|
MacGowan SA, Senge MO. Contribution of bacteriochlorophyll conformation to the distribution of site-energies in the FMO protein. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2016; 1857:427-42. [PMID: 26851682 DOI: 10.1016/j.bbabio.2016.02.001] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/04/2015] [Revised: 01/30/2016] [Accepted: 02/01/2016] [Indexed: 10/22/2022]
Abstract
The structural data for the Fenna-Matthews-Olson (FMO) protein indicate that the bacteriochlorophylls (BChls) display a significant degree of conformational heterogeneity of their peripheral substituents and the protein-induced nonplanar skeletal deformations of the tetrapyrrole macrocycle. As electronic properties of chromophores are altered by such differences, a conformational effect may influence the site-energies of specific pigments and thus play a role in mediating the excitation energy transfer dynamics, but this has not yet been established. The difficulty of assessing this question is shown to be partly the result of the inability of the sequential truncation approach usually employed to account for interactions between the conformations of the macrocycle and its substituents and an alternative approach is suggested. By assigning the BChl atoms to meaningful atom groups and performing all possible permutations of partial optimizations in a full-factorial design, where each group is either frozen in the crystal geometry or optimized in vacuo, followed by excited state calculations on each resulting structure (PM6//ZIndo/S), the specific effects of the conformations of each BChl component as well as mutual interactions between the molecular fragments on the site-energy can be delineated. This factorial relaxation procedure gives different estimates of the macrocycle conformational perturbation than the approach of sequentially truncating the BChl periphery. The results were evaluated in the context of published site-energies for the FMO pigments from three species to identify how conformational effects contribute to their distribution and instances of cross-species conservation and functional divergence of the BChl nonplanarity conformational contribution are described.
Collapse
Affiliation(s)
- Stuart A MacGowan
- School of Chemistry, SFI Tetrapyrrole Laboratory, Trinity Biomedical Sciences Institute, 152-160 Pearse Street, Trinity College Dublin, The University of Dublin, Dublin 2, Ireland
| | - Mathias O Senge
- School of Chemistry, SFI Tetrapyrrole Laboratory, Trinity Biomedical Sciences Institute, 152-160 Pearse Street, Trinity College Dublin, The University of Dublin, Dublin 2, Ireland.
| |
Collapse
|
37
|
Cupellini L, Jurinovich S, Prandi IG, Caprasecca S, Mennucci B. Photoprotection and triplet energy transfer in higher plants: the role of electronic and nuclear fluctuations. Phys Chem Chem Phys 2016; 18:11288-96. [DOI: 10.1039/c6cp01437b] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Photosynthetic organisms employ several photoprotection strategies to avoid damage due to the excess energy in high light conditions.
Collapse
Affiliation(s)
- Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
| | - Sandro Jurinovich
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
| | - Ingrid G. Prandi
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
| | - Stefano Caprasecca
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
| |
Collapse
|
38
|
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
|
39
|
Jia X, Mei Y, Zhang JZH, Mo Y. Hybrid QM/MM study of FMO complex with polarized protein-specific charge. Sci Rep 2015; 5:17096. [PMID: 26611739 PMCID: PMC4661465 DOI: 10.1038/srep17096] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2015] [Accepted: 10/26/2015] [Indexed: 12/20/2022] Open
Abstract
The Fenna-Matthews-Olson (FMO) light-harvesting complex is now one of the primary model systems for the study of excitation energy transfer (EET). However, the mechanism of the EET in this system is still controversial. In this work, molecular dynamics simulations and the electrostatic-embedding quantum-mechanics/molecular-mechanics single-point calculations have been employed to predict the energy transfer pathways utilizing the polarized protein-specific charge (PPC), which provides a more realistic description of Coulomb interaction potential in the protein than conventional mean-field charge scheme. The recently discovered eighth pigment has also been included in this study. Comparing with the conventional mean-field charges, more stable structures of FMO complex were found under PPC scheme during molecular dynamic simulation. Based on the electronic structure calculations, an exciton model was constructed to consider the couplings during excitation. The results show that pigments 3 and 4 dominate the lowest exciton levels whereas the highest exciton level are mainly constituted of pigments 1 and 6. This observation agrees well with the assumption based on the spatial distribution of the pigments. Moreover, the obtained spectral density in this study gives a reliable description of the diverse local environment embedding each pigment.
Collapse
Affiliation(s)
- Xiangyu Jia
- State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - John Z H Zhang
- State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Department of Chemistry, New York University, New York, NY 10003
| | - Yan Mo
- State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
| |
Collapse
|
40
|
Goez A, Neugebauer J. A Local Variant of the Conductor-Like Screening Model for Fragment-Based Electronic-Structure Methods. J Chem Theory Comput 2015; 11:5277-90. [DOI: 10.1021/acs.jctc.5b00832] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Albrecht Goez
- Theoretische Organische Chemie,
Organisch-Chemisches Institut and Center for Multiscale Theory and
Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße
40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie,
Organisch-Chemisches Institut and Center for Multiscale Theory and
Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße
40, 48149 Münster, Germany
| |
Collapse
|
41
|
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
|
42
|
Artiukhin DG, Jacob CR, Neugebauer J. Excitation energies from frozen-density embedding with accurate embedding potentials. J Chem Phys 2015; 142:234101. [DOI: 10.1063/1.4922429] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Denis G. Artiukhin
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Christoph R. Jacob
- Institute of Physical and Theoretical Chemistry, TU Braunschweig, Hans-Sommer-Straße 10, 38106 Braunschweig, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
43
|
Wesolowski TA, Shedge S, Zhou X. Frozen-Density Embedding Strategy for Multilevel Simulations of Electronic Structure. Chem Rev 2015; 115:5891-928. [DOI: 10.1021/cr500502v] [Citation(s) in RCA: 219] [Impact Index Per Article: 24.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Tomasz A. Wesolowski
- Department of Physical Chemistry, University of Geneva, CH-1211 Geneva 4, Switzerland
| | - Sapana Shedge
- Department of Physical Chemistry, University of Geneva, CH-1211 Geneva 4, Switzerland
| | - Xiuwen Zhou
- Department of Physical Chemistry, University of Geneva, CH-1211 Geneva 4, Switzerland
| |
Collapse
|
44
|
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
|
45
|
Fornaro T, Carnimeo I, Biczysko M. Toward Feasible and Comprehensive Computational Protocol for Simulation of the Spectroscopic Properties of Large Molecular Systems: The Anharmonic Infrared Spectrum of Uracil in the Solid State by the Reduced Dimensionality/Hybrid VPT2 Approach. J Phys Chem A 2014; 119:5313-26. [DOI: 10.1021/jp510101y] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Teresa Fornaro
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
| | - Ivan Carnimeo
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
- Compunet, Istituto Italiano di Tecnologia (IIT), Via Morego 30, I-16163 Genova, Italy
| | - Malgorzata Biczysko
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
- Istituto di Chimica
dei Composti OrganoMetallici (ICCOM-CNR), UOS di Pisa, Area della
Ricerca CNR, Consiglio Nazionale delle Ricerche, Via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
46
|
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
|
47
|
Goerigk L, Collyer CA, Reimers JR. Recommending Hartree–Fock Theory with London-Dispersion and Basis-Set-Superposition Corrections for the Optimization or Quantum Refinement of Protein Structures. J Phys Chem B 2014; 118:14612-26. [DOI: 10.1021/jp510148h] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Lars Goerigk
- School
of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Charles A. Collyer
- School
of Molecular Bioscience, The University of Sydney, Sydney, New South Wales 2006, Australia
| | - Jeffrey R. Reimers
- Centre
for Quantum and Molecular Structure, College of Sciences, Shanghai University, Shanghai 200444, China
- School
of Physics and Advanced Materials, The University of Technology, Sydney, New South Wales 2007, Australia
| |
Collapse
|
48
|
Higashi M, Kosugi T, Hayashi S, Saito S. Theoretical study on excited states of bacteriochlorophyll a in solutions with density functional assessment. J Phys Chem B 2014; 118:10906-18. [PMID: 25153487 DOI: 10.1021/jp507259g] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The excited-state properties of bacteriochlorophyll (BChl) a in triethylamine, 1-propanol, and methanol are investigated with the time-dependent density functional theory by using the quantum mechanical and molecular mechanical reweighting free energy self-consistant field method. It is found that no prevalent density functionals can reproduce the experimental excited-state properties, i.e., the absorption and reorganization energies, of BChl a in the solutions. The parameter μ in the range-separated hybrid functional is therefore optimized to reproduce the differences of the absorption energies in the solutions. We examine the origin of the differences of the absorption energies in the solutions and find that sensitive balance between contributions of structural changes and solute-solvent interactions determines the differences. The accurate description of the excitation with the density functional with the adjusted parameter is therefore essential to the understanding of the excited-state properties of BChl a in proteins and also the mechanism of the photosynthetic systems.
Collapse
Affiliation(s)
- Masahiro Higashi
- Department of Chemistry, Biology and Marine Science, University of the Ryukyus , 1 Senbaru, Nishihara, Okinawa 903-0213, Japan
| | | | | | | |
Collapse
|
49
|
Barone V, Biczysko M, Borkowska-Panek M, Bloino J. A Multifrequency Virtual Spectrometer for Complex Bio-Organic Systems: Vibronic and Environmental Effects on the UV/Vis Spectrum of Chlorophylla. Chemphyschem 2014; 15:3355-64. [DOI: 10.1002/cphc.201402300] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2014] [Indexed: 11/07/2022]
|
50
|
Jacob CR, Neugebauer J. Subsystem density-functional theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2014. [DOI: 10.1002/wcms.1175] [Citation(s) in RCA: 229] [Impact Index Per Article: 22.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Christoph R. Jacob
- Center for Functional Nanostructures and Institute of Physical Chemistry; Karlsruhe Institute of Technology (KIT); Karlsruhe Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut; Westfälische Wilhelms-Universität Münster; Münster Germany
| |
Collapse
|