PubMed Health⌕ Search

Biomedical subjects

Michiel Sprik

Publications and source records attributed to Michiel Sprik.

14 recordsLinked to original sources
PubMed2006

Long-range solvent effects on the orbital interaction mechanism of water acidity enhancement in metal ion solutions: a comparative study of the electronic structure of aqueous Mg and Zn dications.

We study the dissociation of water coordinated to a divalent metal ion center, M2+ = Mg2+, Zn2+ using density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. First, the proton affinity of a coordinated OH- group is computed from gas-phase Mg2+(H2O)5(OH-), which yields a relative higher gas-phase acidity for a Zn2+-coordinated as compared to a Mg2+-coordinated water molecule, DeltapKa(gp) = 5.3. We explain this difference on the basis of a gain in stabilization energy of the Zn2+(H2O)5(OH-) system arising from direct orbital interaction between the coordinated OH- and the empty 4s state of the cation. Next, we compute the acidity of coordinated water molecules in solution using free-energy thermodynamic integration with constrained AIMD. This approach yields pKa Mg2+ = 11.2 and pKa Zn2+ = 8.4, which compare favorably to experimental data. Finally, we examine the factors responsible for the apparent decrease in the relative Zn2+-coordinated water acidity in going from the gas-phase (DeltapKa(gp) = 5.3) to the solvated (DeltapKa = 2.8) regime. We propose two simultaneously occurring solvation-induced processes affecting the relative stability of Zn2+(H2O)5(OH-), namely: (a) reduction of the Zn 4s character in solution states near the bottom of the conduction band; (b) hybridization between OH- orbitals and valence-band states of the solvent. Both effects contribute to hindering the OH- --> Zn2+ charge transfer, either by making it energetically unfavorable or by delocalizing the ligand charge density over several water molecules.

Acids↗
PubMed2006

Density functional theory study of tetrathiafulvalene and thianthrene in acetonitrile: structure, dynamics, and redox properties.

The redox potentials of the organic compounds tetrathiafulvalene (TTF) and thianthrene (TH) in an explicit aprotic polar solvent, acetonitrile, have been computed using ab initio molecular dynamics simulation based on a Gaussian basis set methodology. The density functional description of the pure solvent yields a diffuse and mobile liquid, with structural and dynamical properties that are in good agreement with earlier classical models and experiment. Molecular dynamics simulation of both solute species in their neutral and radical cation states combined with free energy difference calculations result in estimates for the redox potentials of the reactions TH*+ + TTF --> TH + TTF*+ and TH2+ + TTF*+ --> TH*+ + TTF2+. The obtained values are 0.95 +/- 0.06 and 1.09 +/- 0.06 V, respectively, in excellent agreement with experimental data of 0.93 and 1.08 V. Our computational approach is based on Marcus theory, assuming quadratic free energy surfaces. We show that this approximation can still be accurate in systems, such as TH, that undergo a significant change in geometry upon oxidation. Furthermore, despite the different localization of the spin density in the radical cations, results based on self-interaction-corrected functionals and on standard generalized gradient approximations are identical to within 10 meV.

Journal Article↗
PubMed2006

Diabatic free energy curves and coordination fluctuations for the aqueous Ag+Ag2+ redox couple: a biased Born-Oppenheimer molecular dynamics investigation.

Biased Born-Oppenheimer molecular dynamics simulations are performed to compute redox potential and free energy curves for the redox half reaction Ag(+)-->Ag(2+)+e(-) in aqueous solution. The potential energy surfaces of reactant and product state are linearly coupled and the system transferred from the reduced state to the oxidized state by variation of the coupling parameter from 0 to 1. The redox potential is obtained by thermodynamic integration of the average ionization energy of Ag(+). Diabatic free energy curves of reduced (R) and oxidized (O) states are obtained to good statistical accuracy by reweighting and combining the set of biased distributions of the ionization energy. The diabatic free energy curves of Ag(+) and Ag(2+) are parabolic over a wide range of the reaction coordinate in agreement with the linear response assumption that underlies Marcus theory. However, we observe deviations from parabolic behavior in the equilibrium region of Ag(+) and find different values for the reorganization free energy of R (1.4 eV) and O (0.9 eV). The computed reorganization free energy of Ag(2+) is in good agreement with the experimental estimate of 0.9-1.2 eV obtained from photoelectron spectroscopy. As suggested by our calculations, the moderate deviation from linear response behavior found for Ag(+) is likely related to the highly fluxional solvation shell of this ion, which exhibits water exchange reactions on the picosecond time scale of the present molecular dynamics simulation.

