Site Selective and Single Complex Laser-Based Spectroscopies: A Window on Excited State Electronic Structure, Excitation Energy Transfer, and Electron–Phonon Coupling of Selected Photosynthetic Complexes

PigmentHunter: A point-and-click application for automated chlorophyll-protein simulations

Chlorophyll proteins (CPs) are the workhorses of biological photosynthesis, working together to absorb solar energy, transfer it to chemically active reaction centers, and control the charge-separation process that drives its storage as chemical energy. Yet predicting CP optical and electronic properties remains a serious challenge, driven by the computational difficulty of treating large, electronically coupled molecular pigments embedded in a dynamically structured protein environment. To address this challenge, we introduce here an analysis tool called PigmentHunter, which automates the process of preparing CP structures for molecular dynamics (MD), running short MD simulations on the nanoHUB.org science gateway, and then using electrostatic and steric analysis routines to predict optical absorption, fluorescence, and circular dichroism spectra within a Frenkel exciton model. Inter-pigment couplings are evaluated using point-dipole or transition-charge coupling models, while site energies can be estimated using both electrostatic and ring-deformation approaches. The package is built in a Jupyter Notebook environment, with a point-and-click interface that can be used either to manually prepare individual structures or to batch-process many structures at once. We illustrate PigmentHunter's capabilities with example simulations on spectral line shapes in the light harvesting 2 complex, site energies in the Fenna-Matthews-Olson protein, and ring deformation in photosystems I and II.

High-Resolution Frequency-Domain Spectroscopic and Modeling Studies of Photosystem I (PSI), PSI Mutants and PSI Supercomplexes

Photosystem I (PSI) is one of the two main pigment-protein complexes where the primary steps of oxygenic photosynthesis take place. This review describes low-temperature frequency-domain experiments (absorption, emission, circular dichroism, resonant and non-resonant hole-burned spectra) and modeling efforts reported for PSI in recent years. In particular, we focus on the spectral hole-burning studies, which are not as common in photosynthesis research as the time-domain spectroscopies. Experimental and modeling data obtained for trimeric cyanobacterial Photosystem I (PSI3), PSI3 mutants, and PSI3-IsiA18 supercomplexes are analyzed to provide a more comprehensive understanding of their excitonic structure and excitation energy transfer (EET) processes. Detailed information on the excitonic structure of photosynthetic complexes is essential to determine the structure-function relationship. We will focus on the so-called "red antenna states" of cyanobacterial PSI, as these states play an important role in photochemical processes and EET pathways. The high-resolution data and modeling studies presented here provide additional information on the energetics of the lowest energy states and their chlorophyll (Chl) compositions, as well as the EET pathways and how they are altered by mutations. We present evidence that the low-energy traps observed in PSI are excitonically coupled states with significant charge-transfer (CT) character. The analysis presented for various optical spectra of PSI3 and PSI3-IsiA18 supercomplexes allowed us to make inferences about EET from the IsiA18 ring to the PSI3 core and demonstrate that the number of entry points varies between sample preparations studied by different groups. In our most recent samples, there most likely are three entry points for EET from the IsiA18 ring per the PSI core monomer, with two of these entry points likely being located next to each other. Therefore, there are nine entry points from the IsiA18 ring to the PSI3 trimer. We anticipate that the data discussed below will stimulate further research in this area, providing even more insight into the structure-based models of these important cyanobacterial photosystems.

Empirical asymmetric phonon sideband due to phonons in the protein matrix present in photosynthetic complexes: Time-domain response theory

The phonon spectral density plays a key role in probing the dynamical and spectral behavior of molecular aggregates. One may utilize the intimate connection between the one-phonon profile and the phonon spectral density to extract a plausible form of the spectral density of media with rich structure using advanced optical spectroscopy. The excitonic transition is normally accompanied by a broad, asymmetric phonon-side band due to the coupling to the phonons in the surrounding protein matrix present in photosynthetic complexes. The asymmetry in the one-phonon profile of a homogeneous absorption spectrum and other experiments performed on photosynthetic bacterial reaction centers (BRCs) led the Small group to employ a half-Gaussian distribution function on the red side and half-Lorentzian distribution function on the blue side of the absorption lineshape to account for the one-phonon profile asymmetrical shape and relaxation effects contributing to spectroscopy and dynamics of BRCs at hand. Different research groups successfully employed the theory of Small to simulate their photosynthetic spectral data so they could calculate the homogeneous absorption and hole-burned spectra of photosynthetic complexes. Although this report does not directly use the formulae of homogeneous absorption, hole-burning, and fluorescence line-narrowed spectra of BRCs, and photosynthetic complexes, developed by Hayes-Small, it builds on their idea of the phonon sideband asymmetric shape in deriving an accurate and computationally efficient linear electronic transition dipole moment time correlation function. Besides the compelling tractability and efficiency of this correlation function, it accounts for excitonic coupling and eliminates all the inconsistencies arising in the Hayes-Small theory.

Chlorophyll triplet states in thylakoid membranes of Acaryochloris marina. Evidence for a triplet state sitting on the photosystem I primary donor populated by intersystem crossing

Photo-induced triplet states in the thylakoid membranes isolated from the cyanobacterium Acaryocholoris marina, that harbours Chlorophyll (Chl) d as its main chromophore, have been investigated by Optically Detected Magnetic Resonance (ODMR) and time-resolved Electron Paramagnetic Resonance (TR-EPR). Thylakoids were subjected to treatments aimed at poising the redox state of the terminal electron transfer acceptors and donors of Photosystem II (PSII) and Photosystem I (PSI), respectively. Under ambient redox conditions, four Chl d triplet populations were detectable, identifiable by their characteristic zero field splitting parameters, after deconvolution of the Fluorescence Detected Magnetic Resonance (FDMR) spectra. Illumination in the presence of the redox mediator N,N,N',N'-Tetramethyl-p-phenylenediamine (TMPD) and sodium ascorbate at room temperature led to a redistribution of the triplet populations, with T3 (|D|= 0.0245 cm-1, |E|= 0.0042 cm-1) becoming dominant and increasing in intensity with respect to untreated samples. A second triplet population (T4, |D|= 0.0248 cm-1, |E|= 0.0040 cm-1) having an intensity ratio of about 1:4 with respect to T3 was also detectable after illumination in the presence of TMPD and ascorbate. The microwave-induced Triplet-minus-Singlet spectrum acquired at the maximum of the |D|-|E| transition (610 MHz) displays a broad minimum at 740 nm, accompanied by a set of complex spectral features that overall resemble, despite showing further fine spectral structure, the previously reported Triplet-minus-Singlet spectrum attributed to the recombination triplet of PSI reaction centre,3 P 740 [Schenderlein M, Çetin M, Barber J, et al. Spectroscopic studies of the chlorophyll d containing photosystem I from the cyanobacterium Acaryochloris marina. Biochim Biophys Acta 1777:1400-1408]. However, TR-EPR experiments indicate that this triplet displays an eaeaea electron spin polarisation pattern which is characteristic of triplet sublevels populated by intersystem crossing rather than recombination, for which an aeeaae polarisation pattern is expected instead. It is proposed that the observed triplet, which leads to the bleaching of the P740 singlet state, sits on the PSI reaction centre.

