A method for calculating solvation structure on a sample surface from a force curve between a probe and the sample: One-dimensional version

Ken-ichi Amano
DOI: https://doi.org/10.48550/arXiv.1209.0303
2013-03-06
Abstract:Recent surface force apparatus (SFA) and atomic force microscopy (AFM) can measure force curves between a probe and a sample surface in solvent. The force curve is thought as the solvation structure in some articles, because its shape is generally oscilltive and pitch of the oscillation is about the same as diameter of the solvent. However, it is not the solvation structure. It is only the force between the probe and the sample surface. Therefore, this brief paper presents a method for calculating the solvation structure from the force curve. The method is constructed by using integral equation theory, a statistical mechanics of liquid (Ornstein-Zernike equation coupled by hypernetted-chain closure). This method is considered to be important for elucidation of the solvation structure on a sample surface.
Chemical Physics,Biological Physics,Computational Physics,Data Analysis, Statistics and Probability
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to calculate the solvation structure from the force curve between the probe and the sample surface. The author points out that although some literatures believe that the shape of the force curve reflects the solvation structure because its oscillation period is approximately equal to the diameter of the solvent molecule, in fact it is only the force between the probe and the sample surface. Therefore, this paper proposes a method based on the integral equation theory (especially the Ornstein - Zernike equation combined with the hypernetted - chain closure approximation) to calculate the solvation structure from the force curve. This method is of great significance for revealing the solvation structure on the sample surface. Specifically, this method is implemented through the following steps: 1. **Calculate the mean - force potential**: Starting from the force curve measured experimentally, use integration to calculate the mean - force potential \(\Phi(r)\). 2. **Extract the solvent - mediated potential**: Subtract the direct interaction potential \(u_{MP}(r)\) between the probe and the sample from the mean - force potential to obtain the solvent - mediated potential \(\varphi(r)\). 3. **Obtain \(w_{MP}\)**: Calculate \(w_{MP}(r)\) according to the solvent - mediated potential. 4. **Calculate CPS**: Use the Ornstein - Zernike equation and the hypernetted - chain closure approximation to calculate CPS. 5. **Introduce a trial function**: Set the trial function form of the total correlation function \(h'_{MS}\) between the sample and the solvent. 6. **Perform a three - dimensional Fourier transform**: Perform a three - dimensional Fourier transform on \(h'_{MS}\) to obtain \(H'_{MS}\). 7. **Calculate \(W'_{MP}\)**: Use \(H'_{MS}\) and CPS to calculate \(W'_{MP}\) through the Ornstein - Zernike equation. 8. **Inverse Fourier transform**: Perform an inverse Fourier transform on \(W'_{MP}\) to obtain \(w'_{MP}\). 9. **Evaluate the difference**: Calculate the difference \(\omega\) between \(w'_{MP}\) and \(w_{MP}\). 10. **Iterative optimization**: If \(\omega\) is small enough, then the current \(h'_{MS}\) trial function is the desired \(h_{MS}\); otherwise, return to step 5, adjust the parameters in the trial function and recalculate. Through the above steps, this method can accurately calculate the solvation structure on the sample surface from the force curve measured experimentally, thus providing a new tool for understanding and studying the solvation phenomenon.