Acidity constants from DFT-based molecular dynamics simulations

Ab initio molecular dynamics free energy study of enhanced copper (II) dimerization on mineral surfaces

Understanding the adsorption of isolated metal cations from water on to mineral surfaces is critical for toxic waste retention and cleanup in the environment. Heterogeneous nucleation of metal oxyhydroxides and other minerals on material surfaces is key to crystal growth and dissolution. The link connecting these two areas, namely cation dimerization and polymerization, is far less understood. In this work we apply ab initio molecular dynamics calculations to examine the coordination structure of hydroxide-bridged Cu(II) dimers, and the free energy changes associated with Cu(II) dimerization on silica surfaces. The dimer dissociation pathway involves sequential breaking of two Cu2+-OH- bonds, yielding three local minima in the free energy profiles associated with 0-2 OH- bridges between the metal cations, and requires the design of a (to our knowledge) novel reaction coordinate for the simulations. Cu(II) adsorbed on silica surfaces are found to exhibit stronger tendency towards dimerization than when residing in water. Cluster-plus-implicit-solvent methods yield incorrect trends if OH- hydration is not correctly depicted. The predicted free energy landscapes are consistent with fast equilibrium times (seconds) among adsorbed structures, and favor Cu2+ dimer formation on silica surfaces over monomer adsorption.

Automated Workflow for Computation of Redox Potentials, Acidity Constants and Solvation Free Energies Accelerated by Machine Learning

Fast evolution of modern society stimulates intense development of new materials with novel functionalities in energy and environmental applications. Due to rapid progress of computer science, computational design of materials with target properties has recently attracted lots of interests. Accurate and efficient calculation of fundamental thermodynamic properties, including redox potentials, acidity constants, and solvation free energies, is of great importance for selection and design of desirable materials. Free energy calculation based on ab initio molecular dynamics (AIMD) can predict these properties with high accuracy at complex environments, however being impeded by high computational costs. To address this issue, this work develops an automated scheme that combines iterative training of machine learning potentials (MLPs) and free energy calculation, and deomonstrates that these thermodynamic properties can be computed by ML accelerated MD with ab initio accuracy and much longer time scale at cheaper costs, improving poor statistics and convergence of numerical integration by AIMD. Our automated scheme lays the foundation for computational chemistry-assisted materials design.

Lower degree of dissociation of pyruvic acid at water surfaces than in bulk

Understanding the acid/base behavior of environmentally relevant organic acids is of key relevance for accurate climate modelling. Here we investigate the effect of pH on the (de)protonation state of pyruvic acid at the air-water interface and in bulk by using the analytical techniques surface-specific vibrational sum frequency generation and attenuated total reflection spectroscopy. To provide a molecular interpretation of the observed behavior, simulations are carried out using a free energy perturbation approach in combination with electronic structure-based molecular dynamics. In both the experimental and theoretical results we observe that the protonated form of pyruvic acid is preferred at the air-water interface. The increased proton affinity is the result of the specific microsolvation at the interface.

How to Predict the pK a of Any Compound in Any Solvent

Acid-base properties of molecules in nonaqueous solvents are of critical importance for almost all areas of chemistry. Despite this very high relevance, our knowledge is still mostly limited to the pK a of rather few compounds in the most common solvents, and a simple yet truly general computational procedure to predict pK a's of any compound in any solvent is still missing. In this contribution, we describe such a procedure. Our method requires only the experimental pK a of a reference compound in water and a few standard quantum-chemical calculations. This method is tested through computing the proton solvation energy in 39 solvents and by comparing the pK a of 142 simple compounds in 12 solvents. Our computations indicate that the method to compute the proton solvation energy is robust with respect to the detailed computational setup and the construction of the solvation model. The unscaled pK a's computed using an implicit solvation model on the other hand differ significantly from the experimental data. These differences are partly associated with the poor quality of the experimental data and the well-known shortcomings of implicit solvation models. General linear scaling relationships to correct this error are suggested for protic and aprotic media. Using these relationships, the deviations between experiment and computations drop to a level comparable to that observed in water, which highlights the efficiency of our method.

