Marialore Sulpizi

Redox potentials and pKa for benzoquinone from density functional theory based molecular dynamics

The density functional theory based molecular dynamics (DFTMD) method for the computation of redox free energies presented in previous publications and the more recent modification for computation of acidity constants are reviewed. The method uses a half reaction scheme based on reversible insertion/removal of electrons and protons. The proton insertion is assisted by restraining potentials acting as chaperones. The procedure for relating the calculated deprotonation free energies to Brønsted acidities (pK(a)) and the oxidation free energies to electrode potentials with respect to the normal hydrogen electrode is discussed in some detail. The method is validated in an application to the reduction of aqueous 1,4-benzoquinone. The conversion of hydroquinone to quinone can take place via a number of alternative pathways consisting of combinations of acid dissociations, oxidations, or dehydrogenations. The free energy changes of all elementary steps (ten in total) are computed. The accuracy of the calculations is assessed by comparing the energies of different pathways for the same reaction (Hesss law) and by comparison to experiment. This two-sided test enables us to separate the errors related with the restrictions on length and time scales accessible to DFTMD from the errors introduced by the DFT approximation. It is found that the DFT approximation is the main source of error for oxidation free energies.

Vibrational Sum Frequency Generation Spectroscopy of the Water Liquid–Vapor Interface from Density Functional Theory-Based Molecular Dynamics Simulations

The vibrational sum frequency generation (VSFG) spectrum of the water liquid-vapor (LV) interface is calculated using density functional theory-based molecular dynamics simulations. The real and imaginary parts of the spectrum are in good agreement with the experimental data, and we provide an assignment of the SFG bands according to the dipole orientation of the interfacial water molecules. We use an instantaneous definition of the surface, which is more adapted to the study of interfacial phenomena than the Gibbs dividing surface. By calculating the vibrational (infrared, Raman) properties for interfaces of varying thickness, we show that the bulk spectra signatures appear after a thin layer of 2-3 Å only. We therefore use this value as a criterion for calculating the VSFG spectrum.

The oxidation of tyrosine and tryptophan studied by a molecular dynamics normal hydrogen electrode

The thermochemical constants for the oxidation of tyrosine and tryptophan through proton coupled electron transfer in aqueous solution have been computed applying a recently developed density functional theory (DFT) based molecular dynamics method for reversible elimination of protons and electrons. This method enables us to estimate the solvation free energy of a proton (H(+)) in a periodic model system from the free energy for the deprotonation of an aqueous hydronium ion (H(3)O(+)). Using the computed solvation free energy of H(+) as reference, the deprotonation and oxidation free energies of an aqueous species can be converted to pK(a) and normal hydrogen electrode (NHE) potentials. This conversion requires certain thermochemical corrections which were first presented in a similar study of the oxidation of hydrobenzoquinone [J. Cheng, M. Sulpizi, and M. Sprik, J. Chem. Phys. 131, 154504 (2009)]. Taking a different view of the thermodynamic status of the hydronium ion, these thermochemical corrections are revised in the present work. The key difference with the previous scheme is that the hydronium is now treated as an intermediate in the transfer of the proton from solution to the gas-phase. The accuracy of the method is assessed by a detailed comparison of the computed pK(a), NHE potentials and dehydrogenation free energies to experiment. As a further application of the technique, we have analyzed the role of the solvent in the oxidation of tyrosine by the tryptophan radical. The free energy change computed for this hydrogen atom transfer reaction is very similar to the gas-phase value, in agreement with experiment. The molecular dynamics results however, show that the minimal solvent effect on the reaction free energy is accompanied by a significant reorganization of the solvent.

Redox potentials and acidity constants from density functional theory based molecular dynamics.

CONSPECTUS: All-atom methods treat solute and solvent at the same level of electronic structure theory and statistical mechanics. All-atom computation of acidity constants (pKa) and redox potentials is still a challenge. In this Account, we review such a method combining density functional theory based molecular dynamics (DFTMD) and free energy perturbation (FEP) methods. The key computational tool is a FEP based method for reversible insertion of a proton or electron in a periodic DFTMD model system. The free energy of insertion (work function) is computed by thermodynamic integration of vertical energy gaps obtained from total energy differences. The problem of the loss of a physical reference for ionization energies under periodic boundary conditions is solved by comparing with the proton work function computed for the same supercell. The scheme acts as a computational hydrogen electrode, and the DFTMD redox energies can be directly compared with experimental redox potentials. Consistent with the closed shell nature of acid dissociation, pKa estimates computed using the proton insertion/removal scheme are found to be significantly more accurate than the redox potential calculations. This enables us to separate the DFT error from other sources of uncertainty such as finite system size and sampling errors. Drawing an analogy with charged defects in solids, we trace the error in redox potentials back to underestimation of the energy gap of the extended states of the solvent. Accordingly the improvement in the redox potential as calculated by hybrid functionals is explained as a consequence of the opening up of the bandgap by the Hartree-Fock exchange component in hybrids. Test calculations for a number of small inorganic and organic molecules show that the hybrid functional implementation of our method can reproduce acidity constants with an uncertainty of 1-2 pKa units (0.1 eV). The error for redox potentials is in the order of 0.2 V.

