Mario Vazdar

The nonadiabatic deactivation paths of pyrrole

Multireference configuration interaction (MRCI) calculations have been performed for pyrrole with the aim of providing an explanation for the experimentally observed photochemical deactivation processes. Potential energy curves and minima on the crossing seam were determined using the analytic MRCI gradient and nonadiabatic coupling features of the COLUMBUS program system. A new deactivation mechanism based on an out-of-plane ring deformation is presented. This mechanism directly couples the charge transfer 1pipi* and ground states. It may be responsible for more than 50% of the observed photofragments of pipi*-excited pyrrole. The ring deformation mechanism should act complementary to the previously proposed NH-stretching mechanism, thus offering a more complete interpretation of the pyrrole photodynamics.

Automerization reaction of cyclobutadiene and its barrier height: an ab initio benchmark multireference average-quadratic coupled cluster study.

The problem of the double bond flipping interconversion of the two equivalent ground state structures of cyclobutadiene (CBD) is addressed at the multireference average-quadratic coupled cluster level of theory, which is capable of optimizing the structural parameters of the ground, transition, and excited states on an equal footing. The barrier height involving both the electronic and zero-point vibrational energy contributions is 6.3 kcal mol(-1), which is higher than the best earlier theoretical estimate of 4.0 kcal mol(-1). This result is confirmed by including into the reference space the orbitals of the CC sigma bonds beyond the standard pi orbital space. It places the present value into the middle of the range of the measured data (1.6-10 kcal mol(-1)). An adiabatic singlet-triplet energy gap of 7.4 kcal mol(-1) between the transition state (1)B(tg) and the first triplet (3)A(2g) state is obtained. A low barrier height for the CBD automerization and a small DeltaE((3)A(2g),(1)B(1g)) gap bear some relevance on the highly pronounced reactivity of CBD, which is briefly discussed.

Excited-state non-adiabatic dynamics simulations of pyrrole

Non-adiabatic on-the-fly-dynamics simulations of the photodynamics of pyrrole were performed at multireference configuration interaction level involving five electronic states with a simulation time of 200 fs. The analysis of the time dependence of the average state occupations shows that the deactivation of pyrrole to the electronic ground state takes place in about 140 fs. This deactivation time agrees very well with the experimentally measured time constant of 110 fs for the formation of fast hydrogen atoms. After excitation into the S4 state, 80% of the trajectories followed the NH-stretching mechanism giving rise to a population of fast H atoms. The computed average kinetic energy is in good accord with the experimentally observed average kinetic energy of the fast hydrogen atoms. It is found that 10% of trajectories followed the ring-puckering mechanism and 3% followed the ring-opening mechanism. This latter mechanism was characterized in pyrrole for the first time and involves the conical intersection of lowest energy of this molecule.

Like-Charge Guanidinium Pairing from Molecular Dynamics and Ab Initio Calculations

Pairing of guanidinium moieties in water is explored by molecular dynamics simulations of short arginine-rich peptides and ab initio calculations of a pair of guanidinium ions in water clusters of increasing size. Molecular dynamics simulations show that, in an aqueous environment, the diarginine guanidinium like-charged ion pairing is sterically hindered, whereas in the Arg-Ala-Arg tripeptide, this pairing is significant. This result is supported by the survey of protein structure databases, where it is found that stacked arginine pairs in dipeptide fragments exist solely as being imposed by the protein structure. In contrast, when two arginines are separated by a single amino acid, their guanidinium groups can freely approach each other and they frequently form stacked pairs. Molecular dynamics simulations results are also supported by ab initio calculations, which show stabilization of stacked guanidinium pairs in sufficiently large water clusters.

Absorption and fluorescence of PRODAN in phospholipid bilayers: a combined quantum mechanics and classical molecular dynamics study.

Absorption and fluorescence spectra of PRODAN (6-propionyl-2-dimethylaminonaphthalene) were studied by means of the time-dependent density functional theory and the algebraic diagrammatic construction method. The influence of environment, a phosphatidylcholine lipid bilayer and water, was taken into account employing a combination of quantum chemical calculations with empirical force-field molecular dynamics simulations. Additionally, experimental absorption and emission spectra of PRODAN were measured in cyclohexane, water, and lipid vesicles. Both planar and twisted configurations of the first excited state of PRODAN were taken into account. The twisted structure is stabilized in both water and a lipid bilayer, and should be considered as an emitting state in polar environments. Orientation of the excited dye in the lipid bilayer significantly depends on configuration. In the bilayer, the fluorescence spectrum can be regarded as a combination of emission from both planar and twisted structures.

