Search

Photoinduced charge transfer (CT) excited states and their relaxation mechanisms can be highly interdependent on the environment effects and the consequent changes in the electronic density. Providing a molecular interpretation of the ultrafast (subpicosecond) interplay between initial photoexcited states in such dense electronic manifolds in condensed phase is crucial for improving and understanding such phenomena. Real-time time-dependent density functional theory is here the method of choice to observe the charge density, explicitly propagated in an ultrafast time domain, along with all time-dependent properties that can be easily extracted from it. A designed protocol of analysis for real-time electronic dynamics to be applied to time evolving electronic density related properties to characterize both in time and in space CT dynamics in complex systems is here introduced and validated, proposing easy to be read cross-correlation maps. As case studies to test such tools, we present the photoinduced charge-transfer electronic dynamics of 5-benzyluracil, a mimic of nucleic acid/protein interactions, and the metal-to-ligand charge-transfer electronic dynamics in water solution of [Ru(dcbpy)2(NCS)2]4-, dcbpy = (4,4'-dicarboxy-2,2'-bipyridine), or "N34-", a dye sensitizer for solar cells. Electrostatic and explicit ab initio treatment of solvent molecules have been compared in the latter case, revealing the importance of the accurate modeling of mutual solute-solvent polarization on CT kinetics. We observed that explicit quantum mechanical treatment of solvent slowed down the charge carriers mobilities with respect to the gas-phase. When all water molecules were modeled instead as simpler embedded point charges, the electronic dynamics appeared enhanced, with a reduced hole-electron distance and higher mean velocities due to the close fixed charges and an artificially increased polarization effect. Such analysis tools and the presented case studies can help to unveil the influence of the electronic manifold, as well as of the finite temperature-induced structural distortions and the environment on the ultrafast charge motions.

The excited state proton transfer (ESPT) reaction from the photoacid 8-hydroxypyrene-1,3,6-trisulfonic acid (HPTS or pyranine) to an acetate molecule has been investigated in explicit aqueous solution via excited state ab initio molecular dynamics simulations based on hybrid quantum/molecular mechanics (QM/MM) potentials. In all the trajectories, the direct proton transfer has been observed in the excited state within 1 ps. We find that the initial structural configuration extracted from the ground state distribution strongly affects the ESPT kinetics. Indeed, the relative orientation of the proton donor-acceptor pair and the presence of a water molecule hydrogen bonded to the phenolic acid group of the pyranine are the key factors to facilitate the ESPT. Furthermore, we analyze the vibrational fingerprints of the ESPT reaction, reproducing the blue shift of the acetate CO stretching (COac), from 1666 to 1763 cm-1 testifying the transformation of acetate to acetic acid. Finally, our findings suggest that the acetate CC stretching (CCac) is also sensitive to the progress of the ESPT reaction. The CCac stretching is indeed ruled by the two vibrational modes (928 and 1426 cm-1), that in the excited state are alternately activated when the proton is shared or bound to the donor/acceptor, respectively.

Fluorescent Proteins are widely studied for their multiple applications in technological and biotechnological fields. Despite this, they continue to represent a challenge in terms of a complete understanding of all the non-equilibrium photo-induced processes that rule their properties. In this context, a theoretical-computational approach can support experimental results in unveiling and understanding the processes taking place after electronic excitation. A non-standard cyan fluorescent protein, psamFP488, is characterized by an absorption maximum that is blue-shifted in comparison to other cyan fluorescent proteins. This protein is characterized by an extended Stokes shift and an ultrafast (170 fs) excited state proton transfer. In this work, a theoretical-computational study, including excited state ab initio dynamics, is performed to help understanding the reaction mechanism and propose new hypotheses on the role of the residues surrounding the chromophore. Our results suggest that the proton transfer could be indirect toward the acceptor (Glu167) and involves other residues surrounding the chromophore, despite the ultrafast kinetics.

