Solubility Prediction of Organic Molecules with Molecular Dynamics Simulations
OPEN ACCESS
Loading...
Author / Producer
Date
2021-09-01
Publication Type
Journal Article
ETH Bibliography
yes
Citations
Altmetric
OPEN ACCESS
Data
Abstract
We present a molecular dynamics simulation method for the computation of the solubility of organic crystals in solution. The solubility is calculated based on the equilibrium free energy difference between the solvated solute and its crystallized state at the crystal surface kink site. To efficiently sample the growth and dissolution process, we have carried out well-tempered metadynamics simulations with a collective variable that captures the slow degrees of freedom, namely, the solute diffusion to and adsorption at the kink site together with the desolvation of the kink site. Simulations were performed at different solution concentrations using constant chemical potential molecular dynamics, and the solubility was identified at the concentration at which the free energy values between the grown and dissolved kink states were equal. The effectiveness of this method is demonstrated by its success in reproducing the experimental trends of solubility of urea and naphthalene in a variety of solvents.
Permanent link
Publication status
published
External links
Editor
Book title
Journal / series
Volume
21 (9)
Pages / Article No.
5198 - 5205
Publisher
American Chemical Society
Event
Edition / version
Methods
Software
Geographic location
Date collected
Date created
Subject
Organisational unit
03484 - Mazzotti, Marco (emeritus) / Mazzotti, Marco (emeritus)