Journal Article↗
PubMed2005

On the position of the highest occupied molecular orbital in aqueous solutions of simple ions.

The energies of the highest occupied molecular orbital (HOMO) of four simple microsolvated aqua ion clusters (Na+, Ag+, Cl-, CN-) are computed for varying numbers of water molecules. Extrapolating to infinite hydration numbers we find that these energies approach a value of -6 eV. This limiting one-electron energy is within a margin of +/-1 eV independent of the character of the ion and is 4 eV lower compared to the estimate obtained for the HOMO energy of the ions in aqueous solution under periodic boundary conditions. We argue that this discrepancy must the attributed to a shift in the reference of the one-electron potential of the periodic solvent model.

Journal Article↗
PubMed2005

Time-dependent density functional theory description of on-site electron repulsion and ligand field effects in the optical spectrum of hexaaquoruthenium(II) in solution.

Time-dependent density functional theory (TDDFT) and density functional-based molecular dynamics were used to simulate the finite temperature t(2g)5 e(g) <-- t(2g)6 absorption band of the Ru2+ hexahydrate coordination complex in aqueous solution. The (1)T1 <-- (1)A1 and (1)T2 <-- (1)A1 molecular term splitting of this transition, which is not accounted for by the Kohn-Sham excitation spectrum, is shown to be satisfactorily reproduced by TDDFT at the BLYP/ALDA level of theory. Comparison to the spectrum of the Ru2+ (H2O)6 complex in vacuo computed by similar density functional classical molecular dynamics methods leads to the observation that bulk solvation has a negligible effect on the position and the shape of the absorption profile.

Journal Article↗
PubMed2005

Density-functional molecular-dynamics study of the redox reactions of two anionic, aqueous transition-metal complexes.

The thermochemistry of the RuO(4)(2-)+MnO(4)(-)-->RuO(4)(-)+MnO(4)(2-) redox reaction in aqueous solution is studied by separate density-functional-based ab initio molecular-dynamics simulations of the component half reactions RuO(4)(2-)-->RuO(4)(-)+e(-) and MnO(4)(2-)-->MnO(4)(-)+e(-). We compare the results of a recently developed grand-canonical method for the computation of oxidation free energies to the predictions by the energy-gap relations of the Marcus theory that can be assumed to apply to these reactions. The calculated redox potentials are in good agreement. The subtraction of the half-reaction free energies gives an estimate of the free energy of the full reaction. The result obtained from the grand-canonical method is -0.4 eV, while the application of the Marcus theory gives -0.3 eV. These should be compared to the experimental value of 0.0 eV. Size effects, in response to increasing the number of water molecules in the periodic model system from 30 to 48, are found to be small ( approximately 0.1 eV). The link to the Marcus theory also has enabled us to compute reorganization free energies for oxidation. For both the MnO(4)(2-) and RuO(4)(2-) redox reactions we find the same reorganization free energy of 0.8 eV (1.0 eV in the larger system). The results for the free energies and further analysis of solvation and electronic structure confirm that these two tetrahedral oxoanions show very similar behavior in solution in spite of the central transition-metal atoms occupying a different row and column in the periodic table.

Journal Article↗
PubMed2005

Ab initio molecular dynamics simulation of the aqueous Ru2+/Ru3+ redox reaction: the Marcus perspective.

A well-behaved (low spin) transition metal aqua ion, Ru(aq)(2+), is used as a model system in an ab initio molecular dynamics study of a redox half reaction to which the Marcus theory of electron transfer is assumed to apply. Using constraint methods, we show that aqueous Ru(2+) can be reversibly transformed to Ru(3+) under the control of the classical solvent electrostatic potential as order parameter. The mean force is found to vary linearly with the order parameter in accordance with the Marcus theory. As can be expected for a half reaction, the slope in the oxidized and reduced states are asymmetric differing by approximately a factor of two. As a further test, we verify that the corresponding quadratic potential of mean force is in excellent agreement with the free energy profile obtained from the Gaussian distribution of potential fluctuations sampled from free (unconstrained) runs of the reduced and oxidized system. Similar to experimental electrochemical methods, our simulation scheme enables us to manipulate the effective thermodynamic driving force and align the free energy minima of product and reactant state. The activation energy and reaction entropy computed under these conditions are discussed and analyzed from the Marcus perspective.