Environmental effects can drastically influence the optical properties and photoreactivity of molecules, particularly in the presence of polar and/or protic solvents. In this work we investigate a negatively charged Ru(ii) complex, [Ru(dcbpy)2(NCS)2]4- [dcbpy = (4,4′-dicarboxy-2,2′-bipyridine)], in water solution, since this system belongs to a broader class of transition-metal compounds undergoing upon photo-excitation rapid and complex charge transfer (CT) dynamics, which can be dictated by structural rearrangement and solvent environment. Ab initio molecular dynamics (AIMD) relying on a hybrid quantum/molecular mechanics scheme is used to probe the equilibrium microsolvation around the metal complex in terms of radial distribution functions of the main solvation sites and solvent effects on the overall equilibrium structure. Then, using our AIMD-based generalized normal mode approach, we investigate how the ligand vibrational spectroscopic features are affected by water solvation, also contributing to the interpretation of experimental Infra-Red spectra. Two solvation sites are found for the ligands: the sulfur and the oxygen sites can interact on average with ∼4 and ∼3 water molecules, respectively, where a stronger interaction of the oxygen sites is highlighted. On average an overall dynamic distortion of the C2 symmetric gas-phase structure was found to be induced by water solvation. Vibrational analysis reproduced experimental values for ligand symmetric and asymmetric stretchings, linking the observed shifts with respect to the gas-phase to a complex solvent distribution around the system. This is the groundwork for future excited-state nuclear and electronic dynamics to monitor non-equilibrium processes of CT excitation in complex environments, such as exciton migration in photovoltaic technologies. This journal is

Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical-computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute-solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural-vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.

Light induces non-equilibrium time evolving molecular phenomena. The computational modeling of photo-induced processes in large systems, embedded in complex environments (i.e., solutions, proteins, materials), demands for a quantum and statistical mechanic treatment to achieve the required accuracy in the description of both the excited-state energy potentials and the choice of the initial conditions for dynamical simulations. On the other hand, the theoretical investigation on the atomistic scale of times and sizes of the ultrafast photo-induced reactivity and non-equilibrium relaxation dynamics right upon excitation requests tailored computational protocols. These methods often exploit hierarchic computation schemes, where a large part of the degrees of freedom are required to be treated explicitly to achieve the right accuracy. Additionally, part of the explicit system needs to be treated at ab initio level, where density functional theory, using hybrid functionals, represents a good compromise between accuracy and computational cost, when proton transfers, non-covalent interactions, and hydrogen bond dynamics play important roles. Thus, the modeling strategies presented in this review stress the importance of hierarchical quantum/molecular mechanics with effective non-periodic boundary conditions and efficient phase-sampling schemes to achieve chemical accuracy in ultrafast time-resolved spectroscopy and photo-induced phenomena. These approaches can allow explicit and accurate treatment of molecule/environment interactions, including also the electrostatic and dispersion forces of the bulk. At the same time, the specificities of the different case studies of photo-induced phenomena in solutions and biological environments are highlighted and discussed, with special attention to the computational and modeling challenges.

By coupling an enhanced sampling algorithm with an orbital-localized variant of Car-Parrinello molecular dynamics, the so-called atomic centered density matrix propagation model, we reconstruct the free energy profiles along reaction pathways using different density functional approximations (DFAs) ranging from locals to hybrids. In particular, we compare the computed free energy barrier height of proton transfer (PT) reactions to those obtained by a more traditional static approach, based on the intrinsic reaction coordinate (IRC), for two case systems, namely malonaldehyde and formic acid dimer. The obtained results show that both the IRC profiles and the potentials of mean force, derived from biased dynamic trajectories, are very sensitive to the density functional approximation applied. More precisely, we observe that, with the notable exception of M06-L, local density functionals always strongly underestimate the reaction barrier heights. More generally, we find that also the shape of the free energy profile is very sensitive to the density functional choice, thus highlighting the effect, often neglected, that the choice of DFA has also in the case of dynamics simulations.