Identification of Residues Potentially Involved in Optical Shifts in the Water-Soluble Chlorophyll a-Binding Protein through Molecular Dynamics Simulations

Reversible light and thermally induced spectral shifts are universally observed in a wide variety of pigment-protein complexes at temperatures ranging from cryogenic to ambient. In this paper, we employed large-scale molecular dynamics (MD) simulations of a prototypical pigment-protein complex to better understand these shifts at a molecular scale. Although multiple mechanisms have been proposed over the years, no verification of these proposals via MD simulations has thus far been performed; our work represents the first step in this direction. From simulations of the water-soluble chlorophyll-binding protein complex, we determined that rearrangements of long hydrogen bonds were unlikely to be the origin of the multiwell landscape features necessary to explain observed spectral shifts. We also assessed small motions of amino acid residues and identified side chain rotations of some of these residues as likely candidates for the origin of relevant multiwell landscape features. The protein free-energy landscapes associated with side chain rotations feature energy barriers of around 1100-1600 cm-1, in agreement with optical spectroscopy results, with the most promising residue type associated with experimental signatures being serine, which possesses a symmetric triple-well landscape and moment of inertia of a relevant magnitude.

Homogeneous Dephasing in Photosynthetic Bacterial Reaction Centers: Time Correlation Function Approach

A new tractable linear electronic transition dipole moment time correlation function (ETDMTCF) that accurately accounts for electronic dephasing, asymmetry, and width of 1-phonon profile, which the zero-phonon line (ZPL) contributes to it, in Rhodopseudomonas viridis bacterial reaction center is derived. This time correlation function proves to be superior to other frequency-domain expressions in case of strong electron-phonon coupling (which is often the case in bacterial RCs and pigment-protein complexes), many vibrational modes involved, and high temperature, whereby more vibronic and electronic (sequence) transitions would arise. The Fourier transform of this ETDMTCF leads to asymmetric multiphonon profiles composed of Lorentzian distribution and Gaussian distribution on the high- and low-energy sides, respectively, whereby the overtone widths fold themselves with that of the one-phonon profile. This ETDMTCF also features expedient computation in large systems using asymmetric phonon profiles to account correctly for dephasing and pigment-protein interaction (electron-phonon coupling). The derived ETDMTCF allows computing all nonlinear optical signals in both time and frequency domains, through the nonlinear dipole moment time correlation functions (as guided by nonlinear optical response theory) in line with the eight Liouville space pathways. The linear transition dipole moment time correlation function is of a central value as the nonlinear transition dipole moment time correlation function is expressed in terms of the linear transition dipole moment time correlation function, derived herein. One of the great advantages of presenting this ETDMTCF is its applicability to nonlinear transition dipole moment time correlation functions in line with the eight Liouville space pathways needed in computing nonlinear signals. As such, there is more to the utility and applicability of the presented ETDMTCF besides computational expediency and efficiency. Results show good agreement with the reported literature. The intimate connection between a one-phonon profile and the corresponding bath spectral density in photosynthetic complexes is discussed.

Bioexcitons by Design: How Do We Get There?

Biological pigment-protein complexes (PPCs) exhibit a remarkable ability to tune the optical properties of biological excitons (bioexcitons) through specific pigment-protein interactions. While such fine-tuning allows natural systems (e.g., photosynthetic proteins) to carry out their native functions with near-optimal performance, native function itself is often suboptimal for applications such as biofuel production or quantum technology development. This perspective offers a look at near-term prospects for the rational reoptimization of PPC bioexcitons for new functions using site-directed mutagenesis. The primary focus is on the "structure-spectrum" challenge of understanding the relationships between structural features and spectroscopic properties. While recent examples demonstrate that site-directed mutagenesis can be used to tune nearly all key bioexciton parameters (e.g., site energies, interpigment couplings, and electronic-vibrational interactions), critical challenges remain before we achieve truly rational design of bioexciton properties.

Modulation of Membrane Microviscosity by Protein-Mediated Carotenoid Delivery as Revealed by Time-Resolved Fluorescence Anisotropy

Carotenoids are potent antioxidants with a wide range of biomedical applications. However, their delivery into human cells is challenging and relatively inefficient. While the use of natural water-soluble carotenoproteins capable to reversibly bind carotenoids and transfer them into membranes is promising, the quantitative estimation of the delivery remains unclear. In the present work, we studied echinenone (ECN) delivery by cyanobacterial carotenoprotein AnaCTDH (C-terminal domain homolog of the Orange Carotenoid Protein from Anabaena), into liposome membranes labelled with BODIPY fluorescent probe. We observed that addition of AnaCTDH-ECN to liposomes led to the significant changes in the fast-kinetic component of the fluorescence decay curve, pointing on the dipole-dipole interactions between the probe and ECN within the membrane. It may serve as an indirect evidence of ECN delivery into membrane. To study the delivery in detail, we carried out molecular dynamics modeling of the localization of ECN within the lipid bilayer and calculate its orientation factor. Next, we exploited FRET to assess concentration of ECN delivered by AnaCTDH. Finally, we used time-resolved fluorescence anisotropy to assess changes in microviscosity of liposomal membranes. Incorporation of liposomes with β-carotene increased membrane microviscosity while the effect of astaxanthin and its mono- and diester forms was less pronounced. At temperatures below 30 °C addition of AnaCTDH-ECN increased membrane microviscosity in a concentration-dependent manner, supporting the protein-mediated carotenoid delivery mechanism. Combining all data, we propose FRET-based analysis and assessment of membrane microviscosity as potent approaches to characterize the efficiency of carotenoids delivery into membranes.

