Four resonance structures elucidate double-bond isomerisation of a biological chromophore.

Photoinduced double-bond isomerisation of the chromophore of photoactive yellow protein (PYP) is highly sensitive to chromophore-protein interactions. On the basis of high-level ab initio calculations, we scrutinise the effect of hydrogen bonds on the photophysical and photochemical properties of the chromophore. We identify four resonance structures - two closed-shell and two biradicaloid - that elucidate the electronic structure of the ground and first excited states involved in the isomerisation process. Changing the relative energies of the resonance structures by hydrogen-bonding interactions tunes all photochemical properties of the chromophore in an interdependent manner. Our study sheds new light on the role of the chromophore electronic structure in tuning in photosensors and fluorescent proteins.

Caged-Electron States in Endohedral Li Fullerenes.

By employing large-scale high-level EA-EOM-CCSD calculations, we have computed and analyzed the low-lying states of neutral Li@C60. Apart from one state, all states are found to be charge-separated states of the type Li+@C60-. The new state is the first reported non-charge-separated state in endohedral alkali fullerenes. This caged-electron state is analyzed in detail. Arguments are given that in larger highly symmetric endohedral fullerenes the caged-electron state can be the electronic ground state of the system. HF and DFT calculations on Li@C180 indeed find that the caged-electron state is the ground state and that in its equilibrium geometry Li sits at the center of the cage. Applications are mentioned.

Charge separated states of endohedral fullerene Li@C20.

We report on high-level coupled-cluster calculations of electronic states of the neutral endohedral fullerene Li@C20. All computed states of neutral Li@C20 are found to be the charge separated states of the Li+@C20 - type. Using the state-of-the-art EA-EOM-CCSD method, we found that neutral Li@C20 (D3d) possesses several valence and superatomic charge separated states with considerable electron binding energies, the strongest bound state of Li+@C20 - being the 12Eu state (6.73 eV). The valence charge separated states correspond to two sets of states of C20 -. The states 12Eu, 12A2u, 22Eu, and 22A2u correspond to the respective bound states of C20 -, and the states 22A2g, 12Eg, 12A1g, and 42Eu correspond to the unbound states of C20 -. There are eight superatomic states with electron binding energy higher than 1.0 eV, being much stronger bound than the single weakly bound superatomic state of the parent fullerene anion. The analysis of the radial density distribution of the excess electron on the carbon cage indicates the important role of the inner part of the superatomic states in forming the charge separated states.

Bound electronic states of the smallest fullerene C20- anion.

We report on high-level coupled-cluster calculations for the anion states of the smallest fullerene C20. Using the state-of-the-art EA-EOM-CCSD method we revealed that the C20- anion has five bound electronic states at the C20 neutral ground-state D3d equilibrium configuration. These are two pairs of 2Eu and 2A2u states and one 2A1g state. The binding energies vary from 2.05 eV for the most bound 2Eu state to <1 meV for the 2A1g state. An analysis in terms of radial and angular density distribution of the excess electron revealed that the two 2Eu/2A2u pairs are valence-like states while 2A1g corresponds to a super-atomic-like (SAMO) state. The valence states of the first 2Eu/2A2u pair were found to be of p-type whereas those of the second pair are of hybrid sp-type. We have also applied a simple model to understand the binding of the excess electron in C20-. The model confirms the existence of only one SAMO state of the s-type for C20-. At the same time, it overestimates the number of the bound valence states of C20-.

How many bound valence states does the C60⁻ anion have?

We report on unprecedentedly large coupled cluster calculations for the C60(-) anion, and on a heuristic model uncovering the valence states of C60(-) that allow the resolution of the headlined question. Our results convincingly demonstrate that C60(-) possesses as many as four bound valence states: (2)T1u, (2)T1g, (2)T2u and (2)Hg. Our findings reconcile previous controversies regarding the existence of the bound (2)T2u and (2)Hg states. For all bound states of C60(-) we present an analysis of the radial and angular distributions of the excess electron, which reveals some unique properties of the valence states. Some interesting features of the introduced model are analyzed and discussed.

Unveiling the mechanism of photoinduced isomerization of the photoactive yellow protein (PYP) chromophore.

