Analysis of Bonding by Quantum Chemistry─Resolving Delocalization Stabilization in a Mechanistic Basis and New Hückel Model.

No general and unique understanding of the mechanism of covalent bonding in physical terms is provided by current computational methods or by a consensus among experts. Bonding is studied by energy decomposition analysis but may also be related to the interatomic motion of valence electrons within the molecule. This dynamical view of the mechanism of bonding is not widely appreciated. The aim here is to make it accessible by translation into a corresponding form of quantum chemical energy analysis. The interatomic electron motion is directly related to the delocalization taking place when atomic basis functions are combined into molecular orbitals. A "tribasis method" is introduced, allowing an atomic basis set to form subsets of (1) strictly localized atomic functions and (2) interatomic bridge functions which allow delocalization. Calculations can then identify ground states without (no bridge functions) and with delocalization. The scheme is based on exact quantum mechanics but demonstrated by a minimal basis treatment of H2+ and H2 in Hartree-Fock and valence bond approximations which show that the bond energy is a sum of repulsive localization and more strongly attractive delocalization energies. The tribasis method is used to reconstruct the Hückel theory of π-electron delocalization in planar hydrocarbon molecules to account for the "overlap problem". In its empirically fitted form, the new theory can accurately resolve both π → π* transition energy and aromatic stabilization energy. The picture of covalent bonding emerging from both hydrogenic and Hückel calculations is that there is a presence of a Pauli repulsion of localization which is overcome by a roughly twice as strong delocalization stabilization to form the bond.

From Electronegativity towards Reactivity-Searching for a Measure of Atomic Reactivity.

Pauling introduced the concept of electronegativity of an atom which has played an important role in understanding the polarity and ionic character of bonds between atoms. We set out to define a related concept of atomic reactivity in such a way that it can be quantified and used to predict the stability of covalent bonds in molecules. Guided by the early definition of electronegativity by Mulliken in terms of first ionization energies and Pauling in terms of bond energies, we propose corresponding definitions of atomic reactivity. The main goal of clearly distinguishing the inert gas atoms as nonreactive is fulfilled by three different proposed measures of atomic reactivity. The measure likely to be found most useful is based on the bond energies in atomic hydrides, which are related to atomic reactivities by a geometric average. The origin of the atomic reactivity is found in the symmetry of the atomic environment and related conservation laws which are also the origin of the shell structure of atoms and the periodic table. The reactive atoms are characterized by degenerate or nearly degenerate (several states of the same or nearly the same energy) ground states, while the inert atoms have nondegenerate ground states and no near-degeneracies. We show how to extend the use of the Aufbau model of atomic structure to qualitatively describe atomic reactivity in terms of ground state degeneracy. The symmetry and related conservation laws of atomic electron structures produce a strain (energy increase) in the structure, which we estimate by use of the Thomas-Fermi form of DFT implemented approximately with and without the symmetry and conservation constraints. This simplified and approximate analysis indicates that the total strain energy of an atom correlates strongly with the corresponding atomic reactivity measures but antibonding mechanisms prevent full conversion of strain relaxation to bonding.

The Basics of Covalent Bonding in Terms of Energy and Dynamics.

We address the paradoxical fact that the concept of a covalent bond, a cornerstone of chemistry which is well resolved computationally by the methods of quantum chemistry, is still the subject of debate, disagreement, and ignorance with respect to its physical origin. Our aim here is to unify two seemingly different explanations: one in terms of energy, the other dynamics. We summarize the mechanistic bonding models and the debate over the last 100 years, with specific applications to the simplest molecules: H2+ and H2. In particular, we focus on the bonding analysis of Hellmann (1933) that was brought into modern form by Ruedenberg (from 1962 on). We and many others have helped verify the validity of the Hellmann-Ruedenberg proposal that a decrease in kinetic energy associated with interatomic delocalization of electron motion is the key to covalent bonding but contrary views, confusion or lack of understanding still abound. In order to resolve this impasse we show that quantum mechanics affords us a complementary dynamical perspective on the bonding mechanism, which agrees with that of Hellmann and Ruedenberg, while providing a direct and unifying view of atomic reactivity, molecule formation and the basic role of the kinetic energy, as well as the important but secondary role of electrostatics, in covalent bonding.

