Large integration time-step in molecular dynamics simulation artificially enhances the strength of hydrophobic interaction
Abstract
Molecular dynamics simulations are used to compute the potential of mean force (PMF) between two united-atom methane molecules at several temperatures and integration time-steps. At a fixed time-step, the contact minimum of the PMF deepens with increasing temperature, as expected for hydrophobicity driven association. Compared with a time-step of 0.5 fs, one that preserves equipartition, larger time-steps alter the PMF and make the contact minimum more favorable. Thus even relative free energy values are sensitive to time-steps that break equipartition. Using quasichemical theory, we partition the free energy of association into hydrophobic and hydrophilic contributions. The hydrophilic contribution opposes association and is insensitive to the time-step. Thus the artificial enhancement of association at larger time-steps comes entirely from the hydrophobic contribution. We explain this behavior through the temperature dependence of the internal pressure of the liquid. The same analysis also accounts for earlier observations of the liquid's p-V behavior under conditions that break equipartition.
Related
- Guest-guest interactions and their impact on molecular mobility: a simulation study of 5-fluorouracil adsorption and diffusion in hydrophobic zeolites
- Hydrophobic Hydration at Fluorinated and Hydrocarbon Self-Assembled Monolayers Investigated by In Situ Soft X-Ray Emission Spectroscopy and Molecular Dynamics Simulation.
- Unveiling the Molecular Driving Forces of Pollutant Extraction by Hydrophobic Eutectic Solvents.
- DBS‐Based Hydrophobic Nonionic Supramolecular Eutectogels: Bridging Molecular Insights to Water Remediation
- Integration of Network Toxicology and Machine Learning to Prioritize Candidate Molecular Features Associated with NNK-Related Lung Cancer
BibTeX
@misc{asthagiri2026large,
title = {{Large integration time-step in molecular dynamics simulation artificially enhances the strength of hydrophobic interaction}},
author = {Dilipkumar N. Asthagiri},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.30109},
eprint = {2609.30109},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.30109},
}
RIS
TY - UNPB
TI - Large integration time-step in molecular dynamics simulation artificially enhances the strength of hydrophobic interaction
AU - Dilipkumar N. Asthagiri
PY - 2026
DA - 2026/09/24
PB - arXiv
DO - 10.48550/arxiv.2609.30109
UR - https://arxiv.org/abs/2609.30109
AB - Molecular dynamics simulations are used to compute the potential of mean force (PMF) between two united-atom methane molecules at several temperatures and integration time-steps. At a fixed time-step, the contact minimum of the PMF deepens with increasing temperature, as expected for hydrophobicity driven association. Compared with a time-step of 0.5 fs, one that preserves equipartition, larger time-steps alter the PMF and make the contact minimum more favorable. Thus even relative free energy values are sensitive to time-steps that break equipartition. Using quasichemical theory, we partition the free energy of association into hydrophobic and hydrophilic contributions. The hydrophilic contribution opposes association and is insensitive to the time-step. Thus the artificial enhancement of association at larger time-steps comes entirely from the hydrophobic contribution. We explain this behavior through the temperature dependence of the internal pressure of the liquid. The same analysis also accounts for earlier observations of the liquid's p-V behavior under conditions that break equipartition.
ER -