DOI:
10.1039/C6SC03833F
(Edge Article)
Chem. Sci., 2017,
8, 1621-1630
ortho and para chromophores of green fluorescent protein: controlling electron emission and internal conversion†
Received
26th August 2016
, Accepted 5th November 2016
First published on 7th November 2016
Abstract
Green fluorescent protein (GFP) continues to play an important role in the biological and biochemical sciences as an efficient fluorescent probe and is also known to undergo light-induced redox transformations. Here, we employ photoelectron spectroscopy and quantum chemistry calculations to investigate how the phenoxide moiety controls the competition between electron emission and internal conversion in the isolated GFP chromophore anion, following photoexcitation with ultraviolet light in the range 400–230 nm. We find that moving the phenoxide group from the para position to the ortho position enhances internal conversion back to the ground electronic state but that adding an additional OH group to the para chromophore, at the ortho position, impedes internal conversion. Guided by quantum chemistry calculations, we interpret these observations in terms of torsions around the C–C–C bridge being enhanced by electrostatic repulsions or impeded by the formation of a hydrogen-bonded seven-membered ring. We also find that moving the phenoxide group from the para position to the ortho position reduces the energy required for detachment processes, whereas adding an additional OH group to the para chromophore at the ortho position increases the energy required for detachment processes. These results have potential applications in tuning light-induced redox processes of this biologically and technologically important fluorescent protein.
1 Introduction
Green fluorescent protein (GFP), together with its family of variants, is the most widely used fluorescent protein for in vivo monitoring of biological and biochemical processes.1–4 The chromophore that lies at the heart of GFP is 4-hydroxybenzylidene-1,2-dimethylimidazolinone, p-HBDI (the relevant anionic deprotonated form is shown in Fig. 1). The chromophore is anchored, both covalently and via a hydrogen-bonded network, to the protein that is wrapped around it in a β-barrel structure. GFP has two absorption bands centered around 395 nm and 480 nm, attributed to the neutral and anionic forms, respectively. Following excitation of the neutral form of the chromophore, excited-state proton transfer takes place via a proton relay between water molecules and amino acid residues, resulting in very intense fluorescence (quantum yield Φ ≈ 0.8) from the deprotonated anionic form of the chromophore.5 This fluorescence is lost when the protein is denatured, although it returns upon renaturation or cooling below the glass transition.6,7 The isolated chromophore is virtually non-fluorescent in solution8 and in vacuo,9 but can become fluorescent when incorporated in rigid scaffolds such as metal–organic frameworks or other non-native protein environments.10–13 The lack of fluorescence from the free chromophores is widely accepted to be due to ultrafast isomerisation around the C2–C3–C4 bridge (Fig. 1) followed by efficient internal conversion taking excited state population back to the electronic ground state.14,15 When the molecular framework of the chromophore is held rigid by the protein or a rigid scaffold, isomerisation around the C–C–C bridge is impeded and fluorescence from the excited state is the primary decay channel.
![image file: c6sc03833f-f1.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-f1.gif) |
| Fig. 1 Structures of the model GFP chromophores employed in this work. | |
Knowing that the environment of the chromophore plays a pivotal role in determining the spectral properties of GFP, which are important for imaging applications, led to numerous investigations to determine how structural modifications to the protein, such as changing single amino acid residues, can be used to tune the absorption or fluorescence properties.4,16 Recently, there has also been a great deal of interest in determining how the structure and environment of the isolated chromophore influences the absorption properties17,18 and dynamics.19 For example, variants of the isolated chromophore in which the hydroxyl group was moved from the para position to the ortho position were found to be weakly fluorescent in cyclohexane (Φ ≈ 0.003).20,21 The fluorescence was attributed to an intramolecular hydrogen-bond between the ortho OH group and the N atom of the imidazolinone group forming a seven-membered ring, which partly locked the molecular framework and restricted isomerisation around the C–C–C bridge. Recently, it was shown that the fluorescence quantum yield could be enhanced further (Φ ≈ 0.2) by introducing a five-membered ring to completely lock the C3–C4 single bond in the C–C–C bridge (Fig. 1).22 Chromophores in which the hydroxyl group has been moved from the para position to the meta position have also been found to be weakly fluorescent (Φ ≈ 0.002);23,24 in these cases, the fluorescence was attributed to changes in the electronic structure introducing a barrier along the isomerisation coordinate and thus decreasing the rate of isomerisation around the C–C–C bridge.
Upon irradiation with UV light, not only does GFP fluoresce, but it also has the ability to act as a light-induced electron donor,25,26 which is important for monitoring and manipulating redox processes in cells. The idea that the position of the hydroxy group and structural rigidity play key roles in controlling the fluorescence emission quantum yield, led us to ask how moving the O− group and partly locking the molecular framework of the deprotonated chromophore anions might affect redox transformations.
A number of recent photoelectron spectroscopy and theoretical studies of the gas-phase isolated deprotonated chromophore anion, p-HBDI−, have shown that the vertical excitation energy of the first bright electronically excited ππ* state lies below the adiabatic threshold for photodetachment.27–30 The most accurate measurement of the adiabatic threshold for photodetachment has been made using vibrationally-resolved photoelectron spectroscopy and is found to be 2.73 eV.31 A separate photoelectron spectroscopy study has shown that adding water molecules to solvate the chromophore anion raises the adiabatic photodetachment threshold,32 supporting earlier computational studies showing that the hydrogen-bonding environment stabilises the first excited state against electron emission.33 Other computational studies have focused on the effect of the protein and specific interactions on the detachment energy.34,35 A femtosecond time-resolved study of the excited state dynamics of this first excited state found two characteristic timescales: 330 fs for isomerisation around the C–C–C bridge and 1.4 ps for subsequent internal conversion to the anion ground state,36 remarkably similar to the excited state dynamics of the isolated chromophore in solution measured by fluorescence up-conversion.14 A subsequent time-resolved study found the lifetime for competing autodetachment from the first excited state to be 30 ps,37 consistent with calculations.38 Just above the first detachment threshold, the second and third excited ππ* states are accessible39 and autodetachment from these states is found to compete with direct detachment.27,30,40
The aim of the work presented in this paper was to investigate how moving the O− group from the para position to the ortho position, or adding an additional OH group at the ortho position to partly lock the molecular framework of the deprotonated chromophore anions, might affect the competition between electron detachment and internal conversion. Specifically, we recorded photoelectron spectra of p-HBDI, o-HBDI and op-DHBDI in their deprotonated anionic forms (Fig. 1) following photoexcitation with ultraviolet light in the range 400–230 nm. Using quantum chemistry calculations to guide our interpretation of the spectra, we find that our results support the idea that partial locking of the C–C–C bridge inhibits internal conversion of the deprotonated chromophore anion and we find that moving the position of the O− group to the ortho position results in changes to the electronic structure that enhance internal conversion. These results have potential applications in tuning light-induced redox processes in biology.
2 Methodology
2.1 Chromophores
4-Hydroxybenzylidene-1,2-dimethylimidazolinone (p-HBDI) was prepared using reported procedures41 and 2-hydroxybenzylidene-1,2-dimethylimidazolinone (o-HBDI) was prepared using a modified version of a literature procedure.42 2,4-Dihydroxybenzylidene-1,2-dimethylimidazolinone22 (op-DHBDI) was synthesised using an analogous route to those employed for p-HBDI and o-HBDI, via the dimethoxy analogue43 and deprotection, which was higher yielding than using 2,4-dihydroxybenzaldehyde directly.†
2.2 Photoelectron spectroscopy
Photoelectron spectra were recorded using our instrument that combines electrospray ionisation (ESI) and photoelectron velocity map imaging (VMI).27,30,44–46 Briefly, deprotonated anions were generated in the gas-phase by ESI of ∼0.1 mM solutions of chromophores in methanol with a few drops of aqueous ammonia. Chromophore anions were mass-selected by a quadrupole before passing into a hexapole ion trap. Anions were pulsed out of the trap and guided to the photodetachment region of a collinear VMI spectrometer where they interacted with laser pulses with wavelengths of 400 nm (3.10 eV), 346.1 nm (3.582 eV) and 230.7 nm (5.374 eV). The 400 nm pulses were generated by frequency-doubling the output of a commercial, amplified, Ti:sapphire femtosecond laser system operating at 1 kHz and the 346.1 nm and 230.7 nm pulses were generated by frequency doubling or tripling the output of a Nd:YAG pumped dye laser operating at 20 Hz. Following the anion–laser interaction, photodetached electrons were accelerated towards a position sensitive detector and imaged using a phosphor screen and CCD camera (Photek). Since anions and scattered laser light generate background counts, anion-only and laser-only images were subtracted from images recorded following the anion–laser interaction. The photoelectron images were inverted using the pBASEX method.47 The wavelength of the laser light was measured using a wavemeter and the energy scale of the detector was calibrated by recording the photoelectron spectrum of I−. The energy resolution of our instrument is <5% and the error in electron kinetic energy (eKE) is around 0.05 eV.
2.3 Calculations
In order to calculate VDEs of each chromophore, a combination of electron propagator theory48 (EPT) and time-dependent density functional theory49 (TD-DFT) calculations were employed using the Gaussian 09 computational chemistry software package.50 The geometries of the chromophore anions were optimised using B3LYP51,52 with the 6-311++G(3df,3pd)53 basis set. Vertical detachment energies from the minimum energy of the ground electronic state of the anion, S0, to the ground electronic state of the corresponding neutral radical, D0, were calculated using EPT applying the outer valence Green's function (OVGF) propagator54 with the 6-311++G(3df,3pd) basis set. This method has been found in earlier work to yield reliable VDEs.30,45,46 To calculate VDEs for higher lying electronically excited states of the neutral radical, D1–D3, TD-DFT calculations were carried out on the ground electronic state of the neutral radical at the minimum energy geometry of the ground electronic state of the anion using the CAM-B3LYP55/6-311++G(3df,3pd) method. These VEEs were added to the previously calculated vertical detachment energies to D0 to obtain the higher lying detachment energies. This combined methodology is referred to as EPT + TD-DFT.
In addition, the equation-of-motion coupled-cluster method with single and double excitations for the calculation of ionisation potentials56 (EOM-IP-CCSD) was employed together with the aug-cc-pVDZ57 basis set to calculate the first four VDEs. These calculations were carried out using the Q-Chem computational chemistry software package.58 The geometries of the deprotonated chromophore anions were minimised with Møller–Plesset perturbation theory of second order59 (MP2) and the aug-cc-pVDZ basis set using the Gaussian computational chemistry package.
VEEs of the excited singlet states of the anion were computed at the MP2 minima with the algebraic diagrammatic construction method to second order60,61 (ADC(2)) using the aug-cc-pVDZ basis set, employing the Turbomole computational chemistry software package.62 ADC(2) was chosen to be able to calculate higher electronically excited states reliably at a reasonable computational cost. Due to the number of diffuse functions in the basis set, a high number of continuum states are being calculated where that number of states depends on the selected basis set; however, the lower continuum states converge with increasing number of calculated states and smaller basis sets that exclude continuum states neglect the interaction of resonance states with the continuum.63
3 Results and discussion
3.1 Photoelectron spectra
Photoelectron spectra of the deprotonated chromophores were recorded as a function of electron kinetic energy (eKE) and are presented as a function of electron binding energy, eBE = hν − eKE in Fig. 2. The electron emission processes contributing the photoelectron spectra are illustrated on a Jablonski diagram in Fig. 3.
![image file: c6sc03833f-f2.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-f2.gif) |
| Fig. 2 Photoelectron spectra of p-HBDI− (left) o-HBDI− (centre) and op-DHBDI− (right) recorded at 400 nm (3.10 eV), 346.1 nm (3.582 eV) and 230.7 nm (5.374 eV). Individual spectra have been normalised to their maximum intensity. Vertical dashed orange and blue lines mark the first three vertical detachment energies, calculated using the EPT/6-311++G(3df,3pd) and EOM-IP-CCSD/aug-cc-pVDZ methods, respectively. Vertical green lines mark the VEEs of the anions calculated using the ADC(2)/aug-cc-pVDZ method (solid green). | |
![image file: c6sc03833f-f3.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-f3.gif) |
| Fig. 3 Jablonski diagram illustrating the electronic relaxation and emission processes following UV photoexcitation of the model GFP chromophores employed in this work. Horizontal black lines represent the vibrational levels of the excited electronic states and the solid green area represents the electron detachment continuum. Vertical blue arrows represent the eKE of direct and indirect photodetachment (PD) processes and the solid grey areas represent the vibrational energy left in the neutral radical after electron detachment (determined by the propensity for conservation of vibrational energy). The horizontal black arrows represent some of the possible internal conversion (IC) processes and the horizontal blue arrow represents electrons with eKE ∼0 eV following thermionic emission (TE) from vibrationally hot S0. | |
3.1.1
p-HBDI−.
To benchmark our experiments, we recorded a photoelectron spectrum of p-HBDI− at 350.0 nm (3.542 eV), which is similar to previously reported spectra recorded around this wavelength.28–31 It should be noted that this spectrum is significantly narrower than our own previously reported spectrum.
Spectra recorded at 400 nm (3.10 eV), 346.1 nm (3.582 eV) and 230.7 nm (5.374 eV) all have steep rising edges at around 2.7 eV, corresponding to direct detachment from S0 to D0 (Fig. 3), as reported by several experimental and theoretical groups.27–31 This is very close to the adiabatic detachment threshold which has been determined experimentally to be 2.73 eV and to be equivalent to the VDE.31 It is worth noting that the peak maxima are blue-shifted by around 0.1 eV from the calculated VDE, in agreement with calculated photoelectron spectra.64
The 400 nm photoelectron spectrum has a sharp feature at high eBE around 3.10 eV, very close to the photon energy of 3.10 eV. This high eBE (low eKE) feature has a profile that decays exponentially with eKE; this statistical eKE distribution is characteristic of thermionic emission from the hot ground electronic state of an anion, S0, following internal conversion from an excited electronic state (Fig. 3). The 346.1 nm photoelectron spectrum is dominated by direct detachment from S0 to D0 and there is no evidence of thermionic emission from the hot ground state. The 230.7 nm photoelectron spectrum is dominated by three, broad, unresolved features centered around 2.8 eV, 4.4 eV and 5.3 eV, which are attributed to detachment from S0 to D0, D1 and D2, respectively (Section 3.2). The broadening of the direct detachment peak on its high eBE edge in the 230.7 nm spectrum, compared to 346.1 nm spectrum, can be attributed to indirect detachment to D0 after population of high lying electronically excited states of the anion around 230.7 nm (5.374 eV) as illustrated in Fig. 3.
3.1.2
o-HBDI−.
Photoelectron spectra of o-HBDI− recorded at 400 nm (3.10 eV), 346.1 nm (3.582 eV) and 230.7 nm (5.374 eV) all have rising edges at around 2.5 eV, corresponding to direct detachment from S0 to D0. This is ∼0.2 eV lower than in p-HBDI−, indicating that moving the O− group from the para position to the ortho position lowers the detachment energy. This can be rationalised in terms of repulsion between the O− group and the N-atom on the imidazolinone group pushing charge from the O− group into the conjugated system, if we assume that the deprotonated chromophore anions are formed in the gas-phase in a cis geometry, i.e. the same as their parent neutral molecules. Calculations also show that in the cis geometry the conjugation over the bridge is weakened by the O− group, which will also destabilise the molecular orbital (MO) from which the electron is removed (see Section 3.2). This is consistent with detailed investigations of fluorescent protein chromophores showing that the trend in detachment energy for the series of deprotonated chromophore anions o-HBDI−, m-HBDI− and p-HBDI−, can be explained in terms of the resonance stabilisation of the anion.65,66 This is significant since it has been shown that the trend in redox potentials of a series of equivalent neutral chromophores is determined by the variation in ionisation energy since the solvent effects for structurally similar chromophores are similar.67
The 400 nm photoelectron spectrum is dominated by a sharp peak with eBE around the photon energy of 3.10 eV and whilst there is some contribution from direct detachment, it is much less significant. This suggests that moving the O− group from the para position to the ortho position results in a more efficient internal conversion process funnelling excited state population back to the ground electronic state. The photoelectron spectrum recorded at 346.1 nm has a broad feature peaking around 2.7 eV and a sharp feature at high eBE (low eKE). Thus, unlike p-HBDI−, in o-HBDI−, internal conversion is still able to compete with direct detachment following excitation at 346.1 nm. The 230.7 nm photoelectron spectrum is dominated by two, broad, unresolved features centered around 2.6 eV and 4.2 eV, which are attributed to detachment from S0 to D0 and D1 (Section 3.2). There is also a third, less intense feature centered around 5.1 eV.
3.1.3
op-DHBDI−.
Photoelectron spectra of op-DHBDI− recorded at 400 nm (3.10 eV), 346.1 nm (3.582 eV) and 230.7 nm (5.374 eV) all have steep rising edges at an eBE of around 2.9 eV, corresponding to direct detachment from S0 to D0. The photoelectron spectra of op-DHBDI− are remarkably similar to those of p-HBDI−, but shifted by ∼0.2 eV to higher eBE. The op-DHBDI chromophore has two deprotonation sites. Since the ortho-OH group stabilises the anion by forming a hydrogen bond with the nearby N atom of the imidazolinone unit, we conclude that the op-DHBDI− is deprotonated predominantly at the para site. From the MOs (see Sec. 3.2), the O-atom of the ortho-OH group does not play a role in the ring conjugation, so it acts as an electron withdrawing group and stabilises the MO from which the electron is removed.
In contrast to p-HBDI− and o-HBDI−, the spectrum of op-HBDI− at 400 nm is dominated by direct detachment, although there is a slight increase in photoelectron counts at high eBE (low eKE) suggesting that there is a small contribution from thermionic emission from the ground electronic state. This seems most likely to be the result of hydrogen bonding between the ortho-OH group and the N atom of the imidazolinone unit restricting torsional motion around the central C–C–C bridge, which is required for efficient internal conversion back to the ground electronic state.14,15,36,38 The 230.7 nm photoelectron spectrum is dominated by two, broad, unresolved features centered around 3.0 eV and 4.4 eV, which are attributed to detachment from S0 to D0 and D1, respectively (Section 3.2). As with p-HBDI−, broadening of the direct detachment peak on its high eBE edge in the 230.7 nm spectrum, compared to the 346.1 nm spectrum, can be attributed to indirect detachment to D0 following population of an electronically excited state of the anion at 230.7 nm (5.374 eV). It is interesting to note that the band associated with detachment from S0 to D1 is more intense than the band associated with detachment from S0 to D0, which is opposite to the relative intensities observed in the photoelectron spectra of the p-HBDI− and o-HBDI− chromophores.
3.2 Computational results
All three molecules are found to have stable planar geometries calculated at the MP2/aug-cc-pVDZ level of theory (Fig. 4).† The effect of adding the substituents at the ortho position is to open up the bridge. In p-HBDI− the C–C–C angle of the bridge joining the two rings is 131.1°, in op-DHBDI− it is 134.4° and in o-HBDI− 137.8°. The planar geometry of cis o-HBDI− is, however, not a minimum which is actually found to be with a C2–C3–C4–C5 torsion angle at the bridge of 23.5° (using the numbering in Fig. 1), which allows the bridge angle to close up to 133.3°. This is a very shallow minimum of 0.017 eV, indicating that the conjugation on the bridge has been weakened making it more susceptible to rotation.
![image file: c6sc03833f-f4.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-f4.gif) |
| Fig. 4 (a) Overlay of the optimised structures of p-HBDI−, op-DHBDI− and cis o-HBDI−. (b)–(d) Geometries of the anions with selected geometrical data. | |
It should be noted that the o-HBDI− chromophore has two stable conformations, with the O− group either cis or trans to the imidazolinone nitrogen (atom 1 in Fig. 1). In the parent neutral molecule, the cis chromophore is the minimum energy geometry, because it is stabilised by a hydrogen bond between the OH group and the imidizolinone nitrogen.20,68 A relaxed scan on the anion around the torsion angle starting in the cis geometry shows a barrier height of approximately 0.4 eV to the isomerisation, which is high compared to the energy available in the formation of the ion. Thus it seems likely that the chromophore is formed in the cis geometry during the ESI process even though the global minimum for this molecule is found to be the trans geometry by approximately 0.4 eV.† Our analysis of the photoelectron spectra supports this (Section 3.3), and for the rest of the discussion all results for o-HBDI refer to the cis conformer, unless stated otherwise.
In op-DHBDI−, the hydrogen bond between the ortho-OH group and the imidizolinone nitrogen stabilises the cis geometry of the anion and this is the minimum energy structure. Relaxed scans around the torsion support the reduced flexibility of the bridge in op-DHBDI− as the barrier to rotation from cis to trans in this molecule is approximately 1.4 eV.†
The VDEs calculated using both the EPT + TD-DFT and EOM-IP-CCSD methods are listed in Table 1, together with the characters of the EOM-IP-CCSD states. VDEs for trans o-HBDI− are very similar to those for the cis conformer.† The orbitals are presented in Fig. 5. The EOM-IP-CCSD calculated values for D0–D2 are indicated in Fig. 2, together with the EPT calculated value for D0 as a comparison. The higher lying ionisations calculated at the EPT + TD-DFT values are similar to the higher level of theory EOM-IP-CCSD results and are omitted from the plots for clarity. The VDEs to D0 for p-HBDI−, o-HBDI− and op-DHBDI− are 2.69 eV (2.78 eV), 2.54 eV (2.64 eV) and 2.90 eV (2.96 eV), calculated using the EOM-IP-CCSD (EPT + TDDFT) method. These calculated values follow the same trend as the experimental photoelectron spectra, with the VDEs increasing in the order, o-DHBDI− < p-HBDI− < op-HBDI−. For all the chromophores, the D0 state corresponds to a hole in a π orbital which is the HOMO.
Table 1 VDEs (in eV) for all anions calculated using the EPT + TD-DFT/6-311++G(3df,3pd) and EOM-IP-CCSD/aug-cc-pVDZ methods for the first four ionisation continua. The EPT value for D0 and all the EOM-IP-CCSD VDEs are plotted on the photoelectron spectra in Fig. 2. For EOM-IP-CCSD VDEs, the state character is given in parentheses to indicate the orbital in which the electron hole is present. The corresponding Hartree–Fock orbitals are presented in Fig. 5
Anion |
State |
EPT + TD-DFT |
EOM-IP-CCSD |
p-HBDI− |
D0 |
2.78 |
2.69 |
(π4) |
D1 |
4.34 |
4.51 |
(nO1) |
D2 |
4.85 |
4.87 |
(π3) |
D3 |
5.40 |
5.39 |
(π2) |
cis o-HBDI− |
D0 |
2.64 |
2.54 |
(π4) |
D1 |
3.97 |
3.92 |
(nO1) |
D2 |
4.31 |
4.20 |
(π3) |
D3 |
5.24 |
5.33 |
(nO2) |
op-DHBDI− |
D0 |
2.96 |
2.90 |
(π4) |
D1 |
4.46 |
4.48 |
(π3) |
D2 |
4.50 |
4.64 |
(nO1) |
D3 |
5.14 |
5.16 |
(π2) |
![image file: c6sc03833f-f5.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-f5.gif) |
| Fig. 5 Hartree–Fock (HF) orbitals of (a) p-HBDI−, (b) o-HBDI− and (c) op-DHBDI− used in the ADC(2)/aug-cc-pVDZ calculations. The HOMO is π4 and the LUMO is . The HF orbitals that were employed in the EOM-IP-CCSD calculations are very similar and therefore not shown. | |
To provide a benchmark of the calculated values, the VDE of p-HBDI− compares well to the experimental value determined by high-resolution vibrationally-resolved photoelectron spectroscopy of 2.73 eV (ref. 31) and the value calculated at the XMCQDPT2/(aug)-cc-pVTZ level of theory of 2.68 eV.29,39 The VDEs of the first two excited radical states D1 and D2 in Table 1 are slightly higher than the XMCQDPT2/(aug)-cc-pVTZ values of 3.97 eV and 4.40 eV; however, the characters are the same, with D1 having a hole in the nO orbital and D2 in a π orbital.39
As in p-HBDI−, D1 of o-HBDI− is of nO character, whereas D1 for op-DHBDI− has a hole in a π orbital. This explains the change in the relative intensities of the D0 and D1 bands observed in the 230.7 nm photoelectron spectrum of op-DHBDI− compared to the other two analogues since the transition probability (defined as the square of the norm of the Dyson orbital) between a continuum state of π character with a ππ* excited state would be greater than that for an nπ* excited state, especially if the same π orbital is involved. For all three chromophores, D1 and D2 are accessible at 230.7 nm; however, the VDEs for D1 are higher in energy than the 346.1 nm and 400 nm photon energies. Thus, at 346.1 nm and 400 nm, direct detachment to D1 and D2 is not possible and the D1 and D2 continua would only be accessible by indirect detachment from high-lying electronically excited states of the anions.
Relevant ADC(2)/aug-cc-pVDZ VEEs and their corresponding oscillator strengths and configurations are presented in Table 2 and plotted in Fig. 2. VEEs for trans o-HBDI− are very similar to those for the cis conformer.† These VEEs are for the electronically excited states with the strongest oscillator strengths and, thus, realistically accessible. All of these states, with one exception, are of ππ* character. The details of all calculated excited states for all three chromophores, including other ππ*, nπ* and continuum states, can be found in the ESI.† It is difficult to directly compare the ADC(2) excited states of p-HBDI to those calculated at the XMCQDPT2/(aug)-cc-pVTZ level69 as these short-lived states are very sensitive to basis set and method used. In both cases, however, a set of accessible (i.e. non-negligible oscillator strength) ππ* states are found at around 4 eV, which have shape resonance character with respect to the D0 continuum.
Table 2 VEEs (in eV) of all anions calculated with ADC(2)/aug-cc-pVDZ with main configurations and corresponding oscillator strengths, f. The orbitals involved in the excitations are shown in Fig. 5
State |
p-HBDI− |
cis o-HBDI− |
op-DHBDI− |
E/eV |
Config. |
f
|
E/eV |
Config. |
f
|
E/eV |
Config. |
f
|
S1 |
2.53 |
π4 → ![image file: c6sc03833f-t1.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t1.gif) |
0.9483 |
2.36 |
π4 → ![image file: c6sc03833f-t2.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t2.gif) |
0.4429 |
2.65 |
π4 → ![image file: c6sc03833f-t3.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t3.gif) |
0.9848 |
S2 |
3.99 |
π4 → ![image file: c6sc03833f-t4.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t4.gif) |
0.0445 |
3.88 |
π3 → ![image file: c6sc03833f-t5.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t5.gif) |
0.1743 |
4.15 |
π4 → ![image file: c6sc03833f-t6.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t6.gif) |
0.0176 |
S3 |
4.74 |
π2 → ![image file: c6sc03833f-t7.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t7.gif) |
0.0229 |
4.16 |
π4 → ![image file: c6sc03833f-t8.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t8.gif) |
0.0261 |
4.31 |
π4 → /π4 → ![image file: c6sc03833f-t10.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t10.gif) |
0.0390 |
S4 |
4.95 |
π4 → ![image file: c6sc03833f-t11.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t11.gif) |
0.0735 |
4.67 |
π4 → ![image file: c6sc03833f-t12.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t12.gif) |
0.0206 |
4.66 |
π4 → /π4 → ![image file: c6sc03833f-t14.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t14.gif) |
0.0666 |
S5 |
5.00 |
nO2 → ![image file: c6sc03833f-t15.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t15.gif) |
0.0114 |
4.93 |
π4 → ![image file: c6sc03833f-t16.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t16.gif) |
0.2487 |
4.78 |
π4 → ![image file: c6sc03833f-t17.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t17.gif) |
0.0236 |
S6 |
|
|
|
5.00 |
π2 → /π4 → ![image file: c6sc03833f-t19.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-t19.gif) |
0.4115 |
|
|
|
The VEEs and their corresponding configurations for p-HBDI− and op-DHBDI− show similar trends. In both chromophores, the first excited state (ππ*) is very bright with oscillator strengths close to unity, whereas all other states have very low oscillator strengths. Additionally, the energy levels of p-HBDI− and op-DHBDI− are very similar to one another. In contrast, the S1 state of o-HBDI− has a lower oscillator strength of 0.4429 but on the other hand, the S2 state has a higher oscillator strength than the S2 states in p-HBDI− and op-DHBDI−. The S2 state of o-HBDI− is described by a transition from the HOMO−1 into a π* orbital delocalised across the whole chromophore, whereas for p-HBDI− and op-DHBDI−, the S2 state corresponds to the promotion of an electron from the HOMO into a π* orbital located at the benzyl moiety. The corresponding MOs are presented in Fig. 6 and in the ESI.† Additionally in o-HBDI−, there are two higher excited states at 4.93 and 5.00 eV with reasonably high oscillator strengths. The state at 4.93 eV corresponds to a transition from a bound π orbital to a π* orbital with continuum character and the state at 5.00 eV has some of the same character. These states could play a role in the electron emission process in the 230.7 nm photoelectron spectra.
![image file: c6sc03833f-f6.tif](/image/article/2017/SC/c6sc03833f/c6sc03833f-f6.gif) |
| Fig. 6 Energy level diagram of (a) p-HBDI− and (b) o-HBDI− in which possible excitations and subsequent detachment processes are indicated (note that it is only the vibrationally excited states of S1 lying in the D0 continuum that have shape resonance character). S1 and S2 refer to the first two excited states selected in Table 2. The HOMO is π4 and the LUMO is . Alongside are the Hartree–Fock molecular orbitals that are involved in the excitations. For op-DHBDI−, a similar scheme applies as for p-HBDI− except the excited state resonance is S3 (Table 2). | |
3.3 Competition between detachment and internal conversion
Fig. 6 shows electronic energy level diagrams for p-HBDI−/op-DHBDI− and o-HBDI−, together with relevant molecular orbitals. Included are the states that are most important in the dynamics of detachment to the D0 continuum, which are the first two excited states with ππ* character, referred to as S1 and S2. Possible excitations and their corresponding detachment processes are indicated.
In the case of p-HBDI− (Fig. 6a), detachment from either electronically excited state requires the removal of one electron to reach D0. Hence, they have shape resonance character and autodetachment from both states is extremely fast. Due to its very large oscillator strength (0.9483), the transition from S0 to S1 is strongly resonant. The photon energy at 400 nm (3.10 eV) is low enough to excite S1 significantly, as illustrated by the action absorption spectrum of this chromophore.39 Population from S1 can then rapidly internally convert to highly excited vibrational levels of S0. The internal vibrational energy in S0 is high enough for thermionic emission to occur, which gives rise to the high eBE (low eKE) peak in the 400 nm spectrum. Such processes have been previously observed for p-HBDI−.28,30,40 and photoactive yellow protein chromophores.45,46,70
The S2 state lies at higher energy (3.99 eV), so it can only be excited significantly at the shortest laser wavelength used in our experiments. For this reason, it is unlikely that resonant processes from S2 to D0 impact the electron emission dynamics in the 400 nm or the shape of the D0 band. However, at 230.7 nm (5.374 eV) the photon energy is high enough to excite S2, whose minimum structure is significantly different from the S0 minimum. Moreover, the S2 state has excited shape resonance character with respect to the D0 continuum, which results in rapid indirect detachment and vibrational autodetachment from S2 and accounts for the broadening of the D0 band on the high eBE side.69 Such broadening is also observed on the high eBE side of the 346.1 nm photoelectron spectrum following excitation of low vibrational states of S2.69
The energy level diagram for op-DHBDI− is essentially identical to that of p-HBDI−, only shifted to slightly higher energies. At 400 nm, the photon energy is closer to the VDE of S1, so at first sight it would seem likely that indirect detachment might play a more important role than direct detachment in electron emission from op-DHBDI− compared to p-HBDI−. Since this is not the case and it appears that there is less thermionic emission from op-DHBDI− than p-HBDI− following excitation at 400 nm (Fig. 2), we conclude that the S1/S0 conical intersection is less accessible in op-DHBDI− than in p-HBDI−. This can be rationalised in terms of the hydrogen-bonding between the OH group in the ortho position and the N atom on the imidazolinone moiety restricting torsional motions around the C–C–C bridge, which are known to be important in directing the excited state population to the S1/S0 conical intersection. Importantly, these observations also support solution phase studies of neutral o-HBDI and op-DHBDI chromophores, in which the formation of a hydrogen-bond between the ortho-OH and the N of the imidazolinone group has been shown to lock the molecular framework sufficiently to enhance the fluorescence quantum yield from the first excited state by impeding internal conversion.21,22,42 Similar restrictions to torsional motions are also responsible for the enhanced fluorescence yield of chromophores locked inside rigid scaffolds.10–13
Due to the differences in the oscillator strengths and state characters for o-HBDI−, we obtain a very different picture for the dynamics, as illustrated in Fig. 6b. The S2 state has Feshbach resonance character with respect to D0 (detachment is a two-electron process) indicating that it is only weakly coupled to the continuum. Hence, autodetachment from S2 will be slow. Since there is the possibility of significantly populating S2 due to its relatively high oscillator strength (0.1743) compared with p-HBDI− and op-DHBDI−, and subsequent autodetachment is slow, we would expect internal conversion to lower lying electronic states to dominate. This explains why the peak at high eBE (low eKE), which can be attributed to thermionic emission from S0 or autodetachment from high lying energy levels of S1, appears in the 346.1 nm photoelectron spectrum of o-HBDI− but not in the 346.1 nm photoelectron spectra of p-HBDI− or op-DHBDI−. Also, at 400 nm, the photoelectron spectrum of o-HBDI− is completely dominated by a high eBE (low eKE) thermionic emission peak. It seems very unlikely that this could be attributed to internal conversion from S2 since at this photon energy it is unlikely that even the lowest vibrational levels of S2 would be accessible. It seems far more likely that excitation of high lying vibrational states of S1 results in more efficient internal conversion to S0 in o-HBDI− compared to p-HBDI− or op-DHBDI−. This can be rationalised in terms of the repulsion between the O− group at the ortho position and the N atom on the imidazolinone group enhancing torsional motions around the C–C–C bridge as discussed above (Section 3.2); and supports our suggestion that the o-HBDI− chromophore is formed predominantly in its cis geometry in the gas phase.
4 Conclusions
From the combination of experimental and computational results, we have advanced our understanding of how the OH group affects the photodynamics of the model GFP chromophore. For longer wavelengths of light closer to the visible range, we have shown that by moving the O− group from the para position to the ortho position, we have an effective way of enhancing internal conversion back to the S0 ground state of the chromophore as a decay pathway which is highly competitive with electron emission. Conversely, introducing an OH group at the ortho position in addition to the para O− group allows us to enhance the electron emission process by impeding internal conversion. This dynamical control is driven by changes to the molecular and nuclear structure that control torsional motions around the C–C–C bridge in the chromophore. We also find that moving the O− group from the para position to the ortho position results in a decrease in photon energy required for electron detachment by ∼0.2 eV whereas introducing an OH group at the ortho position in addition to the para O− group increases the photon energy required for electron detachment by ∼0.2 eV. This type of chemical control allows us to control the wavelength of light required to initiate the redox processes in GFP chromophores. The information obtained for processes at higher photon energies remains more general. Although this region may be of less interest for the fluorescence properties of GFP chromophores, it could be more important in the control of light induced electron driven processes. For the highest energy photons employed in this study (230.7 nm), we found that many additional electronically excited states of the anion and neutral radical become accessible for all chromophores. The photodynamics are similar for all the chromophores with the O− and OH groups simply acting as tuning substituents. The oscillator strengths of the higher lying states of the o-HBDI− chromophore are brighter than in the other two analogues, so may play a more important role in the electron emission processes at shorter wavelengths. The next step is to follow the excited state dynamics in time-resolved experiments, which we plan to do.
In summary, we have shown that the electron emission and competing processes can be controlled by modifying the substitution pattern in the phenoxide moiety. We propose that the ability to control these processes will play an important role in tuning light-induced redox processes of the biologically and technologically important family of fluorescent proteins derived from GFP.
Acknowledgements
This work was made possible by EPSRC grants (EP/L005646/1, EP/L005697/1 and EP/D054508/1). We acknowledge use of the EPSRC UK National Service for Computational Chemistry Software (NSCCS) at Imperial College London and computational support from Dr Jörg Saßmannshausen at UCL. We gratefully acknowledge the UCL mass spectrometry facility and A. E. Aliev in the NMR facility of the Dept. of Chemistry at UCL.
References
- R. Y. Tsien, Annu. Rev. Biochem., 1998, 67, 509–544 CrossRef CAS PubMed
.
- M. Zimmer, Chem. Rev., 2002, 102, 759–782 CrossRef CAS PubMed
.
- S. R. Meech, Chem. Soc. Rev., 2009, 38, 2922–2934 RSC
.
- A. S. Mishin, V. V. Belousov, K. M. Solntsev and K. A. Lukyanov, Curr. Opin. Chem. Biol., 2015, 27, 1–9 CrossRef CAS PubMed
.
- A. Usman, O. F. Mohammed, E. T. J. Nibbering, J. Dong, K. M. Solntsev and L. M. Tolbert, J. Am. Chem. Soc., 2005, 127, 11214–11215 CrossRef CAS PubMed
.
- M. Chattoraj, B. A. King, G. U. Bublitz and S. G. Boxer, Proc. Natl. Acad. Sci. U. S. A., 1996, 93, 8362–8367 CrossRef CAS
.
- N. M. Webber, K. L. Litvinenko and S. R. Meech, J. Phys. Chem. B, 2001, 105, 8036–8039 CrossRef CAS
.
- H. Niwa, S. Inouye, T. Hirano, T. Matsuno, S. Kojima, M. Kubota, M. Ohashi and F. Tsuji, Proc. Natl. Acad. Sci. U. S. A., 1996, 93, 13617–13622 CrossRef CAS
.
- M. W. Forbes and R. A. Jockusch, J. Am. Chem. Soc., 2009, 131, 17038–17039 CrossRef CAS PubMed
.
- A. Baldridge, S. R. Samanta, N. Jayaraj, V. Ramamurthy and L. M. Tolbert, J. Am. Chem. Soc., 2011, 133, 712–715 CrossRef CAS PubMed
.
- A. Baldridge, S. Feng, Y. T. Chang and L. M. Tolbert, ACS Comb. Sci., 2011, 13, 214–217 CrossRef CAS PubMed
.
- K. D. Warner, M. C. Chen, W. Song, R. L. Strack, A. Thorn, S. R. Jaffrey and A. R. Ferré-D'Amaré, Nat. Struct. Mol. Biol., 2014, 21, 658–663 CAS
.
- E. A. Dolgopolova, A. M. Rice, M. D. Smith and N. B. Shustova, Inorg. Chem., 2016, 55, 7257–7264 CrossRef CAS PubMed
.
- D. Mandal, T. Tahara and S. R. Meech, J. Phys. Chem. B, 2004, 108, 1102–1108 CrossRef CAS
.
- M. E. Martin, F. Negri and M. Olivucci, J. Am. Chem. Soc., 2004, 126, 5452–5464 CrossRef CAS PubMed
.
- R. Tsien, Angew. Chem., Int. Ed., 2009, 48, 5612–5626 CrossRef CAS PubMed
.
- J. B. Greenwood, J. Miles, S. de Camillis, P. Mulholland, L. Zhang, M. A. Parkes, H. C. Hailes and H. H. Fielding, J. Phys. Chem. Lett., 2014, 5, 3588–3592 CrossRef CAS PubMed
.
- H. V. Kiefer, E. Lattouf, N. W. Persen, A. V. Bochenkova and L. H. Andersen, Phys. Chem. Chem. Phys., 2015, 17, 20056–20060 RSC
.
- A. Toniolo, S. Olsen, L. Manohar and T. J. Martínez, Faraday Discuss., 2004, 127, 149–163 RSC
.
- K.-Y. Chen, Y.-M. Cheng, C.-H. Lai, C.-C. Hsu, M.-L. Ho, G.-H. Lee and P.-T. Chou, J. Am. Chem. Soc., 2007, 129, 4534–4535 CrossRef CAS PubMed
.
- C. C. Hsieh, P. T. Chou, W.-T. Chuang, M.-W. Chung, J. Lee and T. Joo, J. Am. Chem. Soc., 2011, 133, 2932–2943 CrossRef CAS PubMed
.
- Y.-H. Hsu, Y.-A. Chen, H.-W. Tseng, Z. Zhang, J.-Y. Shen, W.-T. Chuang, T.-C. Lin, C.-S. Lee, W.-Y. Hung and B.-C. Hong, J. Am. Chem. Soc., 2014, 136, 11805–11812 CrossRef CAS PubMed
.
- K. M. Solntsev, O. Poizat, J. Dong, J. Rehault, Y. Lou, C. Burda and L. M. Tolbert, J. Phys. Chem. B, 2008, 112, 2700–2711 CrossRef CAS PubMed
.
- J. Dong, K. M. Solntsev, O. Poizat and L. M. Tolbert, J. Am. Chem. Soc., 2007, 129, 10084–10085 CrossRef CAS PubMed
.
- A. M. Bogdanov, A. S. Mishin, I. V. Yampolsky, V. V. Belousov, D. M. Chudakov, F. V. Subach, V. V. Verkhusha, S. Lukyanov and K. A. Lukyanov, Nat. Chem. Biol., 2009, 5, 459–461 CrossRef CAS PubMed
.
- A. M. Bogdanov, A. Acharya, A. V. Titelmayer, A. V. Mamontova, K. B. Bravaya, A. B. Kolomeisky, K. A. Lukyanov and A. I. Krylov, J. Am. Chem. Soc., 2016, 138, 4807–4817 CrossRef CAS PubMed
.
- C. R. S. Mooney, M. E. Sanz, A. R. McKay, R. J. Fitzmaurice, A. E. Aliev, S. Caddick and H. H. Fielding, J. Phys. Chem. A, 2012, 116, 7943–7949 CrossRef CAS PubMed
.
- D. A. Horke and J. R. R. Verlet, Phys. Chem. Chem. Phys., 2012, 14, 8511–8515 RSC
.
- Y. Toker, D. B. Rahbek, B. Klærke, A. V. Bochenkova and L. H. Andersen, Phys. Rev. Lett., 2012, 109, 128101 CrossRef CAS PubMed
.
- C. R. S. Mooney, M. A. Parkes, L. Zhang, H. C. Hailes, A. Simperler, M. J. Bearpark and H. H. Fielding, J. Chem. Phys., 2014, 140, 205103 CrossRef PubMed
.
- S. H. M. Deng, X.-Y. Kong, G. Zhang, Y. Yang, W.-J. Zheng, Z.-R. Sun, D.-Q. Zhang and X.-B. Wang, J. Phys. Chem. Lett., 2014, 5, 2155–2159 CrossRef CAS PubMed
.
- K. Bhaskaran-Nair, M. Valiev, S. H. M. Deng, W. A. Shelton, K. Kowalski and X.-B. Wang, J. Chem. Phys., 2015, 143, 224301 CrossRef PubMed
.
- D. Zuev, K. B. Bravaya, M. V. Makarova and A. I. Krylov, J. Chem. Phys., 2011, 135, 194304 CrossRef PubMed
.
- K. B. Bravaya, M. G. Khrenova, B. L. Grigorenko, A. V. Nemukhin and A. I. Krylov, J. Phys. Chem. B, 2011, 115, 8296–8303 CrossRef CAS PubMed
.
- K. B. Bravaya, B. L. Grigorenko, A. V. Nemukhin and A. I. Krylov, Acc. Chem. Res., 2012, 45, 265–275 CrossRef CAS PubMed
.
- C. R. S. Mooney, D. A. Horke, A. S. Chatterley, A. Simperler, H. H. Fielding and J. R. R. Verlet, Chem. Sci., 2013, 4, 921–927 RSC
.
- C. W. West, A. S. Hudson, S. L. Cobb and J. R. R. Verlet, J. Chem. Phys., 2013, 139, 071104 CrossRef PubMed
.
- A. V. Bochenkova and L. H. Andersen, Faraday Discuss., 2013, 163, 297–319 RSC
.
- A. V. Bochenkova, B. Klærke, D. B. Rahbek, J. Rajput, Y. Toker and L. H. Andersen, Angew. Chem., Int. Ed., 2014, 53, 9797–9801 CrossRef CAS PubMed
.
- C. W. West, J. N. Bull, A. S. Hudson, S. L. Cobb and J. R. R. Verlet, J. Phys. Chem. B, 2015, 119, 3982–3987 CrossRef CAS PubMed
.
- J. S. Paige, K. Y. Wu and S. R. Jaffrey, Science, 2011, 333, 642–646 CrossRef CAS PubMed
.
- K.-Y. Chen, Y.-M. Cheng, C.-H. Lai, C.-C. Hsu, M.-L. Ho, G.-H. Lee and P.-T. Chou, J. Am. Chem. Soc., 2007, 129, 4535 Search PubMed
.
- T. Chatterjee, M. Mandal, A. Das, K. Bhattacharyya, A. Datta and P. K. Mandal, J. Phys. Chem. B, 2016, 120, 3503–3510 CrossRef CAS PubMed
.
- A. R. McKay, M. E. Sanz, C. R. S. Mooney, R. S. Minns, E. M. Gill and H. H. Fielding, Rev. Sci. Instrum., 2010, 81, 123101 CrossRef CAS PubMed
.
- C. R. S. Mooney, M. A. Parkes, A. Iskra and H. H. Fielding, Angew. Chem., Int. Ed., 2015, 54, 5646–5649 CrossRef CAS PubMed
.
- M. A. Parkes, C. Phillips, M. J. Porter and H. H. Fielding, Phys. Chem. Chem. Phys., 2016, 18, 10329–10336 RSC
.
- G. A. Garcia, L. Nahon and I. Powis, Rev. Sci. Instrum., 2004, 75, 4989–4996 CrossRef CAS
.
-
J. Linderberg and Y. Öhrn, Propagators in Quantum Chemistry, John Wiley and Sons, Hoboken, NJ, 2014 Search PubMed
.
- E. Runge and E. K. U. Gross, Phys. Rev. Lett., 1984, 52, 997–1000 CrossRef CAS
.
-
M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, J. A. Montgomery Jr, J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers, K. N. Kudin, V. N. Staroverov, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth, P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels, Ö. Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski and D. J. Fox, Gaussian 09 Revision E.01, Gaussian Inc., Wallingford CT, 2009 Search PubMed
.
- C. Lee, W. Yang and R. G. Parr, Phys. Rev. B: Condens. Matter Mater. Phys., 1988, 37, 785–789 CrossRef CAS
.
- A. D. Becke, J. Chem. Phys., 1993, 98, 5648–5652 CrossRef CAS
.
- R. Krishnan, J. S. Binkley, R. Seeger and J. A. Pople, J. Chem. Phys., 1980, 72, 650–654 CrossRef CAS
.
- W. von Niessen, J. Schirmer and L. S. Cederbaum, Comput. Phys. Rep., 1984, 1, 57–125 CrossRef CAS
.
- T. Yanai, D. P. Tew and N. C. Handy, Chem. Phys. Lett., 2004, 393, 51–57 CrossRef CAS
.
- A. I. Krylov, Annu. Rev. Phys. Chem., 2008, 59, 433–462 CrossRef CAS PubMed
.
- R. A. Kendall, T. H. Dunning Jr and R. J. Harrison, J. Chem. Phys., 1992, 96, 6796–6806 CrossRef CAS
.
- Y. Shao, Z. Gan, E. Epifanovsky, A. T. Gilbert, M. Wormit, J. Kussmann, A. W. Lange, A. Behn, J. Deng, X. Feng, D. Ghosh, M. Goldey, P. R. Horn, L. D. Jacobson, I. Kaliman, R. Z. Khaliullin, T. Kuś, A. Landau, J. Liu, E. I. Proynov, Y. M. Rhee, R. M. Richard, M. A. Rohrdanz, R. P. Steele, E. J. Sundstrom, H. L. Woodcock III, P. M. Zimmerman, D. Zuev, B. Albrecht, E. Alguire, B. Austin, G. J. O. Beran, Y. A. Bernard, E. Berquist, K. Brandhorst, K. B. Bravaya, S. T. Brown, D. Casanova, C.-M. Chang, Y. Chen, S. H. Chien, K. D. Closser, D. L. Crittenden, M. Diedenhofen, R. A. DiStasio Jr, H. Do, A. D. Dutoi, R. G. Edgar, S. Fatehi, L. Fusti-Molnar, A. Ghysels, A. Golubeva-Zadorozhnaya, J. Gomes, M. W. Hanson-Heine, P. H. Harbach, A. W. Hauser, E. G. Hohenstein, Z. C. Holden, T.-C. Jagau, H. Ji, B. Kaduk, K. Khistyaev, J. Kim, J. Kim, R. A. King, P. Klunzinger, D. Kosenkov, T. Kowalczyk, C. M. Krauter, K. U. Lao, A. D. Laurent, K. V. Lawler, S. V. Levchenko, C. Y. Lin, F. Liu, E. Livshits, R. C. Lochan, A. Luenser, P. Manohar, S. F. Manzer, S.-P. Mao, N. Mardirossian, A. V. Marenich, S. A. Maurer, N. J. Mayhall, E. Neuscamman, C. M. Oana, R. Olivares-Amaya, D. P. O'Neill, J. A. Parkhill, T. M. Perrine, R. Peverati, A. Prociuk, D. R. Rehn, E. Rosta, N. J. Russ, S. M. Sharada, S. Sharma, D. W. Small, A. Sodt, T. Stein, D. Stück, Y.-C. Su, A. J. Thom, T. Tsuchimochi, V. Vanovschi, L. Vogt, O. Vydrov, T. Wang, M. A. Watson, J. Wenzel, A. White, C. F. Williams, J. Yang, S. Yeganeh, S. R. Yost, Z.-Q. You, I. Y. Zhang, X. Zhang, Y. Zhao, B. R. Brooks, G. K. Chan, D. M. Chipman, C. J. Cramer, W. A. Goddard III, M. S. Gordon, W. J. Hehre, A. Klamt, H. F. Schaefer III, M. W. Schmidt, C. D. Sherrill, D. G. Truhlar, A. Warshel, X. Xu, A. Aspuru-Guzik, R. Baer, A. T. Bell, N. A. Besley, J.-D. Chai, A. Dreuw, B. D. Dunietz, T. R. Furlani, S. R. Gwaltney, C.-P. Hsu, Y. Jung, J. Kong, D. S. Lambrecht, W. Liang, C. Ochsenfeld, V. A. Rassolov, L. V. Slipchenko, J. E. Subotnik, T. V. Voorhis, J. M. Herbert, A. I. Krylov, P. M. Gill and M. Head-Gordon, Mol. Phys., 2015, 113, 184–215 CrossRef CAS
.
- C. Møller and M. S. Plesset, Phys. Rep., 1934, 46, 618–622 Search PubMed
.
- J. Schirmer, Phys. Rev. A, 1982, 26, 2395–2416 CrossRef CAS
.
- A. B. Trofimov and J. Schirmer, J. Phys. B: At., Mol. Opt. Phys., 1995, 28, 2299–2324 CrossRef CAS
.
-
TURBOMOLE V6.3 2011, a development of University of Karlsruhe and Forschungszentrum Karlsruhe GmbH, 1989-2007, TURBOMOLE GmbH, since 2007, http://www.turbomole.com Search PubMed
.
- E. Epifanovsky, I. Polyakov, B. Grigorenko, A. Nemukhin and A. I. Krylov, J. Chem. Theory Comput., 2009, 5, 1895–1906 CrossRef CAS PubMed
.
- K. B. Bravaya and A. I. Krylov, J. Phys. Chem. A, 2013, 117, 11815–11822 CrossRef CAS PubMed
.
- J. Lazzari-Dean, A. I. Krylov and K. B. Bravaya, Int. J. Quantum Chem., 2015, 115, 1258–1264 CrossRef CAS
.
- D. Ghosh, A. Acharya, S. C. Tiwari and A. I. Krylov, J. Phys. Chem. B, 2012, 116, 12398–12405 CrossRef CAS PubMed
.
- K. M. Solntsev, D. Ghosh, A. Amador, M. Josowicz and A. I. Krylov, J. Phys. Chem. Lett., 2011, 2, 2593–2597 CrossRef CAS
.
- Y.-H. Hsu, Y.-A. Chen, H.-W. Tseng, Z. Zhang, J.-Y. Shen, W.-T. Chuang, T.-C. Lin, C.-S. Lee, W.-Y. Hung, B.-C. Hong, S.-H. Liu and P.-T. Chou, J. Am. Chem. Soc., 2014, 136, 11805–11812 CrossRef CAS PubMed
.
-
A. V. Bochenkova, private communication.
- I.-R. Lee, W. Lee and A. H. Zewail, Proc. Natl. Acad. Sci. U. S. A., 2006, 103, 258–262 CrossRef CAS PubMed
.
Footnote |
† Electronic supplementary information (ESI) available: Details of chromophore synthesis and NMR data; potential energy scans along bridge and CCOH torsion; geometrical data of optimised structures; ionisation potentials with state characters; excited state energies and characters; coordinates of optimised geometries. See DOI: 10.1039/c6sc03833f |
|
This journal is © The Royal Society of Chemistry 2017 |