A classical density functional theory for the asymmetric restricted primitive model of ionic liquids.

A new three-parameter (valency, ion size, and charge asymmetry) model, the asymmetric restricted primitive model (ARPM) of ionic liquids, has recently been proposed. Given that ionic liquids generally are composed of monovalent species, the ARPM effectively reduces to a two-parameter model. Monte Carlo (MC) simulations have demonstrated that the ARPM is able to reproduce key properties of room temperature ionic liquids (RTILs) in bulk and at charged surfaces. The relatively modest complexity of the model raises the possibility, which is explored here, that a classical density functional theory (DFT) could resolve its properties. This is relevant because it might generate great improvements in terms of both numerical efficiency and understanding in the continued research of RTILs and their applications. In this report, a DFT for rod-like molecules is proposed as an approximate theoretical tool for an ARPM fluid. Borrowing data on the ion pair fraction from a single bulk simulation, the ARPM is modelled as a mixture of dissociated ions and connected ion pairs. We have specifically studied an ARPM where the hard-sphere diameter is 5 Å, with the charge located 1 Å from the hard-sphere centre. We focus on fluid structure and electrochemical behaviour of this ARPM fluid, into which a model electrode is immersed. The latter is modelled as a perfect conductor, and surface polarization is handled by the method of image charges. Approximate methods, which were developed in an earlier study, to take image interactions into account, are also incorporated in the DFT. We make direct numerical comparisons between DFT predictions and corresponding simulation data. The DFT theory is implemented both in the normal mean field form with respect to the electrostatic interactions and in a correlated form based on hole formation by both steric repulsions and ion-ion Coulomb interactions. The results clearly show that ion-ion correlations play a very important role in the screening of the charged surfaces by our ARPM ionic liquid. We have studied electrostatic potentials and ion density profiles as well the differential capacitance. The mean-field DFT fails to reproduce these properties, but the inclusion of ion-ion correlation by a simple approximate treatment yields quite reasonable agreement with the corresponding simulation results. An interesting finding is that there appears to be a surface phase transition at relatively low surface charge which is readily explored by DFT, but seen also in the MC simulations at somewhat higher asymmetry.

The Virial Theorem and Covalent Bonding.

A long-held view of the origin of covalent binding is based on the notion that electrostatic forces determine the stability of a system of charged particles and that, therefore, potential energy changes drive the stabilization of molecules. A key argument advanced for this conjecture is the rigorous validity of the virial theorem. Rigorous in-depth analyses have however shown that the energy lowering of covalent bonding is due to the wave mechanical drive of electrons to lower their kinetic energy through expansion. Since the virial theorem applies only to systems with Coulombic interaction potentials, its relevance as a foundation of the electrostatic view is tested here by calculations on analogues of the molecules H2+ and H2, where all 1/ r interaction potentials are replaced by Gaussian-type potentials that yield one-electron "atoms" with realistic stability ranges. The virial theorem does not hold in these systems, but covalent bonds are found to form nonetheless, and the wave mechanical bonding analysis yields analogous results as in the case of the Coulombic potentials. Notably, the key driving feature is again the electron delocalization that lowers the interatomic kinetic energy component. A detailed discussion of the role of the virial theorem in the context of covalent binding is given.

Ionic liquid interface at an electrode: simulations of electrochemical properties using an asymmetric restricted primitive model.