Charge-transfer dynamics and interligand electron transfer (ILET) phenomena play a pivotal role in dye-sensitizers, mostly represented by the Ru-based polypyridyl complexes, for TiO2 and ZnO-based solar cells. Starting from metal-to-ligand charge-transfer (MLCT) excited states, charge dynamics and ILET can influence the overall device efficiency. In this letter, we focus on N34– dye ( [Ru(dcbpy)2(NCS)2]4–, dcbpy = 4,4′-dicarboxy-2,2′-bipyridine) to provide a first direct observation with high time resolution (

The interplay between light absorption and the molecular environment has a central role in the observed photophysics of a wide range of photoinduced chemical and biological phenomena. The understanding of the interplay between vibrational and electronic transitions is the focus of this work, since it can provide a rationale to tune the optical properties of charge transfer (CT) materials used for technological applications. A clear description of these processes poses a nontrivial challenge from both the theoretical and experimental points of view, where the main issue is how to accurately describe and probe drastic changes in the electronic structure and the ultrafast molecular relaxation and dynamics. In this work we focused on the intermolecular CT reaction that occurs upon photon absorption in a π-stacked model system in dichloromethane solution, in which the 1-chloronaphthalene (1ClN) acts as the electron donor and tetracyanoethylene (TCNE) is the electron acceptor. Density functional theory calculations have been carried out to characterize both the ground-state properties and more importantly the low-lying CT electronic transition, and excellent agreement with recently available experimental results [Mathies, R. A.; et al. J. Phys. Chem. A 2018, 122 (14), 3594] was obtained. The minima of the ground state and first singlet excited state have been accurately characterized in terms of spatial arrangements and vibrational Raman frequencies, and the CT natures of the first two low-lying electronic transitions in the absorption spectra have been addressed and clarified too. Finally, by modeling the possible coordination sites of the TCNE electron acceptor with respect to monovalent ions (Na+, K+) in an implicit solution of acetonitrile, we find that TCNE can accommodate a counterion in two different arrangements, parallel and orthogonal to the C═C axis, leading to the formation of a contact ion pair. The nature of the counterion and its relative position entail structural modifications of the TCNE radical anion, mainly the central C═C and C≡N bonds, compared to the isolated case. An important red shift of the C═C stretching frequency was observed when the counterion is orthogonal to the double bond, to a greater extent for Na+. On the contrary, in the second case, where the counterion ion lies along the internuclear C═C axis, we find that K+ polarizes the electron density of the double bond more, resulting in a greater red shift than with Na+.

Environmental effects can drastically influence the optical properties and photoreactivity of molecules, particularly in the presence of polar and/or protic solvents. In this work we investigate a negatively charged Ru(II) complex, [Ru(dcbpy)

In the search for new therapeutic strategies to contrast SARS-CoV-2, we here studied the interaction of polydatin (PD) and resveratrol (RESV)—two natural stilbene polyphenols with manifold, well known biological activities—with Spike, the viral protein essential for virus entry into host cells, and ACE2, the angiotensin-converting enzyme present on the surface of multiple cell types (including respiratory epithelial cells) which is the main host receptor for Spike binding. Molecular Docking simulations evidenced that both compounds can bind Spike, ACE2 and the ACE2:Spike complex with good affinity, although the interaction of PD appears stronger than that of RESV on all the investigated targets. Preliminary biochemical assays revealed a significant inhibitory activity of the ACE2:Spike recognition with a dose-response effect only in the case of PD.

In this work, excited state ab initio molecular dynamics together with a time resolved vibrational analysis is employed to shed light on the vibrational photoinduced dynamics of a well-known diarylethene molecule experiencing a ring opening reaction upon electronic excitation. The photoreactivity of diarylethenes is recognized to be controlled by a non-adiabatic intersection point between the ground and the first excited state surfaces. The computation of an energy scan, along a suitable reaction coordinate, allows us to identify the region of potential energy surfaces in which the ground (S0) and the first excited (S1) state are well separated. The adiabatic sampling of that region in S1 shows that in the first 3 picoseconds, the central CC bond, which is subject to break, oscillates in an antiphase with respect to the energy gap ΔE(S1 - S0). A multiresolution analysis based on the wavelet transform was then applied to the structural parameters extracted from the excited state dynamics. The wavelet maps show characteristic oscillations of the frequencies, mainly CC stretching and CCC bending localized on the central 4-ring moiety. Moreover, we have identified the main frequency (methyl wagging motion) involved in the modulation of these oscillations. The anharmonic coupling within a group of vibrational modes was therefore highlighted, in good agreement with experimental evidence. For the first time, a quantitative analysis of time resolved signals from a wavelet transform/ab initio molecular dynamics approach was performed.

We present electronic structure methods to unveil the non-radiative pathways of photoinduced charge transfer (CT) reactions that play a main role in photophysics and light harvesting technologies. A prototypical π-stacked molecular complex consisting of an electron donor (1-chloronaphthalene, 1ClN) and an electron acceptor (tetracyanoethylene, TCNE) was investigated in dichloromethane solution for this purpose. The characterization of TCNE:π:1ClN in both its equilibrium ground and photoinduced low-lying CT electronic states was performed by using a reliable and accurate theoretical–computational methodology exploiting ab initio molecular dynamics simulations. The structural and vibrational time evolution of key vibrational modes is found to be in excellent agreement with femtosecond stimulated Raman spectroscopy experiments [R. A. Mathies et al., J. Phys. Chem. A, 2018, 122, 14, 3594], unveiling a correlation between vibrational fingerprints and electronic properties. The evaluation of nonadiabatic coupling matrix elements along generalized normal modes has made possible the interpretation on the molecular scale of the activation of nonradiative relaxation pathways towards the ground electronic state. In particular, two low frequency vibrational modes such as the out of plane bending and dimer breathing and the TCNE central C Created by potrace 1.16, written by Peter Selinger 2001-2019 C stretching play a prominent role in relaxation phenomena from the electronic CT state to the ground state one., We present electronic structure methods to unveil the non-radiative pathways of photoinduced charge transfer (CT) reactions that play a main role in photophysics and light harvesting technologies.

In this work, we simulate the excited state proton transfer (ESPT) reaction involving the pyranine photoacid and an acetate molecule as proton acceptor, connected by a bridge water molecule. We employ ab initio molecular dynamics combined with an hybrid quantum/molecular mechanics (QM/MM) framework. Furthermore, a time-resolved vibrational analysis based on the wavelet-transform allows one to identify two low frequency vibrational modes that are fingerprints of the ESPT event: a ring wagging and ring breathing. Their composition suggests their key role in optimizing the structure of the proton donor-acceptor couple and promoting the ESPT event. We find that the choice of the QM/MM partition dramatically affects the photoinduced reactivity of the system. The QM subspace was gradually extended including the water molecules directly interacting with the pyranine-water-acetate system. Indeed, the ESPT reaction takes place when the hydrogen bond network around the reactive system is taken into account at full QM level.

Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical-computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute-solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural-vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.

We present a novel time-resolved vibrational analysis for studying photoinduced nuclear relaxation. Generalized modes velocities are defined from ab initio molecular dynamics and wavelet transformed, providing the time localization of vibrational signals in the electronic excited state. The photoexcited pyranine in aqueous solution is presented as a case study. The transient and sequential activation of the simulated vibrational signals is in good agreement with vibrational dynamics obtained from femtosecond stimulated Raman spectroscopy data.

The development of efficient artificial leaves relies on the subtle combination of molecular assemblies able to absorb sunlight, converting light energy into electrochemical potential energy and finally transducing it into accessible chemical energy. The electronic design of these charge transfer molecular machines is crucial to build a complex supramolecular architecture for the light energy conversion. Here, we present an ab initio simulation of the whole decay pathways of a recently proposed artificial molecular reaction center. A complete structural and energetic characterization has been carried out with methods based on density functional theory, its time-dependent version, and a broken-symmetry approach. On the basis of our findings we provide a revision of the pathway only indirectly postulated from an experimental point of view, along with unprecedented and significant insights on the electronic and nuclear structure of intramolecular charge-separated states, which are fundamental for the application of this molecular assembly in photoelectrochemical cells. Importantly, we unravel the molecular driving forces of the various charge transfer steps, in particular those leading to the proton-coupled electron transfer final product, highlighting key elements for the future design strategies of such molecular assays.

Advances in time resolved vibrational spectroscopy techniques provided a new stimulus for understanding the transient molecular dynamics triggered by the electronic excitation. The detailed interpretation of such time dependent spectroscopic signals is a challenging task from both experimental and theoretical points of view. We simulated and analyzed the transient photorelaxation of the pyranine photoacid in aqueous solution, with special focus on structural parameters and low frequency skeleton modes that are possibly preparatory for the photoreaction occurring at later time, as suggested by experimental spectroscopic studies. To this aim, we adopted an accurate computational protocol that combines excited state ab-initio molecular dynamics within a hybrid quantum mechanical/molecular mechanics framework and a time resolved vibrational analysis based on the wavelet transform. According to our results, the main nuclear relaxation on the excited potential energy surface is completed in about 500 fs, in agreement with experimental data. The rearrangement of C-C bonds occurs according to a complex vibrational dynamics, showing oscillatory patterns that are out of phase and modulated by modes below 200 cm-1. We also analyzed in both the ground and the excited state the evolution of some structural parameters involved in excited state proton transfer reaction, namely those involving the pyranine and the water molecule hydrogen bonded to the phenolic O-H group. Both the hydrogen bond distance and the intermolecular orientation are optimized in the excited state, resulting in a tighter proton donor-acceptor couple. Indeed, we found evidence that collective low frequency skeleton modes, such as the out of plane wagging at 108 cm-1and the deformation at 280 cm-1are photoactivated since the ultrafast part of the relaxation and modulate the pyranine-water molecule rearrangement, favouring the step preparatory for the photoreactivity.

Synthetic receptors for biomacromolecules lack the supramolecular self-assembly behavior typical of biological systems. Here we propose a new method for the preparation of protein imprinted polymers based on the specific interaction of a peptide multi-functional block with a protein target. This peptide block contains a protein-binding peptide domain, a polymerizable moiety at the C-terminus and an environment-sensitive fluorescent molecule at the N-terminus. The method relies on a preliminary step consisting of peptide/protein supramolecular assembly, followed by copolymerization with the most common acrylate monomers (acrylamide, acrylic acid and bis-acrylamide) to produce a protein imprinted hydrogel polymer. Such a peptide block can function as an active assistant recognition element to improve affinity, and guarantees its effective polymerization at the protein/cavity interface, allowing for proper placement of a dye. As a proof of concept, we chose Bovine Serum Albumin (BSA) as the protein target and built the peptide block around a BSA binding dodecapeptide, with an allyl group as the polymerizable moiety and a dansyl molecule as the responsive dye. Compared to conventional approaches these hydrogels showed higher affinity (more than 45%) and imprinted sensitivity (about twenty fold) to the target, with a great BSA selectivity with respect to ovalbumin (alpha = 1.25) and lysozyme (alpha = 6.02). Upon protein binding, computational and experimental observations showed a blue shift of the emission peak (down to 440 nm) and an increase of fluorescence emission (twofold) and average lifetime (Delta(tau) = 4.3 ns). Such an approach generates recognition cavities with controlled chemical information and represents an a priori method for self-responsive materials. Provided a specific peptide and minimal optimization conditions are used, such a method could be easily implemented for any protein target.

We simulated an excited state proton transfer in green fluorescent protein by excited state ab initio dynamics, and examined the reaction mechanism in both the time and the frequency domain through a multi resolution wavelet analysis. This original approach allowed us, for the first time, to directly compare the trends of photoactivated vibrations to femtosecond stimulated Raman spectroscopy results, and to give an unequivocal interpretation of the role played by low frequency modes in promoting the reaction. We could attribute the main driving force of the reaction to an important photoinduced softening of the ring-ring orientational motion of the chromophore, thus permitting the tightening of the hydrogen bond network and the opening of the reaction pathway. We also found that both the chromophore (in terms of its inter-ring dihedral angle and phenolic C-O and imidazolinone C-N bond distances) and its pocket (in terms of the inter-molecular oxygen's dihedral angle of the chromophore pocket) relaxations are modulated by low frequency (about 120 cm-1) modes involving the oxygen atoms of the network. This is in agreement with the femtosecond Raman spectroscopy findings in the time-frequency domain. Moreover, the rate in proximity to the Franck Condon region involves a picosecond time scale, with a significant influence from fluctuations of nearby hydrogen bonded residues such as His148. This approach opens a new scenario with ab initio simulations as routinely used tools to understand photoreactivity and the results of advanced time resolved spectroscopy techniques.

Among polyphenols, trans-resveratrol (tRES) and trans-polydatin (tPD) exert multiple biological effects, particularly antioxidant and antiproliferative. In this work, we have investigated the interaction of tPD with three cancer-related DNA sequences able to form G-quadruplex (G4) structures, as well as with a model duplex, and compared its behaviour with tRES. Interestingly, fluorescence analysis evidenced the ability of tPD to bind all the studied DNA systems, similarly to tRES, with tRES displaying a higher ability to discriminate G4 over duplex with respect to tPD. However, neither tRES nor tPD produced significant conformational changes of the analyzed DNA upon binding, as determined by CD-titration analysis. Computational analysis and biological data confirmed the biophysical results: indeed, molecular docking evidenced the stronger interaction of tRES with the promoter of c-myc oncogene, and immunoblotting assays revealed a reduction of c-myc expression, more effective for tRES than tPD. Furthermore, in vitro assays on melanoma cells proved that tPD was able to significantly reduce telomerase activity, and inhibit cell proliferation, with tRES producing higher effects than tPD.

We simulated for the first time the oxygen-oxygen bond formation in a synthetic calcium-tetra manganese complex recently developed by Zhang and co-workers. In spite of promising structural similarities to the native oxygen evolving complex (OEC) in Photosystem II, several uncertainties on the mimic stability in water and on its potential catalytic activity still persist. Here, we characterized at density functional theory level the electronic and structural features of the Sn states of the complex, along with the oxygen-oxygen bond formation reaction, proposing a reasonable model for the hydrate complex. As a main finding, both the synthetic compound and the natural OEC show very close energetic barriers for the oxo-oxyl coupling process, suggesting that key electronic features of the natural OEC reactivity are well reproduced. This result strongly encourages the use of this synthetic complex in combination with other molecular assemblies for the design of successful artificial leaves.

We simulated the intrinsic reaction path of the Green Fluorescent Protein (GFP) proton shuttle in both the ground state (S0) and first singlet excited state (S1), accounting for the main energetic and steric effects of the protein in a convenient model including the chromophore, the crystallographic water, and the residues directly involved in the proton transfer event. We adopted density functional theory (DFT) and time-dependent density functional theory (TD-DFT) levels to define the potential energy surfaces of the two electronic states, and we compared results obtained by the Damped Velocity Verlet and the Hessian-based Predictor-Corrector integrators of the intrinsic reaction coordinate, which gave a comparable and consistent picture of the mechanism. We show that, at S1, the GFP proton transfer becomes favored, with respect to S0, as suggested by the experimental evidence. As an important finding, this change is strictly related to the rearrangement of the hydrogen bond network composing the reaction path, which, in S1, relaxes to a tighter and planar configuration, as a consequence of the photoinduced relaxation in the GFP chromophore structure, thus prompting more effectively for the proton shuttle. Therefore, we give an unprecedented direct proof of the key role played by the photoinduced structural relaxation of the GFP on the chromophore photoacidity, validating, in particular, the hypothesis of Fang and co-workers [Nature 2009, 462, 200].

The formation of polaron pairs is one of the important photophysical processes that take place after the excitation in semiconducting organic polymers. First-principles Ehrenfest excited-state dynamics is a unique tool to investigate ultrafast photoinduced charge carrier dynamics and related nonequilibrium processes involving correlated electron-nuclear dynamics. In this work the formation of polaron pairs and their dynamical evolution in an oligomer of seven thiophene units is investigated with a combined approach of first-principles exciton-nuclear dynamics and wavelet analysis. The real-time formation of a polaron pair can be observed in the dipole evolution during the excited-state dynamics. The possible driving force of the polaron pair formation is investigated through qualitative correlation between the structural dynamics and the dipole evolution. The time-dependent characteristics and spectroscopic consequences of the polaron pair formation are probed using the wavelet analysis.

The growing production of biodiesel as a promising alternative and renewable fuel led as the main problem the dramatic increase of its by-product: glycerol. Different strategies for glycerol derivatization have been reported so far, some more efficient or sustainable than others. Herein, we report a very promising and eco-friendly transformation of glycerol in nontoxic solvents and chemicals (i.e., solketal, ketals), proposing three new families of Fe(III) compounds capable of catalysing glycerol acetalization with unpublished turn over frequencies (TOFs), and adhering most of the principles of green chemistry. The comparison between the activity of complexes of formula [FeCl3(1-R)] (1-R = substituted pyridinimine), [FeCl(2-R,R′)] (2-R,R′ = substituted O,O′-deprotonated salens) and their corresponding simple salts reveals that the former are extremely convenient because they are able to promote solketal formation with excellent TOFs, up to 105 h–1. Satisfactory performances were shown with respect to the entire range of substrates, with results being competitive to those reported in the literature so far. Moreover, the experimental activity was supported by an accurate and complete ab initio study, which disclosed the fundamental role of iron(III) as Lewis acid in promoting the catalytic activity. The unprecedented high activity and the low loading of the catalyst, combined with the great availability and the good eco-toxicological profile of iron, foster future applications of this catalytic process for the sustainable transformation of an abundant by-product in a variety of chemicals.

The detailed knowledge of excited state proton transfer mechanisms in complex environments is of paramount importance in chemistry. However, the definition of an effective reaction coordinate and the understanding of the driving force of the reaction can be difficult from both the experimental and the theoretical points of view. Here we analyzed by theoretical approaches the mechanism and the driving forces of the excited state proton transfer reaction occurring between the 7-hydroxy-4-(trifluoromethyl)coumarin photoacid and the 1-methylimidazole base molecules in toluene solution. In particular, we compared the intrinsic and the dynamical reaction pathways, obtained by integrating the reaction coordinate, and by performing ab initio simulations of molecular dynamics, respectively. Time-dependent density functional theory and polarizable solvation continuum models were adopted to define the excited state potential energy surface. Results were analyzed by means of the D(CT) electronic density based index. Our findings suggest that the reaction coordinate is mainly composed of several intra- and intermolecular modes of the reactants. An analysis of both the intrinsic coordinate and the dynamical results shows that the charge transfer induced by electronic excitation of the coumarin molecule is the main proton transfer driving force. With regards to the methodological validation, the combination of ab initio molecular dynamics with time-dependent density functional theory appears to be feasible and reliable to study excited state proton transfer reactions in the condensed phase.

A new time-resolved vibrational analysis unveils the mechanism of an excited state proton shuttle in green fluorescent protein., We simulated an excited state proton transfer in green fluorescent protein by excited state ab initio dynamics, and examined the reaction mechanism in both the time and the frequency domain through a multi resolution wavelet analysis. This original approach allowed us, for the first time, to directly compare the trends of photoactivated vibrations to femtosecond stimulated Raman spectroscopy results, and to give an unequivocal interpretation of the role played by low frequency modes in promoting the reaction. We could attribute the main driving force of the reaction to an important photoinduced softening of the ring–ring orientational motion of the chromophore, thus permitting the tightening of the hydrogen bond network and the opening of the reaction pathway. We also found that both the chromophore (in terms of its inter-ring dihedral angle and phenolic C–O and imidazolinone C–N bond distances) and its pocket (in terms of the inter-molecular oxygen’s dihedral angle of the chromophore pocket) relaxations are modulated by low frequency (about 120 cm–1) modes involving the oxygen atoms of the network. This is in agreement with the femtosecond Raman spectroscopy findings in the time-frequency domain. Moreover, the rate in proximity to the Franck Condon region involves a picosecond time scale, with a significant influence from fluctuations of nearby hydrogen bonded residues such as His148. This approach opens a new scenario with ab initio simulations as routinely used tools to understand photoreactivity and the results of advanced time resolved spectroscopy techniques.

The mechanism of base to base intermolecular proton shuttling occurring in the excited state proton transfer reaction between 7-hydroxy-4-(trifluoromethyl)coumarin (CouOH) and concentrated 1-methylimidazole base (1-MeId) in toluene solution is disclosed here by means of a computational approach based on Density Functional Theory (DFT) and Time Dependent DFT (TD-DFT). These methods allow us to characterize both the ground and excited state potential energy surfaces along the proton shuttling coordinate, and to assess the nature of the emitting species in the presence of an excess of 1-MeId. As a result, the tautomerism of CouOH is found to be photo-activated and, from a mechanistic point of view, the calculations clearly show that the overall driving force of the entire shuttling is the coumarin photoacidity, which is responsible for both the first proton transfer event and the strengthening of the following chain mechanism of base to base proton hopping.

In-chain donor/acceptor block copolymers comprised of alternating electron rich/poor moieties are emerging as promising semiconducting chromophores for use in organic photovoltaic devices. The mobilities of charge carriers in these materials are experimentally probed using gated organic field-effect transistors to quantify electron and hole mobilities, but a mechanistic understanding of the relevant charge diffusion pathways is lacking. To elucidate the mechanisms of electron and hole transport following excitation to optically accessible low-lying valence states, we utilize mean-field quantum electronic dynamics in the TDDFT formalism to explicitly track the evolution of these photo-accessible states. From the orbital pathway traversed in the dynamics, p- and n-type conductivities can be distinguished. The electronic dynamics of the studied polymers show the time-resolved transitions between the initial photoexcited state, a tightly-bound excitonic state that is dark to the ground state, and a partially charge separated state indicated by long-lived, out-of-phase charge oscillations along the polymer backbone. The frequency of these charge oscillations yields an insight into the characteristic mobilities of charge carriers in these materials. When the barycenters of the electron and hole densities are followed during the dynamics, a pseudo-classical picture for the translation of charge carrier densities along the polymer backbone emerges that clarifies a crucial aspect in the design of efficient organic photovoltaic materials.

The excited state proton transfer (ESPT) reaction taking place between 7-hydroxy-4-(trifluorometyl)coumarin and 1-methylimidazole, recently experimentally characterized, has been here considered as a case study to illustrate the possibility of using theoretical approaches rooted in density functional theory (DFT) and time-dependent DFT (TD-DFT) for the description of complex reactions occurring at the excited state. In particular, beside identifying all stable species occurring at the ground and excited state during the ESPT reaction, a quantitative characterization of their photophysical properties, such as absorption and emission, is obtained by properly including solvent effects. More interestingly, a computational protocol enabling one to locate possible reaction pathways for the ESPT is here proposed. This protocol is based on the use of density based indices purposely developed to characterize the properties of vertical and relaxed excited states which allow one to discriminate the most favorable reaction paths on potential energy surfaces that are in the case of ESPT intrinsically very flat and difficult to characterize based on sole energy criteria, thus opening a new scenario for the description of photoinduced proton transfer reactions.

A new descriptor tool enabling to qualitatively identify excited-state potential energy regions with high decay probability is here disclosed and applied to analyze a photoinduced ring-opening reaction already well characterized from the experimental and theoretical point of view. The analysis based on such descriptor allows one to highlight a high probability of excited-state deactivation within the Franck–Condon region, in agreement with experiments, and to qualitatively indentify the main mechanisms providing efficient pathways of photoreactivity.

The photophysical signature of the tautomeric species of the asymmetric (N,N-dimethylanilino)-1,3-diketone molecule are investigated using approaches rooted in density functional theory (DFT) and time-dependent DFT (TD-DFT). In particular, since this molecule, in the excited state, can undergo proton transfer reactions coupled to intramolecular charge transfer events, the different radiative and nonradiative channels are investigated by making use of different density-based indexes. The use of these tools, together with the analysis of both singlet and triplet potential energy surfaces, provide new insights into excited-state reactivity allowing one to rationalize the experimental findings including different behavior of the molecule as a function of solvent polarity.

In spite of the detailed information provided by advanced time-resolved spectroscopy, the understanding of the excited-state proton transfer (ESPT) reactivity remains difficult to obtain at molecular level. In this work we studied three photoacids showing different strength: the 8-hydroxy-1,3,6-pyrenetrisulfonate weak photoacid, the N-methyl-6-hydroxyquinolinium strong photoacid and the phenol-carboxyether dipicolinium cyanine (QCy9) superphotoacid, focusing on the intermolecular ESPT toward a solvent molecule or a base molecule in aqueous solution. To this aim, the ground and the first singlet excited-state potential energy surfaces of the three systems were characterized by means of the time-dependent density functional theory and a hybrid implicit/explicit model of the solvent. Main structural and photophysical features of the photoacids were assessed and satisfactorily compared with the experimental data. Energy profiles along the PT coordinate were analyzed in both the electronic states. We reproduced many important features of the photoacidity experimentally observed. The results suggest that the relative strength is mainly due to the different extent of charge transfer caused by the electronic transition in proximity of the acid group. Remarkably, we found that even in the case of the strongest photoacid (QCy9), showing a ESPT rate as rapid as to escape the solvent dynamics control, the PT is modulated and supported by the first solvation shell of the proton-accepting molecule. However, a complete understanding of this fascinating field needs the full account for the electronic and the molecular dynamics in play at different timescales.

In the present work, the general liquid optimized boundary (GLOB) model, which is a discrete/continuum approach for molecular simulations of liquids and solutions, has been extended to include a pressure coupling algorithm based on an extended phase-space scheme. To this end, a definition of the instantaneous pressure for a microscopic system has been derived from the minimum energy principle for the Helmholtz free energy. Applications to a pure liquid, such as water, and an aqueous solution of myoglobin are presented and the results are compared to those obtained using standard periodic boundary conditions (PBC).

In this work, the GLOB model, an effective and reliable computational approach well suited for ab initio and QM/MM molecular dynamics simulations of complex molecular systems in solution, has been applied to study two representative open-shell systems, the cobalt(II) ion and the glycine radical in aqueous solution, with special reference to their structural and magnetic properties. The main structural features of the solvent cage around the cobalt ion and the hydrogen bonding patterns around the neutral and zwitterionic forms of the glycine radical have been investigated in some detail. The general good agreement with experiments supports the use of the present model to investigate more challenging and biological/technological relevant open-shell systems.

Semi-empirical correction to density functional theory for dispersion (DFT-D) has been implemented for energies, analytical gradients, and Hessians in order to explore potential energy surfaces by means of a complete set of first-principle methods. The impact of non-bonding interactions on structures, binding energies and zero-point energy contributions as well as on ab initio molecular dynamics trajectories have been investigated for the well known case of benzene dimer. While the static results are in remarkable agreement with the most sophisticated post-Hartree–Fock approaches, the low cost of DFT-D allows to unravel dynamical aspects too, which are mandatory for situations ruled by weak interactions.