Arbitrarily accurate quantum alchemy

Doping compounds can be considered a perturbation to the nuclear charges in a molecular Hamiltonian. Expansions of this perturbation in a Taylor series, i.e., quantum alchemy, have been used in the literature to assess millions of derivative compounds at once rather than enumerating them in costly quantum chemistry calculations. So far, it was unclear whether this series even converges for small molecules, whether it can be used for geometry relaxation, and how strong this perturbation may be to still obtain convergent numbers. This work provides numerical evidence that this expansion converges and recovers the self-consistent energy of Hartree-Fock calculations. The convergence radius of this expansion is quantified for dimer examples and systematically evaluated for different basis sets, allowing for estimates of the chemical space that can be covered by perturbing one reference calculation alone. Besides electronic energy, convergence is shown for density matrix elements, molecular orbital energies, and density profiles, even for large changes in electronic structure, e.g., transforming He₃ into H₆. Subsequently, mixed alchemical and spatial derivatives are used to relax H₂ from the electronic structure of He alone, highlighting a path to spatially relaxed quantum alchemy. Finally, the underlying code that allows for arbitrarily accurate evaluation of restricted Hartree-Fock energies and arbitrary order derivatives is made available to support future method development.

Computational Prediction of All Lanthanide Aqua Ion Acidity Constants

The protonation state of lanthanide-ligand complexes, or lanthanide-containing porous materials, with many Brønsted acid sites can change due to proton loss/gain reactions with water or other heteroatom-containing compounds. Consequently, variations in the protonation state of lanthanide-containing species affect their molecular structure and desired properties. Lanthanide(III) aqua ions undergo hydrolysis and form hydroxides; they are the best characterized lanthanide-containing species with multiple Brønsted acid sites. We employed constrained ab initio molecular dynamics simulations and electronic structure calculations to determine all acidity constants of the lanthanide(III) aqua ions solely from computation. The first, second, and third acidity constants of lanthanide(III) aqua ions were predicted, on average, within 1.2, 2.5, and 4.7 absolute pK_a units from experiment, respectively. A table includes our predicted pK_a values alongside most experimentally measured pK_a values known to date. The approach presented is particularly suitable to determine the Brønsted acidity of lanthanide-containing systems with multiple acidic sites, including those whose measured acidity constants cannot be linked to specific acid sites.

An Atomic-Scale Understanding of the Solution Chemistry of Antimony(V): Insights from First-Principles Molecular Dynamics Simulation

In this study, the structure, hydrolysis, and complexation of Sb(V) in aqueous solution has been elucidated by using first-principles molecular dynamics (FPMD) simulations. The results show that both antimonic acid and its deprotonated form have an octahedral configuration, with average Sb-OH₂ and Sb-OH distances of 2.25 and 2.05 Å, respectively. The computed pK_a of [Sb(OH)₅(OH₂)] is 1.8, while [Sb(OH)₆]^- has an extremely high pK_a. Consequently, [Sb(OH)₆]^- is the most dominant species of Sb(V) under common environmental conditions. A stable aqueous complex can form between [Sb(OH)₆]^- and common cations, and an Sb-Al bidentate complex has the largest dissociation free energy, followed by a Sb-Mg bidentate complex, indicating that they have significantly higher stabilities. For Na⁺ and Ca²⁺, their respective monodentate and bidentate complexes have similar dissociation free energies, indicating very close possibilities. These findings provide a comprehensive understanding of the solution chemistry of Sb(V) from a quantitative and microscopic perspective.

Computing pKa Shifts Using Traditional Molecular Dynamics: Example of Acid-Base Indicator Dyes in Organized Solutions