We use Monte Carlo simulations of a coarse-grained model to investigate structure and electrochemical behaviours at an electrode immersed in room temperature ionic liquids (RTILs). The simple RTIL model, which we denote the asymmetric restricted primitive model (ARPM), is composed of monovalent hard-sphere ions, all of the same size, in which the charge is asymmetrically placed. Not only the hard-sphere size (d), but also the charge displacement (b), is identical for all species, i.e. the monovalent RTIL ions are fully described by only two parameters (d, b). In earlier work, it was demonstrated that the ARPM can capture typical static RTIL properties in bulk solutions with remarkable accuracy. Here, we investigate its behaviour at an electrode surface. The electrode is assumed to be a perfect conductor and image charge methods are utilized to handle polarization effects. We find that the ARPM of the ionic liquid reproduces typical (static) electrochemical properties of RTILs. Our model predicts a declining differential capacitance with increasing temperature, which is expected from simple physical arguments. We also compare our ARPM, with the corresponding RPM description, at an elevated temperature (1000 K). We conclude that, even though ion pairing occurs in the ARPM system, reducing the concentration of 'free' ions, it is still better able to screen charge than a corresponding RPM melt. Finally, we evaluate the option to coarse-grain the model even further, by treating the fraction of the ions that form ion pairs implicitly, only through the contribution to the dielectric constant of the corresponding dipolar (ion pair) fluid. We conclude that this primitive representation of ion pairing is not able to reproduce the structures and differential capacitances of the system with explicit ion pairs. The main problem seems to be due to a limited dielectric screening in a layer near the electrode surface, resulting from a combination of orientational restrictions and a depleted dipole density.

Covalent Bonding in the Hydrogen Molecule.

This work addresses the continuing disagreement between two schools of thought concerning the mechanism of covalent bonding. According to Hellmann, Ruedenberg, and Kutzelnigg, covalent bonding is a quantum mechanical phenomenon whereby lowering of the kinetic energy associated with electron sharing, i.e., delocalization, is the key stabilization mechanism. The opposing view of Slater, Feynman, and Bader has maintained that the source of stabilization is electrostatic potential energy lowering due to electron density redistribution to binding regions between nuclei. Following our study of H2+ we present an analogous detailed study of H2 where bonding involves an electron pair with repulsion and correlation playing a significant role in its properties. We use a range of different computational approaches to study and reveal the relevant contributions to bonding as seen in the electron density and corresponding kinetic and potential energy distributions. The energetics associated with the more complex electronic structure of H2, when examined in detail, clearly agrees with the analysis of Ruedenberg; i.e., covalent bonding is due to a decrease in the interatomic kinetic energy resulting from electronic delocalization. Our results support the view that covalent bonding is a quantum dynamical phenomenon requiring a properly quantized kinetic energy to be used in its description.

Ion pairing and phase behaviour of an asymmetric restricted primitive model of ionic liquids.

An asymmetric restricted primitive model (ARPM) of electrolytes is proposed as a simple three parameter (charge q, diameter d, and charge displacement b) model of ionic liquids and solutions. Charge displacement allows electrostatic and steric interactions to operate between different centres, so that orientational correlations arise in ion-ion interactions. In this way the ionic system may have partly the character of a simple ionic fluid/solid and of a polar fluid formed from ion pairs. The present exploration of the system focuses on the ion pair formation mechanism, the relative concentration of paired and free ions and the consequences for the cohesive energy, and the tendency to form fluid or solid phase. In contrast to studies of similar (though not identical) models in the past, we focus on behaviours at room temperature. By MC and MD simulations of such systems composed of monovalent ions of hard-sphere (or essentially hard-sphere) diameter equal to 5 Å and a charge displacement ranging from 0 to 2 Å from the hard-sphere origin, we find that ion pairing dominates for b larger than 1 Å. When b exceeds about 1.5 Å, the system is essentially a liquid of dipolar ion pairs with a small presence of free ions. We also investigate dielectric behaviours of corresponding liquids, composed of purely dipolar species. Many basic features of ionic liquids appear to be remarkably consistent with those of our ARPM at ambient conditions, when b is around 1 Å. However, the rate of self-diffusion and, to a lesser extent, conductivity is overestimated, presumably due to the simple spherical shape of our ions in the ARPM. The relative simplicity of our ARPM in relation to the rich variety of new mechanisms and properties it introduces, and to the numerical simplicity of its exploration by theory or simulation, makes it an essential step on the way towards representation of the full complexity of ionic liquids.