Frequency-Domain Spectroscopic Study of the Photosystem I Supercomplexes, Isolated IsiA Monomers, and the Intact IsiA Ring

The PSI₃-IsiA₁₈ supercomplex is one of the largest and most complicated assemblies in photosynthesis. The IsiA ring, composed of 18 IsiA monomers (IsiA₁₈) surrounding the PSI trimer (PSI₃), forms under iron-deficient conditions in cyanobacteria and acts as a peripheral antenna. Based on the supercomplex structure recently determined via cryo-EM imaging, we model various optical spectra of the IsiA monomers and IsiA₁₈ ring. Comparison of the absorption and emission spectra of the isolated IsiA monomers and the full ring reveals that about 2.7 chlorophylls (Chls) are lost in the isolated IsiA monomers. The best fits for isolated monomers spectra are obtained assuming the absence of Chl 508 and Chl 517 and 70% loss of Chl 511. The best model describing all three hexamers and the entire ring suggests that the lowest energy pigments are Chls 511, 514, and 517. Based on the modeling results presented in this work, we conclude that there are most likely three entry points for EET from the IsiA₆ hexamer to the PSI core monomer, with two of these entry points likely being located next to each other (i.e., nine entry points from IsiA₁₈ to the PSI₃ trimer). Finally, we show that excitation energy transfer inside individual monomers is fast (<2 ps at T = 5 K) and at least 20 times faster than intermonomer energy transfer.

Light-Harvesting Complex II Adopts Different Quaternary Structures in Solution as Observed Using Small-Angle Scattering

The high-resolution crystal structure of the trimeric major light-harvesting complex of photosystem II (LHCII) is often perceived as the basis for understanding its light-harvesting and photoprotective functions. However, the LHCII solution structure and its oligomerization or aggregation state may generally differ from the crystal structure and, moreover, also depend on its functional state. In this regard, small-angle scattering experiments provide the missing link by offering structural information in aqueous solution at physiological temperatures. Herein, we use small-angle scattering to investigate the solution structures of two different preparations of solubilized LHCII employing the nonionic detergents n-octyl-β-d-glucoside (OG) and n-dodecyl-β-D-maltoside (β-DM). The data reveal that the LHCII-OG complex is equivalent to the trimeric crystal structure. Remarkably, however, we observe─for the first time─a stable oligomer composed of three LHCII trimers in the case of the LHCII-β-DM preparation, implying additional pigment-pigment interactions. The latter complex is assumed to mimic trimer-trimer interactions which play an important role in the context of photoprotective nonphotochemical quenching.

Molecular Photocatalytic Water Splitting by Mimicking Photosystems I and II

In nature, water is oxidized by plastoquinone to evolve O₂ and form plastoquinol in Photosystem II (PSII), whereas NADP⁺ is reduced by plastoquinol to produce NADPH and regenerate plastoquinone in Photosystem I (PSI), using homogeneous molecular photocatalysts. However, water splitting to evolve H₂ and O₂ in a 2:1 stoichiometric ratio has yet to be achieved using homogeneous molecular photocatalysts, remaining as one of the biggest challenges in science. Herein, we demonstrate overall water splitting to evolve H₂ and O₂ in a 2:1 ratio using a two liquid membranes system composed of two toluene phases, which are separated by a solvent mixture of water and trifluoroethanol (H₂O/TFE, 3:1 v/v), with a glass membrane to combine PSI and PSII molecular models. A PSII model contains plastoquinone analogs [p-benzoquinone derivatives (X-Q)] in toluene and an iron(II) complex as a molecular oxidation catalyst in H₂O/TFE (3:1 v/v), which evolves a stoichiometric amount of O₂ and forms plastoquinol analogs (X-QH₂) under photoirradiation. On the other hand, a PSI model contains nothing in toluene but contains X-QH₂, 9-mesityl-10-methylacridinium ion (Acr⁺-Mes) as a photocatalyst, and a cobalt(III) complex as an H₂ evolution catalyst in H₂O/TFE (3:1 v/v), which evolves a stoichiometric amount of H₂ and forms X-Q under photoirradiation. When a PSII model system is combined with a PSI model system with two glass membranes and two liquid membranes, photocatalytic water splitting with homogeneous molecular photocatalysts is achieved to evolve hydrogen and oxygen with the turnover number (TON) of >100.

Accurate prediction of mutation-induced frequency shifts in chlorophyll proteins with a simple electrostatic model

Photosynthetic pigment-protein complexes control local chlorophyll (Chl) transition frequencies through a variety of electrostatic and steric forces. Site-directed mutations can modify this local spectroscopic tuning, providing critical insight into native photosynthetic functions and offering the tantalizing prospect of creating rationally designed Chl proteins with customized optical properties. Unfortunately, at present, no proven methods exist for reliably predicting mutation-induced frequency shifts in advance, limiting the method's utility for quantitative applications. Here, we address this challenge by constructing a series of point mutants in the water-soluble chlorophyll protein of Lepidium virginicum and using them to test the reliability of a simple computational protocol for mutation-induced site energy shifts. The protocol uses molecular dynamics to prepare mutant protein structures and the charge density coupling model of Adolphs et al. [Photosynth. Res. 95, 197-209 (2008)] for site energy prediction; a graphical interface that implements the protocol automatically is published online at http://nanohub.org/tools/pigmenthunter. With the exception of a single outlier (presumably due to unexpected structural changes), we find that the calculated frequency shifts match the experiment remarkably well, with an average error of 1.6 nm over a 9 nm spread in wavelengths. We anticipate that the accuracy of the method can be improved in the future with more advanced sampling of mutant protein structures.

Vibration-mediated energy transport in bacterial reaction center: Simulation study

Exciton energy relaxation in a bacterial Reaction Center (bRC) pigment-protein aggregate presumably involves emission of high energy vibrational quanta to cover wide energy gaps between excitons. Here, we assess this hypothesis utilizing vibronic two-particle theory in modeling of the excitation relaxation process in bRC. Specific high frequency molecular vibrational modes are included explicitly one at a time in order to check which high frequency vibrations are involved in the excitation relaxation process. The low frequency bath modes are treated perturbatively within Redfield relaxation theory. The analysis of the population relaxation rate data indicates energy flow pathways in bRC and suggests that specific vibrations may be responsible for the excitation relaxation process.

