Resumen: Hemos llevado a cabo una intervención con futuros maestros de Educación Primaria con la intención de trabajar el concepto de sostenibilidad.La secuencia desarrollada incluye la elaboración de un proyecto y una unidad didáctica por investigación escolar.A través de un portafolio los alumnos expresaron sus reflexiones sobre su idea de sostenibilidad.Para analizarlas se diseñó una tabla de niveles con indicadores.Más

A comparison between the free-energy surfaces of the all- trans-retinal protonated Schiff base (RPSB) and its 10-methylated derivative in gas phase and methanol solution is performed at CASSCF//CASSCF and CASPT2//CASSCF levels. Solvent effects were included using the average solvent electrostatic potential from molecular dynamics method. This is a QM/MM (quantum mechanics/molecular mechanics) method that makes use of the mean field approximation. It is found that the methyl group bonded to C10 produces noticeable changes in the solution free-energy profile of the S1 excited state, mainly in the relative stability of the minimum energy conical intersections (MECIs) with respect to the Franck-Condon (FC) point. The conical intersections yielding the 9- cis and 11- cis isomers are stabilized while that yielding the 13- cis isomer is destabilized; in fact, it becomes inaccessible by excitation to S1. Furthermore, the planar S1 minimum is not present in the methylated compound. The solvent notably stabilizes the S2 excited state at the FC geometry. Therefore, if the S2 state has an effect on the photoisomerization dynamics, it must be because it permits the RPSB population to branch around the FC point. All these changes combine to speed up the photoisomerization in the 10-methylated compound with respect to the native compound.

The absorption spectra of several systems containing an adenine (A) cation are computed in water by TD-DFT calculations, using three different functionals and including solvent effects by a mixed discrete–continuum approach. Our calculations well reproduce the experimental absorption spectrum of deoxyadenosine nucleoside cation (dA+), which provides an intense peak at 350 nm and a weaker one at 600 nm, allowing its assignment. M052X and CAM-B3LYP predict that the hole is essentially localized on a single A base also in the cationic forms of dApdA dinucleotide, both in single and in double strand, and of (dApdT) duplex, exhibiting absorption spectra similar to that of dA+. For all these compounds, B3LYP predicts instead the delocalization of the hole over two A bases. On this ground, it is possible to propose which would be the spectral signature of partial hole delocalization. Vertical and adiabatic ionization potentials for the compounds under investigation are also computed, providing values in good agreement with the available experimental results.

Substituent and solvent effects on the excited state dynamics of the Photoactive Yellow Protein chromophore are studied using the average solvent electrostatic potential from molecular dynamics (ASEP/MD) method. Four molecular models were considered: the ester and thioester derivatives of the p-coumaric acid anion and their methylated derivatives. We found that the solvent produces dramatic modifications on the free energy profile of the S1 state: 1) Two twisted structures that are minima in the gas phase could not be located in aqueous solution. 2) Conical intersections (CIs) associated with the rotation of the single bond adjacent to the phenyl group are found for the four derivatives in water solution but only for thio derivatives in the gas phase. 3) The relative stability of minima and CIs is reverted with respect to the gas phase values, affecting the prevalent de-excitation paths. As a consequence of these changes, three competitive de-excitation channels are open in aqueous solution: the fluorescence emission from a planar minimum on S1, the trans-cis photoisomerization through a CI that involves the rotation of the vinyl double bond, and the nonradiative, nonreactive, de-excitation through the CI associated with the rotation of the single bond adjacent to the phenyl group. In the gas phase, the minima are the structures with the lower energy, while in solution these are the conical intersections. In solution, the de-excitation prevalent path seems to be the photoisomerization for oxo compounds, while thio compounds return to the initial trans ground state without emission.

We present a quantum-chemical investigation of the excited states of the complex [Re(CO)3(Im)(Phen)]+ (Im = imidazole; Phen = 1,10-phenanthroline) in solution including spin-orbit couplings and vibrational sampling. To this aim, we implemented electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) in the Amsterdam Density Functional program suite, suitable for time-dependent density functional calculations including spin-orbit couplings. The new implementation is employed to simulate the absorption spectrum of the complex, which is compared to the results of implicit continuum solvation and frozen-density embedding. Molecular dynamics simulations are used to sample the ground state conformations in solution. The results demonstrate that any study of the excited states of [Re(CO)3(Im)(Phen)]+ in solution and their dynamics should include extensive sampling of vibrational motion and spin-orbit couplings.