Covalent bonding: the fundamental role of the kinetic energy.

This work addresses the continuing disagreement between two prevalent schools of thought concerning the mechanism of covalent bonding. According to Hellmann, Ruedenberg, and Kutzelnigg, a lowering of the kinetic energy associated with electron delocalization is the key stabilization mechanism. The opposing view of Slater, Feynman, and Bader has maintained that the source of stabilization is electrostatic potential energy lowering due to electron density redistribution to binding regions between nuclei. Despite the large body of accurate quantum chemical work on a range of molecules, the debate concerning the origin of bonding continues unabated, even for H2(+), the simplest of covalently bound molecules. We therefore present here a detailed study of H2(+), including its formation, that uses a sequence of computational methods designed to reveal the relevant contributing mechanisms as well as the spatial density distributions of the kinetic and potential energy contributions. We find that the electrostatic mechanism fails to provide real insight or explanation of bonding, while the kinetic energy mechanism is sound and accurate but complex or even paradoxical to those preferring the apparent simplicity of the electrostatic model. We further argue that the underlying mechanism of bonding is in fact of dynamical character, and analyses that focus on energy do not reveal the origin of covalent bonding in full clarity.

Polyelectrolyte mediated interactions in colloidal dispersions: hierarchical screening, simulations, and a new classical density functional theory.

The pair interaction between two charged colloidal particles, in the presence of a polyelectrolyte as well as simple salt, is analyzed theoretically. Of particular interest is the way in which such a combination of salts can be used to induce a strong, long-range attraction, with at most a minor free energy barrier. We show that the nature of the simple salt is highly relevant, i.e., 2:1, 1:1, and 1:2 salts generate quite different particle interaction free energies at the same overall ionic strength. We adopt several different theoretical levels of description. Defining simulations at the primitive model level with explicit simple salt as our reference, we invoke stepwise coarse-graining with careful evaluations of each approximation. Representing monovalent simple ions by the ionic screening they generate is one such simplification. In order to proceed further, with additional computational savings, we also develop a correlation-corrected classical density functional theory. We analyze the performance of this theory with explicit spherical particles as well as in a flat surface geometry, utilizing Derjaguin's approximation. The calculations are particularly fast in the latter case, facilitating computational savings of many (typically 5-7) orders of magnitude, compared to corresponding simulations with explicit salt. Yet, the predictions are remarkably accurate, and considering the crudeness of the model itself, the density functional theory is a very attractive alternative to simulations.

Monte Carlo studies of drug nucleation 1: formation of crystalline clusters of bicalutamide in water.

A computational method of predicting the effects of the metastability of drug solutions is sought. A simple extension of our in silicio approach to thermodynamic drug solubility is tested on the drug bicalutamide for which we performed vapor pressure measurements complementing earlier measurements of aqueous solubility and crystal-water interfacial tension. The free energy of formation of an N-cluster of the drug molecule is estimated semiempirically by use of an Einstein model of the crystal in which experiment supplies the crystal structure, enthalpy of sublimation, and Einstein frequency of vibration. The rigid drug clusters with N from 2 to 14 are extracted from the bulk crystal by minimization of either cluster energy or radius of gyration. The free energy of hydration is estimated by Monte Carlo simulation combined with simplified response theory based on the OPLS-AA/COMPASS force field for the drug-water interaction and the TIP4P water model. The results have been interpreted in terms of an apparent crystal-water interfacial tension according to classical nucleation theory. The energy-minimal and radius of gyration-minimal clusters seem to give very similar crystal-water interfacial tensions for both the monoclinic and the triclinic polymorph. The interfacial tension of the monoclinic polymorph is significantly higher (by around 20%) than that of the triclinic polymorph in accordance with experiment. For the triclinic polymorph a substantial overestimation of the interfacial tension compared to estimates from crystal nucleation experiments is found, mitigated somewhat by an empirical scaling of the simulated binding energies and free energies of hydration.