Orientational Dependence of the Affinity of Guanidinium Ions to the Water Surface

The behavior of guanidinium chloride at the surface of aqueous solutions is investigated using classical molecular dynamics (MD) simulations. It is found that the population of guanidinium ions oriented parallel to the interface is greater in the surface region than in bulk. The opposite is true for ions in other orientations. Overall, guanidinium chloride is depleted in the surface region, in agreement with the fact that the addition of guanidinium chloride increases the surface tension of water. The orientational dependence of the surface affinity of the guanidinium cation is related to its anisotropic hydration. To bring the ion to the surface in the parallel orientation does not require hydrogen bonds to be broken, in contrast to other orientations. The surface enrichment of parallel-oriented guanidinium indicates that its solvation is more favorable near the surface than in bulk solution for this orientation. The dependence of the bulk and surface properties of guanidinium on the force field parameters is also investigated. Despite significant quantitative differences between the force fields, the surface behavior is qualitatively robust. The implications for the orientations of the guanidinium groups of arginine side chains on protein surfaces are also outlined.

Electrodeposition of nano-sized nuclei of magnetic Co-Ni alloys onto n-Si (100)

Nanometer-sized nuclei of Co–Ni alloys were electrodeposited onto n-Si (1 0 0) electrodes from buffered solutions of the metallic ions. The energy levels of the interface Si/(Co(II),Ni(II))aq were determined from Mott–Schottky plots. Electrodeposition of the alloys occurs by an anomalous mechanism and preferential deposition of Co took place at all the Co(II)/Ni(II) ratios investigated. The morphology and magnetic properties of the nuclei were studied by AFM and MFM, respectively. It is showed that nuclei of height <30 nm and low aspect ratio present a single-domain magnetic behaviour. Larger nuclei depict complex magnetic characteristics, with a progressive built-up of multi-domains.

Cholesterol under oxidative stress-How lipid membranes sense oxidation as cholesterol is being replaced by oxysterols.

The behavior of oxysterols in phospholipid membranes and their effects on membrane properties were investigated by means of dynamic light scattering, fluorescence spectroscopy, NMR, and extensive atomistic simulations. Two families of oxysterols were scrutinized-tail-oxidized sterols, which are mostly produced by enzymatic processes, and ring-oxidized sterols, formed mostly via reactions with free radicals. The former family of sterols was found to behave similar to cholesterol in terms of molecular orientation, roughly parallel to the bilayer normal, leading to increasing membrane stiffness and suppression of its membrane permeability. In contrast, ring-oxidized sterols behave quantitatively differently from cholesterol. They acquire tilted orientations and therefore disrupt the bilayer structure with potential implications for signaling and other biochemical processes in the membranes.

Aqueous Guanidinium–Carbonate Interactions by Molecular Dynamics and Neutron Scattering: Relevance to Ion–Protein Interactions

Guanidinium carbonate was used in this study as a simple proxy for the biologically relevant arginine-carbonate interactions in water. Molecular dynamics (MD) simulations of guanidinium carbonate were performed with nonpolarizible water using two implementations of the ion force fields. In the first, the ions had full charges, while in the second, the ions had reduced charges in order to effectively account for electronic polarization effects of water. The results from the simulations were then compared to data from previous neutron scattering experiments. It was found that there were significant discrepancies between the full charge force field MD simulations and the experimental results due to excessive ion pairing and clustering in the former. In contrast, reducing the ionic charges yields a more regular solution with a simulated structure, which fits well the experimental data.

Ab initio quantum dynamical study of the multi-state nonadiabatic photodissociation of pyrrole

There has been a substantial amount of theoretical investigations on the photodynamics of pyrrole, often relying on surface hopping techniques or, if fully quantal, confining the study to the lowest two or three singlet states. In this study we extend ab initio based quantum dynamical investigations to cover simultaneously the lowest five singlet states, two π-σ∗ and two π-π∗ excited states. The underlying potential energy surfaces are obtained from large-scale MRCI ab initio computations. These are used to extract linear and quadratic vibronic coupling constants employing the corresponding coupling models. For the N-H stretching mode Q(24) an anharmonic treatment is necessary and also adopted. The results reveal a sub-picosecond internal conversion from the S(4) (π-π∗) state, corresponding to the strongly dipole-allowed transition, to the S(1) and S(2) (π-σ∗) states and, hence, to the ground state of pyrrole. The significance of the various vibrational modes and coupling terms is assessed. Results are also presented for the dissociation probabilities on the three lowest electronic states.