A detailed theoretical analysis, based on extensive ab initio second-order approximate coupled cluster calculations, has been performed for the S1 potential energy surface (PES) of four photoactive yellow protein (PYP) chromophore derivatives that are hydrogen bonded with two water molecules and differ merely in the carbonyl substituent. The main focus is put on contrasting the isomerization properties of these four species in the S1 excited state, related to torsion around the chromophore's single and double carbon-carbon bonds. The analysis provides evidence of the different isomerization behavior of these four chromophore complexes, which relates to the difference in their carbonyl substituents. While a stable double-bond torsion pathway exists on the S1 PES of the chromophores bearing the -O-CH3 and -NH2 substituents, this is not the case for the -S-CH3 and -CH3 substituted species. The presence of the -S-CH3 group leads to a strong instability of the chromophore with respect to the single-bond twist, whereas in the case of the -CH3 substituent a crossing of the S1 and S2 PESs occurs, which perturbs the pathway. Based on this analysis, the key factors that support the double-bond torsion have been identified. These are (i) the hydrogen bonds at the phenolic oxygen of the chromophore, (ii) the weak electron-acceptor character of the carbonyl group, and (iii) the ethylene-like pattern of the torsion in the beginning of the process. Our results suggest that the interplay between these factors determines the chromophore's isomerization in the solvent environment and in the native PYP environment.

Ionization satellites of the ArHe dimer.

Ionization satellites are key ingredients in the control of post ionization processes such as molecular dissociation and interatomic Coulombic decay. Here, using the high-level ab initio method of multi-reference configuration interaction up to triple excitations, we study the potential energy curves (PECs) of the ionization satellites of the ArHe dimer. With this model system, we demonstrate that the simple model used in alkaline earth metal and rare gas complexes to describe the satellites as a Rydberg electron moving on top of a dicationic core does not fully hold for the rare gas clusters. The more complex valence structure in the rare gas atom leads to the mixing of different electronic configurations of the dimer. This prevents one from assigning a single dicationic parent state to some of the ionization satellites. We further analyze the structure of the different PECs, demonstrating how the density of the Rydberg electron is reflected in the structure of the PEC wherever the simple model is applicable.

All for one and one for all: accommodating an extra electron in C60.

Much like the neutral C60 fullerene, the C60(-) anion possesses certain unique properties which have attracted a great deal of research. One of these special properties, only recently fully uncovered, is that the C60(-) anion supports a substantial number of electronically stable excited states in contrast to other molecular anions with comparable electron affinity. In this work, we clarify how the C60(-) anion can support so many stable states by analyzing the radial and angular distributions of the excess electron bound to the anion. The analysis is based on ab initio calculations which are by far the most accurate on the C60(-) anion to date. Surprisingly, the radial distributions are highly similar for states of very different binding energies and the analysis stresses the importance of angular correlation in binding the excess electron. We further analyze the effect of the single excess electron on the electrons of the underlying neutral molecule. We demonstrate how this substantially modifies the actual distribution of the excess charge by shifting the underlying electron density. Implications of these findings are discussed.

Photoinduced isomerization of the photoactive yellow protein (PYP) chromophore: interplay of two torsions, a HOOP mode and hydrogen bonding.

We report on a detailed theoretical analysis, based on extensive ab initio calculations at the CC2 level, of the S(1) potential energy surface (PES) of the photoactive yellow protein (PYP) chromophore. The chromophore's photoisomerization pathway is shown to be fairly complex, involving an intimate coupling between single-bond and double-bond torsions. Furthermore, these torsional modes are shown to couple to a third coordinate of hydrogen out-of-plane (HOOP) type whose role in the isomerization is here identified for the first time. In addition, it is demonstrated that hydrogen bonding at the phenolate moiety of the chromophore can hinder the single-bond torsion and thus facilitates double-bond isomerization. These results suggest that the interplay between intramolecular factors and H-bonding determines the isomerization in native PYP.

An experimental and theoretical core-level study of tautomerism in guanine.

The core level photoemission and near edge X-ray photoabsorption spectra of guanine in the gas phase have been measured and the results interpreted with the aid of high level ab initio calculations. Tautomers are clearly identified spectroscopically, and their relative free energies and Boltzmann populations at the temperature of the experiment (600 K) have been calculated and compared with the experimental results and with previous calculations. We obtain good agreement between experiment and the Boltzmann weighted theoretical photoemission spectra, which allows a quantitative determination of the ratio of oxo to hydroxy tautomer populations. For the photoabsorption spectra, good agreement is found for the C 1s and O 1s spectra but only fair agreement for the N 1s edge.

Tautomerism in cytosine and uracil: an experimental and theoretical core level spectroscopic study.