A compound's acidity constant (K_a) in a given medium determines its protonation state and, thus, its behavior and physicochemical properties. Therefore, it is among the key characteristics considered during the design of new compounds for the needs of advanced technology, medicine, and biological research, a notable example being pH sensors. The computational prediction of K_a for weak acids and bases in homogeneous solvents is presently rather well developed. However, it is not the case for more complex media, such as microheterogeneous solutions. The constant-pH molecular dynamics (MD) method is a notable contribution to the solution of the problem, but it is not commonly used. Here, we develop an approach for predicting K_a changes of weak small-molecule acids upon transfer from water to colloid solutions by means of traditional classical molecular dynamics. The approach is based on free energy (ΔG) computations and requires limited experiment data input during calibration. It was successfully tested on a series of pH-sensitive acid-base indicator dyes in micellar solutions of surfactants. The difficulty of finite-size effects affecting ΔG computation between states with different total charges is taken into account by evaluating relevant corrections; their impact on the results is discussed, and it is found non-negligible (0.1-0.4 pK_a units). A marked bias is found in the ΔG values of acid deprotonation, as computed from MD, which is apparently caused by force-field issues. It is hypothesized to affect the constant-pH MD and reaction ensemble MD methods as well. Consequently, for these methods, a preliminary calibration is suggested.

Computing Surface Acidity Constants of Proton Hopping Groups from Density Functional Theory-Based Molecular Dynamics: Application to the SnO2(110)/H2O Interface

Proton transfer at metal oxide/water interfaces plays an important role in electrochemistry, geochemistry, and environmental science. The key thermodynamic quantity to characterize this process is the surface acidity constant. An ab initio method that combines density functional theory-based molecular dynamics (DFTMD) and free energy perturbation theory has been established for computing surface acidity constants. However, it involves a reversible proton insertion procedure in which frequent proton hopping, e.g., for strong bases and some oxide surfaces (e.g., SnO₂), can cause instability issues in electronic structure calculation. In the original implementation, harmonic restraining potentials are imposed on all O-H bonds (denoted by the V_rH scheme) to prevent proton hopping and thus may not be applicable for systems involving spontaneous proton hopping. In this work, we introduce an improved restraining scheme with a repulsive potential V_rep to compute the surface acidities of systems in which proton hopping is spontaneous and fast. In this V_rep scheme, a Buckingham-type repulsive potential V_rep is applied between the deprotonation site and all other protons in DFTMD simulations. We first verify the V_rep scheme by calculating the pK_a values of H₂O and aqueous HS^- solution (i.e., strong conjugate bases) and then apply it to the SnO₂(110)/H₂O interface. It is found that the V_rep scheme leads to a prediction of the point of zero charge (PZC) of 4.6, which agrees well with experiment. The intrinsic individual pK_a values of the terminal five-coordinated Sn site (Sn_5cOH₂) and bridge oxygen site (Sn₂O_brH⁺) are 4.4 and 4.7, respectively, both being almost the same as the PZC. The similarity of the two pK_a values indicates that dissociation of terminal water has almost zero free energy at this proton hopping interface (i.e., partial water dissociation), as expected from the acid-base equilibrium on SnO₂.

Simulating Solvation and Acidity in Complex Mixtures with First-Principles Accuracy: The Case of CH3SO3H and H2O2 in Phenol

We present a generally applicable computational framework for the efficient and accurate characterization of molecular structural patterns and acid properties in an explicit solvent using H₂O₂ and CH₃SO₃H in phenol as an example. To address the challenges posed by the complexity of the problem, we resort to a set of data-driven methods and enhanced sampling algorithms. The synergistic application of these techniques makes the first-principle estimation of the chemical properties feasible without renouncing to the use of explicit solvation, involving extensive statistical sampling. Ensembles of neural network (NN) potentials are trained on a set of configurations carefully selected out of preliminary simulations performed at a low-cost density functional tight-binding (DFTB) level. The energy and forces of these configurations are then recomputed at the hybrid density functional theory (DFT) level and used to train the neural networks. The stability of the NN model is enhanced by using DFTB energetics as a baseline, but the efficiency of the direct NN (i.e., baseline-free) is exploited via a multiple-time-step integrator. The neural network potentials are combined with enhanced sampling techniques, such as replica exchange and metadynamics, and used to characterize the relevant protonated species and dominant noncovalent interactions in the mixture, also considering nuclear quantum effects.