In an attempt to shed light on the environmental effects on the deactivation channels of the PYP chromophore, radiative and non-radiative deactivation mechanisms of the anionic p-coumaric acid methyl ester (pCE-) in the gas phase and water solution are compared at the CASPT2//CASSCF/cc-pVDZ level and, when necessary, at the CASPT2//CASPT2/cc-pVDZ level. We find that the solvent produces dramatic modifications on the free energy profile of the S1 state. Two twisted structures that are minima in the gas phase could not be localized in aqueous solution. Furthermore, the relative stability of minima and conical intersections (CIs) is reverted with respect to the gas phase values, affecting the prevalent de-excitation paths. As a consequence of these changes, three competitive de-excitation channels are open in aqueous solution: the fluorescence emission from a planar minimum on S1, the trans-cis photoisomerization through a CI that involves the rotation of the vinyl double bond and the non-radiative, non-reactive, de-excitation through the CI associated with the rotation of the single bond adjacent to the phenyl group. In the gas phase, the minima are the structures with lower energy, while in solution the CIβ structure, characterized by a large charge separation, is strongly stabilized by interactions with water molecules and becomes the structure with the lowest energy on S1. These facts explain the low fluorescence signal of pCE- in aqueous solution and the presence of partial trans-cis photoisomerization in this system.

Adenyl radicals generated in DNA single and double strands, (dA)20 and (dA)20·(dT)20, by one- and two-photon ionization by 266 nm laser pulses decay at 600 nm with half-times of 1.0 ± 0.1 and 4 ± 1 ms, respectively. Though ionization initially forms the cation radical, the radicals detected for (dA)20 are quantitatively identified as N6-deprotonated adenyl radicals by their absorption spectrum, which is computed quantum mechanically employing TD-DFT. Theoretical calculations show that deprotonation of the cation radical induces only weak spectral changes, in line with the spectra of the adenyl radical cation and the deprotonated radical trapped in low temperature glasses.

The reaction intermediate in the path leading to UV-induced formation of adenine dimers A═A and AA* is identified for the first time quantum mechanically, using PCM/TD-DFT calculations on (dA)2 (dA: 2'deoxyadenosine). In parallel, its fingerprint is detected in the absorption spectra recorded on the millisecond time-scale for the single strand (dA)20 (dA: 2'deoxyadenosine).

It is well known that solvents can modify the frequency and intensity of the solute spectral bands, the thermodynamics and kinetics of chemical reactions, the strength of molecular interactions or the fate of solute excited states. The theoretical study of solvent effects is quite complicated since the presence of the solvent introduces additional difficulties with respect to the study of analogous problems in gas phase. The mean field approximation (MFA) is used for many of the most employed solvent effect theories as it permits to reduce the computational cost associated to the study of processes in solution. In this chapter we revise the performance of ASEP/MD, a quantum mechanics/molecular mechanics method developed in our laboratory that makes use of this approximation. It permits to combine state of the art calculations of the solute electron distribution with a detailed, microscopic, description of the solvent. As examples of application of the method we study solvent effects on the absorption spectra of some molecules involved in photoisomerization processes of biological systems.

It is well known that solvents can modify the frequency and intensity of the solute spectral bands, the thermodynamics and kinetics of chemical reactions, the strength of molecular interactions or the fate of solute excited states. The theoretical study of solvent effectsSolvent effects is quite complicated since the presence of the solvent introduces additional difficulties with respect to the study of analogous problems in gas phase. The mean fieldMean field theory approximation (MFA) is used for many of the most employed solvent effect theories as it permits to reduce the computational cost associated to the study of processes in solution. In this chapter we revise the performance of ASEP/MD, a quantum mechanics/molecular mechanics method developed in our laboratory that makes use of this approximation. It permits to combine state of the art calculations of the solute electron distribution with a detailed, microscopic, description of the solvent. As examples of application of the method we study solvent effects on the absorption spectra of some molecules involved in photoisomerization processes of biological systems.

A State-Specific PCM–DFT method taking into account dynamical solvent effects in the calculation of ionization energies (IEs) is reported. Mono-methylated nucleobases, representative of DNA blocks, are used as a test case. The accuracy of several functionals is tested, the M05-2X functional providing the most consistent performance. The inclusion of some molecules of the first hydration shell does not significantly affect the computed IE. The obtained vertical and adiabatic IEs are in good agreement with the available experimental and computational results, supporting the reliability of our method, enabling the calculation of first and upper ionization energies by using a single approach.