Simulation of Nonlinear Femtosecond Signals at Finite Temperature via a Thermo Field Dynamics-Tensor Train Method: General Theory and Application to Time- and Frequency-Resolved Fluorescence of the Fenna-Matthews-Olson Complex

Addressing needs of contemporary nonlinear femtosecond optical spectroscopy, we have developed a fully quantum, numerically accurate wave function-based approach for the calculation of third-order spectroscopic signals of polyatomic molecules and molecular aggregates at finite temperature. The systems are described by multimode nonadiabatic vibronic-coupling Hamiltonians, in which diagonal terms are treated in harmonic approximation, while off-diagonal interstate couplings are assumed to be coordinate independent. The approach is based on the Thermo Field Dynamics (TFD) representation of quantum mechanics and tensor-train (TT) machinery for efficient numerical simulation of quantum evolution of systems with many degrees of freedom. The developed TFD-TT approach is applied to the calculation of time- and frequency-resolved fluorescence spectra of the Fenna-Matthews-Olson (FMO) antenna complex at room temperature taking into account finite time-frequency resolution in fluorescence detection, orientational averaging, and static disorder.

Role of Spatially Correlated Fluctuations in Photosynthetic Excitation Energy Transfer with an Equilibrium and a Nonequilibrium Initial Bath

The transfer of excitation energy in photosynthetic light-harvesting complexes has inspired growing interest for its scientific and engineering significance. Recent experimental findings have suggested that spatially correlated environmental fluctuations may account for the existence of long-lived quantum coherent energy transfer observed even at physiological temperature. In this paper, we investigate the effects of spatial correlations on the excitation energy transfer dynamics by including a nonequilibrium initial bath in a simulated donor-acceptor model. The initial bath state, which is assumed to be either equilibrium or nonequilibrium, is expanded in powers of coupling strength within the polaron formalism of a quantum master equation. The spatial correlations of bath fluctuations strongly influence the decay of coherence in the dynamics. The role of a nonequilibrium initial bath is also influenced by spatial correlations and becomes the most conspicuous for certain degrees of spatial correlations from which we propose a picture that the spatial correlations of bath fluctuations open up new energy transfer pathways, playing a role of protecting coherence. Besides, we apply the polaron master equation approach to study the dynamics in a two-site subsystem of the FMO complex and provide a practical example that shows the versatility of this approach.

Effects of Chlorophyll Triplet States on the Kinetics of Spectral Hole Growth

Spectral hole burning has been employed for decades to study various amorphous solids and proteins. Triplet states and respective transient holes were incorporated into theoretical models and software simulating nonphotochemical spectral hole burning (NPHB) and including all relevant distributions, in particular the distribution of the angle between the electric field of light E and transient dipole moment of the chromophore μ. The presence of a chlorophyll a triplet state with a lifetime of several milliseconds explains the slowdown of NPHB (on the depth vs illumination dose scale) with the increase of the light intensity, as well as larger hole depths observed in weak probe beam experiments, compared to those deduced from the hole growth kinetics (HGK) measurements (signal collected at a fixed wavelength while a stronger burning beam is on) in cytochrome b₆f and chemically modified LH2. We also considered the solvent deuteration effects on triplet lifetime and concluded that both triplet states and local heating likely play a role in slowing down the HGK with increasing burn intensity.

Exciton Lifetime Distributions and Population Dynamics in the FMO Protein Complex from Prosthecochloris aestuarii

Significant protein rearrangement upon excitation and energy transfer in Fenna-Matthews-Olson protein of Prosthecochloris aestuarii results in a modified energy landscape, which induces more changes in pigment site energies than predicted by the "standard" hole-burning theory. The energy changes are elucidated by simulations while investigating the effects of site-dependent disorder, both static (site-energy distribution widths) and dynamic (spectral density shapes). The resulting optimized site energies and their fluctuations are consistent with relative differences observed in inhomogeneous widths calculated by recent molecular dynamic simulations. Two sets of different spectral densities reveal how their shapes affect the population dynamics and distribution of exciton lifetimes. Calculations revealed the wavelength-dependent distributions of exciton lifetimes (T ₁) in the femtosecond to picosecond time frame. We suggest that the calculated multimodal and asymmetric wavelength-dependent T ₁ distributions offer more insight into the interpretation of resonant hole-burned (HB) spectra, kinetic traces in two-dimensional (2D) electronic spectroscopy experiments, and widely used global analyses in fitting data from transient absorption experiments.

On wavelength-dependent exciton lifetime distributions in reconstituted CP29 antenna of the photosystem II and its site-directed mutants

To provide more insight into the excitonic structure and exciton lifetimes of the wild type (WT) CP29 complex of photosystem II, we measured high-resolution (low temperature) absorption, emission, and hole burned spectra for the A2 and B3 mutants, which lack chlorophylls a612 and b614 (Chls), respectively. Experimental and modeling results obtained for the WT CP29 and A2/B3 mutants provide new insight on the mutation-induced changes at the molecular level and shed more light on energy transfer dynamics. Simulations of the A2 and B3 optical spectra, using the second-order non-Markovian theory, and comparison with improved fits of WT CP29 optical spectra provide more insight into their excitonic structure, mutation induced changes, and frequency-dependent distributions of exciton lifetimes (T₁). A new Hamiltonian obtained for WT CP29 reveals that deletion of Chls a612 or b614 induces changes in the site energies of all remaining Chls. Hamiltonians obtained for A2 and B3 mutants are discussed in the context of the energy landscape of chlorophylls, excitonic structure, and transfer kinetics. Our data suggest that the lowest exciton states in A2 and B3 mutants are contributed by a611(57%), a610(17%), a615(15%) and a615(58%), a611(20%), a612(15%) Chls, respectively, although other compositions of lowest energy states are also discussed. Finally, we argue that the calculated exciton decay times are consistent with both the hole-burning and recent transient absorption measurements. Wavelength-dependent T₁ distributions offer more insight into the interpretation of kinetic traces commonly described by discrete exponentials in global analysis/global fitting of transient absorption experiments.

Improvement of luminescence properties of NaYF4 :Yb3+ /Er3+ upconversion materials by a cross-relaxation mechanism based on co-doped Ho3+ ion concentrations

