About this Journal Submit a Manuscript Table of Contents
ISRN Thermodynamics
Volume 2012 (2012), Article ID 204104, 3 pages
doi:10.5402/2012/204104
Research Article

The Born Formula Describes Enthalpy of Ions Solvation

Institute of Chemical Kinetics and Combustion, Novosibirsk State University, Novosibirsk 630090, Russia

Received 3 February 2012; Accepted 26 February 2012

Academic Editor: M. Manciu

Copyright © 2012 Nikolai Bazhin. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The process of ion solvation has been studied in the reversible system using the Van't Hoff equilibrium box. It is shown that the Born formula for solvation energy describes a change in enthalpy rather than in Gibbs energy.

1. Introduction

Solvation energy is usually taken as the difference in work which results from ion charging in both vacuum and liquid phase due to a mental process [1]. For an approximate, simple estimate of solvation energy, one can use the Born formula [1] according to which the solvation energy ( 𝐸 s o l v ) is the difference in the energy of charging the sphere of radius 𝑎 to charge 𝑍 𝑒 in both vacuum and medium with a dielectric constant 𝜀 𝐸 s o l v = 𝑍 2 𝑒 2 8 𝜋 𝜀 0 𝑎 1 1 𝜀 𝑁 A , ( 1 ) where 𝜀 0 is an electric constant, and 𝑁 A is the Avogadro number. According to Born, (1) is the work obtained by transferring a charged sphere from vacuum to condensed medium. The work manifests itself in solvation heat [1]. Since the system temperature is constant, the heat emitted passes to a thermostat. It is assumed then that Born has taken solvation energy as the change in either internal energy or enthalpy upon solvation. However, discussing Born’s work, the authors [2] assert that without any doubts, value (1) calculated by Born is the maximum work performed by electric forces during ion dissolving. Therefore, in [2], (1) is equated with the change in the Gibbs energy which allowed the authors [2] to calculate solvation enthalpy as Δ 𝐻 s o l v = 𝑍 2 𝑒 2 8 𝜋 𝜀 0 𝑎 1 1 𝜀 𝑇 𝜀 2 𝑑 𝜀 𝑁 𝑑 𝑇 A . ( 2 ) Equation (1) fails to quantitatively describe ion solvation in numerical form. The advantage is that it can be used to qualitatively estimate the dependence of solvation energy on ion size, medium’s dielectric constant, and temperature. It is of importance then to ascertain whether (1) is the Gibbs function or the solvation enthalpy.

We will omit here the problems of the applicability of (1), considered, for example, in [310].

Since the energy of charging is the work, formula (1) is commonly accepted as the change in the Gibbs function upon solvation which causes, however, some doubts for the following reasons.(1)To calculate the change in the Gibbs function, one should accurately determine the initial and final states of the system. To derive formula (1), the charging of only one particle is considered, and the result obtained is multiplied by the Avogadro number to calculate the charging energy per one mole. As compared with enthalpy, in the initial state, the Gibbs function is pressure dependent, and, in the final state, it depends on ion concentration. Thus, it is impossible to describe the change in the Gibbs function in terms of (1).(2)There is no evidence that formula (1) is sure to describe the change in the Gibbs function rather than the change in other values, for example, in enthalpy.

The goal of the present work is to solve this problem by means of the Van’t Hoff equilibrium box (VHEB) proposed by Van’t Hoff to derive an equilibrium constant. To this end, the manuals on chemical thermodynamics [11, 12] consider the VHEB as an old alternative of the method of chemical potentials. At present, the VHEB is very useful to solve qualitative problems. Thus, for example, in [13], the VHEB was used to demonstrate that the reversible work at 𝑇 , 𝑝 = C o n s t in the systems, including chemical reactions, arises only due to the thermal energy of the environment, and the energy evolved by reactions is not employed to produce useful work.

2. Solvation

The change in the Gibbs function between two states of the system is independent of the type of transition between these states. This transition can be either reversible or irreversible. However, to calculate the change in the Gibbs function, the transition should be reversible because, in this case, the calculation of the change in the Gibbs function becomes more available. The reversible solvation process is performed in this paper by analogy as described in [13] for coal oxidation. The system under study is considered ideal.