Monte Carlo simulations of salt solutions: exploring the validity of primitive models.

An extensive series of Monte Carlo (MC) simulations were performed in order to explore the validity of simple primitive models of electrolyte solutions and in particular the effect of ion size asymmetry on the bulk thermodynamic properties of real salt solutions. Ionic activity and osmotic coefficients were calculated for 1:1, 2:1, and 3:1 electrolytes by using the unrestricted primitive model (UPM); i.e., ions are considered as charged hard spheres of different sizes dissolved in a dielectric continuum. Mean ionic activity and osmotic coefficients calculated by the MC simulations were fitted simultaneously to the experimental data by adjusting only the cation radius while keeping the anion radius fixed at its crystallographic value. Ionic radii were further optimized by systematically varying the cation and anion radii at a fixed sum of ionic radii. The success of this approach is found to be highly salt specific. For example, experimental data (mean ionic activity and osmotic coefficients) of salts which are usually considered as dissociated such as HCl, HBr, LiCl, LiBr, LiClO(4), and KOH were successfully fitted up to 1.9, 2.5, 1.9, 3, 2.5, and 4.5 M concentrations, respectively. In the case of partially dissociated salts such as NaCl, the successful fits were only obtained in a more restricted concentration range. Consistent sets of the best fitted cation radii were obtained for acids, alkali, and alkaline earth halides. A list of recommended ionic radii is also provided. The reliability of the optimized ionic radii was further tested in simulations of the osmotic coefficients of LiCl-NaCl-KCl salt mixtures. A very good agreement between the simulated and experimental data was obtained up to ionic strength of 4.5 M.

In silico prediction of drug solubility: 4. Will simple potentials suffice?

In view of the extreme importance of reliable computational prediction of aqueous drug solubility, we have established a Monte Carlo simulation procedure which appears, in principle, to yield reliable solubilities even for complex drug molecules. A theory based on judicious application of linear response and mean field approximations has been found to reproduce the computationally demanding free energy determinations by simulation while at the same time offering mechanistic insight. The focus here is on the suitability of the model of both drug and solvent, i.e., the force fields. The optimized potentials for liquid simulations all atom (OPLS-AA) force field, either intact or combined with partial charges determined either by semiempirical AM1/CM1A calculations or taken from the condensed-phase optimized molecular potentials for atomistic simulation studies (COMPASS) force field has been used. The results illustrate the crucial role of the force field in determining drug solubilities. The errors in interaction energies obtained by the simple force fields tested here are still found to be too large for our purpose but if a component of this error is systematic and readily removed by empirical adjustment the results are significantly improved. In fact, consistent use of the OPLS-AA Lennard-Jones force field parameters with partial charges from the COMPASS force field will in this way produce good predictions of amorphous drug solubility within 1 day on a standard desktop PC. This is shown here by the results of extensive new simulations for a total of 47 drug molecules which were also improved by increasing the water box in the hydration simulations from 500 to 2000 water molecules.

In silico prediction of drug solubility. 3. Free energy of solvation in pure amorphous matter.

The solubility of drugs in water is investigated in a series of papers. In this work, we address the process of bringing a drug molecule from the vapor into a pure drug amorphous phase. This step enables us to actually calculate the solubility of amorphous drugs in water. In our general approach, we, on one hand, perform rigorous free energy simulations using a combination of the free energy perturbation and thermodynamic integration methods. On the other hand, we develop an approximate theory containing parameters that are easily accessible from conventional Monte Carlo simulations, thereby reducing the computation time significantly. In the theory for solvation, we assume that DeltaG* = DeltaGcav + ELJ + EC/2, where the free energy of cavity formation, DeltaGcav, in pure drug systems is obtained using a theory for hard-oblate spheroids, and ELJ and EC are the Lennard-Jones and Coulomb interaction energies between the chosen molecule and the others in the fluid. The theoretical predictions for the free energy of solvation in pure amorphous matter are in good agreement with free energy simulation data for 46 different drug molecules. These results together with our previous studies support our theoretical approach. By using our previous data for the free energy of hydration, we compute the total free energy change of bringing a molecule from the amorphous phase into water. We obtain good agreement between the theory and simulations. It should be noted that to obtain accurate results for the total process, high precision data are needed for the individual subprocesses. Finally, for eight different substances, we compare the experimental amorphous and crystalline solubility in water with the results obtained by the proposed theory with reasonable success.