NaYF₄ :Yb³⁺ /Er³⁺ /Ho³⁺ nanophosphors were successfully synthesized using a solvothermal method and with various concentrations of Ho³⁺ ions. The crystal structure, grain size, morphology, and luminescence properties were analyzed by X-ray diffraction, field-emission scanning electron microscopy, and photoluminescence measurements. All samples were crystallized as a cubic structure; it was confirmed that all samples exhibited strong green emission and weak red emission generated at a particular level of the activated ions. The strongest upconversion fluorescence intensity was observed in the Ho³⁺ and Er³⁺ ions co-doped NaYF₄ materials with a Ho³⁺ ion concentration of 0.005 mol. Only the green fluorescence intensity at the 542 nm centre increased strongly due to the ⁴ S_3/2 →⁴ I_15/2 energy transfer. This increase in upconversion fluorescence intensity at a selected wavelength was described as a cross-relaxation mechanism due to the addition of Ho³⁺ ions.

Asymmetry in the Q y Fluorescence and Absorption Spectra of Chlorophyll a Pertaining to Exciton Dynamics

Significant asymmetry found between the high-resolution Q_y emission and absorption spectra of chlorophyll-a is herein explained, providing basic information needed to understand photosynthetic exciton transport and photochemical reactions. The Q_y spectral asymmetry in chlorophyll has previously been masked by interference in absorption from the nearby Q_x transition, but this effect has recently been removed using extensive quantum spectral simulations or else by analytical inversion of absorption and magnetic circular dichroism data, allowing high-resolution absorption information to be accurately determined from fluorescence-excitation spectra. To compliment this, here, we measure and thoroughly analyze the high-resolution differential fluorescence line narrowing spectra of chlorophyll-a in trimethylamine and in 1-propanol. The results show that vibrational frequencies often change little between absorption and emission, yet large changes in line intensities are found, this effect also being strongly solvent dependent. Among other effects, the analysis in terms of four basic patterns of Duschinsky-rotation matrix elements, obtained using CAM-B3LYP calculations, predicts that a chlorophyll-a molecule excited into a specific vibrational level, may, without phase loss or energy relaxation, reemit the light over a spectral bandwidth exceeding 1,000 cm⁻¹ (0.13 eV) to influence exciton-transport dynamics.

Static Disorder in Excitation Energies of the Fenna-Matthews-Olson Protein: Structure-Based Theory Meets Experiment

Inhomogeneous broadening of optical lines of the Fenna-Matthews-Olson (FMO) light-harvesting protein is investigated by combining a Monte Carlo sampling of low-energy conformational substates of the protein with a quantum chemical/electrostatic calculation of local transition energies (site energies) of the pigments. The good agreement between the optical spectra calculated for the inhomogeneous ensemble and the experimental data demonstrates that electrostatics is the dominant contributor to static disorder in site energies. Rotamers of polar amino acid side chains are found to cause bimodal distribution functions of site energy shifts, which can be probed by hole burning and single-molecule spectroscopy. When summing over the large number of contributions, the resulting distribution functions of the site energies become Gaussians, and the correlations in site energy fluctuations at different sites practically average to zero. These results demonstrate that static disorder in the FMO protein is in the realm of the central limit theorem of statistics.

Nature of low-energy exciton levels in light-harvesting complex II of green plants as revealed by satellite hole structure

Persistent non-photochemical hole burning at 4.2 K is an efficient experimental tool to unravel position and nature of low-energy excitonic states in pigment-protein complexes. This is demonstrated here for the case of the trimeric chlorophyll (Chl) a/b light-harvesting complexes of Photosystem II (LHC II) of green plants, where previous work (Pieper et al. J Phys Chem B 103:2412, 1999a) reported a highly localized lowest energy state at 680 nm. At that time, this finding appeared to be consistent with the contemporary knowledge about the LHC II structure, which mainly suggested the presence of weakly coupled Chl heterodimers. Currently, however, it is widely accepted that the lowest state is associated with an excitonically coupled trimer of Chl molecules at physiological temperatures. This raises the question, why an excitonically coupled state has not been identified by spectral hole burning. A re-inspection of the hole burning data reveals a remarkable dependence of satellite hole structure on burn fluence, which is indicative of the excitonic coupling of the low-energy states of trimeric LHC II. At low fluence, the satellite hole structure of the lowest/fluorescing  ~ 680 nm state is weak with only one shallow satellite hole at 649 nm in the Chl b spectral range. These findings suggest that the lowest energy state at  ~ 680 nm is essentially localized on a Chl a molecule, which may belong to a Chl a/b heterodimer. At high fluence, however, the lowest energy hole shifts blue to  ~ 677 nm and is accompanied by two satellite holes at  ~ 673 and 663 nm, respectively, indicating that this state is excitonically coupled to other Chl a molecules. In conclusion, LHC II seems to possess two different, but very closely spaced lowest energy states at cryogenic temperatures of 4.2 K.

Delocalization Effects in Chlorophyll Fluorescence: Nonperturbative Line Shape Analysis of a Vibronically Coupled Dimer

Non-adiabatic vibrational/electronic (vibronic) interactions in photosynthetic pigment/protein complexes (PPCs) have recently attracted considerable interest as a potential source for long-lived dynamic coherence and optimized light harvesting. The analysis of such effects is limited, however, by the complexity of the vibrational spectrum of biological pigments such as chlorophyll (Chl) molecules, which often makes numerical calculations prohibitively expensive and complicates the interpretation of experimental spectroscopic data. This work contributes to both challenges by using numerically exact computational methods to systematically examine vibronic mixing effects in the low-temperature fluorescence spectra of a Chl dimer possessing a full complement of local vibrations, using parameters extracted from experimental data. The results highlight the varying roles local vibrations can play in energy-transfer dynamics, both enhancing delocalization through vibronic resonance and, conversely, inducing dynamic localization by acting as a "self-bath" for local electronic transitions. In the specific context of line-narrowed fluorescence, the results indicate that, while low-frequency features are strongly suppressed by delocalization, high-frequency modes are likely to be dynamically localized in the parameter regime relevant to most photosynthetic complexes.

Photocatalytic Hydrogen Evolution from Plastoquinol Analogues as a Potential Functional Model of Photosystem I