Determination of pKa Values via ab initio Molecular Dynamics and its Application to Transition Metal-Based Water Oxidation Catalysts

The p K a values are important for the in-depth elucidation of catalytic processes, the computational determination of which has been challenging. The first simulation protocols employing ab initio molecular dynamics simulations to calculate p K a values appeared almost two decades ago. Since then several slightly different methods have been proposed. We compare the performance of various evaluation methods in order to determine the most reliable protocol when it comes to simulate p K a values of transition metal-based complexes, such as the here investigated Ru-based water oxidation catalysts. The latter are of high interest for sustainable solar-light driven water splitting, and understanding of the underlying reaction mechanism is crucial for their further development.

Acidity Constants of the Hematite-Liquid Water Interface from Ab Initio Molecular Dynamics

The interface between transition metal oxides (TMO) and liquid water plays a crucial role in environmental chemistry, catalysis, and energy science. Yet, the mechanism and energetics of chemical transformations at solvated TMO surfaces is often unclear, largely because of the difficulty to characterize the active surface species experimentally. The hematite (α-Fe₂O₃)-liquid water interface is a case in point. Here we demonstrate that ab initio molecular dynamics is a viable tool for determining the protonation states of complex interfaces. The p K_a values of the oxygen-terminated (001) surface group of hematite, ≡OH, and half-layer terminated (012) surface groups, ≡²OH and ≡¹OH₂, are predicted to be (18.5 ± 0.3), (18.9 ± 0.6), and (10.3 ± 0.5) p K_a units, respectively. These are in good agreement with recent bond-valence theory based estimates, and suggest that the deprotonation of these surfaces require significantly more free energy input than previously thought.

Mechanistic Understanding of Uranyl Ion Complexation on Montmorillonite Edges: A Combined First-Principles Molecular Dynamics-Surface Complexation Modeling Approach

Systematic first-principles molecular dynamics (FPMD) simulations were carried out to study the structures, free energies, and acidity constants of UO₂²⁺ surface complexes on montmorillonite in order to elucidate the surface complexation mechanisms of the uranyl ion (UO₂²⁺) on clay mineral edges at the atomic scale. Four representative complexing sites were investigated, that is, ≡Al(OH)₂ and ≡AlOHSiO on the (010) surface and ≡AlOHOa and ≡SiOOa on the (110) surface. The results show that uranyl ions form bidentate complexes on these sites. All calculated binding free energies for these complexes are very similar. These bidentate complexes can be hydrolyzed, and their corresponding derived p K_a values (around 5.0 and 9.0 for p K_a1 and p K_a2, respectively) indicate that UO₂(OH)⁺ and UO₂(OH)₂ surface groups are the dominant surface species in the environmental pH range. The OH groups of UO₂(OH)₂ surface complexes can act as complexing sites for subsequent metals. Additional simulations showed that such multinuclear adsorption is feasible and can be important at high pH. Furthermore, FPMD simulation results served as input parameters for an electrostatic thermodynamic surface complexation model (SCM) that adequately reproduced adsorption data from the literature. Overall, this study provides an improved understanding of UO₂²⁺ complexation on clay mineral edge surfaces.

Atypical titration curves for GaAl12 Keggin-ions explained by a joint experimental and simulation approach