Absolute pKa Values and Solvation Structure of Amino Acids from Density Functional Based Molecular Dynamics Simulation

Absolute pKa values of the amino acid side chains of arginine, aspartate, cysteine, histidine, and tyrosine; the C- and N-terminal group of tyrosine; and the tryptophan radical cation are calculated using a revised density functional based molecular dynamics simulation technique introduced previously [ Cheng , J. ; Sulpizi , M. ; Sprik , M. J. Chem. Phys. 2009 , 131 , 154504 ]. In the revised scheme, acid deprotonation is considered as a dissociation rather than a proton transfer reaction, and a correction term for treating the proton as a hydronium ion is suggested. The acidity constants of the amino acids are obtained from the vertical energy gaps for removal or insertion of the acidic proton and the computed solvation free energy of the proton. The unsigned mean error relative to experimental results is 2.1 pKa units with a maximum error of 4.0 pKa units. The estimated mean statistical uncertainty due to the finite length of the trajectories is ±1.1 pKa units. The solvation structures of the protonated and deprotonated amino acids are analyzed in terms of radial distribution functions, which can serve as reference data for future force field developments.

Hybrid QM/MM Car-Parrinello simulations of catalytic and enzymatic reactions

A review. First-principles mol. dynamics (Car-Parrinello) simulations based on d. functional theory have emerged as a powerful tool for the study of phys., chem. and biol. systems. At present, using parallel computers, systems of a few hundreds of atoms can be routinely investigated. By extending this method to a mixed quantum mech. - mol. mech. (QM/MM) hybrid scheme, the system size can be enlarged further. Such an approach is esp. attractive for the in situ investigation of chem. reactions that occur in a complex and heterogeneous environment. Here, we review some recent applications of hybrid Car-Parrinello simulations of chem. and biol. systems as illustrative examples of the current potential and limitations of this promising novel technique. [on SciFinder (R)]

Aqueous Redox Chemistry and the Electronic Band Structure of Liquid Water

The electronic states of aqueous species can mix with the extended states of the solvent if they are close in energy to the band edges of water. Using density functional theory-based molecular dynamics simulation, we show that this is the case for OH(-) and Cl(-). The effect is, however, badly exaggerated by the generalized gradient approximation leading to systematic underestimation of redox potentials and spurious nonlinearity in the solvent reorganization. Drawing a parallel to charged defects in wide gap solid oxides, we conclude that misalignment of the valence band of water is the main source of error turning the redox levels of OH(-) and Cl(-) in resonant impurity states. On the other hand, the accuracy of energies of levels corresponding to strongly negative redox potentials is acceptable. We therefore predict that mixing of the vertical attachment level of CO2 and the unoccupied states of water is a real effect.

Acidity constants from DFT-based molecular dynamics simulations

In this contribution we review our recently developed method for the calculation of acidity constants from density functional theory based molecular dynamics simulations. The method is based on a half reaction scheme in which protons are formally transferred from solution to the gas phase. The corresponding deprotonation free energies are computed from the vertical energy gaps for insertion or removal of protons. Combined to full proton transfer reactions, the deprotonation energies can be used to estimate relative acidity constants and also the Brønsted pK(a) when the deprotonation free energy of a hydronium ion is used as a reference. We verified the method by investigating a series of organic and inorganic acids and bases spanning a wide range of pK(a) values (20 units). The thermochemical corrections for the biasing potentials assisting and directing the insertion are discussed in some detail.

The electron attachment energy of the aqueous hydroxyl radical predicted from the detachment energy of the aqueous hydroxide anion.

Combining photoemission and electrochemical data from the literature we argue that the difference between the vertical and adiabatic ionization energy of the aqueous hydroxide anion is 2.9 eV. We then use density functional theory based molecular dynamics to show that the solvent response to ionization is nonlinear. Adding this to the experimental data we predict a 4.1 eV difference between the energy for vertical attachment of an electron to the aqueous hydroxyl radical and the corresponding adiabatic electron affinity. This places the state accepting the electron only 2.2 eV below vacuum or 7.7 eV above the edge of the valence band of water.

Understanding the Microscopic Origin of Gold Nanoparticle Anisotropic Growth from Molecular Dynamics Simulations

We use molecular dynamics simulations in order to understand the microscopic origin of the asymmetric growth mechanism in gold nanorods. We provide the first atomistic model of different surfaces on gold nanoparticles in a growing electrolyte solution, and we describe the interaction of the metal with the surfactants, namely, cetyltrimethylammonium bromide (CTAB) and the ions. An innovative aspect is the inclusion of the role of the surfactants, which are explicitly modeled. We find that on all the investigated surfaces, namely, (111), (110), and (100), CTAB forms a layer of distorted cylindrical micelles where channels among micelles provide direct ion access to the surface. In particular, we show how AuCl2(-) ions, which are found in the growth solution, can freely diffuse from the bulk solution to the gold surface. We also find that the (111) surface exhibits a higher CTAB packing density and a higher electrostatic potential. Both elements would favor the growth of gold nanoparticles along the (111) direction. These findings are in agreement with the growth mechanisms proposed by the experimental groups of Murphy and Mulvaney.