Ion Pair Dissociation Dynamics in an Aqueous Premelting Layer

Samuel P. Niblett,David T. Limmer

DOI: https://doi.org/10.48550/arXiv.2012.09881

2020-12-18

Abstract:Using molecular dynamics simulations and methods of importance sampling, we study the thermodynamics and dynamics of sodium chloride in the aqueous premelting layer formed spontaneously at the interface between ice and its vapor. We uncover a hierarchy of timescales that characterize the relaxation dynamics of this system, spanning the picoseconds of ionic motion to the 10s-100s of nanoseconds associated with fluctuations of the liquid-crystal interface in their presence. We find that ions distort both local interfaces, incurring restoring forces that result in the ions preferentially residing in the middle of the layer. While ion pair dissociation is thermodynamically favorable, these structural and dynamic effects cause its rate to vary by over an order of magnitude through the layer, with a maximum rate significantly depressed from the corresponding bulk value. The solvation environment of ions in the premelting layer is distinct from that in a bulk liquid, being dominated by slow reorganization of water molecules and a water structure intermediate between ice and its melt.

Statistical Mechanics,Chemical Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to understand the kinetic and thermodynamic characteristics of the dissociation of ion pairs (such as sodium chloride) in the premelting layer of ice. Specifically, through molecular dynamics simulations and importance sampling methods, the researchers explored the dissociation behavior of sodium chloride in the premelting layer of water that spontaneously forms at the ice - vapor interface.

### Main problems
1. **Changes in dissociation rate**: The paper reveals that in the premelting layer of water, the dissociation rate of ion pairs is significantly different from that in the bulk liquid. Although dissociation is thermodynamically favorable, structural and kinetic effects cause the dissociation rate to vary by more than an order of magnitude at different positions within the layer.
2. **Influence of solvent environment**: The solvent environment in the premelting layer of water is different from that in the bulk liquid, manifested as the slow reorganization of water molecules and the water structure between ice and liquid water. These factors lead to the complexity of the dissociation behavior of ion pairs.
3. **Influence of interface effects on dissociation**: The study found that the dissociation of ion pairs is affected not only by the solvent environment but also by interface effects. For example, when ions approach the liquid - gas interface, the local interface will deform, generating a restoring force, causing the ions to tend to stay in the middle part of the layer.

### Research background
Chemical reactions at the ice - air interface have an important impact on atmospheric chemistry, especially in polar atmospheric chemistry. Understanding the influence of these interfaces on chemical reactions helps to explain some important atmospheric phenomena, such as the photolysis of halogen molecules and their destruction of the ozone layer. In addition, chemical reactions on the ice surface can also affect the concentration of trace gases in the atmosphere, such as the oxidation of sulfur dioxide to sulfuric acid, which in turn leads to the acidification of ice and snow.

### Methods and models
The researchers used molecular dynamics simulations and classical theories to study this problem. Specific methods include:
- Using the TIP4P - 2005 water model and the OPLS - aa force field to simulate water, sodium ions, and chloride ions.
- Defining coarse - grained order parameters to identify the liquid - crystal interface and the liquid - gas interface.
- Calculating the free - energy surface to describe the thermodynamic properties of ion - pair dissociation.
- Calculating the dissociation rate through the Bennett - Chandler method and comparing the results of transition - state theory (TST), Kramers theory, and Bennett - Chandler theory.

### Conclusions
The research results show that in the premelting layer of water, the dissociation rate of ion pairs is significantly lower than that in the bulk liquid, and shows significant differences at different positions within the layer. These differences are mainly caused by the particularity of the solvent environment and interface effects. Understanding these phenomena is of great significance for further exploring the chemical reaction mechanisms at the ice - air interface.

### Formula summary
- **Coarse - grained order parameter**:
\[
\phi_{\alpha}(r; \xi) = \int dr' \phi_{\alpha}(r') e^{-\frac{(r - r')^2}{2\xi^2}} (2\pi\xi)^{-3/2}
\]
- **Height autocorrelation function**:
\[
C_h(t) = \langle \delta h_{\alpha}(x, 0) \delta h_{\alpha}(x, t) \rangle
\]
- **Transition - state theory estimate**:
\[
k_{\text{TST}}(\bar{z}) = \frac{k_BT}{2\mu} e^{-\beta F(R^\ddagger, \bar{z})} \int_{-\infty}^{R^\ddagger} dR e^{-\beta F(R, \bar{z})}
\]
- **Transport coefficient**:
\[
\kappa(\bar{z}) = \lim_{t \to \infty} \frac{\langle \dot{R}(0) \theta[R(t) - R^\ddagger(\bar{z})] \rangle^\ddagger}{\langle |\dot{R}(0)| \rangle^\ddagger / 2}
\]

These formulas help to quantify the behavior of ion pairs in the premelting layer of water and explain the changes in their dissociation rates.