The recent development of a functional model of photosystem II (PSII) has paved a new way to connect the PSII model with a functional model of photosystem I (PSI). However, PSI functional models have yet to be reported. We report herein the first potential functional model of PSI, in which plastoquinol (PQH₂) analogues were oxidized to plastoquinone (PQ) analogues, accompanied by hydrogen (H₂) evolution. Photoirradiation of a deaerated acetonitrile (MeCN) solution containing hydroquinone derivatives (X-QH₂) as a hydrogen source, 9-mesityl-10-methylacridinium ion (Acr⁺-Mes) as a photoredox catalyst, and a cobalt(III) complex, Co^III(dmgH)₂pyCl (dmgH = dimethylglyoximate monoanion; py = pyridine) as a redox catalyst resulted in the evolution of H₂ and formation of the corresponding p-benzoquinone derivatives (X-Q) quantitatively. The maximum quantum yield for photocatalytic H₂ evolution from tetrachlorohydroquinone (Cl₄QH₂) with Acr⁺-Mes and Co^III(dmgH)₂pyCl and H₂O in deaerated MeCN was determined to be 10%. Photocatalytic H₂ evolution is started by electron transfer (ET) from Cl₄QH₂ to the triplet ET state of Acr⁺-Mes to produce Cl₄QH₂^•+ and Acr^•-Mes with a rate constant of 7.2 × 10⁷ M^-1 s^-1, followed by ET from Acr^•-Mes to Co^III(dmgH)₂pyCl to produce [Co^II(dmgH)₂pyCl]^-, accompanied by the regeneration of Acr⁺-Mes. On the other hand, Cl₄QH₂^•+ is deprotonated to produce Cl₄QH^•, which transfers either a hydrogen-atom transfer or a proton-coupled electron transfer to [Co^II(dmgH)₂pyCl]^- to produce a cobalt(III) hydride complex, [Co^III(H)(dmgH)₂pyCl]^-, which reacts with H⁺ to evolve H₂, accompanied by the regeneration of Co^III(dmgH)₂pyCl. The formation of [Co^II(dmgH)₂pyCl]^- was detected by electron paramagnetic resonance measurements.

Analysis of Photosynthetic Systems and Their Applications with Mathematical and Computational Models

In biological and life science applications, photosynthesis is an important process that involves the absorption and transformation of sunlight into chemical energy. During the photosynthesis process, the light photons are captured by the green chlorophyll pigments in their photosynthetic antennae and further funneled to the reaction center. One of the most important light harvesting complexes that are highly important in the study of photosynthesis is the membrane-attached Fenna–Matthews–Olson (FMO) complex found in the green sulfur bacteria. In this review, we discuss the mathematical formulations and computational modeling of some of the light harvesting complexes including FMO. The most recent research developments in the photosynthetic light harvesting complexes are thoroughly discussed. The theoretical background related to the spectral density, quantum coherence and density functional theory has been elaborated. Furthermore, details about the transfer and excitation of energy in different sites of the FMO complex along with other vital photosynthetic light harvesting complexes have also been provided. Finally, we conclude this review by providing the current and potential applications in environmental science, energy, health and medicine, where such mathematical and computational studies of the photosynthesis and the light harvesting complexes can be readily integrated.

Establishment of the Qy Absorption Spectrum of Chlorophyll a Extending to Near-Infrared

A weak absorption tail related to the Q_y singlet electronic transition of solvated chlorophyll a is discovered using sensitive anti-Stokes fluorescence excitation spectroscopy. The quasi-exponentially decreasing tail was, at ambient temperature, readily observable as far as -2400 cm^-1 from the absorption peak and at relative intensity of 10^-7. The tail also weakened rapidly upon cooling the sample, implying its basic thermally activated nature. The shape of the spectrum as well as its temperature dependence were qualitatively well reproduced by quantum chemical calculations involving the pigment intramolecular vibrational modes, their overtones, and pairwise combination modes, but no intermolecular/solvent modes. A similar tail was observed earlier in the case of bacteriochlorophyll a, suggesting generality of this phenomenon. Long vibronic red tails are, thus, expected to exist in all pigments of light-harvesting relevance at physiological temperatures.

On uncorrelated inter-monomer Förster energy transfer in Fenna-Matthews-Olson complexes

The Fenna-Matthews-Olson (FMO) light-harvesting antenna protein of green sulfur bacteria is a long-studied pigment-protein complex which funnels energy from the chlorosome to the reaction centre where photochemistry takes place. The structure of the FMO protein from Chlorobaculum tepidum is known as a homotrimeric complex containing eight bacteriochlorophyll a per monomer. Owing to this structure FMO has strong intra-monomer and weak inter-monomer electronic coupling constants. While long-lived (sub-picosecond) coherences within a monomer have been a prevalent topic of study over the past decade, various experimental evidence supports the presence of subsequent inter-monomer energy transfer on a picosecond time scale. The latter has been neglected by most authors in recent years by considering only sub-picosecond time scales or assuming that the inter-monomer coupling between low-energy states is too weak to warrant consideration of the entire trimer. However, Förster theory predicts that energy transfer of the order of picoseconds is possible even for very weak (less than 5 cm^-1) electronic coupling between chromophores. This work reviews experimental data (with a focus on emission and hole-burned spectra) and simulations of exciton dynamics which demonstrate inter-monomer energy transfer. It is shown that the lowest energy 825 nm absorbance band cannot be properly described by a single excitonic state. The energy transfer through FMO is modelled by generalized Förster theory using a non-Markovian, reduced density matrix approach to describe the electronic structure. The disorder-averaged inter-monomer transfer time across the 825 nm band is about 27 ps. While only isolated FMO proteins are presented, the presence of inter-monomer energy transfer in the context of the overall photosystem is also briefly discussed.

Structural basis of light-harvesting in the photosystem II core complex

Photosystem II (PSII) is a membrane-spanning, multi-subunit pigment-protein complex responsible for the oxidation of water and the reduction of plastoquinone in oxygenic photosynthesis. In the present review, the recent explosive increase in available structural information about the PSII core complex based on X-ray crystallography and cryo-electron microscopy is described at a level of detail that is suitable for a future structure-based analysis of light-harvesting processes. This description includes a proposal for a consistent numbering scheme of protein-bound pigment cofactors across species. The structural survey is complemented by an overview of the state of affairs in structure-based modeling of excitation energy transfer in the PSII core complex with emphasis on electrostatic computations, optical properties of the reaction center, the assignment of long-wavelength chlorophylls, and energy trapping mechanisms.

On Excitation Energy Transfer within the Baseplate BChl a-CsmA Complex of Chloroflexus aurantiacus

