1
|
Kim Y, Mitchell Z, Lawrence J, Morozov D, Savikhin S, Slipchenko LV. Predicting Mutation-Induced Changes in the Electronic Properties of Photosynthetic Proteins from First Principles: The Fenna-Matthews-Olson Complex Example. J Phys Chem Lett 2023; 14:7038-7044. [PMID: 37524046 DOI: 10.1021/acs.jpclett.3c01461] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2023]
Abstract
Multiscale molecular modeling is utilized to predict optical absorption and circular dichroism spectra of two single-point mutants of the Fenna-Matthews-Olson photosynthetic pigment-protein complex. The modeling approach combines classical molecular dynamics simulations with structural refinement of photosynthetic pigments and calculations of their excited states in a polarizable protein environment. The only experimental input to the modeling protocol is the X-ray structure of the wild-type protein. The first-principles modeling reproduces changes in the experimental optical spectra of the considered mutants, Y16F and Q198V. Interestingly, the Q198V mutation has a negligible effect on the electronic properties of the targeted bacteriochlorophyll a pigment. Instead, the electronic properties of several other pigments respond to this mutation. The molecular modeling demonstrates that a single-point mutation can induce long-range effects on the protein structure, while extensive structural changes near a pigment do not necessarily lead to significant changes in the electronic properties of that pigment.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Zach Mitchell
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Jack Lawrence
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Dmitry Morozov
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, P.O. Box 35, 40014 Jyväskylä, Finland
| | - Sergei Savikhin
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| |
Collapse
|
2
|
Kaiser A, Daoud RE, Aquilante F, Kühn O, De Vico L, Bokarev SI. A Multiconfigurational Wave Function Implementation of the Frenkel Exciton Model for Molecular Aggregates. J Chem Theory Comput 2023; 19:2918-2928. [PMID: 37115036 DOI: 10.1021/acs.jctc.3c00185] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
We present an implementation of the Frenkel exciton model into the OpenMolcas program package enabling calculations of collective electronic excited states of molecular aggregates based on a multiconfigurational wave function description of the individual monomers. The computational protocol avoids using diabatization schemes and, thus, supermolecule calculations. Additionally, the use of the Cholesky decomposition of the two-electron integrals entering pair interactions enhances the efficiency of the computational scheme. The application of the method is exemplified for two test systems, that is, a formaldehyde oxime and a bacteriochlorophyll-like dimer. For the sake of comparison with the dipole approximation, we restrict our considerations to situations where intermonomer exchange can be neglected. The protocol is expected to be beneficial for aggregates composed of molecules with extended π systems, unpaired electrons such as radicals or transition metal centers, where it should outperform widely used methods based on time-dependent density functional theory.
Collapse
Affiliation(s)
- Andy Kaiser
- Institut für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock. Germany
| | - Razan E Daoud
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Francesco Aquilante
- Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Oliver Kühn
- Institut für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock. Germany
| | - Luca De Vico
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Sergey I Bokarev
- Institut für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock. Germany
- Chemistry Department, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany
| |
Collapse
|
3
|
Maity S, Kleinekathöfer U. Recent progress in atomistic modeling of light-harvesting complexes: a mini review. PHOTOSYNTHESIS RESEARCH 2023; 156:147-162. [PMID: 36207489 PMCID: PMC10070314 DOI: 10.1007/s11120-022-00969-w] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Accepted: 09/22/2022] [Indexed: 06/16/2023]
Abstract
In this mini review, we focus on recent advances in the atomistic modeling of biological light-harvesting (LH) complexes. Because of their size and sophisticated electronic structures, multiscale methods are required to investigate the dynamical and spectroscopic properties of such complexes. The excitation energies, in this context also known as site energies, excitonic couplings, and spectral densities are key quantities which usually need to be extracted to be able to determine the exciton dynamics and spectroscopic properties. The recently developed multiscale approach based on the numerically efficient density functional tight-binding framework followed by excited state calculations has been shown to be superior to the scheme based on pure classical molecular dynamics simulations. The enhanced approach, which improves the description of the internal vibrational dynamics of the pigment molecules, yields spectral densities in good agreement with the experimental counterparts for various bacterial and plant LH systems. Here, we provide a brief overview of those results and described the theoretical foundation of the multiscale protocol.
Collapse
Affiliation(s)
- Sayan Maity
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759, Bremen, Germany
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759, Bremen, Germany.
| |
Collapse
|
4
|
Reiter S, Kiss FL, Hauer J, de Vivie-Riedle R. Thermal site energy fluctuations in photosystem I: new insights from MD/QM/MM calculations. Chem Sci 2023; 14:3117-3131. [PMID: 36970098 PMCID: PMC10034153 DOI: 10.1039/d2sc06160k] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2022] [Accepted: 02/03/2023] [Indexed: 02/09/2023] Open
Abstract
Cyanobacterial photosystem I (PSI) is one of the most efficient photosynthetic machineries found in nature. Due to the large scale and complexity of the system, the energy transfer mechanism from the antenna complex to the reaction center is still not fully understood. A central element is the accurate evaluation of the individual chlorophyll excitation energies (site energies). Such an evaluation must include a detailed treatment of site specific environmental influences on structural and electrostatic properties, but also their evolution in the temporal domain, because of the dynamic nature of the energy transfer process. In this work, we calculate the site energies of all 96 chlorophylls in a membrane-embedded model of PSI. The employed hybrid QM/MM approach using the multireference DFT/MRCI method in the QM region allows to obtain accurate site energies under explicit consideration of the natural environment. We identify energy traps and barriers in the antenna complex and discuss their implications for energy transfer to the reaction center. Going beyond previous studies, our model also accounts for the molecular dynamics of the full trimeric PSI complex. Via statistical analysis we show that the thermal fluctuations of single chlorophylls prevent the formation of a single prominent energy funnel within the antenna complex. These findings are also supported by a dipole exciton model. We conclude that energy transfer pathways may form only transiently at physiological temperatures, as thermal fluctuations overcome energy barriers. The set of site energies provided in this work sets the stage for theoretical and experimental studies on the highly efficient energy transfer mechanisms in PSI.
Collapse
Affiliation(s)
- Sebastian Reiter
- Department of Chemistry, Ludwig-Maximilians-Universität München Butenandtstr. 11 81377 Munich Germany
| | - Ferdinand L Kiss
- Department of Chemistry, Ludwig-Maximilians-Universität München Butenandtstr. 11 81377 Munich Germany
| | - Jürgen Hauer
- Department of Chemistry, Technical University of Munich Lichtenbergstr. 4, Garching 85747 Germany
| | - Regina de Vivie-Riedle
- Department of Chemistry, Ludwig-Maximilians-Universität München Butenandtstr. 11 81377 Munich Germany
| |
Collapse
|
5
|
Feighan O, Manby FR, Bourne-Worster S. An efficient protocol for excited states of large biochromophores. J Chem Phys 2023; 158:024107. [PMID: 36641400 DOI: 10.1063/5.0132417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Qy transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory's accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies.
Collapse
Affiliation(s)
- Oliver Feighan
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Frederick R Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Susannah Bourne-Worster
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
6
|
Förster A, Visscher L. Quasiparticle Self-Consistent GW-Bethe-Salpeter Equation Calculations for Large Chromophoric Systems. J Chem Theory Comput 2022; 18:6779-6793. [PMID: 36201788 PMCID: PMC9648197 DOI: 10.1021/acs.jctc.2c00531] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The GW-Bethe–Salpeter equation
(BSE) method
is promising for calculating the low-lying excitonic states of molecular
systems. However, so far it has only been applied to rather small
molecules and in the commonly implemented diagonal approximations
to the electronic self-energy, it depends on a mean-field starting
point. We describe here an implementation of the self-consistent and
starting-point-independent quasiparticle self-consistent (qsGW)-BSE approach, which is suitable for calculations on
large molecules. We herein show that eigenvalue-only self-consistency
can lead to an unfaithful description of some excitonic states for
chlorophyll dimers while the qsGW-BSE vertical excitation
energies (VEEs) are in excellent agreement with spectroscopic experiments
for chlorophyll monomers and dimers measured in the gas phase. Furthermore,
VEEs from time-dependent density functional theory calculations tend
to disagree with experimental values and using different range-separated
hybrid (RSH) kernels does change the VEEs by up to 0.5 eV. We use
the new qsGW-BSE implementation to calculate the
lowest excitation energies of the six chromophores of the photosystem
II (PSII) reaction center (RC) with nearly 2000 correlated electrons.
Using more than 11,000 (6000) basis functions, the calculation could
be completed in less than 5 (2) days on a single modern compute node.
In agreement with previous TD-DFT calculations using RSH kernels on
models that also do not include environmental effects, our qsGW-BSE calculations only yield states with local characters
in the low-energy spectrum of the hexameric complex. Earlier works
with RSH kernels have demonstrated that the protein environment facilitates
the experimentally observed interchromophoric charge transfer. Therefore,
future research will need to combine correlation effects beyond TD-DFT
with an explicit treatment of environmental electrostatics.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HVAmsterdam, The Netherlands
| | - Lucas Visscher
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HVAmsterdam, The Netherlands
| |
Collapse
|
7
|
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
|
8
|
Reiter S, Bäuml L, Hauer J, de Vivie-Riedle R. Q-Band relaxation in chlorophyll: new insights from multireference quantum dynamics. Phys Chem Chem Phys 2022; 24:27212-27223. [DOI: 10.1039/d2cp02914f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The ultrafast relaxation within the Q-bands of chlorophyll plays a crucial role in photosynthetic light-harvesting. We investigate this process via nuclear and electronic quantum dynamics on multireference potential energy surfaces.
Collapse
Affiliation(s)
- Sebastian Reiter
- Department of Chemistry, Ludwig-Maximilians-Universität München, Butenandtstr. 11, 81377 Munich, Germany
| | - Lena Bäuml
- Department of Chemistry, Ludwig-Maximilians-Universität München, Butenandtstr. 11, 81377 Munich, Germany
| | - Jürgen Hauer
- Department of Chemistry, Technical University of Munich, Lichtenbergstr. 4, 85747 Garching, Germany
| | - Regina de Vivie-Riedle
- Department of Chemistry, Ludwig-Maximilians-Universität München, Butenandtstr. 11, 81377 Munich, Germany
| |
Collapse
|
9
|
|
10
|
Anda A, Cole JH. Two-dimensional spectroscopy beyond the perturbative limit: The influence of finite pulses and detection modes. J Chem Phys 2021; 154:114113. [PMID: 33752354 DOI: 10.1063/5.0038550] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
Ultra-fast and multi-dimensional spectroscopy gives a powerful looking glass into the dynamics of molecular systems. In particular, two-dimensional electronic spectroscopy (2DES) provides a probe of coherence and the flow of energy within quantum systems, which is not possible with more conventional techniques. While heterodyne-detected (HD) 2DES is increasingly common, more recently fluorescence-detected (FD) 2DES offers new opportunities, including single-molecule experiments. However, in both techniques, it can be difficult to unambiguously identify the pathways that dominate the signal. Therefore, the use of numerically modeling of 2DES is vitally important, which, in turn, requires approximating the pulsing scheme to some degree. Here, we employ non-perturbative time evolution to investigate the effects of finite pulse width and amplitude on 2DES signals. In doing so, we identify key differences in the response of HD and FD detection schemes, as well as the regions of parameter space where the signal is obscured by unwanted artifacts in either technique. Mapping out parameter space in this way provides a guide to choosing experimental conditions and also shows in which limits the usual theoretical approximations work well and in which limits more sophisticated approaches are required.
Collapse
Affiliation(s)
- André Anda
- ARC Centre of Excellence in Exciton Science and Chemical and Quantum Physics, School of Science, RMIT University, Melbourne, Australia
| | - Jared H Cole
- ARC Centre of Excellence in Exciton Science and Chemical and Quantum Physics, School of Science, RMIT University, Melbourne, Australia
| |
Collapse
|
11
|
Hashemi Z, Leppert L. Assessment of the Ab Initio Bethe-Salpeter Equation Approach for the Low-Lying Excitation Energies of Bacteriochlorophylls and Chlorophylls. J Phys Chem A 2021; 125:2163-2172. [PMID: 33656894 PMCID: PMC8028335 DOI: 10.1021/acs.jpca.1c01240] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Bacteriochlorophyll and chlorophyll molecules are crucial building blocks of the photosynthetic apparatus in bacteria, algae, and plants. Embedded in transmembrane protein complexes, they are responsible for the primary processes of photosynthesis: excitation energy and charge transfer. Here, we use ab initio many-body perturbation theory within the GW approximation and Bethe-Salpeter equation (BSE) approach to calculate the electronic structure and optical excitations of bacteriochlorophylls a, b, c, d, and e and chlorophylls a and b. We systematically study the effects of the structure, basis set size, partial self-consistency in GW, and the underlying exchange-correlation approximation and compare our calculations with results from time-dependent density functional theory, multireference RASPT2, and experimental literature results. We find that optical excitations calculated with GW+BSE are in excellent agreement with experimental data, with an average deviation of less than 100 meV for the first three bright excitations of the entire family of (bacterio)chlorophylls. Contrary to state-of-the-art time-dependent density functional theory (TDDFT) with an optimally tuned range-separated hybrid functional, this accuracy is achieved in a parameter-free approach. Moreover, GW+BSE predicts the energy differences between the low-energy excitations correctly and eliminates spurious charge transfer states that TDDFT with (semi)local approximations is known to produce. Our study provides accurate reference results and highlights the potential of the GW+BSE approach for the simulation of larger pigment complexes.
Collapse
Affiliation(s)
- Zohreh Hashemi
- Institute of Physics, University of Bayreuth, Bayreuth 95440, Germany
| | - Linn Leppert
- Institute of Physics, University of Bayreuth, Bayreuth 95440, Germany.,MESA+ Institute for Nanotechnology, University of Twente, 7500 AE Enschede, The Netherlands
| |
Collapse
|
12
|
Sirohiwal A, Berraud-Pache R, Neese F, Izsák R, Pantazis DA. Accurate Computation of the Absorption Spectrum of Chlorophyll a with Pair Natural Orbital Coupled Cluster Methods. J Phys Chem B 2020; 124:8761-8771. [PMID: 32930590 PMCID: PMC7584356 DOI: 10.1021/acs.jpcb.0c05761] [Citation(s) in RCA: 55] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
![]()
The
ability to accurately compute low-energy excited states of
chlorophylls is critically important for understanding the vital roles
they play in light harvesting, energy transfer, and photosynthetic
charge separation. The challenge for quantum chemical methods arises
both from the intrinsic complexity of the electronic structure problem
and, in the case of biological models, from the need to account for
protein–pigment interactions. In this work, we report electronic
structure calculations of unprecedented accuracy for the low-energy
excited states in the Q and B bands of chlorophyll a. This is achieved by using the newly developed domain-based local
pair natural orbital (DLPNO) implementation of the similarity transformed
equation of motion coupled cluster theory with single and double excitations
(STEOM-CCSD) in combination with sufficiently large and flexible basis
sets. The results of our DLPNO–STEOM-CCSD calculations are
compared with more approximate approaches. The results demonstrate
that, in contrast to time-dependent density functional theory, the
DLPNO–STEOM-CCSD method provides a balanced performance for
both absorption bands. In addition to vertical excitation energies,
we have calculated the vibronic spectrum for the Q and B bands through
a combination of DLPNO–STEOM-CCSD and ground-state density
functional theory frequency calculations. These results serve as a
basis for comparison with gas-phase experiments.
Collapse
Affiliation(s)
- Abhishek Sirohiwal
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.,Fakultät für Chemie und Biochemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Romain Berraud-Pache
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Róbert Izsák
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Dimitrios A Pantazis
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
13
|
Kjær C, Gruber E, Nielsen SB, Andersen LH. Color tuning of chlorophyll a and b pigments revealed from gas-phase spectroscopy. Phys Chem Chem Phys 2020; 22:20331-20336. [PMID: 32895686 DOI: 10.1039/d0cp03210g] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Chlorophyll (Chl) pigments are responsible for vital mechanisms in photosynthetic proteins: light harvesting, energy transfer and charge separation. A complex interplay between the Chl molecule and its microenvironment determines its transition energy. Interactions such as excitonic coupling with one or more pigments (Chls or carotenoids), axial ligation to the magnesium center, or electrostatic interactions between Chl and nearby amino-acid residues all influence the photophysical properties. Here we use time-resolved photodissociation action spectroscopy to determine transition energies of Chla/b complexes in vacuo to directly compare the impact of a negatively charged axial ligand (formate) to that of exciton coupling between two Chls. Experiments carried out at the electrostatic ion storage ring ELISA allow dissociation to be sampled on hundreds of milliseconds time scale. Absorption-band maxima of Chla-formate complexes are found at 433 ± 4 nm/2.86 ± 0.03 eV (Soret band) and in the region 654-675 nm/1.84-1.90 eV (Q band) and those of Chla dimers tagged by a quaternary ammonium ion at 419 ± 5 nm/2.96 ± 0.04 eV (Soret band) and 647 nm/1.92 eV (Q band). The axial ligand strongly affects the Chla transition energies causing redshifts of 0.21 eV of the Soret band and 0.04-0.1 eV of the Q band compared to Chla tagged by a quaternary ammonium. Slightly smaller shifts were found in case of Chlb. The redshifts are approximately twice that induced by excitonic coupling between two Chlas, also tagged by a quaternary ammonium ion. Axial ligation brings the absorption by isolated Chls very close to that of photosynthetic proteins.
Collapse
Affiliation(s)
- Christina Kjær
- Department of Physics and Astronomy, Aarhus University, Denmark.
| | - Elisabeth Gruber
- Department of Physics and Astronomy, Aarhus University, Denmark.
| | | | - Lars H Andersen
- Department of Physics and Astronomy, Aarhus University, Denmark.
| |
Collapse
|
14
|
Aquilante F, Autschbach J, Baiardi A, Battaglia S, Borin VA, Chibotaru LF, Conti I, De Vico L, Delcey M, Fdez Galván I, Ferré N, Freitag L, Garavelli M, Gong X, Knecht S, Larsson ED, Lindh R, Lundberg M, Malmqvist PÅ, Nenov A, Norell J, Odelius M, Olivucci M, Pedersen TB, Pedraza-González L, Phung QM, Pierloot K, Reiher M, Schapiro I, Segarra-Martí J, Segatta F, Seijo L, Sen S, Sergentu DC, Stein CJ, Ungur L, Vacher M, Valentini A, Veryazov V. Modern quantum chemistry with [Open]Molcas. J Chem Phys 2020; 152:214117. [PMID: 32505150 DOI: 10.1063/5.0004835] [Citation(s) in RCA: 220] [Impact Index Per Article: 55.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.
Collapse
Affiliation(s)
- Francesco Aquilante
- Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, Buffalo, New York 14260-3000, USA
| | - Alberto Baiardi
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Stefano Battaglia
- Department of Chemistry - BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
| | - Veniamin A Borin
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Liviu F Chibotaru
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Irene Conti
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, Bologna I-40136, Italy
| | - Luca De Vico
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, via Aldo Moro 2, 53100 Siena, Italy
| | - Mickaël Delcey
- Department of Chemistry - Ångström Laboratory, Uppsala University, SE-751 21 Uppsala, Sweden
| | - Ignacio Fdez Galván
- Department of Chemistry - BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
| | - Nicolas Ferré
- Aix-Marseille University, CNRS, Institut Chimie Radicalaire, Marseille, France
| | - Leon Freitag
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, Bologna I-40136, Italy
| | - Xuejun Gong
- Department of Chemistry, University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Stefan Knecht
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Ernst D Larsson
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, Lund 22100, Sweden
| | - Roland Lindh
- Department of Chemistry - BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
| | - Marcus Lundberg
- Department of Chemistry - Ångström Laboratory, Uppsala University, SE-751 21 Uppsala, Sweden
| | - Per Åke Malmqvist
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, Lund 22100, Sweden
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, Bologna I-40136, Italy
| | - Jesper Norell
- Department of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Michael Odelius
- Department of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Massimo Olivucci
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, via Aldo Moro 2, 53100 Siena, Italy
| | - Thomas B Pedersen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Laura Pedraza-González
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, via Aldo Moro 2, 53100 Siena, Italy
| | - Quan M Phung
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Chikusa, Nagoya 464-8602, Japan
| | - Kristine Pierloot
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Igor Schapiro
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Javier Segarra-Martí
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, 80 Wood Lane, London W12 0BZ, United Kingdom
| | - Francesco Segatta
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, Bologna I-40136, Italy
| | - Luis Seijo
- Departamento de Química, Instituto Universitario de Ciencia de Materiales Nicolás Cabrera, and Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Saumik Sen
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | | | - Christopher J Stein
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Liviu Ungur
- Department of Chemistry, University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Morgane Vacher
- Laboratoire CEISAM - UMR CNRS 6230, Université de Nantes, 44300 Nantes, France
| | - Alessio Valentini
- Theoretical Physical Chemistry, Research Unit MolSys, Université de Liège, Allée du 6 Août, 11, 4000 Liège, Belgium
| | - Valera Veryazov
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, Lund 22100, Sweden
| |
Collapse
|