Let us describe first the standard initial and the final states. The final state is the electrolyte solution. Due to Alberty [14], the standard state of an electrolyte solution is “the hypothetically ideal solution at unit activity in which ions have the properties they do at infinite dilution.” We assume then that a standard state of ions in gas is a hypothetical gas at 1 bar, in which the ions have the properties they do at infinite expansion. Thus, in these hypothetical standard states, the electrostatic interaction of ions is cut off which alleviates the calculation of thermodynamic values, describing the process of solvation.

Consider now the solvation process. Gaseous ions are assumed to occupy a very large volume. In addition, the volume of the liquid-phase system is also assumed very large. The system is located in a thermostat. Let us reversibly convert the part of the gaseous system from the initial state into the final one. We calculate the change in the thermodynamic functions upon transition of a hypothetical system from the standard state in the gas phase to that in the liquid state in the course of a conceptual process, consisting of two stages.

Step 1. At step 1 we take 1 mole of ions in the gaseous state at 𝑝 s t = 1 bar and expand the ion gas reversible to pressure 𝑝 n , at which the chemical potential of ions in gas becomes equal to that in solution. The reversible process is accompanied obligatory by useful work production ( 𝑤 u s e f u l ) , and the change in the Gibbs function ( Δ 𝐺 1 ) of the ions is equal to 𝑤 u s e f u l Δ 𝐺 1 = 𝑤 u s e f u l 𝑝 = 𝑅 𝑇 l n n 𝑝 s t . ( 3 ) The change in enthalpy ( Δ 𝐻 1 ) of ideal ion gas expansion is of the form Δ 𝐻 1 = 0 . ( 4 ) The value of 𝑝 n is unknown. Value (3) cannot be given in terms of system parameters in the framework of the approach used. Upon expansion, the reversible work is performed which can be converted into the potential energy of the load. The work is performed only by the thermal thermostat energy because, in the system, there are no other energy sources. Thus, the change in the thermostat enthalpy ( Δ 𝐻 1 t h r e r m o s t a t ) during gas expansion at the first step is Δ 𝐻 1 t h r e r m o s t a t 𝑝 = 𝑅 𝑇 l n n 𝑝 s t . ( 5 )

Step 2. At the final pressure 𝑝 n , we slowly introduce ions from the gas phase to solution through, for example, a semipermeable membrane. In this case, the Gibbs function does not change because of the equality of the chemical potentials of ions in the gas and liquid phases, and, thus, Δ 𝐺 2 = 0 . ( 6 ) Hence, the full change in the Gibbs function upon solvation process ( Δ 𝐺 s o l v ) is equal to useful work of the ion gas expansion during first step Δ 𝐺 s o l v = Δ 𝐺 1 + Δ 𝐺 2 𝑝 = 𝑅 𝑇 l n n 𝑝 s t . ( 7 ) It is worth noting that no relation was established between Δ 𝐺 s o l v and formula (1).
The process of solvation results in heat release. The reaction of ion dissolution is similar to that of coal oxidation [13]. The value of the change in the ions enthalpy at the second step ( Δ 𝐻 2 ) is calculated from (1) because, in terms of the law of energy conservation, the difference in work on sphere charging in both vacuum and liquid should convert into another type of energy, represented in our case, by the solvation reaction heat Δ 𝐻 2 𝑍 = 2 𝑒 2 8 𝜋 𝜀 0 𝑎 1 1 𝜀 𝑁 A . ( 8 ) Hence, the full enthalpy of solvation process ( Δ 𝐻 s o l v ) is calculated via the Born formula (1) Δ 𝐻 s o l v = Δ 𝐻 1 + Δ 𝐻 2 𝑍 = 2 𝑒 2 8 𝜋 𝜀 0 𝑎 1 1 𝜀 𝑁 A . ( 9 ) The change in the thermostat enthalpy at the second step ( Δ 𝐻 2 t h r e r m o s t a t ) during gas expansion is Δ 𝐻 2 t h r e r m o s t a t = Δ 𝐻 2 . ( 1 0 ) The full thermostat enthalpy change in the reversible process is Δ 𝐻 t h r e r m o s t a t = Δ 𝐻 1 t h r e r m o s t a t + Δ 𝐻 2 t h r e r m o s t a t 𝑝 = 𝑅 𝑇 l n n 𝑝 s t + 𝑍 2 𝑒 2 8 𝜋 𝜀 0 𝑎 1 1 𝜀 𝑁 A . ( 1 1 ) In the case of the irreversible process ( 𝑤 u s e f u l = 0 ) , the full thermostat enthalpy change is equal to Δ 𝐻 t h r e r m o s t a t = 𝑍 2 𝑒 2 8 𝜋 𝜀 0 𝑎 1 1 𝜀 𝑁 A , Δ 𝐻 t h r e r m o s t a t = Δ 𝐻 s o l v . ( 1 2 ) The change in the Gibbs function (7) cannot be equated to the negative value of solvation energy (1) because the change in the Gibbs function in the reversible process determines the value of useful, reversible work rather than of the heat emitted.
Thus, the Born function (1) can be used to calculate only the solvation enthalpy.