Recently, a hybrid approach combining solid-state NMR spectroscopy and cryo-electron microscopy showed that the baseplate in green sulfur bacterium Chlorobaculum tepidum is a 2D lattice of BChl a-CsmA dimers [Nielsen, J. T.; et al., Nat. Commun. 2016, 7, 12454-12465]. While the existence of the BChl a-CsmA subunit was previously known, the proposed orientations of the BChl a pigments had only been elucidated from spectral data up to this point. Regarding the electronic structure of the baseplate, two models have been proposed. 2D electronic spectroscopy data were interpreted as revealing that at least four excitonically coupled BChl a might be in close contact. Conversely, spectral hole burning data suggested that the lowest energy state was localized, yet additional states are sometimes observed because of the presence of the Fenna-Matthews-Olson (FMO) antenna protein. To solve this conundrum, this work studies the chlorosome-baseplate complex from Chloroflexus aurantiacus, which does not contain the FMO protein. The results confirm that in both C. tepidum and C. aurantiacus, excitation energy is transferred to a localized low-energy trap state near 818 nm with similar rates, most likely via exciton hopping.

Influence of Hydrogen Bonds on the Electron-Phonon Coupling Strength/Marker Mode Structure and Charge Separation Rates in Reaction Centers from Rhodobacter sphaeroides

Low-temperature persistent and transient hole-burning (HB) spectra are presented for the triple hydrogen-bonded L131LH + M160LH + M197FH mutant of Rhodobacter sphaeroides. These spectra expose the heterogeneous nature of the P-, B-, and H-bands, consistent with a distribution of electron transfer (ET) times and excitation energy transfer (EET) rates. Transient P⁺Q_A^- holes are observed for fast (tens of picoseconds or faster) ET times and reveal strong coupling to phonons and marker mode(s), while the persistent holes are bleached in a fraction of reaction centers with long-lived excited states characterized by much weaker electron-phonon coupling. Exposed differences in electron-phonon coupling strength, as well as a different coupling to the marker mode(s), appear to affect the ET times. Both resonantly and nonresonantly burned persistent HB spectra show weak blue- (∼150 cm^-1) and large, red-shifted (∼300 cm^-1) antiholes of the P band. Slower EET times from the H- and B-bands to the special pair dimer provide new insight on the influence of hydrogen bonds on mutation-induced heterogeneity.

Higher Order Vibronic Sidebands of Chlorophyll a and Bacteriochlorophyll a for Enhanced Excitation Energy Transfer and Light Harvesting

Optical absorption and fluorescence spectra of molecules in condensed phases often show extensive sidebands. Originating from electron-vibrational and electron-phonon couplings, these spectral tails bear important information on the dynamics of electronic states and processes the molecules are involved in. The vibronic sidebands observed in conjugate Q_y absorption and fluorescence spectra of chlorophyll a and bacteriochlorophyll a are relatively weak, characterized by the total Huang-Rhys factor which is less than one. Therefore, it is widely considered that only fundamental intramolecular modes are responsible for their formation. Here, we provide evidence for extra-long and structured fluorescence tails of chlorophyll a and bacteriochlorophyll a as far as 4000 cm^-1 from respective spectral origins, far beyond the frequency range of fundamental modes. According to quantum chemical simulations, these sidebands extending to ∼960 nm in chlorophyll a and ∼1140 nm in bacteriochlorophyll a into the infrared part of the optical spectrum are mainly contributed to by vibrational overtones of the fundamental modes. Because energy transfer and relaxation processes generally depend on vibronic overlap integrals, these findings potentially contribute to better understanding of many vital photo-induced phenomena, including photosynthetic light harvesting.

How Well Does the Hole-Burning Action Spectrum Represent the Site-Distribution Function of the Lowest-Energy State in Photosynthetic Pigment-Protein Complexes?

For the first time, we combined Monte Carlo and nonphotochemical hole burning (NPHB) master equation approaches to allow for ultrahigh-resolution (<0.005 cm^-1, smaller than the typical homogeneous line widths at 5 K) simulations of the NPHB spectra of dimers and trimers of interacting pigments. These simulations reveal significant differences between the zero-phonon hole (ZPH) action spectrum and the site-distribution function (SDF) of the lowest-energy state. The NPHB of the lowest-energy pigment, following the excitation energy transfer (EET) from the higher-energy pigments which are excited directly, results in the shifts of all excited states. These shifts affect the ZPH action spectra and EET times derived from the widths of the spectral holes burned in the donor-dominated regions. The effect is present for a broad variety of realistic antihole functions, and it is maximal at relatively low values of interpigment coupling (V ≤ 5 cm^-1) where the use of the Förster approximation is justified. These findings need to be considered in interpreting various optical spectra of photosynthetic pigment-protein complexes for which SDFs (describing the inhomogeneous broadening) are often obtained directly from the ZPH action spectra. Water-soluble chlorophyll-binding protein (WSCP) was considered as an example.

Robust and Fragile Quantum Effects in the Transfer Kinetics of Delocalized Excitons between B850 Units of LH2 Complexes

Aggregates of light harvesting 2 (LH2) complexes form the major exciton-relaying domain in the photosynthetic unit of purple bacteria. Application of a generalized master equation to pairs of the B850 units of LH2 complexes, where excitons predominantly reside, provides quantitative information on how the inter-LH2 exciton transfer depends on the distance, relative rotational angle, and the relative energies of the two LH2s. The distance dependence demonstrates significant enhancement of the rate due to quantum delocalization of excitons, the qualitative nature of which remains robust against the disorder. The angle dependence reflects isotropic nature of exciton transfer, which remains similar for the ensemble of disorder. The variation of the rate on relative excitation energies of LH2 exhibits resonance peaks, which, however, is fragile as the disorder becomes significant. Overall, the average transfer times between two LH2s are estimated to be in the range of 4-25 ps for physically plausible inter-LH2 distances.

Wavelength-Dependent Exciton-Vibrational Coupling in the Water-Soluble Chlorophyll Binding Protein Revealed by Multilevel Theory of Difference Fluorescence Line-Narrowing

One of the most powerful line-narrowing techniques used to unravel the homogeneous lineshapes of inhomogeneously broadened systems is difference fluorescence line-narrowing spectroscopy. When this spectroscopy was applied to multichromophoric systems so far, the spectra were analyzed by an effective two-level system approach, composed of the electronic ground state and the lowest exciton state. An effective Huang-Rhys factor was assigned for the coupling of this state to the vibrations. Here, we extend this approach by including a multilevel line shape theory, which takes into account the excitonic coupling between pigments and thereby the effect of the delocalization of the excited states explicitly. In this way, it becomes possible to extract the spectral density of the local exciton-vibrational coupling. The theory is applied to the recombinant water-soluble chlorophyll binding protein reconstituted with chlorophyll a or b and reveals a significant decrease of the Huang-Rhys factor of the local exciton-vibrational coupling with decreasing transition energy of the chlorophylls. This decrease could be due to the increase in steric interactions reducing the flexibility of the environment and red-shifting the site energy of the pigments.

