Baochi Doan,Patrick K. Schelling
Abstract:Dissipation and adhesion are important in many areas of materials science, including friction and lubrication, cold spray deposition, and micro-electromechanical systems (MEMS). Another interesting problem is the adhesion of mineral grains during the early stages of planetesimal formation in the early solar system. Molecular-dynamics (MD) simulation has often been used to elucidate dissipative properties, most often in the simulation of sliding friction. In this paper, we demonstrate how the reversible and irreversible work associated with interactions between planar surfaces can be calculated using the dynamical contact simulation approach based on MD and empirical potentials. Moreover, it is demonstrated how the approach can obtain the free-energy $\Delta A(z)$ as a function of separation between two slabs using the Jarzynksi equality applied to an ensemble of trajectories which deviate significantly from equilibrium. Furthermore, the dissipative work can also be obtained using this method without the need to compute an entire cycle from approach to retraction. It is expected that this technique might be used to efficiently compute dissipative properties which might enable the use of more accurate approaches including density-functional theory. In this paper, we present results obtained for forsterite surfaces both with and without MgO-vacancy surface defects. It is shown that strong dissipation is possible when MgO-vacancy defects are present. The mechanism for strong dissipation is connected to the tendency of less strongly-bound surface units to undergo large displacements including mass transfer between the two surfaces. Systems with strong dissipation tend to exhibit a long-tailed distribution rather than the Gaussian distribution often anticipated in near-equilibrium applications of the JE.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to understand the dissipation and adhesion hysteresis phenomena between (010) olivine surfaces. Specifically, the authors use molecular dynamics simulations and the Jarzynski equation to calculate the reversible and irreversible work between two planar surfaces and explore the changes of these properties under different conditions.
### Research Background and Problem Description
1. **Importance of Dissipation and Adhesion**
- In many fields of materials science, such as friction and lubrication, cold - spray deposition, micro - electro - mechanical systems (MEMS), etc., dissipation and adhesion are very important phenomena.
- Another interesting example is the adhesion problem of mineral particles in the early stage of planetesimal formation in the early solar system, which is not very clear at present.
2. **Limitations of Existing Methods**
- Molecular dynamics (MD) simulations are usually used to study the dissipation characteristics in sliding friction, but these methods are often unable to directly calculate the free - energy difference \(\Delta A(z)\) at the separation distance \(z\), because these processes are usually irreversible.
3. **Research Objectives**
- By combining molecular dynamics simulations and the Jarzynski equation, calculate the free - energy difference \(\Delta A(z)\) and the dissipation work \(W_{\text{diss}}\) between two planar surfaces.
- Explore the influence of MgO vacancy defects on dissipation and adhesion, especially the strong dissipation mechanism and its relationship with large - displacement of surface units and mass transfer.
### Research Methods
1. **Molecular Dynamics Simulations**
- Use the LAMMPS simulation code for calculations, taking into account Coulomb interactions and other pairwise potential - energy terms.
- Simulate the interactions of (010) olivine surfaces under different conditions, including the cases with and without MgO vacancy defects.
2. **Application of the Jarzynski Equation**
- Through the Jarzynski equation (JE), calculate the free - energy difference \(\Delta A(z)\) from a set of trajectories far from equilibrium.
- Calculate the dissipation work \(W_{\text{diss}}(z)\) and explore its changes at different positions \(z\).
### Main Results
1. **Surfaces without Defects**
- For the (010) olivine surfaces without defects, the results show that there is an obvious adhesion hysteresis phenomenon during the approaching and retracting processes.
- The changes of the free - energy difference \(\Delta A(z)\) and the dissipation work \(W_{\text{diss}}\) with the separation distance \(z\) are as expected.
2. **Surfaces with MgO Vacancy Defects**
- When there are MgO vacancy defects, the strong dissipation phenomenon is more obvious, indicating that the surface units are more likely to have large displacements and mass transfers.
- The results show that the defective surfaces exhibit a longer - tailed distribution rather than a Gaussian distribution in the retracting stage, which further verifies the existence of strong dissipation.
### Conclusion
This study successfully calculates the free - energy difference and the dissipation work between (010) olivine surfaces through molecular dynamics simulations and the Jarzynski equation, and reveals the influence of MgO vacancy defects on dissipation and adhesion. These results are helpful for better understanding the microscopic interaction mechanisms at material interfaces and provide new methods and technical means for future research.
### Formula Summary
- **Dissipation Work Formula**
\[
W_{\text{diss}}=\langle W\rangle-\Delta A
\]
- **Jarzynski Equation**
\[
e^{-\beta\Delta A}=\langle e^{-\beta W}\rangle
\]
where \(\beta = \frac{1}{k_B T}\), \(\Delta A\) is the free - energy difference between two equilibrium states, and \(W\) is the work done by external forces.
- **Free - Energy Difference Calculation**
\[
\Delta A(z)=-\frac{1}{\beta}\ln\left\langle e^{-\beta[W_0(z)-K(z)]}\right\rangle