|Home | About | Journals | Submit | Contact Us | Français|
We review the optical properties of the FMO complex as found by spectroscopic studies of the Qy band over the last two decades. This article emphasizes the different methods used, both experimental and theoretical, to elucidate the excitonic structure and dynamics of this pigment–protein complex.
The online version of this article (doi:10.1007/s11120-010-9540-1) contains supplementary material, which is available to authorized users.
In 1975, Fenna was the first to resolve the X-ray structure of the Fenna–Matthews–Olson (FMO) complex of Prosthecochloris aestuarii. In photosynthetic membranes of green sulfur bacteria, this protein channels the excitations from the chlorosomes to the reaction center. Since it was the first photosynthetic antenna complex of which the X-ray structure became available, it triggered a wide variety of studies of spectroscopic and theoretical nature, and it therefore has become one of the most widely studied and well-characterized pigment–protein complexes. Owing to its relatively simple structure amongst the light-harvesting complexes, with only seven interacting bacteriochlorophyll a (BChl a) molecules, and with the level of sophistication at which the optical properties are known, it comes as no surprise that the FMO complex serves as a guinea pig for new and ever-improving simulation methods as well as new optical techniques. Remarkably, FMO is still a subject of active investigation and new insights continue to emerge. Even fundamental properties, such as the pigment–protein ratio, remain controversial. The goal of this article is to guide the reader through the mass of information that has appeared over the last ~20 years on the optical properties of the FMO complex. We attempt to provide an objective view of the experimental data and the parameters and methods used in simulations. Also, where applicable, it is indicated which data and parameter sets have become most favored and for which reasons. In order to keep this article insightful and focused, it is restricted to a discussion of the spectral structure of the Qy transition band of a BChl a molecule at 800 nm. This article will specifically address optical properties of the FMO protein from the most thoroughly characterized green sulfur bacterium Prosthecochloris aestuarii. Similar data on the FMO protein from Chlorobium tepidum can be found in the electronic supplementary material. We hope that this article will enable researchers applying novel, complex techniques to gain an overview of the current understanding of the optical properties and photosynthetic processes in the FMO protein.
In 1962, John Olson isolated a water-soluble bacteriochlorophyll (BChl a) protein (150 kDa) from green sulfur bacteria (Olson and Romano 1962). This specific protein is part of the light-harvesting system in green sulfur bacteria where it acts as a subantenna to collect sunlight and transfer excitation energy from the light-harvesting antennas to the reaction center. Absorption spectroscopy on extracts of strains of Chlorobium showed that the newly discovered protein contained only BChl a chromophores, non-covalently bound to a protein envelope (Fig. 1). In 1975, Roger Fenna and Brian Matthews resolved the X-ray structure of the FMO protein from Prosthecochloris aestuarii at 2.8 Å resolution and found that the complex consists of three identical subunits related by C3 symmetry, each containing seven BChl a pigments (Fenna and Matthews 1975). It showed a protein shell in which the BChl a molecules were enclosed. The major part of the outside of the protein shell exposed to the solvent is composed of 15 strands of β-sheet. The side of the shell that is in contact with one of the other subunits in the trimer consists of four short strands of α-helix alternated by regions of the protein without a clear structure. The average distance between BChl a molecules within one subunit of the trimer is 12 Å while the nearest molecule in the neighboring subunit is found at a distance of 24 Å. Analysis of the X-ray data showed no evidence for interactions—whether these be covalent or noncovalent—between neighboring BChl a molecules; however, the same analysis predicted the presence of extensive interactions between the chlorophyll molecules and the protein shell. Besides hydrophobic interactions, hydrogen bonding and coordination to the Mg ion in the BChl a molecule occurs. Over the years, the structure of the FMO protein from Prosthecochloris aestuarii has been refined (Matthews et al. 1979; Tronrud et al. 1986) and recently a 1.3 Å diffraction dataset of the structure has been obtained (Tronrud et al. 2009).
In 1997, the crystal structure of FMO from Chlorobium tepidum was determined at a 2.2 Å resolution (Li et al. 1997). Similar to Prosthecochloris aestuarii (Fig. 1a), the distance between the pigments within a subunit range from 4 to 11 Å while the nearest neighbor in a different subunit is over 20 Å away. The β sheet is folded in such a way that the strands at the front and the back of the shell are roughly perpendicular to each other (Fig. 1b). The opening in the shell is situated toward the center of the trimer, forming the shape of a shell. The six α-helices are located at the open end of the shell and mainly connect the separated β-strands. BChl a molecules 1 and 2 are situated at the outside of the protein complex, while BChl a 3–7 are located in the center (Fig. 1c). Polar interactions and salt bridges between amino acids insure the formation of a stable trimer. The magnesium ion is a five-coordinate in all the BChl a molecules, although the fifth ligand varies between the pigments. For BChl a 1, 3, 4, 6, and 7, it is a histidine residue, for BChl a 5, it is an oxygen atom from a leucine residue, and for BChl a 2, the electron density suggests a water molecule as the fifth ligand. The structures of the FMO protein present in the two species Prosthecochloris aestuarii and Chlorobium tepidum show a high degree of similarity (the amino acid sequences are identical to one another within 77%). The residues that are not conserved do not alter the interaction between the protein and the BChl a molecules. Besides that, the relative positions of each of the BChl a molecules in the two species match almost perfectly. The main difference is in the planarity of the tetrapyrrole ring of the BChl a molecules. For a more detailed description of the comparison between the two species, see Li et al. (1997) and the discussion at the end of this section.
Various spectroscopic investigations using linear absorption spectroscopy, circular dichroism (CD) and linear dichroism (LD) on samples of the isolated FMO protein and the protein associated with membrane vesicles have revealed the orientation of the proteins with respect to the membrane (Melkozernov et al. 1998). The three subunits of the FMO protein are related by C3 symmetry and can be modeled as disks, with the axis of the disks parallel to the C3 axis (Fig. 2a). The spectroscopic studies show that the C3 symmetry axis of the three subunits of the FMO protein is perpendicular to the membrane plane. This implies that the flat sides of the discs is embedded in the membrane (Fig. 2a).
In two recent studies, the presence of an additional BChl a molecule per monomer was proposed. This observation is based on careful studies of high resolution X-ray data. Ben-Shem et al. noticed additional electron density at the interface between the monomers in their newly crystalized and solved structure. They attributed this to the presence of an eighth BChl a molecule possibly located at the interface between the FMO protein and RC (Ben-Shem et al. 2004). In an attempt to clarify matters, Tronrud et al. (2009) decided to revisit the structure of Chlorobium tepidum as well as collect a new diffraction dataset at 1.3 Å of the FMO protein from Prosthecochloris aestuarii. Their comparison indicated the presence of an eighth BChl a molecule at the same location in both variants, however, with a different local protein structure that could account for the difference in the optical spectra (see “Linear spectra”). The nature of the eighth BChl a molecule is different from the other seven: its occupancy is not unity and it is located in the region of the protein that is directed towards the chlorosome. Its location and the orientation of its transition dipole moment, i.e., parallel to the BChl a in the baseplate, might facilitate energy transfer. In both variants, a carbonyl oxygen binds to the central magnesium atom on one side of the BChl a ring while an α-helix covers the other side. It was shown that between the two variants there are three critical differences concerning the amino acid sequence in this helix, close to the additional BChl a molecule. In Prosthecochloris aestuarii at residue 165, threonine is changed into phenylalanine and at residue 168, alanine is changed into serine. In addition, in the loop that directs the helix back to the protein, an alanine is inserted. These three mutations have the following effect in Prosthecochloris aestuarii: on binding of the eighth BChl a molecule, the side chain of the Phenylalanine has to move out of the binding pocket. As a result, the α-helix moves sufficiently close to the Mg atom to make an additional link, creating a bidentate interaction between protein and BChl a. However, in Chlorobium tepidum, the smaller Threonine does not move on binding of the BChl; on top of that, the shorter loop of the α-helix restricts motion preventing bidentate binding. The differences in binding of this extra BChl a molecule is expected to have a considerable influence on the optical spectra, especially on the CD spectra (vide infra).
This section describes the various aspects that come into play on describing and simulating the optical spectra of the FMO complex. First, the differences between the low-temperature absorption spectra of Prosthecochloris aestuarii and Chlorobium tepidum are discussed. This is followed by an account on the site energies of the BChl a molecules. These values cannot be deduced from optical experiments directly and are usually obtained by fits to optical spectra; however, recent attempts to calculate the site energies have emerged. Simulations of the optical spectra are extremely sensitive to the exact choice of site energies, and hence, a detailed overview of the results of different research groups is provided. Subsequently, a third important optical property of the FMO complex is discussed: the pigment with the lowest site energy. This pigment is thought to be the exit pigment of the complex, i.e., the pigment that transfers the excitation energy to the reaction center. As the individual BChl a molecules interact within the FMO complex, the exciton nature of their excitation is treated and exciton simulations, used to generate various linear spectra, are described. Important parameters in these simulations are the dipolar coupling strength and the linewidth of the transitions. The section ends with a discussion of the controversial nature of the lowest energy absorption band at 825 nm.
Over the years, simulations of the linear spectra have become increasingly sophisticated. Whereas early on, almost all optical properties were hotly debated, in recent times, the tendency is to use parameter sets and methods as obtained and developed by Louwe et al. The validity of their study also extends into the nonlinear regime, as is the topic of the next section.
The linear absorption spectrum of the FMO complex shows several bands in the wavelength range of 200–900 nm (Olson 2004). The Qy (S1) absorption band around 800 nm is the most well-characterized band and the focus of the current study. In membrane factions of Chlorobium tepidum, this band appears in the spectral region between the absorption band of BChl c in the chlorosomes (720–750 nm) and the Qy band of the BChl a in the reaction center at ~834 nm (Melkozernov et al. 1998). The Qy (S1) absorption band has a temperature-dependent shape. At cryogenic temperatures, in a mixture of Tris buffer and glycerol, the absorption band consists of at least three distinct peaks (Johnson and Small 1991; Gulbinas et al. 1996) (Fig. 3). At elevated temperatures, the fine structure disappears, and the absorption spectrum appears as a broad featureless band.
Low-temperature absorption spectra of the Qy (S1) band show a clear difference between the FMO complex of Prosthecochloris aestuarii and Chlorobium tepidum; the former has a strong absorption band at 815 nm, while for the latter, the strongest absorption band is at 809 nm. Comparison between the two species with 97% homology (Chlorobium limicola and Chlorobium tepidum) shows a nearly identical absorption spectrum at 6 K. This indicates that the local protein environment has a limited but observable influence in the spectral differences between the FMO complexes (Francke and Amesz 1997). Li et al. (1997) showed that most ligands of the BChl a pigments are conserved in both species and propose that the spectral differences find their origin in the non-planarity of the ring in the BChl a molecules. Recently, Tronrud et al. showed that the difference in absorption spectra of the FMO complex of various green sulfur bacteria can be explained by the structure. As described in the previous section, an additional BChl a molecule has been observed. Three mutations in the α-helix, covering this molecule, lead to a bidentate binding between pigment and protein in the FMO complex from Prosthecochloris aestuarii. As the other seven BChl a molecules are nearly identical, Tronrud et al. ascribe the differences in the spectra to the presence or absence of the additional link to the eighth BChl a molecule. To support this point, a sequence alignment of the FMO protein of several species was performed. This showed that the three mutations, described above, tend to appear together. However, on top of that, the mutations correlate with the type of spectra, i.e., similar to Prosthecochloris aestuarii in the presence of the mutations, and similar to Chlorobium tepidum in the absence of the mutations.
One of the most debated properties of the FMO complex concerns the site energies of the seven BChl a molecules in the complex. These values are needed for exciton calculations of the linear spectra and simulations of dynamics. They are defined as the transition energy of a pigment in the absence of coupling between the pigments. It does, however, depend on local interactions between the BChl a molecule and the protein envelope, and includes electrostatic interactions and ligation. Since the interactions are difficult to identify and even harder to quantify, the site energies are usually treated as independent parameters that are obtained from a simultaneous fit to several optical spectra. Table 1 gives an overview of the different site energies determined by various research groups, using a range of methods described in this section. One of the main differences between the approaches, to obtain the site energies by simulating the spectra, is whether they restrict the interactions to BChl a molecules within a subunits or wether they include interactions in the whole trimer. These two approaches are labeled in Table 1 with M (only include interactions within a monomer) and T (allow interactions between BChl a molecules in the whole trimer).
At the beginning of the 1990s, the optical spectra were fit, assuming interactions between the BChl a pigments from different subunits in one trimer (Johnson and Small 1991; Van Mourik et al. 1994; Rätsep and Freiberg 2007). Although previous efforts to model the system using the full trimer geometry had not been very successful, Pearlstein still expected the C3 symmetry of the system to amplify the coupling effect between the intersubunit BChl a molecules (Pearlstein 1992). In contrast to earlier simulations, in his later studies, different site energies were assigned to the 21 transitions. Instead of a single transitions at 802.6 nm, 21 site energies were used as fitting parameters, and the best fit was judged by eye. A mixed approach was employed by Lu et al. and Gülen et al.; the full trimer was taken into account while simultaneously fitting linear optical spectra. However, the same site energies were assigned to the symmetry related BChl a pigments, resulting in seven adjustable site energies (Lu and Pearlstein 1993; Gülen 1996). This approach implies that, although there are only seven different site energies assigned, all the 21 possible exciton transitions in the trimer will be included in the fits (vide infra). Lu and Pearlstein (1993) restricted the interactions to a single subunit and improved the fits from Pearlstein, making use of an algorithm to minimize the difference between the measured and the simulated spectra with various adjustable parameters, amongst which are the seven site energies of the monomer. Their fits were based on two sets of absorption and CD spectra at 77 K, obtained by two different groups (referred to as 1 and 2 in Table Table1).1). A similar approach was used by Gülen et al. In contrast to the earlier fits by Pearlstein and Lu et al., CD spectra were excluded from the fits, since they tend to be very sensitive to the experimental conditions like the choice of solvent. Figure 2b shows directions of the individual (not excitonic) transition dipole moments with respect to the C3 axis: BChl a pigments 7, 1, and 4 lie almost parallel to the C3 axis, while the orientation of the dipole moments of BChl a 6, 2, 5, and 3 is almost perpendicular. Gülen used the spatial organization of the individual dipole moments to help restrict and direct the fit. As a start of the fit, the energy of BChl a 6 was fixed between 815 and 820 nm. In order to reproduce the spectra, the site energies of the remaining BChl a pigments were found to be around 805 nm. Fine tuning of the fits was done by the naked eye.
In contrast to the trimer approach, a different group of researchers fitted the optical spectra, only allowing for interactions within one subunit, the monomer approach. Louwe et al. were among the first to use the monomer approach. Similar to Pearlstein, the site energies were obtained by means of adjusting the parameters manually in the simulations of the spectra, starting from a common site energy at 809.7 nm (Louwe et al. 1997b). Four possible parameter sets were obtained based on the orientation of the transition dipole moments, as shown previously by Gülen et al. Three of these improved the other existing simulations. However, only one of the basis sets, containing the seven site energies, produced simulations resembling the shape of the spectra (see Table Table1).1). Vulto et al. attempted to simulate the excited state dynamics using the site energies as proposed by Louwe et al. For a satisfactory fit, the site energies needed to be adapted slightly (see Table 1; Vulto et al. 1999). Simulations of both time-resolved and steady-state spectra were the aim of Iseri et al. The site energies were used as free parameters in a manual-fitting routine (Iseri and Gülen 1999). As reported in a previos study by Gülen et al., the signs of the bands in the LD spectra limits the choice of site energies as they impose a restriction on the direction of the dipole moments with respect to the C3 symmetry axis (see Fig. 2b). An improved fit of absorption and LD spectra was obtained using the site energies as proposed by Louwe et al. and included spectral broadening (vide infra) (Wendling et al. 2002). Further improvements were instigated by a global fit of absorption, CD, and LD spectra. The site energies that were found in these fits are stated in 1, and they are obtained assuming two different types of broadening, denoted by the numbers 1* and 2*.
Adolphs and Renger (2006) used a different approach by calculating the “electrochromic shifts” of the site energies by taking into account the interaction between charged amino acids and the pigments. The individual electrochromic shifts were calculated using the Coulomb coupling between the charged amino acids, approximated by point charges, and the difference between the permanent dipole moments of the BChl a ground and excited state, estimated from Stark experiments. Remarkable is that the red shift of BChl a 3 and the blue shift of BChl a 6 are caused by charged amino acids that are conserved in the structures of Prosthecochloris aestuarii and Chlorobium tepidum. Adolphs et al. show that the fits of the seven site energies for the monomeric and the trimeric structure give similar results. The current method of calculating site energies only succeeded partially in reproducing the site energies obtained from fits to the linear spectra. Therefore, a more elaborate model was needed for better agreement.
Further elaboration, and as a result a good agreement with the experimental spectra, is provided by calculations of the site energies including a detailed description of the charge distribution of both the pigment and the protein states (Müh et al. 2007). Ab initio methods were used to describe the pigments, while a classical electrostatic method was used to describe the whole complex on the atomic level. As a result of the low dielectric constant of water/glycerol below the freezing point, the standard protonation pattern of the amino acids was no longer valid and half of the usually acidic and basic groups turned out to be neutral. This complex method was simplified, without losing the main results by assuming a standard protonation pattern and by the introduction of an effective dielectric constant for screening effects (Adolphs et al. 2008). There exists an earlier account of similar quantum calculations where, amongst others, the effect of the charged amino acids was included (Gudowksa-Nowak et al. 1990). However, the resulting site energies are spread over a range (~770–840 nm) much larger than what is observed in spectra, hence, and these results are not used for exciton calculations.
While, for some of the earlier calculations and fits, the range of site energies only spans 10 nm, the more recent ones seem to converge to a difference between the highest and lowest site energy of almost 30 nm, which is comparable to the total width of the absorption spectrum. The most widely accepted values of the site energies for Prosthecochloris aestuarii are given by Louwe et al. (see Table 1). Nevertheless, improvements have been obtained using more and more elaborate models and by calculations of the site energies rather than fitting them. In general, only seven different site energies are included as parameters in the fits, however wether or not to include interaction between the monomers remains controversial.
The pigment with the lowest site energy is the most likely candidate for an exit pigment, which transfers the excitation energy from the FMO complex to the reaction center. The position of this pigment within the FMO complex cannot be detected optically because this would require a resolution below the diffraction limit, and, therefore, it can only be assigned from the outcome of exciton simulations. However, since photosynthesis occurs at 300 K, at room temperature, none of the exciton states should be excluded from, a transition dipole-weighted, energy transfer to the reaction core complex. Table 2 shows the different “exit pigments” that have been proposed, with consensus now leaning toward pigment 3. A detailed account on the nature of the electronic state of the exit pigment will be given in “Nature of the lowest energy band”.
Pearlstein pioneered the approach of finding the best site energies by looking at absorption, CD, and hole-burning spectra, rather than just at spectra from one experimental technique. His fits showed that BChl a 7 has the lowest site energy (Pearlstein 1992) Subsequently the lowest energy pigment was assigned to be BChl a 7 or 3 depending on the fitted dataset (Lu and Pearlstein 1993). The site energies simultaneously fitted to absorption, LD, and singlet–triplet spectra (Gülen 1996) brought BChl a 6 forward as the pigment with the lowest site energy. It is the best interconnected pigment. Simulations by Buck et al. favored BChl a 7 for that role. They obtained the best fit using parameters deduced from optical spectra of Olson et al. (1976), in which BChl a 7 has the lowest site energy (Buck et al. 1997). By means of fitting new LD and CD data, Louwe et al. (1997b) concluded that the exciton states are mainly localized on one BChl a and that the lowest energy pigment was BChl a 3. This agrees with the results from Stark hole-burning experiments (Rätsep et al. 1998). Since Since then, different theoretical and experimental approaches agree on BChl a 3 being the pigment with the lowest site energy (Vulto et al. 1999; Iseri and Gülen 1999; Wendling et al. 2002; Adolphs and Renger 2006; Müh et al. 2007; Adolphs et al. 2008). Electron microscopy showed the arrangement of the FMO complex with respect to the reaction center (RC) (Rémigy et al. 1999). The technique lacks the resolution to distinguish between the top and the bottom of the FMO complex. However, from the shape of the FMO complex, it can be deduced that either BChl a 1 and 6 or 3 and 4 form the exit pigments from FMO protein to RC. Wen et al. used mass spectrometry to infer the orientation of the FMO complex with respect to the RC, which is embedded in the cytoplastic membrane (Wen et al. 2009). Their results, in agreement with the theoretical predictions, showed that the BChl a 3 side of the FMO complex interacts with the membrane. Hence, pigment number 3 is the closest to the RC and, therefore, likely to be the exit pigment.
By taking a closer look at the environment of BChl a 3, which is generally assumed to have the largest electrochromic shift to lower site energy, a curious arrangement of α-helices was observed (Müh et al. 2007). The dipoles of the two helices can be represented by two partial charges on the ends of the helix. The positive and negative partial charges of helix 5 lie in the negative and positive regions, respectively, of the calculated difference (S0 − S1) electrostatic potential. This results in a red shift of the site energies of about 200 cm−1. Similarly, helix 6 gives an additional, however somewhat reduced, redshift of ~100 cm−1. Both helices have a definite effect on the site energies without being dependent on protonation states.
The close proximity of the BChl a molecules (~10 Å) leads to electronic coupling between them that exceeds the electron-vibrational coupling in the FMO complex. Therefore, the system is usually described by a superposition of the seven molecular BChl a states forming seven exciton states (Van Amerongen et al. 2000), and the electron-vibrational coupling is treated perturbationally. These excitonic interactions Vij between two chromophores i and j are dominated by the relative orientation of the transition dipole moments and the inverse cube of the distance between the BChl a molecules. The exciton levels have different cross sections and linewidths which together with a close energy spacing results in a dense and complex spectrum. This can be seen in the low-temperature absorption spectra in which only three peaks out of seven are clearly visible.
In order to describe the excitonic wavefunctions in the FMO protein, the following electronic Hamiltonian is used:
in which Ej represents the site energies of the uncoupled BChl a molecules and |j the corresponding localized excitations. The exciton wavefunctions |α are obtained by diagonalizing the Hamiltonian as
using the exciton expansion coefficients Cα(i) which represent the contribution of the individual BChl a molecules to an excitonic transition. The results of such calculations, as performed on the same system by a variety of research groups, are shown in the Tables 3, ,4,4, and and5,5, where α runs vertically and i horizontally. Pearlstein used a point-monopole approach to describe the interaction between the individual BChl a molecules. The transition charge density is calculated for each molecule, represented by point charges at the position of individual atoms, and the interactions of all point charges with those of the other chromophore are considered. All the 21 BChl a of the trimeric FMO complex were included in the model, and the parameters for all the 21 degenerate and non-degenerate exciton transitions are displayed in the original article (Pearlstein 1992). It can be concluded that in each case only one or two of the BChl a pigments contribute significantly to the squared amplitude of the eigenvectors of the transitions. This means that none of the exciton states is delocalized over the complete subunit, let alone the trimer. Later this was verified by using a similar approach to model the absorption spectra (Gülen 1996). However, the main difference is that now the simpler point-dipole approximation was used to model the interactions between the excitonic transitions of two chromophores by representing the electronic transition of each molecule by a single dipole moment. Later, an approach by the same group followed, which restricted the dipole interactions to one monomer only (Iseri and Gülen 1999).
Results from linear–dichroic absorbance-detected magnetic resonance experiments on FMO at 1.2 K exhibited similar results as monomeric BChl a molecules in organic solvents. This technique is sensitive to the triplet state of the complex and, therefore, it was concluded that in FMO, the triplet state is localized on a single BChl a pigment and not on its delocalized trimeric counterpart (Louwe et al. 1997a). Simultaneous simulation of the spectra obtained from this technique together with CD spectra were performed considering a single subunit only (Louwe et al. 1997b). This approach was justified by the fact that the simulations predict exciton states that are mainly dominated by a single BChl a, implying that the degree of exciton delocalization is limited in the FMO complex.
For exciton simulations of the various spectra (e.g., absorption, LD, CD) of the FMO protein there are three basic ingredients: the site energies, the dipolar coupling (coupling strength), and the optical linewidth. The first is treated in “Site energies”, while the latter two will be discussed in this section. The transition dipole strength of BChl a d(D2) in the FMO complex, has been reported several times in literature. It is an important parameter in simulations of the optical spectra. The values of this dipole strength vary widely and range between 20 and 60 D2. Simulations by Pearlstein revealed a dipole coupling strength with a value of 51.6 D2 (Pearlstein 1992). This value is similar to the one he used in previous calculations and corresponds to the value of 50.8 D2 used by Fenna. Further successful simulations of steady-state and time-resolved experiments were obtained using values of 51 D2 (Renger and May 1998) and 30-40 D2 (Iseri and Gülen 1999; Wendling et al. 2002). This value was verified by calculations, which resulted in a value of the effective dipole strength of 30 D2 (Adolphs and Renger 2006) obtained by reducing the dipole strength in vacuum by a factor of 1.25.
Broadening in optical spectra has two distinct origins, both of which are of importance in the spectroscopic studies of the FMO complex (May and Kühn 2000). The first phenomenon that causes line broadening is static disorder. The seven pigments in the FMO complex all have a slightly different local environment, since the protein envelope that surrounds them differs from pigment to pigment. As a result, there is a different mean energy, center absorption frequency, for each BChl a. Owing to the differences between, for example, the solvation of all BChl a 1 pigments in the sample, the center absorption frequency of this pigment is broadened. This effect is referred to as inhomogeneous broadening and can lead to a broad band in the linear absorption spectrum. Inhomogeneous broadening is included in the description of optical spectra in two ways: by including a variable linewidth or by introducing one linewidth for all transitions. An example of the first is given by Pearlstein, who employed widths in the range of ~80 to ~170 cm−1 although there was no physical justification for this large difference (Pearlstein 1992). Exciton simulations by Buck et al. (1997) were performed using ~150 cm−1 for all the transitions in the complex and, therefore, discarded the effect of inhomogeneous broadening shown by Pearlstein to be effective in simulation. Around the same time, linewidths obtained from hole-burning experiments, ~70–80 cm−1, were employed by two sets of authors (Gülen 1996; Wendling et al. 2000) to simulate absorption, linear dichroism, singlet–triplet and low-temperature absorption and fluorescence line-narrowing measurements, respectively. Several successful simulations of both steady-state and time-resolved spectra were performed using an inhomogeneous linewidth of ~80 cm−1 (Louwe et al. 1997b; Vulto et al. 1998a, b, 1999).
Besides inhomogeneous broadening, a second physical process that is thought to contribute to broadening of the linewidths is important in the FMO complex. If the changes in the molecular properties are fast compared to the duration of the measurements, then dynamic disorder occurs. This causes homogeneous broadening, because it is equal for the whole system. Initially, hole-burning spectra provided a way to obtain the homogeneous linewidths and revealed values of ~70−80 cm−1 (Johnson and Small 1991). A better description of the spectra was subsequently obtained by fully including the effects of different types of broadening to an existing model proposed earlier by the same authors (Wendling et al. 2002). The two types of broadening were included in simulations of new LD and CD spectra at low temperatures describing the whole trimer. Inhomogeneous broadening due to the variation in site energies in between subunits and complexes could especially influence the simulations of the polarized spectra. Subsequently, the authors added homogeneous broadening due to dephasing, the lifetimes of the exciton states were calculated using their exciton model. Even without changing the site energies and coupling strengths from reference (Louwe et al. 1997b), the absorption spectra were reproduced better taking broadening into account in the system. The simulations of the LD and CD spectra were further improved by fitting the site energies and the coupling strengths to the experiments using a global fit.
In order to determine the different exciton states and the accompanying transition energy, several approaches were used. To begin with, in the reference (Johnson and Small 1991) exciton energies are determined by simultaneous analysis of different hole-burning spectra. In this case, eight exciton components were observed of which the latter two were assigned to contribute to one band around 825 nm (vide infra). Pearlstein followed a similar procedure and fitted 21 exciton energies (of which 14 degenerate, see Table 6) to absorption and CD spectra (Pearlstein 1992). There are two more reports on the exciton levels in the trimer, both based on the method described by Pearlstein (Lu and Pearlstein 1993; Gülen 1996). Improvements were made using algorithms to fit the spectra and changing the site wavelengths, which are used to determine the exciton levels, respectively.
In further attempts to model the spectra, only monomers containing seven BChl a molecules are taken into account (see Table 7). This results in a structure with seven interconnected exciton levels. These simulations require the site energies of the BChl a molecules as input parameters. Louwe et al. (1997b) set a standard for simulations of the spectra. The interaction energies are calculated in the point-dipole approximation assuming a common linewidth for all transitions of ~80 cm−1. Screening by the protein is taken into account by a dielectric constant that was used a global free-fit parameter. The initial calculated dipole strength of 68.9 D2 is thus reduced by a factor 2.4 leading to an effective dipole strength of 28.7 D2, a value that is lower than that proposed by Pearlstein (1992). This value is close to a physically relevant value of the reduced dipole strength in the range of 25–40 D2. In order to simulate the spectra, a minimum of free parameters was used to fit the essential features of the spectra. The authors proposed that the model can be improved by inclusion of vibrations, lifetime broadening of the highest energy exciton states, and by allowing for different dipole strengths for the individual BChl a molecules and a variation of the dielectric constant over the protein. Simulations based on the same exciton model were performed by the following research groups: Vulto et al. (1998a, 1999), Wendling et al. (2000, 2002), and Iseri and Gülen (1999).
The assignment of the bands in the absorption spectrum, especially of the band, the lowest in energy at 825 nm, has proven to be difficult. The number of excitonic states and their respective energies have been the subject of intense debate. Johnson and Small (1991) concluded that lower and higher spectral energy features flanking the hole-burning line can only be explained when excitonic interactions between the BChls are taken into account. Furthermore, the results of spectral hole burning show the presence of eight states. Two of those eight identified exciton states, which have perpendicular symmetry, contribute to this lowest exciton band at 825 nm. Models excluding the interactions between the subunits of the trimer are not successful in describing this experimental data (Johnson and Small 1991). Therefore, Johnson and Small (1991) have developed a model in which this interaction is included leading to a maximum of 14 delocalized states (21 states in total, of which 14 are degenerate). This implies that the 825-nm band comprises of three, slightly shifted, bands of the subunits, of which two are degenerate. For the space group C3, the states having E symmetry are degenerate while the states with A symmetry are not (Atkins 1995). Owing to the varying local environment, there is a certain degree of static disorder in the system leading to an inhomogeneously broadened 825-nm band. The nature of the 825-nm band was confirmed to have a double origin seven years later by means of Stark hole-burning studies (Rätsep et al. 1998). However, in this case, the nature of these states was assumed to be much more localized, with the excitons mainly spread over one BChl a molecule. Structural heterogeneity in the complex leads to a variation in the excitation energy of the lowest energy state in the subunits of the trimer. This view was tested by temperature-dependent hole-burning experiments on the FMO protein from Chlorobium tepidum (Rätsep et al. 1999). The 825-nm absorption band was fitted with three Gaussian bands of ~55cm−1 at 823.0, 825.0, and 827.0 nm, respectively. The dependence of hole width and hole growing kinetics on the burning frequency confirms that there are three bands contributing to the 825-nm band. Triplet minus singlet (T − S) spectra measured by Louwe et al. (1997a) shows that the triplet state is localized on a single BChl a since it demonstrates the same properties as monomeric BChl a a in organic solvents. The orientations of excitonic transitions in the Qy band were determined relative to the triplet-carrying molecule. In contrast to earlier measurements, fluorescence line narrowing experiments showed that the 825-nm absorption band can be accounted for by a single transition in the range of temperature from 4 K to room temperature (Wendling et al. 2000). This transition is coupled to protein phonons and vibrations in the chromophore. The effect of disorder on the lowest energy band in the trimer was further studied by Monte Carlo simulations (Hayes et al. 2002). The lowest energy band could be fitted with three nearly Gaussian bands of almost identical intensity. One of those band was centered at the absorption maximum of the 825 nm band, while the maxima of the other two bands where shifted by ~−17 and ~+26cm−1, respectively. Summarizing, the outcome of different experimental techniques do not agree on the nature of the 825 nm band. While some state that this band is due to a single transition, others include a distribution of the lowest exciton energy in the different subunits of the trimer to account for the observations.
This section will discuss both the experimental and theoretical aspects of the time-resolved spectra of the FMO protein. Previously, in “Exciton nature of the BChl a excitations in the FMO protein” and “Coupling strengths, linewidth, and exciton energies”, the excitonic structure and simulations of the linear optical spectra were reviewed. Starting from this knowledge, it is a small, yet complex step to simulate the time-dependent behavior of the exciton states. After optical excitation, the population in the exciton states eventually decays back to the ground state. The relaxation pathways and the accompanying time scales have been studied intensely by several nonlinear optical techniques and are often accompanied by theoretical investigations. This section will discuss hole-burning experiments, followed by pump-probe and photon-echo experiments, 2D electronic experiments, and finally new theoretical approaches.
Modeling of the exciton dynamics in the BChl a chromophore complex of the FMO protein has been done using two approaches. The first describes energy transfer between chromophores by the incoherent Förster hopping rate equation, which is valid for weak coupling between the chromophores and a strong coupling of the electronic transition to vibrational states, precluding the formation of exciton levels. Excitation energy will hop from one molecule to the other along the energy gradient. However, since the existence of exciton levels in the FMO complex is well established, the Förster hopping rate equation seems not to be the most appropriate way to describe dynamics in the FMO complex. This problem was partially overcome by Iseri et al. who approximated the energy transfer rate between excitons through a linear combination of the Förster rates between the BChl a pigments that dominate the exciton states (Iseri and Gülen 1999). The second approach is to describe the light-induced dissipative dynamics within the framework of the multi-exciton density matrix theory. Often, the Redfield approach for the description of dissipation is used. This theory combines the time-dependent Schrödinger equation for the excitonic transitions with a linear coupling to a classical bath, given by all the vibrational modes of the chromophore complex (Renger and May 1998; Vulto et al. 1999; Brüggemann and May 2004; Brüggemann et al. 2006). Finally, a modified Redfield approach valid for intermediate coupling regimes has been applied by Read et al. (2008).
After all the light-induced coherences have vanished, the time evolution of the excitonic state populations Pα, where for the FMO protein α runs from 1 to 7, can be described by the Master equation (Van Amerongen et al. 2000).
using the rate constants kα→β, which eventually lead to a thermal equilibrium within the singly excited states.
The proposed pathways of downward energy transfer are shown schematically in Fig. 4, as drawn by the respective authors. Although they show little agreement, a few general conclusions can be drawn from these results. The energy transfer from the highest to the lowest exciton level occurs on a very fast time scale; within 5 ps, mainly the lowest exciton state P1 is populated. The population can be transferred downward either by a few big steps or by small steps including all the exciton levels.
Spectral dynamics, in terms of hole widths, obtained from hole-burning experiments follow a temperature dependence power law Tα, with the temperature exponent for glasses α ~1.3 (Matsuzaki et al. 2000). Such a power law is typical for dephasing of the excitons in a pigment coupled to a two-level system, TLS) (Yamaguchi et al. 2002). Low-frequency excitations in a glass are often described by a TLS, modeled by a double-well potential. These excitations can contribute to the dephasing of a pigment, and hence determine the hole widths. Three states were found to contribute to the absorption band at 825 nm. Taking into account the dephasing due to the glasslike protein, the energy transfer between the three levels within the 825-nm band occurs with 99 and 26 ps, respectively (Matsuzaki et al. 2000). Similar reasoning holds for analysis of low-energy states in the FMO protein from Chlorobium tepidum (Rätsep et al. 1999). In order to bridge the gap between steady-state and time-resolved spectroscopy an elaborate hole-burning experiment was performed (Franken et al. 1998). On top of broad (800–820 nm) uncorrelated signals, sharp holes were detected. The observed hole widths are for an inhomogeneously broadened band twice the homogeneous linewidth, from which it is straightforward to calculate the excited state lifetimes (see Table 8). The lifetimes of the exciton states that were obtained from hole-burning studies were fast, (sub)picosecond, and similar to those obtained from other methods (vide infra).
When researchers started to study the excitation energy transfer within the FMO complex in the early 1990s, they soon realized that the dynamics occur on very fast, subpicosecond, timescales. By studying the bleach spectrum at 2 and 10 ps after excitation, it was shown that even at those short delay times, the spectrum does not exhibit a uniform bleach (Lyle and Struve 1990). In this study, the anisotropy decay was 2–4 ps. As was known from the linewidths of hole burning, the relaxation between exciton levels is complete within several hundreds of femtoseconds (Johnson and Small 1991) and does not contribute to one color anisotropy decay. Therefore, the longer, picosecond, time constant obtained from anisotropic decay traces was attributed to hopping of excitation energy between neighboring subunits and not to lifetimes of the higher exciton states.
The obtained dephasing times from hole-burning experiments are considerably faster than values that were obtained from accumulated photon-echo experiments by Louwe and Aartsma (1994). They demonstrated that the photon-echo decay is wavelength and temperature-dependent, ranging from a decay <1 ps at 790 nm to 100 ps at 824 nm at 10 K and exhibits multiexponential decay. They also observed a decrease of the decay times with increasing temperatures. The wavelength-dependent decay rates from the photon-echo experiments are explained on the basis of phonon-assisted dephasing, where the number of lower lying states determine the dephasing time. Initially, it was thought that the relaxation was governed by scattering within the exciton manifold. It was concluded from pump-probe measurements that energy transfer was favored between exciton levels that lie within an energy spacing of 10 nm (120 cm−1) (Vulto et al. 1997). At this energy, the density of acoustic phonons might be high, so that electron–phonon coupling might be the underlying mechanism of downward energy transfer. Pump-probe transients indicated a sequential relaxation of the exciton energy along a ladder of states, as was also seen in exciton simulations (Vulto et al. 1999, 1997; Buck et al. 1997; Iseri and Gülen 1999; Brüggemann and May 2004) (see Tables 9, ,10,10, ,11,11, ,12).12). Figure 4 shows a couple of examples of this type of decay. Only at very low temperatures, the dephasing might be governed by downward coherent exciton transfer. The origin of the disagreement between the dephasing times from both measurements are unclear but might have to do with the distinct experimental conditions tuning into different mechanisms underlying the energy transfer in the complex.
In a more elaborate study, Louwe and Aartsma (1997) decided to take another look at the possible coherent nature of exciton transport by studying the FMO complex at 1.4 K with accumulated photon echoes and transient absorption (see Table 13). Owing to the broad exciton levels, they probed several excitonic transitions at the same time resulting in traces with multiple time constants. At long wavelengths, (815–830 nm) processes with exciton decay times of 5, 30, 110, and 385 ps were found, while at shorter wavelengths (795 nm), the decay was in the order of 100 fs. The FMO complex was described as pigments coupled to a TLS (double-well potential, in between which tunneling occurs), which is, in turn, coupled to a phonon bath. Fluctuations in the interactions between pigments due to transitions in the TLS is the main dephasing pathway in glasses below 10 K. The TLS transitions can both influence the dipole interactions between the pigments (low frequency transitions in TLS corresponding to large displacements in the protein) as well as the site energies (high frequency, smaller displacement). At low temperatures, the coherent energy transfer is mainly limited by this coupling. Above 10 K, the contribution of the TLS tunneling is of minor significance to the dephasing mechanism that are dominated by other processes. With these measurements, the earlier results from a preliminary study by Louwe and Aartsma (1994) were confirmed.
Several years later, interesting features were seen in low-temperature two-photon-echo (2PE) signals of both Chlorobium tepidum and Prosthecochloris aestuarii (Prokhorenko et al. 2002). At 1.27 K, the 2PE signals show oscillations that increase in intensity when the excitation is tuned to the red edge of the absorption spectrum (up to 40% of the total amplitude for excitation at 832 nm). These oscillations last up to 300 ps and are ascribed to vibrational states of the BChl a molecule in the ground state. Fourier transforms of the 2PE traces show that the obtained frequencies match those from previous studies (Savikhin et al. 1997). In the same study, it was shown that the general theory to describe the results of photon-echo experiments did not account for the current results. The typical δ shape for dynamics in the Markov limit at initial time delays was not observed. Therefore, the dynamics were described beyond the Markov limit where system–bath memory effects occur which, among others, result in the delayed growing in of coherence in the system. At that time, it was unclear whether this had a specific function in light harvesting.
Vulto et al. used a similar approach as was used previously by Louwe et al. in the simulation of the static spectra (see “Exciton nature of the BChl a excitations in the FMO protein” and “Coupling strengths, linewidth and exciton energies”); however, to introduce dynamics, coupling of the electronic excitations to the vibrational modes in the system was included (Vulto et al. 1999). Homogeneous broadening within the system was not incorporated in the model. Owing to the weak coupling, the exciton-vibrational coupling can be treated as a perturbative term in the Hamiltonian. The distinct time scales of dephasing and the bath correlation times leads to the decoupling of the dynamics of population and coherence. The excited state dynamics, therefore, is governed by population relaxation. Similarly, in the simulations of Renger and May, the frequency-dependent coupling of the electronic states in the systems to the surroundings is needed. In order to describe this, the phonon-side band in a fluorescence spectrum is fitted. Using this analytical description for the spectral density, the time-resolved spectra can be fitted. As was shown before, the exciton relaxation occurs mainly between adjacent levels. The number of states lower in energy determine the relaxation rate of an exciton level. However, important additional factors are also the energy difference between the two levels and the overlap between the excitation probability densities on a single pigment j (i.e., |Cα(j)|2|Cβ(j)|2). The authors noted that the spectra of Chlorobium tepidum fitted remarkably better than those of Prosthecochloris aestuarii, in particular an experimental decay time of 1.7 ps was not reproduced. This could be partially overcome by adjusting the site energies of especially BChl a 1 and BChl a 4. The energetic order, of these pigments which are the main contributors to the second lowest exciton states (E2), seems of importance for the dynamics in the system. This was further tested by introducing inhomogeneous broadening in the system by a Monte Carlo simulation of the spectra and the dynamics. In addition to the decay time constants, distributions of time constants centered around the originally simulated values were found. At the exciton level E2, this distribution showed a clear distinction between two time domains; one of several hundreds of femtoseconds and another of several picoseconds, the latter is in the same order as the experimentally observed time scale. The spectra resulting from the Monte Carlo simulations are very similar to the dressed stick spectra calculated earlier (Vulto et al. 1998a). Vulto et al. showed that the method of Renger et al. does not reproduce the T − S and LD spectra at all, and concluded that their description of the electronic structure of the FMO complex was not completely correct. However, the ingenious way of describing the spectral broadening of the transitions by Renger et al. could be used to improve future simulations.
The decay time for energy transfer from the lowest exciton state to the ground state varies widely between different techniques and research groups. Table 14 gives a clear indication that there are two timescales concerned with the lowest exciton lifetime; one of about 100 ps and a longer one of several ns. A more elaborate description of this lifetime for Chlorobium tepidum is found in the electronic supplementary material. The discussion therein indicates that the lifetime of the lowest exciton state is influenced by the preparation method of the samples and in particular by the addition of oxidizing or reducing agents. However, a satisfying explanation for the latter is lacking to date.
In the last 5 years an additional technique was used to study exciton dynamics in the FMO complex: 2D spectroscopy. This technique directly shows the frequency correlation between excited states. When there is coupling between the different states, as is the case in the FMO complex, excitation of one state influences the others. 2D electronic spectroscopy on the FMO complex is mainly used to elucidate the time-dependent couplings between exciton states. This does not provide a direct way of measuring the site energies of the individual pigments. However, in 2D electronic spectroscopy, the coupling between the exciton states will appear in the spectra directly as the so-called cross peaks (Brixner et al. 2005). In the FMO complex, the cross peaks in the 2D spectra overlap with broad and strong diagonal peaks, due to the high spectral density. In order to overcome this problem, a technique in which the diagonal peaks can be eliminated and the cross peaks are brought out was developed (Read et al. 2007). The technique is based on a scheme known from 2D-vibrational spectroscopy and uses polarization of the first two pulses to select the cross peaks. Since most of this study has been done using Chlorobium tepidum, more on this topic can be found in the electronic supplementary material.
In order to extract the contributions of the various energy decay processes in a congested 2D spectrum, polarization-dependent 2D spectroscopy was used (Read et al. 2008). In contrast to the previous study (Read et al. 2007) this was a measurement of both the rephasing and non-rephasing spectra. In the non-rephasing spectra, the diagonal linewidths of the exciton transitions are narrower and, therefore, a higher resolution can be obtained. Furthermore, the authors made use of two polarization combinations for separate 2D experiments. Theoretically, it is possible to obtain the projection angle ϕ between a pair of exciton states from the ratio between these two polarization combinations. In the nonrephasing spectra, a strong cross peak at 804 and 814 nm appears while changing the polarization from one to the other polarization combination. By calculating the amplitude factor of the cross peaks depending on ϕ for the two polarization cases, it was shown that an angle of 40° reproduced the measured 2D data. This implies that without previous knowledge about structural properties of the system, a tentative view of the orientation of transition dipoles can be obtained. The current models of the FMO complex predict that excitons 2 and 4 have a high dipole strength and are the main contributors to the peaks in the spectrum at 804 and 814 nm. Furthermore, the corresponding exciton states have a projection angle of 38° verifying the results from the polarization dependent 2D experiments. Using the same idea of polarized fields in a theoretical study, contributions of coherent evolution and incoherent energy relaxation to a 2D spectrum could be separated due to a specific choice of the polarizations of the incoming pulses (Abramavicius et al. 2008b).
Currently, the best simulations of exciton dynamics are based on a method initiated by Vulto et al. (1999). An important parameter in their simulations is the coupling of an exciton state to a phonon bath. This vibronic coupling can account for energy relaxation in the FMO complex and is therefore an important factor in simulations of the exciton dynamics. In order to model the phonon-side band that mediates the coupling, they used an empirical approximation. The electron–phonon coupling was set to be equal for all states. Results of their simulations were that the exciton states preferably decay stepwise downhill along an energy gradient, as energy transfer mainly occurs between two adjacent levels. The rate of relaxation can be enhanced by the high value of the (linear) electron–phonon coupling. Cho et al. (2005) also showed that the rate of exciton transfer depends on the amplitude of the spectral density at the frequency of the transition. Using the coupling constants between the BChls of Vulto et al., except for a reduced coupling between BChl a 5 and 6, the exciton dynamics were simulated using a modified Förster/Redfield theory. Rates calculated using conventional Redfield theory turned out to be too slow in the presence of weakly coupled pigments. Therefore, the weak couplings are not taken into account into the diagonalization of the Hamiltonian, but are used to calculate the rate matrix using Förster theory. Simulations of 2D electronic spectra showed a better agreement with the experiment when the modified theory was used. Adolphs et al. use an elaborate model for the spectral density by also taking into account vibrational sidebands (Adolphs and Renger 2006). In order to simulate exciton relaxation, Redfield theory was compared to the more elaborate modified theory. The latter assumed that there are possible nuclear rearrangement effects that accompany exciton relaxation. Only minor differences between the two methods were observed, where modified Redfield theory predicts slightly lower rates. Two interesting observations from their simulations are that the spectral density of the electron–phonon coupling seems optimized to dissipate excess energy during relaxation. Also, simulations revealed two different exciton relaxation branches, a slow and a fast one, which are used for energy transfer from the chlorosomes to the RC.
As the exciton dynamics in the FMO complex is well studied and understood, a possible next step is to try and influence this dynamics. Following that route, specially shaped laser pulses are calculated using an optimal control algorithm and applied to the FMO complex (Brüggemann and May 2004). These pulses lead to a superposition of excitonic states, an excitonic wavepacket, with the target to populate just a single chromophore at a given time. The theoretical framework is given by the multi-exciton density matrix, and although the dissipation is damping the wavepacket at low temperatures, the target can be reached quite well. In a follow-up article, the additional effects of inhomogeneous broadening and orientational averaging were included (Brüggemann et al. 2006). Again, the target could be reached although to a lesser extend. The introduction of a laser field, shaped in both polarization directions, led to a larger target state population, partially working against the energetic and oriental averaging.
Under conditions encountered by the FMO complex in vivo it is very likely that multiple excitations occur within one complex. These double-excited states are more complicated than its single counterpart and are less well studied. Often 2D spectra are obscured by overlapping contributions of single and double exciton resonances. By looking at a smart representation of the 2D spectra using a particular set of pulses, the correlated dynamics of the double excited states can be probed (Abramavicius et al. 2008a). Strong peaks are observed for double exciton states 1, 7, and 18 that also happen to be the most delocalized states in the system. In addition, weaker signals of exciton states 9, 16, and 17 are observed. Instead of calculating the wavefunctions of the different exciton states, an alternative method can be used to describe the behavior of excitons in aggregates. In the quasiparticle approach, all the properties of the system are described in terms of scattering and double exciton energies are simply given by a sum of single exciton energies. Comparing the spectra resulting from the full calculation with that of the quasiparticle approach shows that the energies at which the peaks appear in the spectra agree, while the fine structure in the spectra of the quasiparticle approach is distorted. In order to approximate the spectra, the quasiparticle approach can be used, however, because the exciton coupling is strong, which is neglected in this approach, and the nonbosonic nature of the excitons a full calculation of the spectra is necessary for detailed analysis.
New types of 2D techniques can be developed by introducing pulse polarizations as variables into standard 2D schemes, as described in the previous section. This, amongst others enables the dissection of the congested 2D spectra into incoherent and coherent contributions and provides interesting perspective for new control strategies (Abramavicius et al. 2008b; Voronine et al. 2008).
Slowly the choice of parameters used to simulate the results obtained from various optical techniques is converging. Around 15–20 years ago, little agreement existed on whether to include the whole trimer or just the monomer in simulations of the spectra, on the values of the site energies and how to obtain them, and on the location of the lowest energy pigment. Nowadays, these issues seem more or less resolved: Only the monomer is taken into account in simulations, as is inhomogeneous broadening due to structural changes, BChl a 3 is principally assigned to have the lowest site energy. The parameter set from Louwe et al., including the site energies, is widely used in increasingly complex simulations. The latest addition to this is a new approach to calculate site energies instead of fitting them, using amongst others quantum chemical methods. The possible influence of the recently proposed eighth BChl a molecule on the variety of optical spectra could invoke new studies. It is conceivable that new detailed simulations including this pigment can lift the remaining discrepancies between experimental and theoretical spectra.
While the exact energy transfer timescales within the exciton manifold vary between techniques, it is commonly agreed that decay to the lowest exciton state occurs within several picoseconds. Despite this rapid decay, an interesting observation is the prolonged presence of coherence in the complex. This coherence with its potential role in mediating efficient energy transfer, is the topic of current research using advanced techniques such as 2D electronic spectroscopy and coherent control strategies with shaped excitation pulses.
Below is the link to the electronic supplementary material.
This study is part of the research program of the Stichting voor Fundamenteel Onderzoek der Materie (FOM), which is supported financially by the Nederlandse organisatie voor Wetenschappelijk Onderzoek (NWO).
Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
Maaike T. W. Milder, Email: ln.floma@redlim.
Jennifer L. Herek, Phone: +31-53-4893805, Fax: +31-53-4893511, Email: firstname.lastname@example.org.