Solvent effects on the absorption and emission spectra of the p-coumaric methyl ester, a model for the chromophore of the photoactive yellow protein, are analyzed using a QM/MM method known as Average Solvent Electrostatic Potential from Molecular Dynamics (ASEP/MD), which combines state of the art calculations in the description of the solute with a detailed description of the solvent through molecular dynamics simulations. Geometries for the ground and the bright excited state are optimized by using multiconfigurational (CASSCF and CASPT2) and DFT methods. Three different functionals (B3LYP, CAM-B3LYP and PBE0) were used with DFT. The effect of specific interactions and of the interplay between static and dynamic electron correlation on the spectra are analyzed. In absorption there is a flux of charge from the phenolic toward the carboxylic moiety, and the opposite is found in emission. These fluxes are similar in gas phase; however, in water solution, the flux increases with respect to the gas phase value during the absorption but not in the emission. As a consequence of these fluxes the solvent structure around the chromophore changes after the excitation. The average number of water molecules decreases around the phenolic oxygen and increases around the carbonyl oxygen. These changes in the solvent structure together the reduction of the charge on the phenolic oxygen atom during the excitation permit to explain the large Stokes shift, 0.84 eV, displayed by this system. (C) 2014 Elsevier B.V. All rights reserved.

The ground and low-lying excited free energy surfaces of 4-amino-4'-cyano azobenzene, a molecule that has been proposed as building block for chiroptical switches, are studied in gas phase and a variety of solvents (benzene, chloroform, acetone, and water). Solvent effects on the absorption and emission spectra and on the cis-trans thermal and photo isomerizations are analyzed using two levels of calculation: TD-DFT and CASPT2/CASSCF. The solvent effects are introduced using a polarizable continuum model and a QM/MM method, which permits one to highlight the role played by specific interactions. We found that, in gas phase and in agreement with the results found for other azobenzenes, the thermal cis-trans isomerization follows a rotation-assisted inversion mechanism where the inversion angle must reach values close to 180° but where the rotation angle can take almost any value. On the contrary, in polar solvents the mechanism is controlled by the rotation of the CN═NC angle. The change in the mechanism is mainly related to a better solvation of the nitrogen atoms of the azo group in the rotational transition state. The photoisomerization follows a rotational pathway both in gas phase and in polar and nonpolar solvents. The solvent introduces only small modifications in the nπ* free energy surface (S1), but it has a larger effect on the ππ* surface (S2) that, in polar solvents, gets closer to S1. In fact, the S2 band of the absorption spectrum is red-shifted 0.27 eV for the trans isomer and 0.17 eV for the cis. In the emission spectrum the trend is similar: only S2 is appreciably affected by the solvent, but in this case a blue shift is found.

The effects of the solvent and protonation state on the electronic absorption spectrum of the para-coumaric acid (pCA), a model of the photoactive yellow protein (PYP), have been studied using the ASEP/MD (averaged solvent electrostatic potential from molecular dynamics) method. Even though, in the protein, the chromophore is assumed to be in its phenolate monoanionic form, when it is found in water solution pH control can favor neutral, monoanionic, and dianionic species. As the pCA has two hydrogens susceptible of deprotonation, both carboxylate and phenolate monoanions are possible. Their relative stabilities are strongly dependent on the medium. In gas phase, the most stable isomer is the phenolate while in aqueous solution it is the carboxylate, although the population of the phenolate form is not negligible. The s-cis, s-trans, syn, and anti conformers have also been included in the study. Electronic excited states of the chromophore have been characterized by SA-CAS(14,12)-PT2/cc-pVDZ level of theory. The bright state corresponds, in all the cases, to a π → π* transition involving a charge displacement in the system. The magnitude and direction of this displacement depends on the protonation state and on the environment (gas phase or solution). In the same way, the calculated solvatochromic shift of the absorption maximum depends on the studied form, being a red shift for the neutral, carboxylate monoanion, and dianionic chromophores and a blue shift for the phenolate monoanion. Finally, the contribution that the solvent electronic polarizability has on the solvent shift was analyzed. It represents a very important part of the total solvent shift in the neutral form, but its contribution is completly negligible in the mono- and dianionic forms.