In silico prediction of drug solubility: 1. Free energy of hydration.

As a first step in the computational prediction of drug solubility the free energy of hydration, DeltaG*(vw) in TIP4P water has been computed for a data set of 48 drug molecules using the free energy of perturbation method and the optimized potential for liquid simulations all-atom force field. The simulations were performed in two steps, where first the Coulomb and then the Lennard-Jones interactions between the solute and the water molecules were scaled down from full to zero strength to provide physical understanding and simpler predictive models. The results have been interpreted using a theory assuming DeltaG*(vw) = A(MS)gamma + E(LJ) + E(C)/2 where A(MS) is the molecular surface area, gamma is the water-vapor surface tension, and E(LJ) and E(C) are the solute-water Lennard-Jones and Coulomb interaction energies, respectively. It was found that by a proper definition of the molecular surface area our results as well as several results from the literature were found to be in quantitative agreement using the macroscopic surface tension of TIP4P water. This is in contrast to the surface tension for water around a spherical cavity that previously has been shown to be dependent on the size of the cavity up to a radius of approximately 1 nm. The step of scaling down the electrostatic interaction can be represented by linear response theory.

In silico prediction of drug solubility: 2. Free energy of solvation in pure melts.

The solubility of drugs in water is investigated in a series of papers and in the current work. The free energy of solvation, DeltaG*(vl), of a drug molecule in its pure drug melt at 673.15 K (400 degrees C) has been obtained for 46 drug molecules using the free energy perturbation method. The simulations were performed in two steps where first the Coulomb and then the Lennard-Jones interactions were scaled down from full to no interaction. The results have been interpreted using a theory assuming that DeltaG*(vl) = DeltaG(cav) + E(LJ) + E(C)/2 where the free energy of cavity formation, DeltaG(cav), in these pure drug systems was obtained using hard body theories, and E(LJ) and E(C) are the Lennard-Jones and Coulomb interaction energies, respectively, of one molecule with the other ones. Since the main parameter in hard body theories is the volume fraction, an equation of state approach was used to estimate the molecular volume. Promising results were obtained using a theory for hard oblates, in which the oblate axial ratio was calculated from the molecular surface area and volume obtained from simulations. The Coulomb term, E(C)/2, is half of the Coulomb energy in accord with linear response, which showed good agreement with our simulation results. In comparison with our previous results on free energy of hydration, the Coulomb interactions in pure drug systems are weaker, and the van der Waals interactions play a more important role.

Statistical theory of collisional energy transfer in molecular collisions. trans-stilbene deactivation by argon, carbon dioxide, and n-heptane.