The O, N, and C 1s core level photoemission spectra of the nucleobases cytosine and uracil have been measured in the vapor phase, and the results have been interpreted via theoretical calculations. Our calculations accurately predict the relative binding energies of the core level features observed in the experimental photoemission results and provide a full assignment. In agreement with previous work, a single tautomer of uracil is populated at 405 K, giving rise to relatively simple spectra. At 450 K, three tautomers of cytosine, one of which may consist of two rotamers, are identified, and their populations are determined. This resolves inconsistencies between recent laser studies of this molecule in which the rare imino-oxo tautomer was not observed and older microwave spectra in which it was reported.

High-resolution spectroscopy of methyl 4-hydroxycinnamate and its hydrogen-bonded water complex.

The electronic structure and spectroscopic properties of the lower excited singlet states of methyl 4-hydroxycinnamate, a model for the chromophore of the photoactive yellow protein in neutral form, have been investigated using various high-resolution gas-phase spectroscopic techniques and quantum-chemical calculations. The experiments show that under our experimental conditions the molecule can adopt four conformations with similar spectroscopic properties. From the detailed assignment of the vibrationally active modes in excitation and emission spectra, it is concluded that the S(1) and S(2) states should be assigned to the V' and V pipi* states that are characterized by, respectively, small and large contributions of the HOMO --> LUMO excitation. We find that complexation with a single water molecule affects the spectroscopic properties of methyl 4-hydroxycinnamate considerably in terms of stabilization of the lowest excited singlet state but in particular with respect to the transition intensities. The latter observation is tentatively interpreted as being caused by an increase in the oscillator strength of the respective electronic transition as well as by a rise/removal of conical intersections with the pisigma* state.

Electronic structure of the PYP chromophore in its native protein environment.

We report on supermolecular ab initio calculations which clarify the role of the local amino acid environment in determining the unique electronic structure properties of the photoactive yellow protein (PYP) chromophore. The extensive ab initio calculations, at the level of the CC2 and EOM-CCSD methods, allow us to explicitly address how the interactions between the deprotonated p-coumaric thio-methyl ester (pCTM-) chromophore and the surrounding amino acids act together to create a specifically stabilized pCTM- species. Particularly noteworthy is the role of the Arg52 amino acid in stabilizing the chromophore against autoionization, and the role of the Tyr42 and Glu46 amino acids in determining the hydrogen-bonding properties that carry the dominant energetic effects.

Combined experimental-theoretical study of the lower excited singlet states of paravinyl phenol, an analog of the paracoumaric acid chromophore.

The low-lying excited singlet states of paravinyl phenol (pVP) are investigated experimentally and theoretically paying attention to their similarity to excited states of paracoumaric acid, the chromophore of the photoactive yellow protein (PYP). Resonance enhanced multiphoton ionization and laser induced fluorescence spectroscopic techniques are employed to obtain supersonically cooled, vibrationally resolved excitation and emission spectra related to the lowest (1)A'(V') excited state of pVP. Comprehensive analyses of the spectral structures are carried out by means of the equation-of-motion coupled cluster singles and doubles and time dependent density functional theory methods in combination with the linear vibronic coupling model and Franck-Condon calculations. The assignments of the spectral patterns are given, mostly in terms of excitations of totally symmetric modes. Weak activity of the non-totally-symmetric modes indicates low probability of photochemical processes in the Franck-Condon region of the (1)A'(V') state. The second (1)A'(V) and third (1)A" (Ryd) excited states of pVP are characterized with regard to their electronic structure, properties, and effects of geometry relaxations. The lengthening of the double bond relevant to the trans-cis isomerization of the PYP chromophore is found for the (1)A'(V) state. A possibility of photochemical processes and strong vibronic interactions in this state can be expected. The theoretical results for the (1)A"(Ryd) state predict that dissociation with respect to the O-H bond is possible.

Impact of sulfur vs oxygen on the low-lying excited states of trans-p-coumaric acid and trans-p-coumaric thio acid.

The low-lying excited singlet states of trans-p-coumaric acid (CA) and trans-p-coumaric thio acid (CTA) are investigated in view of characterizing the chromophore of the photoactive yellow protein (PYP), with particular regard to the impact of sulfur on the chromophore's electronic structure. The comparative ab initio study, performed with the highly accurate EOM-CCSD method, shows that the electronic state ordering upon vertical excitation and following in-plane geometry relaxation indeed depends in a very sensitive fashion on the presence of either sulfur or oxygen. The study identifies three relevant excited singlet states, two of which are of pi-pi type while the third state is of n-pi character. The study highlights the role of the latter n-pi state which is shown to be the lowest-lying excited state of CTA at all in-plane geometries under consideration, whereas this is not the case for CA.