The ASEP/MD (averaged solvent electrostatic potential from molecular dynamics) method was employed in studying the environment effects (solvent and counterion) on the absorption spectrum of a model of the 11-cis-retinal protonated Schiff base. Experimental studies of the absorption spectra of the rhodopsin chromophore show anomalously large solvent shifts in apolar solvents. In order to clarify their origin, we study the role of the counterion and of the solute-solvent interactions. We compare the absorption spectra in the gas phase, cyclohexane, dichloromethane, and methanol. The counterion effect was described from both a classical and quantum point of view. In the latter case, the contribution of the chromophore-counterion charge transfer to the solvent shift could be analyzed. To the best of our knowledge, this is the first time that counterion and solvent effects on the absorption properties of the 11-cis-retinal chromophore have been simultaneously examined. We conclude that the counterion-solute ionic pair in the gas phase is not a good model to represent the solvent shift in nonpolar solvents, as it does not account for the effect that the thermal agitation of the solvent has on the geometry of the ionic pair. In contrast to nonpolar solvents, the experimental solvent shift values in methanol can be exclusively explained by the polarity of the medium. In dichloromethane, the presence of the counterion does not modify the solvent shift of the first absorption band, but it affects the position of the second excited state. In the three solvents considered, the first two excited states become almost degenerate.

The conformational equilibrium of 1,4-dioxane in the gas phase, in the pure liquid, and in aqueous solution has been studied by means of the Average Solvent Electrostatic Potential from Molecular Dynamics (ASEP/MD) method and the Integral Equation Formalism for the Polarizable Continuum Model (IEF-PCM). The dioxane molecule was described at the DFT(B3LYP)/aug-cc-pVTZ level. In the three phases, the equilibrium is almost completely shifted toward the chair conformer, with populations of the twist-boat conformers lower than 0.01 %. The equilibrium is dominated by the internal energy of the molecule, as the solute–solvent interaction free energies are very similar in the three conformers considered (chair, 1,4 twist-boat, and 2,5 twist-boat). In the pure liquid, where the dioxane–dioxane interaction is dominated by the Lennard-Jones term, the structure is characteristic of a van der Waals liquid. However, the decrease in the C–H distance from gas phase to solution, the increase in the C–H vibrational frequencies, and the presence of a shoulder in the O–Haxial pair radial distribution function point to the presence of a weak C–H–O hydrogen bond. The analysis of the occupancy maps of water oxygen and hydrogen atoms around the 1,4-dioxane molecule confirms this conclusion. Contrary to what is found in small water-dioxane clusters, in the liquid, there is a preference for oxygen atoms to interact with axial hydrogen atoms to form C–H–O hydrogen bonds. Comparison of ASEP/MD and IEF-PCM results indicates that including specific interactions is very important for an adequate description of the solute–solvent interaction; however, the influence of these interactions does not translate in changes in the relative stability of the conformers because it cancels out when energy differences are calculated.

We present a new method for simulating fluorescence observables, particularly those related to bulk and single-molecule fluorescence-detected resonance energy transfer (FRET) experiments. In this method, a molecular dynamics (MD) simulation is used to sample configuration space and quantum mechanics calculations are used to estimate the electronic coupling between the donor and acceptor probes for snapshots along the MD trajectory. A Markov chain method is used to sample the resulting electronic coupling trajectory allowing accurate simulation of any desired fluorescence observables, such as single-molecule FRET efficiency histograms or time-resolved donor fluorescence decays. The Markov chain results will be compared with the results of simple histogram and averaging schemes showing that the Markov chain is the only one that yields realistic results in well known examples such as the rapid diffusion limit. This combination of computational methods also avoids some pitfalls of traditional FRET analysis such as the kappa-squared and the ideal dipole approximations. Because the simulation results can be compared directly with experimental observables, this method may allow more detail to be derived from experiment than is traditionally possible.

Ab initio molecular dynamics simulations of the acetone–CO2 complex (MP2/6-31G(d) level) were performed to investigate the effect of dynamics at finite temperature on the weak electron donor–acceptor intermolecular interactions. In addition, we carried out a study of the free energy of formation of the complex by means of umbrella sampling technique at the MP2 level with a perturbative CCSD(T) correction. The potential of mean force was obtained along a reaction coordinate describing the acetone–CO2 interaction. The results obtained here support some hypothesis that we already explored in past works using static electronic calculations. In particular, when interacting with a molecule having a carbonyl function, carbon dioxide displays both Lewis acid and Lewis base behaviour. This property can be exploited to design molecular systems that are easily solubilised in supercritical CO2.