Recent advances in experimental techniques have made it possible to measure the full conditional probability density P(E, E') of the energy transfer between two colliding molecules in the gas phase, one of which is highly energized and the other in thermal equilibrium at a given temperature. Data have now become available for trans-stilbene deactivation by the three bath gas molecules Ar, CO2, and n-heptane (C7H16). The initial energies of trans-stilbene are set to 10 000, 20 000, 30 000, and 40 000 cm (-1). The results show that exceptionally large amounts of energy are transferred in each collision. By application of our partially ergodic collision theory (PECT), we find that the energy transfer efficiency betaE ranges from a rather normal value of 0.15 for n-heptane at the highest excitation energy to 0.93-nearly in the ergodic collision limit-for the argon bath gas at high excitation energy. Generally, the PECT produces a good fit of the data except for the nearly elastic peak in the case of n-heptane, where PECT produces a rounded and downshifted peak in contrast to a sharply defined elastic maximum of the monoexponential functional fit produced from the original experimental data obtained by kinetically controlled selective ionization in the work of the group of Luther in Göttingen. This problem is analyzed and found to be related partly to the lack of treatment of glancing collisions in the theory with a remaining uncertainty due to the weak dependence of energy transfer efficiency on nearly elastic collisions. A summary of the present state of understanding shows that collisional activation and deactivation of reactant molecules is more efficient and more statistical than has been previously realized.

PECT model analysis and predictions of experimental collisional energy transfer probabilities P(E',E) and moments for azulene and biphenylene.

Experimental collisional energy transfer data from kinetically controlled selective ionization (KCSI) and ultraviolet absorption (UVA) experiments are analyzed in the framework of the partially ergodic collision theory (PECT). Collisions of azulene and biphenylene with different colliders are investigated as case studies. The downward wings of the P(E',E) energy transfer distributions obtained from the PECT model are fitted to the recently introduced "variable-shape"-exponential 3-parameter functional form of P(E',E) obtained from KCSI experiments, P(E',E) proportional, variant exp[-{(E - E')/(C(0) + C(1)E)}(Y)]. The PECT model is able to reproduce the characteristic dependence of the KCSI "shape parameter" Y on the choice of collider, the energy dependent width of the KCSI P(E',E) distributions, described by alpha(E) = C(0) + C(1)E, and the temperature dependence of the UVA data above room temperature. The statistical approach of PECT obviously captures the essence of large molecule energy transfer at chemically significant energies without the need of knowing specific features of the detailed collision dynamics. It therefore shows promise for predicting the shape of P(E',E) in master equation kernels for larger molecules.

Corrected Debye-Hückel analysis of surface complexation; III. Spherical particle charging including ion condensation.

Statistical mechanics has been used to derive a model for the charging of a spherical particle in a salt solution to complement our experimental studies and gain a deeper understanding of the processes involved in surface complexation. Our chosen model goes beyond the equilibrium constants and the Gouy-Chapmann theory currently used in surface complexation models. The proton adsorption is taken to occur at a harmonic potential well on the surface characterized by a frequency v and a well depth u(0). Outside the particle surface there is a capacitor layer of width w(c) which is impenetrable to the salt ions. The diffuse screening of the charged particle is described by a corrected Debye-Hückel analysis accounting for ion size in the ion-ion interactions. To account also for nonlinear electrostatic response a layer of condensed counterions has been introduced. The criterion for the onset of ion condensation is that the electrostatic field exceeds a linear response criterion. Ion size effects are accounted for in terms of hole-corrected electrostatic energies and excluded volume. The model has been applied to titrated surface charge data on goethite (alpha-FeOOH) at various background concentrations and good agreement between the experimental data and the model was obtained. Both the size of the screening ions and the central particle size were shown to be of importance for the surface charge.

Corrected Debye-Hückel analysis of surface complexation. II. A theory of surface charging.

A theory of surface charging of colloidal particles suspended in an electrolyte solution is presented. The charging at the particle surface is assumed to originate from the adsorption and desorption of protons and is therefore strongly dependent on the acidity of the solution. The surface binding of protons occurs locally at sites of occupancy zero or one that are described by a binding energy u(0) and a three-dimensional vibration of frequency nu. The diffuse screening of ions at the surface is described by the corrected Debye-Hückel analysis assuming linear response. The model contains a capacitor layer close to the charged surface and the finite size of the electrolyte ions is taken into account. The theory has been applied to titrated surface charge data on goethite (alpha-FeOOH) at NaClO(4) background concentrations ranging from 0.01 to 1.0 M. The protonation mechanism used in the modeling of these data corresponds to the 1-pK approach. A very good description of the experimental data was obtained at the highest ionic strength. Close to the pH(pzc) the theory also gave a good description at lower ionic strengths. However, at low salt concentrations and pH values far away from the pH(pzc) the electrostatic potential outside the capacitor layer becomes so high that nonlinear electrostatic effects become important and the theory therefore underestimates the surface charge. These results were compared with model calculations obtained using existing surface complexation models.