References

  1. M. Born, “Volumen und Hydratationswarme der Ionen,” Zeitschrift für Physik, vol. 1, pp. 45–48, 1920.
  2. N. Bjerrum and E. Larson, “Studien uber Ionenverteilungskoeffizienten I.,” Zeitschrift für Physikalische Chemie, vol. 127, pp. 358–384, 1927.
  3. P. W. Atkins and A. J. MacDermott, “The born equation and ionic solvation,” Journal of Chemical Education, vol. 59, no. 4, pp. 359–360, 1982. View at Scopus
  4. P. M. Qureshi, R. K. Varshney, and S. B. Singh, “Empirical correlations for the evaluation of the free energies of solvation of some gaseous monovalent ions,” Journal of Chemical Education, vol. 66, no. 11, pp. 903–906, 1989. View at Scopus
  5. D. Bashford and D. A. Case, “Generalized born models of macromolecular solvation effects,” Annual Review of Physical Chemistry, vol. 51, pp. 129–152, 2000. View at Scopus
  6. P. K. Yang and C. Lim, “The importance of excluded solvent volume effects in computing hydration free energies,” Journal of Physical Chemistry B, vol. 112, no. 47, pp. 14863–14868, 2008. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  7. B. Roux, H. A. Yu, and M. Karplus, “Molecular basis for the born model of ion solvation,” Journal of Physical Chemistry, vol. 94, no. 11, pp. 4683–4688, 1990. View at Scopus
  8. G. F. de Lima, H. A. Duarte, and J. R. Pliego, “Dynamical discrete/continuum linear response shells theory of solvation: convergence test for NH4+ and OH ions in water solution using DFT and DFTB methods,” Journal of Physical Chemistry B, vol. 114, no. 48, pp. 15941–15947, 2010. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  9. C. S. Babu and C. Lim, “Theory of ionic hydration: insights from molecular dynamics simulations and experiment,” Journal of Physical Chemistry B, vol. 103, no. 37, pp. 7958–7968, 1999. View at Scopus
  10. F. Hirata, P. Redfern, and R. M. Levy, “Viewing the born model for ion hydration through a microscope,” International Journal of Quantum Chemistry, vol. 15, pp. 179–190, 1988.
  11. Ya. I. Gerasimov, Ed., Physical Chemistry, Chimiya, Moscow, Russia, 1980.
  12. L. E. Steiner, Introduction to Chemical Thermodynamic, McGraw-Hill, New York, NY, USA, 2nd edition, 1948.
  13. N. M. Bazhin and V. N. Parmon, “Conversion of the chemical reaction energy into useful work in the Van't Hoff equilibrium box,” Journal of Chemical Education, vol. 84, no. 6, pp. 1053–1055, 2007. View at Scopus
  14. R. A. Alberty, Physical Chemislry, John Wiley & Sons, New York, NY, USA, 7th edition, 1987.