Although they have been widely used as models for oxide surfaces, the deprotonation behaviors of the Keggin-ions (MeAl₁₂⁷⁺) and typical oxide surfaces are very different. On Keggin-ions, the deprotonation occurs over a very narrow pH range at odds with the broad charging curve of larger oxide surfaces. Depending on the Me concentration, the deprotonation curve levels off sooner (high Me concentration) or later (for low Me concentration). The leveling off shows the onset of aggregation before which the Keggin-ions are present as individual units. We show that the atypical titration data previously observed for some GaAl₁₂ solutions in comparison to the originally reported data can be explained by the presence of Ga₂Al₁₁ ions. The pK_a value of aquo-groups bound to octahedral Ga was determined from ab initio molecular dynamics simulations relative to the pure GaAl₁₂ ions. Using these results within a surface complexation model, the onset of deprotonation of the crude solution is surprisingly well predicted and the ratio between the different species is estimated to be in the proportion 20 (Ga₂Al₁₁) : 20 (Al₁₃) : 60 (GaAl₁₂).

Uranyl Arsenate Complexes in Aqueous Solution: Insights from First-Principles Molecular Dynamics Simulations

In this study, the structures and acidity constants (p K_a's) of uranyl arsenate complexes in solutions have been revealed by using the first principle molecular dynamics technique. The results show that uranyl and arsenate form stable complexes with the U/As ratios of 1:1 and 1:2, and the bidentate complexation between U and As is highly favored. Speciation-pH distributions are derived based on free energy and p K_a calculations, which indicate that for the 1:1 species, UO₂(H₂AsO₄)(H₂O)₃⁺ is the major species at pH < 7, while UO₂(HAsO₄)(H₂O)₃⁰ and UO₂(AsO₄)(H₂O)₃^- dominate in acid-to-alkaline and extreme alkaline pH ranges. For the 1:2 species, UO₂(H₂AsO₄)₂(H₂O)⁰ is dominant under acid-to-neutral pH conditions, while UO₂(HAsO₄)(H₂AsO₄)(H₂O)^-, UO₂(HAsO₄)(HAsO₄)(H₂O)^2-, and UO₂(AsO₄)(HAsO₄)(H₂O)^3- become the major forms in the pH range of 7.2-10.7, 10.7-12.1, and >12.1, respectively.

Computational Modeling of Cobalt-Based Water Oxidation: Current Status and Future Challenges

A lot of effort is nowadays put into the development of novel water oxidation catalysts. In this context, mechanistic studies are crucial in order to elucidate the reaction mechanisms governing this complex process, new design paradigms and strategies how to improve the stability and efficiency of those catalysts. This review is focused on recent theoretical mechanistic studies in the field of homogeneous cobalt-based water oxidation catalysts. In the first part, computational methodologies and protocols are summarized and evaluated on the basis of their applicability toward real catalytic or smaller model systems, whereby special emphasis is laid on the choice of an appropriate model system. In the second part, an overview of mechanistic studies is presented, from which conceptual guidelines are drawn on how to approach novel studies of catalysts and how to further develop the field of computational modeling of water oxidation reactions.

Lead and selenite adsorption at water-goethite interfaces from first principles

The complexation of toxic and/or radioactive ions on to mineral surfaces is an important topic in geochemistry. We apply periodic-boundary-conditions density functional theory (DFT) molecular dynamics simulations to examine the coordination of Pb(II), [Formula: see text], and their contact ion pairs to goethite (1 0 1) and (2 1 0) surfaces. The multitude of Pb(II) adsorption sites and possibility of Pb(II)-induced FeOH deprotonation make this a complex problem. At surface sites where Pb(II) is coordinated to three FeO and/or FeOH groups, and with judicious choices of FeOH surface group protonation states, the predicted Fe-Pb distances are in good agreement with EXAFS measurements. Trajectories where Pb(II) is in part coordinated to only two surface Fe-O groups exhibit larger fluctuations in Pb-O distances. Pb(II)/[Formula: see text] contact ion pairs are at least metastable on goethite (2 1 0) surfaces if the [Formula: see text] has a monodentate Se-O-Fe bond. Our DFT-based molecular dynamics calculations are a prerequisite for calculations of finite temperature equilibrium binding constants of Pb(II) and Pb(II)/[Formula: see text] ion pairs to goethite adsorption sites.

Microscopic dynamics of charge separation at the aqueous electrochemical interface

We have used molecular simulation and methods of importance sampling to study the thermodynamics and kinetics of ionic charge separation at a liquid water-metal interface. We have considered this process using canonical examples of two different classes of ions: a simple alkali-halide pair, Na⁺I^-, or classical ions, and the products of water autoionization, H₃O⁺OH^-, or water ions. We find that for both ion classes, the microscopic mechanism of charge separation, including water's collective role in the process, is conserved between the bulk liquid and the electrode interface. However, the thermodynamic and kinetic details of the process differ between these two environments in a way that depends on ion type. In the case of the classical ion pairs, a higher free-energy barrier to charge separation and a smaller flux over that barrier at the interface result in a rate of dissociation that is 40 times slower relative to the bulk. For water ions, a slightly higher free-energy barrier is offset by a higher flux over the barrier from longer lived hydrogen-bonding patterns at the interface, resulting in a rate of association that is similar both at and away from the interface. We find that these differences in rates and stabilities of charge separation are due to the altered ability of water to solvate and reorganize in the vicinity of the metal interface.

Dehydrogenation Free Energy of Co2+(aq) from Density Functional Theory-Based Molecular Dynamics

Electron and proton transfers are important steps occurring in chemical reactions. The often used approach of calculating the energy differences of those steps using methods based on geometry optimizations neglects the influence of dynamic effects. To further investigate this issue and inspired by research in water oxidation, we calculate in the present study the dehydrogenation free energy of aqueous Co²⁺, which is the free energy change associated with the first step of the water oxidation reaction mechanism of recently investigated model Co(II)-aqua catalysts. We employ a method based on a thermodynamic integration scheme with strong ties to Marcus theory to obtain free energy differences, solvent reorganization free energies, and dynamic structural information on the systems from density functional theory-based molecular dynamics. While this method is computationally orders of magnitude more expensive than a static approach, it potentially allows for predicting the validity of the approximation of neglecting dynamic effects.

Structures and Acidity Constants of Silver-Sulfide Complexes in Hydrothermal Fluids: A First-Principles Molecular Dynamics Study

In order to quantify the speciation and structures of silver-sulfide complexes in aqueous solutions, we have carried out systematic first-principles molecular dynamics (FPMD) simulations at three temperatures (25, 200, and 300 °C). It is found that monosulfide (i.e., Ag(HS)) and disulfide species (i.e., Ag(HS)₂^-) are the major silver-sulfide species over a wide T-P range, while Ag(HS)₃^2- can hold stably only at ambient temperatures, and Ag(HS)₄^3- does not exist even at the ambient conditions. Ag(H₂S)⁺ has a tetrahedral structure up to 300 °C (i.e., Ag(H₂S)(H₂O)₃⁺). Ag(H₂S)₂⁺ remains 4-coordinated to 200 °C (i.e., Ag(H₂S)₂(H₂O)₂⁺), but it transforms to 3-coordinated at 300 °C (i.e., Ag(H₂S)₂(H₂O)⁺). All of the other mono- and disulfide species (Ag(HS)(H₂O)⁰, Ag(HS)(OH)^-, Ag(HS)(H₂S)⁰, Ag(HS)₂^-, and AgS(HS)^2-) have 2-fold linear structures. For their solvation structures, the H₂S ligands donate weak H-bonds to water O; the HS^- ligands accept weak H-bonds from water H; the dangling S^2- form strong H-bonds with H of water molecules, and the OH^- ligands can form strong H-bonds as donors and weak H-bonds as acceptors. We further calculated the acidity constants (i.e., pK_as) of Ag(H₂S)⁺ and Ag(H₂S)₂⁺ complexes using FPMD based vertical energy gap method. Based on the calculated pK_as, the mono- and disulfide species distributions versus pH have been derived. We found that for monosulfide species, Ag(HS)(H₂O)⁰, is the major species in near neutral pH, while Ag(H₂S)(H₂O)₃⁺ and Ag(HS)(OH)^- exist in the acid and alkaline pH range at T ≤ 200 °C, respectively. At 300 °C, both Ag(HS)(OH)^- and Ag(HS)(H₂O)⁰ are dominant in the neutral pH range, and Ag(H₂S)(H₂O)₂⁺ only exists in acidic solutions. For disulfide species, Ag(HS)₂^- is dominative in near neutral pH condition at the three temperatures; Ag(HS)(H₂S)⁰ stays in mild acidic pH range only at 25 °C; AgS(HS)^2- and Ag(H₂S)₂(H₂O)₂⁺ (Ag(H₂S)₂(H₂O)⁺ at 300 °C) are trivial at the three conditions. The results of structures and acidity constants provide quantitative and microscopic basis for understanding the behavior of silver complexes in hydrothermal fluids.

Acidity Constant (pKa ) Calculation of Large Solvated Dye Molecules: Evaluation of Two Advanced Molecular Dynamics Methods

pH‐Sensitive dyes are increasingly applied on polymer substrates for the creation of novel sensor materials. Recently, these dye molecules were modified to form a covalent bond with the polymer host. This had a large influence on the pH‐sensitive properties, in particular on the acidity constant (pK_a). Obtaining molecular control over the factors that influence the pK_a value is mandatory for the future intelligent design of sensor materials. Herein, we show that advanced molecular dynamics (MD) methods have reached the level at which the pK_a values of large solvated dye molecules can be predicted with high accuracy. Two MD methods were used in this work: steered or restrained MD and the insertion/deletion scheme. Both were first calibrated on a set of phenol derivatives and afterwards applied to the dye molecule bromothymol blue. Excellent agreement with experimental values was obtained, which opens perspectives for using these methods for designing dye molecules.

Dispersion corrected RPBE studies of liquid water

The structure of liquid water has been addressed by ab initio molecular dynamics simulations based on density functional theory. Exchange-correlation effects have been described by the popular PBE and RPBE functionals within the generalized gradient approximation as these functionals also yield satisfactory results for metals which is important to model electrochemical interfaces from first principles. In addition, dispersive interactions are included by using dispersion-corrected schemes. It turns out that the dispersion-corrected RPBE functional reproduces liquid water properties quite well in contrast to the PBE functional. This is caused by the replacement of the over-estimated directional hydrogen-bonding in the PBE functional by non-directional dispersive interactions.

The amorphous silica-liquid water interface studied by ab initio molecular dynamics (AIMD): local organization in global disorder

The structural organization of water at a model of amorphous silica-liquid water interface is investigated by ab initio molecular dynamics (AIMD) simulations at room temperature. The amorphous surface is constructed with isolated, H-bonded vicinal and geminal silanols. In the absence of water, the silanols have orientations that depend on the local surface topology (i.e. presence of concave and convex zones). However, in the presence of liquid water, only the strong inter-silanol H-bonds are maintained, whereas the weaker ones are replaced by H-bonds formed with interfacial water molecules. All silanols are found to act as H-bond donors to water. The vicinal silanols are simultaneously found to be H-bond acceptors from water. The geminal pairs are also characterized by the formation of water H-bonded rings, which could provide special pathways for proton transfer(s) at the interface. The first water layer above the surface is overall rather disordered, with three main domains of orientations of the water molecules. We discuss the similarities and differences in the structural organization of the interfacial water layer at the surface of the amorphous silica and at the surface of the crystalline (0 0 0 1) quartz surface.

DRASTIC INFLUENCE OF BORON ATOM ON THE ACIDITY OF ALCOHOL IN BOTH GAS PHASE AND SOLUTION PHASE, A DFT STUDY

In this study, the drastic influence of the boron atom on the acidity of alcohol has been considered. The calculated ΔH acid (320.9–338.1 kcal/mol) and pKa range of boron containing alcohol (-0.1–9.4) indicate that the boronation of alcohol leads to considerable enhancement of its acidity. For instance, we have obtained the ΔH acid values 338.1, 335.2 kcal/mol and the pKa values 4.12, 2.81 for BH2CH2OH , BF2CH2OH alcohols, respectively, which are much smaller than that of CH3OH (with ΔH acid = 374.9 kcal/mol and pKa = 15). The increase in the acidity of boronated alcohol can be related to the stabilization of alkoxy ion due to overlap of unoccupied orbital of boron atom with the electron pairs of negative oxygen. All gas phase computations were performed at MP2/6-311++G(d,p)//(B3LYP/6-31+G(d)) level. The primary results indicate that the presence of boron atom in an alcohol might make it as acidic as nitric acid. The geometry optimization of studied structures was performed with DFT computation and optimized structures were used to carry out natural bond orbital (NBO) analysis. NBO analysis revealed that the increase in the acidity of boron-containing alcohols is due to the charge transfer from the negative oxygen (in deprotonated structure) to the empty orbital of - BH2 and - BF2 . Quantum theory of atoms in molecules (QTAIM) was also applied to determine the nature of bonds formed in the deprotonated structure.

Predicting the acidity constant of a goethite hydroxyl group from first principles

Accurate predictions of the acid-base behavior of hydroxyl groups at mineral surfaces are critical for understanding the trapping of toxic and radioactive ions in soil samples. In this work, we apply ab initio molecular dynamics (AIMD) simulations and potential-of-mean-force techniques to calculate the pK(a) of a doubly protonated oxygen atom bonded to a single Fe atom (Fe(I)OH(2)) on the goethite (101) surface. Using formic acid as a reference system, pK(a) = 7.0 is predicted, suggesting that isolated, positively charged groups of this type are marginally stable at neutral pH. Similarities and differences between AIMD and the more empirical multi-site complexation methodology are highlighted, particularly with respect to the treatment of hydrogen bonding with water and proton sharing among surface hydroxyl groups. We also highlight the importance of an electronic structure method that can accurately predict transition metal ion properties for goethite pK(a) calculations.

Oxide/water interfaces: how the surface chemistry modifies interfacial water properties

The organization of water at the interface with silica and alumina oxides is analysed using density functional theory-based molecular dynamics simulation (DFT-MD). The interfacial hydrogen bonding is investigated in detail and related to the chemistry of the oxide surfaces by computing the surface charge density and acidity. We find that water molecules hydrogen-bonded to the surface have different orientations depending on the strength of the hydrogen bonds and use this observation to explain the features in the surface vibrational spectra measured by sum frequency generation spectroscopy. In particular, 'ice-like' and 'liquid-like' features in these spectra are interpreted as the result of hydrogen bonds of different strengths between surface silanols/aluminols and water.

Thermodynamics of chemical reactions with COSMO-RS: the extreme case of charge separation or recombination

Many technically relevant chemical processes in the condensed phase involve as elementary reactive steps the formation of ions from neutral species or, as the opposite, recombination of ions. Such reactions that generate or annihilate charge defy the standard gas phase quantum chemical treatment, and also continuum solvation models are only partially able to account for the right amount of stabilization in solution. In this work, for such types of reaction, a solvation treatment involving the COSMO-RS method is assessed, which leads to improved results, i.e., errors of only around 10 kJ/mol for both protic and aprotic solvents. The examples discussed here comprise protolysis reactions and organo halide heterolysis, for both of which a comparison with reliable experimental data is possible. It is observed that for protolysis, the quality of results does not strongly depend on the quantum chemical method used for energy calculation. In contrast, in the case of heterolytic carbon-chlorine bond cleavage, clearly better results are obtained for higher correlated (coupled cluster) methods or the density functional M06-2X, which is well known for its accuracy if applied to organic chemistry. This hints at least that the right answer is obtained for the right reason and not due to a compensation of errors from gas phase thermodynamics with those from the solvation treatment. Problems encountered with certain critical solvents or upon decomposing Gibbs free energies into heats or entropies of reaction are found to relate mostly to the parameterization of the H-bonding term within COSMO-RS.