Impact of Single-Point Mutations on the Excitonic Structure and Dynamics in a Fenna-Matthews-Olson Complex

Hole burning (HB) spectroscopy and modeling studies reveal significant changes in the excitonic structure and dynamics in several mutants of the FMO trimer from the Chlorobaculum tepidum. The excited-state decay times ( T₁) of the high-energy excitons are significantly modified when mutation occurs near bacteriochlorophyll (BChl) 1 (V152N mutant) or BChl 6 (W184F). Longer (averaged) T₁ times of highest-energy excitons in V152N and W184F mutants suggest that site energies of BChls 1 and 6, believed to play an important role in receiving excitation from the baseplate BChls, likely play a critical role to ensure the femtosecond (fs) energy relaxation observed in wild-type FMO. HB spectroscopy reveals preferentially slower T₁ times (about 1 ps on average) because fs times prohibit HB due to an extremely low HB quantum yield. Uncorrelated (incoherent) excitation energy transfer times between monomers, the composition of exciton states, and average, frequency-dependent, excited-state decay times ( T₁) are discussed.

Excitonic Energy Landscape of the Y16F Mutant of the Chlorobium tepidum Fenna-Matthews-Olson (FMO) Complex: High Resolution Spectroscopic and Modeling Studies

We report high-resolution (low-temperature) absorption, emission, and nonresonant/resonant hole-burned (HB) spectra and results of excitonic calculations using a non-Markovian reduced density matrix theory (with an improved algorithm for parameter optimization in heterogeneous samples) obtained for the Y16F mutant of the Fenna-Matthews-Olson (FMO) trimer from the green sulfur bacterium Chlorobium tepidum. We show that the Y16F mutant is a mixture of FMO complexes with three independent low-energy traps (located near 817, 821, and 826 nm), in agreement with measured composite emission and HB spectra. Two of these traps belong to mutated FMO subpopulations characterized by significantly modified low-energy excitonic states. Hamiltonians for the two major subpopulations (Sub₈₂₁ and Sub₈₁₇) provide new insight into extensive changes induced by the single-point mutation in the vicinity of BChl 3 (where tyrosine Y16 was replaced with phenylalanine F16). The average decay time(s) from the higher exciton state(s) in the Y16F mutant depends on frequency and occurs on a picosecond time scale.

Contribution of low-temperature single-molecule techniques to structural issues of pigment-protein complexes from photosynthetic purple bacteria

As the electronic energies of the chromophores in a pigment-protein complex are imposed by the geometrical structure of the protein, this allows the spectral information obtained to be compared with predictions derived from structural models. Thereby, the single-molecule approach is particularly suited for the elucidation of specific, distinctive spectral features that are key for a particular model structure, and that would not be observable in ensemble-averaged spectra due to the heterogeneity of the biological objects. In this concise review, we illustrate with the example of the light-harvesting complexes from photosynthetic purple bacteria how results from low-temperature single-molecule spectroscopy can be used to discriminate between different structural models. Thereby the low-temperature approach provides two advantages: (i) owing to the negligible photobleaching, very long observation times become possible, and more importantly, (ii) at cryogenic temperatures, vibrational degrees of freedom are frozen out, leading to sharper spectral features and in turn to better resolved spectra.

Tracing feed-back driven exciton dynamics in molecular aggregates

Excitation, exciton transport, dephasing and energy relaxation, and finally detection processes shift molecular systems into a specific superposition of quantum states causing localization, local heating and finally excitonic polaronic effects.

Spectrally Selective Spectroscopy of Native Ca-Containing and Ba-Substituted LH1-RC Core Complexes from Thermochromatium tepidum

The LH1-RC core complex from the thermophilic photosynthetic purple sulfur bacterium Thermochromatium tepidum has recently attracted interest of many researchers because of its several unique properties, such as increased robustness against environmental hardships and the much red-shifted near-infrared absorption spectrum of the LH1 antenna exciton polarons. The known near-atomic-resolution crystal structure of the complex well supported this attention. Yet several mechanistic aspects of the complex prominence remained to be understood. In this work, samples of the native, Ca²⁺-containing core complexes were investigated along with those destabilized by Ba²⁺ substitution, using various spectrally selective steady-state and picosecond time-resolved spectroscopic techniques at physiological and cryogenic temperatures. As a result, the current interpretation of exciton spectra of the complex was significantly clarified. Specifically, by evaluating the homogeneous and inhomogeneous compositions of the spectra, we showed that there is little to no effect of cation substitution on the dynamic or kinetic properties of antenna excitons. Reasons of the extra red shift of absorption/fluorescence spectra observed in the Ca-LH1-RC and not in the Ba-LH1-RC complex should thus be searched in subtle structural differences following the inclusion of different cations into the core complex scaffold.

Excitation energy transfer in phycobiliproteins of the cyanobacterium Acaryochloris marina investigated by spectral hole burning

The cyanobacterium Acaryochloris marina developed two types of antenna complexes, which contain chlorophyll-d (Chl d) and phycocyanobilin (PCB) as light-harvesting pigment molecules, respectively. The latter membrane-extrinsic complexes are denoted as phycobiliproteins (PBPs). Spectral hole burning was employed to study excitation energy transfer and electron-phonon coupling in PBPs. The data reveal a rich spectral substructure with a total of four low-energy electronic states whose absorption bands peak at 633, 644, 654, and at about 673 nm. The electronic states at ~633 and 644 nm can be tentatively attributed to phycocyanin (PC) and allophycocyanin (APC), respectively. The remaining low-energy electronic states including the terminal emitter at 673 nm may be associated with different isoforms of PC, APC, or the linker protein. Furthermore, the hole burning data reveal a large number of excited state vibrational frequencies, which are characteristic for the chromophore PCB. In summary, the results are in good agreement with the low-energy level structure of PBPs and electron-phonon coupling parameters reported by Gryliuk et al. (BBA 1837:1490-1499, 2014) based on difference fluorescence line-narrowing experiments.