Journal Article↗
PubMed2005

The influence of temperature and density functional models in ab initio molecular dynamics simulation of liquid water.

The performance of density functional theory methods for the modeling of condensed aqueous systems is hard to predict and validation by ab initio molecular simulation of liquid water is absolutely necessary. In order to assess the reliability of these tests, the effect of temperature on the structure and dynamics of liquid water has been characterized with 16 simulations of 20 ps in the temperature range of 280-380 K. We find a pronounced influence of temperature on the pair correlation functions and on the diffusion constant including nonergodic behavior on the time scale of the simulation in the lower temperature range (which includes ambient temperature). These observations were taken into account in a consistent comparison of a series of density functionals (BLYP, PBE, TPSS, OLYP, HCTH120, HCTH407). All simulations were carried out using an ab initio molecular dynamics approach in which wave functions are represented using Gaussians and the density is expanded in an auxiliary basis of plane waves. Whereas the first three functionals show similar behavior, it is found that the latter three functionals yield more diffusive dynamics and less structure.

Journal Article↗
PubMed2004

Density functional calculation of the electronic absorption spectrum of Cu+ and Ag+ aqua ions.

The UV absorption of aqueous Cu+ and Ag+ has been studied using Time Dependent Density Functional Theory (TDDFT) response techniques. The TDDFT electronic spectrum was computed from finite temperature dynamical trajectories in solution generated using the Density Functional Theory (DFT) based Ab Initio Molecular Dynamics (AIMD) method. The absorption of the two ions is shown to arise from similar excitation mechanisms, namely transitions from d orbitals localized on the metal center to a rather delocalized state originating from hybridization of the metal s orbital to the conduction band edge of the solvent. The ions differ in the way the spectral profile builds up as a consequence of solvent thermal motion. The Cu+ absorption is widely modulated, both in transition energies and intensities by fluctuations in the coordination environment which is characterized by the formation of strong coordination bonds to two water molecules in an approximately linear geometry. Though, on average, absorption intensities are typical of symmetry forbidden transitions of metal ions in the solid state, occasionally very short (<100 fs) bursts in intensity are observed, associated with anomalous Cu-H interactions. Absorption by the Ag+ complex is in comparison relatively stable in time, and can be interpreted in terms of the energy splitting of the metal 4d manifold in an average crystal field corresponding to a fourfold coordination in a distorted tetrahedral arrangement. Whereas the spectral profile of the Ag+ aqua ion is in good agreement with experiment, the overall position of the band is underestimated by 2 eV in the BLYP approximation to DFT. The discrepancy with experiment is reduced to 1.3 eV when a hybrid functional (PBE0) is used. The remaining inaccuracy of TDDFT in this situation is related to the delocalized character of the target state in d-->s transitions.

Journal Article↗
PubMed2004

Electronic structure and solvation of copper and silver ions: a theoretical picture of a model aqueous redox reaction.

Electronic states and solvation of Cu and Ag aqua ions are investigated by comparing the Cu(2+) + e(-)--> Cu(+) and Ag(2+) + e(-)--> Ag(+) redox reactions using density functional-based computational methods. The coordination number of aqueous Cu(2+) is found to fluctuate between 5 and 6 and reduces to 2 for Cu(+), which forms a tightly bound linear dihydrate. Reduction of Ag(2+) changes the coordination number from 5 to 4. The energetics of the oxidation reactions is analyzed by comparing vertical ionization potentials, relaxation energies, and vertical electron affinities. The model is validated by a computation of the free energy of the full redox reaction Ag(2+) + Cu(+) --> Ag(+) + Cu(2+). Investigation of the one-electron states shows that the redox active frontier orbitals are confined to the energy gap between occupied and empty states of the pure solvent and localized on the metal ion hydration complex. The effect of solvent fluctuations on the electronic states is highlighted in a computation of the UV absorption spectrum of Cu(+) and Ag(+).

Journal Article↗
PubMed2002

Ab initio molecular dynamics for molecules with variable numbers of electrons.

The ab initio molecular dynamics method is extended to treat exchange of electrons between molecules and a reservoir at fixed chemical potential. The method is based on a rigorously grand-canonical density functional approach using separate potential energy surfaces for each oxidation state. It is shown that the resulting discontinuous dependency of excess charge on chemical potential is consistent with the statistical thermodynamics of equilibrium gas-phase reactions. The method is illustrated by an application to the adiabatic redox dynamics of an aniline molecule.

Journal Article↗